Low formation energy and kinetic barrier of Stone-Wales defect in infinite and finite silicene

Stone-Wales (SW) defects in materials having hexagonal lattice are the most common topological defects that affect the electronic and mechanical properties. Using first principles density functional theory based calculations, we study the formation energy and kinetic barrier of SW-defect in infinite and finite sheets of silicene. The formation energies as well as the barriers in both the cases are significantly lower than those of graphene. Furthermore, compared with the infinite sheets, the energy barriers and formation energies are lower for finite sheets. However, due to low barriers these defects are expected to heal out of the finite sheets. (C) 2013 Elsevier B.V. All rights reserved.

Prediction of the stiffness of nanoclay-polypropylene composites using a monte carlo finite element analysis approach

The present work deals with the prediction of stiffness of an Indian nanoclay-reinforced polypropylene composite (that can be termed as a nanocomposite) using a Monte Carlo finite element analysis (FEA) technique. Nanocomposite samples are at first prepared in the laboratory using a torque rheometer for achieving desirable dispersion of nanoclay during master batch preparation followed up with extrusion for the fabrication of tensile test dog-bone specimens. It has been observed through SEM (scanning electron microscopy) images of the prepared nanocomposite containing a given percentage (3–9% by weight) of the considered nanoclay that nanoclay platelets tend to remain in clusters. By ascertaining the average size of these nanoclay clusters from the images mentioned, a planar finite element model is created in which nanoclay groups and polymer matrix are modeled as separate entities assuming a given homogeneous distribution of the nanoclay clusters. Using a Monte Carlo simulation procedure, the distribution of nanoclay is varied randomly in an automated manner in a commercial FEA code, and virtual tensile tests are performed for computing the linear stiffness for each case. Values of computed stiffness modulus of highest frequency for nanocomposites with different nanoclay contents correspond well with the experimentally obtained measures of stiffness establishing the effectiveness of the present approach for further applications.

Vertical uplift resistance of two closely spaced horizontal strip anchors embedded in cohesive-frictional weightless medium

The vertical uplift resistance of two closely spaced horizontal strip plate anchors has been investigated by using lower and upper bound theorems of the limit analysis in combination with finite elements and linear optimization. The interference effect on uplift resistance of the two anchors is evaluated in terms of a nondimensional efficiency factor (eta(c)). The variation of eta(c) with changes in the clear spacing (S) between the two anchors has been established for different combinations of embedment ratio (H/B) and angle of internal friction of the soil (phi). An interference of the anchors leads to a continuous reduction in uplift resistance with a decrease in spacing between the anchors. The uplift resistance becomes a minimum when the two anchors are placed next to each other without any gap. The critical spacing (S-cr) between the two anchors required to eliminate the interference effect increases with an increase in the values of both H/B and phi. The value of S-cr was found to lie approximately in the range 0.65B-1.5B with H/B = 1 and 11B-14B with H/B = 7 for phi varying from 0 degrees to 30 degrees.

Evolution of Nanocrystalline Ba2NaNb5O15 in 2BaO-0.5Na(2)O-2.5Nb(2)O(5)-4.5B(2)O(3) Glass System and Its Refractive Index and Band Gap Tunability

Monophasic Ba2NaNb5O15 was crystallized at nanometer scale (12-36 nm) in 2BaO-0.5Na(2)O-2.5Nb(2)O(5)- 4.5B(2)O(3) glass system. To begin with, optically transparent glasses, in this system, were fabricated via the conventional melt. quenching technique. The amorphous and glassy characteristics of the as-quenched samples were respectively confirmed by X-ray powder diffraction and differential thermal analyses. Nearly homogeneous distribution of Ba2NaNb5O15 (BNN) nanocrystals associated with tungsten bronze structure akin to their bulk parent structure was accomplished by subjecting the as-fabricated glasses to appropriate heat-treatment temperatures. Indeed transmission electron microscopy (TEM) carried out on these samples corroborated the presence of Ba2NaNb5O15 nanocrystals dispersed in a continuous glass matrix. The as-quenched glasses were similar to 75% transparent in the visible range of the electromagnetic spectrum. The optical band gap and refractive index were found to have crystallite size (at nanoscale) dependence. The optical band gap increased with the decrease in crystallite size. The refractive indices of the glass nanocrystal composites as determined by Brewster angle method were rationalized using different empirical models. The refractive index dispersion with wavelength of light was analyzed on the basis of the Sellmeier relations. At room temperature under UV excitation (355 nm) these glass nanocrystal composites displayed violet-blue emission which was ascribed to the defects states.

Upper bound solution for pullout capacity of vertical anchors in sand using finite elements and limit analysis

The horizontal pullout capacity of vertical anchors embedded in sand has been determined by using an upper bound theorem of the limit analysis in combination with finite elements. The numerical results are presented in nondimensional form to determine the pullout resistance for various combinations of embedment ratio of the anchor (H/B), internal friction angle (ϕ) of sand, and the anchor-soil interface friction angle (δ). The pullout resistance increases with increases in the values of embedment ratio, friction angle of sand and anchor-soil interface friction angle. As compared to earlier reported solutions in literature, the present solution provides a better upper bound on the ultimate collapse load.

