Abstract: The effect of different stacking order of graphene multilayers on the electric field induced band gap is investigated. We considered a positively charged top and a negatively charged back gate in order to independently tune the band gap and the Fermi energy of three and four layer graphene systems. A tight-binding approach within a self-consistent Hartree approximation is used to calculate the induced charges on the different graphene layers. We found that the gap for trilayer graphene with the ABC stacking is much larger than the corresponding gap for the ABA trilayer. Also we predict that for four layers of graphene the energy gap strongly depends on the choice of stacking, and we found that the gap for the different types of stacking is much larger as compared to the case of Bernal stacking. Trigonal warping changes the size of the induced electronic gap by approximately 30% for intermediate and large values of the induced electron density.

Abstract: During the synthesis of ultrathin materials with hexagonal lattice structure Stone-Wales (SW) type of defects are quite likely to be formed and the existence of such topological defects in the graphenelike structures results in dramatic changes of their electronic and mechanical properties. Here we investigate the formation and reactivity of such SW defects in silicene. We report the energy barrier for the formation of SW defects in freestanding (similar to 2.4 eV) and Ag(111)-supported (similar to 2.8 eV) silicene and found it to be significantly lower than in graphene (similar to 9.2 eV). Moreover, the buckled nature of silicene provides a large energy barrier for the healing of the SW defect and therefore defective silicene is stable even at high temperatures. Silicene with SW defects is semiconducting with a direct band gap of 0.02 eV and this value depends on the concentration of defects. Furthermore, nitrogen substitution in SW-defected silicene shows that the defect lattice sites are the least preferable substitution locations for the N atoms. Our findings show the easy formation of SW defects in silicene and also provide a guideline for band gap engineering in silicene-based materials through such defects.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: We present a tight-binding investigation of strained bilayer graphene within linear elasticity theory, focusing on the different environments experienced by the A and B carbon atoms of the different sublattices. We find that the inequivalence of the A and B atoms is enhanced by the application of perpendicular strain epsilon(zz), which provides a physical mechanism for opening a band gap, most effectively obtained when pulling the two graphene layers apart. In addition, perpendicular strain introduces electron-hole asymmetry and can result in linear electronic dispersion near the K point. Our findings suggest experimental means for strain-engineered band gaps in bilayer graphene.

Abstract: In the limit of strong electron-phonon coupling, we use either a Pekar-type or an oscillator wave function for the center-of-mass coordinate and either a Coulomb or an oscillator wave function for the relative coordinate, and are able to reproduce all the results from the literature for the large-bipolaron binding energy. Lower bounds are constructed for the critical ratio eta(c) of dielectric constants below which bipolarons can exist. It is found that, in the strong-coupling limit, the stability region for bipolaron formation is much larger in two dimensions (2D) than in 3D. We introduce a model that combines the averaging of the relative coordinate over the asymptotically best wave function with a path-integral treatment of the center-of-mass motion. The stability region for bipolaron formation is increased compared with the full path-integral treatment at large values of the coupling constant alpha. The critical values are alpha(c) almost-equal-to 9.3 in 3D and alpha(c) almost-equal-to 4.5 in 2D. Phase diagrams for the presented models are also obtained in both 2D and 3D.

Abstract: We study the influence of the inelastic relaxation time τ̃E of the quasiparticle distribution function f(E) on the phase slip process in quasi-one-dimensional superconducting rings at a temperature close to the critical temperature Tc. We find that the initial time of growth of the order parameter |Δ| in the phase slip core after the phase slip is a nonmonotonic function of τ̃E which has a maximum at τ̃E≃τ̃GL=πℏ/8kB(Tc−T) and has a tendency to saturate for large τ̃E⪢τ̃GL. The effective heating of the electron subsystem due to the increase in |Δ| in the phase slip center together with the above effect result in a nonmonotonic dependence of the number of subsequent phase slips on τ̃E in rings of relatively large radius (in which each phase slip reduces the current density to a small fraction of its initial value). During the phase slip process the order parameter distribution has two peaks near the phase slip core due to the diffusion of the nonequilibrium quasiparticles from that region.

Abstract: (La2/3Ca1/3MnO3)n/(BaTiO3)m (LCMOn/BTOm) superlattices on MgO and SrTiO3 substrates with different layer thicknesses (n = 10, 38, 40 and m = 5, 18, 20) have been grown by metal organic aerosol deposition (MAD) and have been fully characterized down to the atomic scale to study the interface characteristics. Scanning transmission electron microscopy combined with spatially resolved electron energy-loss spectroscopy provides clear evidence for the existence of atomically sharp interfaces in MAD grown films, which exhibit epitaxial growth conditions, a uniform normal strain, and a fully oxidized state. Below a critical layer thickness the LCMO structure is found to change from the bulk Pnma symmetry to a pseudocubic R3̅ c symmetry. An atomically flat interface reconstruction consisting of a single Ca-rich atomic layer is observed on the compressively strained BTO on LCMO interface, which is thought to partially neutralize the total charge from the alternating polar atomic layers in LCMO as well as relieving strain at the interface. No interface reconstruction is observed at the tensile strained LCMO on BTO interface.

Abstract: Structural phase transitions in two vertically or horizontally coupled channels of strongly interacting particles are investigated. The particles are free to move in the x direction but are confined by a parabolic potential in the y direction. They interact with each other through a screened power-law potential (r(-n)e(-r/lambda)). In vertically coupled systems, the channels are stacked above each other in the direction perpendicular to the (x, y) plane, while in horizontally coupled systems both channels are aligned in the confinement direction. Using Monte Carlo (MC) simulations we obtain the ground-state configurations and the structural transitions as a function of the linear particle density and the separation between the channels. At zero temperature, the vertically coupled system exhibits a rich phase diagram with continuous and discontinuous transitions. On the other hand, the horizontally coupled system exhibits only a very limited number of phase transitions due to its symmetry. Further, we calculated the normal modes for the Wigner crystals in both cases. From MC simulations, we found that in the case of vertically coupled systems, the zigzag transition is only possible for low densities. A Ginzburg-Landau theory for the zigzag transition is presented, which predicts correctly the behavior of this transition from which we interpret the structural phase transition of the Wigner crystal through the reduction of the Brillouin zone.