Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: Using first principles density functional theory in combination with the nonequilibrium Green's function formalism, we study the effect of derivatization on the electronic and transport properties of C-60 fullerene. As a typical example, we consider [6,6]-phenyl-C-61-butyric acid methyl ester (PCBM), which forms one of the most efficient organic photovoltaic materials in combination with electron donating polymers. Extra peaks are observed in the density of states (DOS) due to the formation of new electronic states localized at/near the attached molecule. Despite such peculiar behavior in the DOS of an isolated molecule, derivatization does not have a pronounced effect on the electronic transport properties of the fullerene molecular junctions. Both C-60 and PCBM show the same response to finite voltage biasing with new features in the transmission spectrum due to voltage induced delocalization of some electronic states. We also study the diffusive motion of molecular fullerenes in ethanol solvent and inside poly(3-hexylthiophene) lamella using reactive molecular dynamics simulations. We found that the mobility of the fullerene reduces considerably due to derivatization; the diffusion coefficient of C-60 is an order of magnitude larger than the one for PCBM. (c) 2015 AIP Publishing LLC.

Abstract: Using first principles calculations we showed that the electronic and optical properties of single-layer black phosphorus (BP) depend strongly on the applied strain. Due to the strong anisotropic atomic structure of BP, its electronic conductivity and optical response are sensitive to the magnitude and the orientation of the applied strain. We found that the inclusion of many body effects is essential for the correct description of the electronic properties of monolayer BP; for example, while the electronic gap of strainless BP is found to be 0.90 eV by using semilocal functionals, it becomes 2.31 eV when many-body effects are taken into account within the G(0)W(0) scheme. Applied tensile strain was shown to significantly enhance electron transport along zigzag direction of BP. Furthermore, biaxial strain is able to tune the optical band gap of monolayer BP from 0.38 eV (at -8% strain) to 2.07 eV (at 5.5%). The exciton binding energy is also sensitive to the magnitude of the applied strain. It is found to be 0.40 eV for compressive biaxial strain of -8%, and it becomes 0.83 eV for tensile strain of 4%. Our calculations demonstrate that the optical response of BP can be significantly tuned using strain engineering which appears as a promising way to design novel photovoltaic devices that capture a broad range of solar spectrum.

Abstract: Using evolutionary algorithm for crystal structure prediction, we present a new stable two-dimensional (2D) carbon allotrope composed of polymerized as-indacenes (PAI) in a zigzag pattern, namely PAI-graphene whose energy is lower than most of the reported 2D allotropes of graphene. Crucially, the crystal structure realizes a nonsymmorphic layer group that enforces a nontrivial global topology of the band structure with two Dirac cones lying perfectly at the Fermi level. The absence of electron/hole pockets makes PAI-graphene a pristine crystalline topological semimetal having anisotropic Fermi velocities with a high value of 7.0 x 10(5) m/s. We show that while the semimetallic property of the allotrope is robust against the application of strain, the positions of the Dirac cone and the Fermi velocities can be modified significantly with strain. Moreover, by combining strain along both the x- and y-directions, two band inversions take place at G leading to the annihilation of the Dirac nodes demonstrating the possibility of strain-controlled conversion of a topological semimetal into a semiconductor. Finally we formulate the bulk-boundary correspondence of the topological nodal phase in the form of a generalized Zak-phase argument finding a perfect agreement with the topological edge states computed for different edge-terminations. (C) 2020 The Author(s). Published by Elsevier Ltd.

Abstract: Using density-functional theory in combination with the nonequilibrium Green's function formalism, we study the effect of the crystal lattice structure of organometallic perovskite CH3NH3PbI3 on its electronic transport properties. Both dispersive interactions and spin-orbit coupling are taken into account in describing structural and electronic properties of the system. We consider two different phases of the material, namely the orthorhombic and cubic lattice structures, which are energetically stable at low (< 160 K) and high (> 330 K) temperatures, respectively. The sizable geometrical differences between the two structures in term of lattice parameters, PbI6 octahedral tilts, rotation and deformations, have considerable impact on the transport properties of the material. For example, at zero bias and for all considered electron energies, the cubic phase has a larger transmission than the orthorhombic one, although both show similar electronic densities of states. Depending on the applied voltage, the current in the cubic system can be several orders of magnitude larger as compared to the one obtained for the orthorhombic sample. We attribute this enhancement in the transmission to the presence of extended states in the cubic phase due to the symmetrically shaped and ordered PbI6 octaherdra. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: Using density functional theory-based first-principles calculations, we investigate the stacking order dependence of the electronic and vibrational properties of HfS2-Ca(OH)(2) heterobilayer structures. It is shown that while the different stacking types exhibit similar electronic and optical properties, they are distinguishable from each other in terms of their vibrational properties. Our findings on the vibrational properties are the following: (i) from the interlayer shear (SM) and layer breathing (LBM) modes we are able to deduce the AB' stacking order, (ii) in addition, the AB' stacking type can also be identified via the phonon softening of E-g(I) and A(g)(III) modes which harden in the other two stacking types, and (iii) importantly, the ultrahigh frequency regime possesses distinctive properties from which we can distinguish between all stacking types. Moreover, the differences in optical and vibrational properties of various stacking types are driven by two physical effects, induced biaxial strain on the layers and the layer-layer interaction. Our results reveal that with both the phonon frequencies and corresponding activities, the Raman spectrum possesses distinctive properties for monitoring the stacking type in novel vertical heterostructures constructed by alkaline-earth-metal hydroxides.

Abstract: Using density functional theory, we obtain the mechanical and thermal properties of MX2 monolayers (where M = Cr, Mo, W and X = O, S, Se, Te). The C-centered phonon frequencies (i.e., A(1), A(2)'', E ', and E ''), relative frequency values of A(1), and E ' modes, and mechanical properties (i.e., elastic constants, Young modulus, and Poisson's ratio) display a strong dependence on the type of metal and chalcogenide atoms. In each chalcogenide (metal) group, transition-metal dichalcogenides (TMDCs) with W (O) atom are found to be much stiffer. Consistent with their stability, the thermal expansion of lattice constants for TMDCs with O (Te) is much slower (faster). Furthermore, in a heterostructure of these materials, the difference of the thermal expansion of lattice constants between the individual components becomes quite tiny over the whole temperature range. The calculated mechanical and thermal properties show that TMDCs are promising materials for heterostructures. (C) 2014 AIP Publishing LLC.

Abstract: Using density functional theory based first-principles calculations, we predict the dynamically stable 1H phase of a Janus single layer composed of S-Mo-O atomic layers. It is an indirect band gap semiconductor exhibiting strong polarization arising from the charge difference on the two surfaces. In contrast to 1H phases of MoS2 and MoO2, Janus MoSO is found to possess four Raman active phonon modes and a large out-of-plane piezoelectric coefficient which is absent in fully symmetric single layers of MoS2 and MoO2. We investigated the electronic and phononic properties under applied biaxial strain and found an electronic phase transition with tensile strain while the conduction band edge displays a shift when under compressive strain. Furthermore, single-layer MoSO exhibits phononic stability up to 5% of compressive and 11% of tensile strain with significant phonon shifts. The phonon instability is shown to arise from the soft in-plane and out-of-plane acoustic modes at finite wave vector. The large strain tolerance of Janus MoSO is important for nanoelastic applications. In view of the dynamical stability even under moderate strain, we expect that Janus MoSO can be fabricated in the common 1H phase with a strong out-of-plane piezoelectric coefficient.

Abstract: Using Brownian dynamics simulations, we investigate the dynamics of colloids confined in two-dimensional narrow channels driven by a nonuniform force Fdr(y). We considered linear-gradient, parabolic, and deltalike driving-force profiles. This driving force induces melting of the colloidal solid (i.e., shear-induced melting), and the colloidal motion experiences a transition from elastic to plastic regime with increasing Fdr. For intermediate Fdr (i.e., in the transition region) the response of the system, i.e., the distribution of the velocities of the colloidal chains õi(y), in general does not coincide with the profile of the driving force Fdr(y), and depends on the magnitude of Fdr, the width of the channel, and the density of colloids. For example, we show that the onset of plasticity is first observed near the boundaries while the motion in the central region is elastic. This is explained by: (i) (in)commensurability between the chains due to the larger density of colloids near the boundaries, and (ii) the gradient in Fdr. Our study provides a deeper understanding of the dynamics of colloids in channels and could be accessed in experiments on colloids (or in dusty plasma) with, e.g., asymmetric channels or in the presence of a gradient potential field.

Abstract: Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean-square height fluctuations < h(2)> and the height-height correlation function H(q) we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN (i) is about 16% smaller than the one of GE at room temperature (300 K), and (ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal-expansion coefficient, and Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

Abstract: Using an extension of the Gutzwiller approximation for an inhomogeneous system, we study the two-band Hubbard model with unequal band widths for a slab geometry. The aim is to investigate the mutual effect of individual bands on the spatial distribution of quasi-particle weight and charge density, especially near the surface of the slab. The main effect of the difference in band width is the presence of two different length scales corresponding to the quasi-particle profile of each band. This is enhanced in the vicinity of the critical interaction of the narrow band where an orbitally selective Mott transition occurs and a surface dead layer forms for the narrow band. For the doped case, two different regimes of charge transfer between the surface and the bulk of the slab are revealed. The charge transfer from surface/ center to center/ surface depends on both the doping level and the average relative charge accumulated in each band. Such effects could also be of importance when describing the accumulation of charges at the interface between structures made of multi-band strongly correlated materials.

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: Using advanced Ginzburg-Landau simulations, we study the superconducting state of a thin superconducting film under a ferromagnetic layer, separated by an insulating oxide, in applied external magnetic field and electric current. The taken uniaxial ferromagnet is organized into a series of parallel domains with alternating polarization of out-of-plane magnetization, sufficiently strong to induce vortex-antivortex pairs in the underlying superconductor in absence of other magnetic field. We show the organization of such vortex-antivortex matter into rich configurations, some of which are not matching the periodicity of the ferromagnetic film. The variety of possible configurations is enhanced by applied homogeneous magnetic field, where additional vortices in the superconductor may lower the energy of the system by either annihilating the present antivortices under negative ferromagnetic domains or by lowering their own energy after positioning under positive ferromagnetic domains. As a consequence, both the vortex-antivortex reordering in increasing external field and the evolution of the energy of the system are highly nontrivial. Finally, we reveal the very interesting effects of applied dc electric current on the vortex-antivortex configurations, since resulting Lorentzian force has opposite direction for vortices and antivortices, while direction of the applied current with respect to ferromagnetic domains is of crucial importance for the interaction of the applied and the Meissner current, as well as the consequent vortex-antivortex dynamics-both of which are reflected in the anisotropic critical current of the system. (C) 2016 AIP Publishing LLC.

Abstract: Using a numerical tight-binding approach based on the Chebyshev–Bogoliubov–de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and nondimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in nondimer sites is enhanced when compared to single-layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the topmost layers while for rhombohedral stacking the current flows throughout the whole structure.

Abstract: Using a complementary combination of x-ray diffraction and atomically resolved imaging we investigated the lattice structure of epitaxial LaNiO3/LaAlO3 superlattices grown on a compressive-strain inducing LaSrAlO4 (001) substrate. A refinement of the structure obtained from the x-ray data revealed the monoclinic I 2/c 1 1 space group. The (Ni/Al)O6 octahedral rotation angle perpendicular to the superlattice plane is enhanced, and the one parallel to the plane is reduced with respect to the corresponding bulk values. High-angle annular dark field imaging was used to determine the lattice parameters within the superlattice unit cell. High-resolution electron microscopy images of the oxygen atoms are consistent with the x-ray results.