Finite-temperature dynamics of vortices in Bose-Einstein condensates

We study the dynamics of a single vortex and a pair of vortices in quasi two-dimensional Bose-Einstein condensates at finite temperatures. To this end, we use the stochastic Gross-Pitaevskii equation, which is the Langevin equation for the Bose-Einstein condensate. For a pair of vortices, we study the dynamics of both the vortex-vortex and vortex-antivortex pairs, which are generated by rotating the trap and moving the Gaussian obstacle potential, respectively. Due to thermal fluctuations, the constituent vortices are not symmetrically generated with respect to each other at finite temperatures. This initial asymmetry coupled with the presence of random thermal fluctuations in the system can lead to different decay rates for the component vortices of the pair, especially in the case of two corotating vortices.

Goldstone modes and bifurcations in phase-separated binary condensates at finite temperature

We show that the third Goldstone mode, which emerges in binary condensates at phase separation, persists to higher interspecies interaction for density profiles where one component is surrounded on both sides by the other component. This is not the case with symmetry-broken density profiles where one species is entirely to the left and the other is entirely to the right. We, then, use Hartree-Fock-Bogoliubov theory with Popov approximation to examine the mode evolution at T not equal 0 and demonstrate the existence of mode bifurcation near the critical temperature. The Kohn mode, however, exhibits deviation from the natural frequency at finite temperatures after the phase separation. This is due to the exclusion of the noncondensate atoms in the dynamics.

A finite element method for the Saint-Venant torsion and bending problems for prismatic beams

In this work, we present a finite element formulation for the Saint-Venant torsion and bending problems for prismatic beams. The torsion problem formulation is based on the warping function, and can handle multiply-connected regions (including thin-walled structures), compound and anisotropic bars. Similarly, the bending formulation, which is based on linearized elasticity theory, can handle multiply-connected domains including thin-walled sections. The torsional rigidity and shear centers can be found as special cases of these formulations. Numerical results are presented to show the good coarse-mesh accuracy of both the formulations for both the displacement and stress fields. The stiffness matrices and load vectors (which are similar to those for a variable body force in a conventional structural mechanics problem) in both formulations involve only domain integrals, which makes them simple to implement and computationally efficient. (C) 2014 Elsevier Ltd. All rights reserved.

REDUCING COMPUTATIONAL EFFORT IN SOLVING GEOMECHANICS PROBLEMS WITH UPPER BOUND FINITE ELEMENTS LIMIT ANALYSIS AND LINEAR OPTIMIZATION

This paper presents a simple technique for reducing the computational effort while solving any geotechnical stability problem by using the upper bound finite element limit analysis and linear optimization. In the proposed method, the problem domain is discretized into a number of different regions in which a particular order (number of sides) of the polygon is chosen to linearize the Mohr-Coulomb yield criterion. A greater order of the polygon needs to be selected only in that region wherein the rate of the plastic strains becomes higher. The computational effort required to solve the problem with this implementation reduces considerably. By using the proposed method, the bearing capacity has been computed for smooth and rough strip footings and the results are found to be quite satisfactory.

Efficiency limitations in a low band-gap diketopyrrolopyrrole-based polymer solar cell

We report the performance and photophysics of a low band-gap diketopyrrolopyrrole-based copolymer used in bulk heterojunction devices in combination with PC71BM. We show that the short lifetime of photogenerated excitons in the polymer constitutes an obstacle towards device efficiency by limiting the diffusion range of the exciton to the donor-acceptor heterojunction. We employ ultrafast transient-probe and fluorescence spectroscopy techniques to examine the excited state loss channels inside the devices. We use the high boiling point solvent additive 1,8-diiodooctane (DIO) to study the photoexcited state losses in different blend morphologies. The solvent additive acts as a compatibiliser between the donor and the acceptor material and leads to smaller domain sizes, higher charge formation yields and increased device efficiency.

Solving Axisymmetric Stability Problems by Using Upper Bound Finite Elements, Limit Analysis, and Linear Optimization

A numerical formulation has been proposed for solving an axisymmetric stability problem in geomechanics with upper bound limit analysis, finite elements, and linear optimization. The Drucker-Prager yield criterion is linearized by simulating a sphere with a circumscribed truncated icosahedron. The analysis considers only the velocities and plastic multiplier rates, not the stresses, as the basic unknowns. The formulation is simple to implement, and it has been employed for finding the collapse loads of a circular footing placed over the surface of a cohesive-frictional material. The formulation can be used to solve any general axisymmetric geomechanics stability problem.

Spectral finite element based on an efficient layerwise theory for wave propagation analysis of composite and sandwich beams

