Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

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: (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.