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: Motivated by the recent successful synthesis of double-M carbides, we investigate structural and electronic properties of WCrC and WCrCO2 monolayers and the effects of biaxial and out-of-plane strain and electric field using density functional theory. WCrC and WCrCO2 monolayers are found to be dynamically stable. WCrC is metallic and WCrCO2 display semi-metallic character with narrow band gap, which can be controlled by strain engineering and electric field. WCrCO2 monolayer exhibits a dual band gap which is preserved in the presence of an electric field. The band gap of WCrCO2 monolayer increases under uniaxial strain while it becomes metallic under tensile strain, resulting in an exotic 2D double semi-metallic behavior. Our results demonstrate that WCrCO2 is a new platform for the study of novel physical properties in two-dimensional Dirac materials and which may provide new opportunities to realize high-speed low-dissipation devices.

Abstract: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 1000 T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier, and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudomagnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation changes the sixfold symmetry of the induced pseudomagnetic field.

Abstract: The effects of strain, induced by a Gaussian bump, on the magnetic field dependent transport properties of a graphene Hall bar are investigated. The numerical simulations are performed using both classical and quantum mechanical transport theory and we found that both approaches exhibit similar characteristic features. The effects of the Gaussian bump are manifested by a decrease of the bend resistance, RB, around zero-magnetic field and the occurrence of side-peaks in RB. These features are explained as a consequence of bump-assisted scattering of electrons towards different terminals of the Hall bar. Using these features we are able to give an estimate of the size of the bump. Additional oscillations in RB are found in the quantum description that are due to the population/depopulation of Landau levels. The bump has a minor influence on the Hall resistance even for very high values of the pseudo-magnetic field. When the bump is placed outside the center of the Hall bar valley polarized electrons can be collected in the leads.

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: Recent experimental studies showed that formation of a two-dimensional crystal structure of nitrogenated holey graphene (NHG) is possible. Similar to graphene, NHGs have an atomically thin and strong crystal structure. Using first-principles calculations, we investigate the structural, phononic, and thermal properties of monolayer NHG crystal. Our charge analysis reveals that the charged holey sites of NHG provide a reactive ground for further functionalization by adatoms or molecules. We also found that similar to graphene, the NHG structure has quite high-frequency phonon modes and the presence of nitrogen atoms leads to the emergence of additional vibrational modes. Our phonon analysis reveals the presence of three characteristic Raman-active modes of NHG. Furthermore, the analysis of constant-volume heat capacity showed that the NHG structure has a linear temperature dependence in the low-temperature region. The strong lattice structure and unique thermal properties of the NHG crystal structure are desirable in nanoscale device applications.

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.