In this paper, we present a spectral finite element model (SFEM) using an efficient and accurate layerwise (zigzag) theory, which is applicable for wave propagation analysis of highly inhomogeneous laminated composite and sandwich beams. The theory assumes a layerwise linear variation superimposed with a global third-order variation across the thickness for the axial displacement. The conditions of zero transverse shear stress at the top and bottom and its continuity at the layer interfaces are subsequently enforced to make the number of primary unknowns independent of the number of layers, thereby making the theory as efficient as the first-order shear deformation theory (FSDT). The spectral element developed is validated by comparing the present results with those available in the literature. A comparison of the natural frequencies of simply supported composite and sandwich beams obtained by the present spectral element with the exact two-dimensional elasticity and FSDT solutions reveals that the FSDT yields highly inaccurate results for the inhomogeneous sandwich beams and thick composite beams, whereas the present element based on the zigzag theory agrees very well with the exact elasticity solution for both thick and thin, composite and sandwich beams. A significant deviation in the dispersion relations obtained using the accurate zigzag theory and the FSDT is also observed for composite beams at high frequencies. It is shown that the pure shear rotation mode remains always evanescent, contrary to what has been reported earlier. The SFEM is subsequently used to study wavenumber dispersion, free vibration and wave propagation time history in soft-core sandwich beams with composite faces for the first time in the literature. (C) 2014 Elsevier Ltd. All rights reserved.

Globally coupled stochastic two-state oscillators: Fluctuations due to finite numbers

Infinite arrays of coupled two-state stochastic oscillators exhibit well-defined steady states. We study the fluctuations that occur when the number N of oscillators in the array is finite. We choose a particular form of global coupling that in the infinite array leads to a pitchfork bifurcation from a monostable to a bistable steady state, the latter with two equally probable stationary states. The control parameter for this bifurcation is the coupling strength. In finite arrays these states become metastable: The fluctuations lead to distributions around the most probable states, with one maximum in the monostable regime and two maxima in the bistable regime. In the latter regime, the fluctuations lead to transitions between the two peak regions of the distribution. Also, we find that the fluctuations break the symmetry in the bimodal regime, that is, one metastable state becomes more probable than the other, increasingly so with increasing array size. To arrive at these results, we start from microscopic dynamical evolution equations from which we derive a Langevin equation that exhibits an interesting multiplicative noise structure. We also present a master equation description of the dynamics. Both of these equations lead to the same Fokker-Planck equation, the master equation via a 1/N expansion and the Langevin equation via standard methods of Ito calculus for multiplicative noise. From the Fokker-Planck equation we obtain an effective potential that reflects the transition from the monomodal to the bimodal distribution as a function of a control parameter. We present a variety of numerical and analytic results that illustrate the strong effects of the fluctuations. We also show that the limits N -> infinity and t -> infinity(t is the time) do not commute. In fact, the two orders of implementation lead to drastically different results.

Design and synthesis of new low band gap organic semiconductors for photovoltaic applications

Donor-acceptor (D-A) conjugated polymers have attracted a good deal of attention in recent years. In D-A systems, the introduction of electron withdrawing groups reduces E-g by lowering the LUMO levels whereas, the introduction of electron donating groups reduces E-g by raising the HOMO levels. Also, conjugated polymers with desired HOMO and LUMO energy levels could be obtained by the proper selection of donor and acceptor units. Because of this reason, D-A conjugated polymers are emerging as promising materials particularly for polymer light emitting diodes (PLEDs) and polymer solar cells (PSCs). We report the design and synthesis of four new narrow band gap donor-acceptor (D-A) conjugated polymers, PTCNN, PTCNF, PTCNV and PTCNO, containing electron donating 3,4-didodecyloxythiophene and electron accepting cyanovinylene units. The effects of further addition of electron donating and electron withdrawing groups to the repeating unit of a D-A conjugated polymer (PTCNN) on its optical and electrochemical properties are discussed. The studies revealed that the nature of D and A units as well as the extent of alternate D-A structure influences the optical and the electrochemical properties of the polymers. All the polymers are thermally stable up to a temperature of 300 degrees C under nitrogen atmosphere. The electrochemical studies revealed that the polymers possess low-lying HOMO energy levels and low-lying LUMO energy levels. In the UV-Vis absorption study, the polymer films displayed broad absorption in the wavelength region of 400-700 nm. The polymers exhibited low optical band gaps in the range 1.70 - 1.77 eV.

An extended finite element model coupled with level set method for stable pitting corrosion

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution), which requires only three inputs, namely the solid metal concentration, saturation concentration of the dissolved metal ions and diffusion coefficient. A combined eXtended Finite Element Model (XFEM) and level set method is developed in this paper. The extended finite element model handles the jump discontinuity in the metal concentrations at the interface, by using discontinuous-derivative enrichment formulation for concentration discontinuity at the interface. This eliminates the requirement of using front conforming mesh and re-meshing after each time step as in conventional finite element method. A numerical technique known as level set method tracks the position of the moving interface and updates it over time. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed method is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions.