Abstract: Using first-principles calculations, we evaluate the electrochemical performance of heterostructures made up of Ti2CO2 and chemically modified graphene for Li batteries. We find that heteroatom doping and molecule intercalation have a significant impact on the storage capacity and Li migration barrier energies. While N and S doping do not improve the storage capacity, B doping together with molecule interaction make it possible to intercalate two layers of Li, which stick separately to the surface of Ti2CO2 and B-doped graphene. The calculated diffusion-barrier energies (E-diff), which are between 0.3 and 0.4 eV depending on Li concentration, are quite promising for fast charge and discharge rates. Besides, the predicted E-diff as much as 2 eV for the diffusion of the Li atom from the Ti2CO2 surface to the B-doped graphene surface significantly suppresses the interlayer Li migration, which diminishes the charge and discharge rates. The calculated volume and lattice parameter changes indicate that Ti2CO2/graphene hybrid structures exhibit cyclic stability against Li loading and unloading. Consequently, first-principles calculations we perform evidently highlight the favorable effect of molecular intercalation on the capacity improvement of ion batteries.

Abstract: The wide measured range of thermal conductivities (k) for monolayer MoS2 and the corresponding incongruent calculated values in the literature all suggest that extrinsic defect thermal resistance is significant and varied in synthesized samples of this material. Here we present defect-mediated thermal transport calculations of MoS2 using interatomic forces derived from density functional theory combined with Green's function methods to describe phonon-point-defect interactions and a Peierls-Boltzmann formalism for transport. Conductivity calculations for bulk and monolayer MoS2 using different density functional formalisms are compared. Nonperturbative first-principles methods are used to describe defect-mediated spectral functions, scattering rates, and phonon k, particularly from sulfur vacancies (VS), and in the context of the plethora of measured and calculated literature values. We find that k of monolayer MoS2 is sensitive to phonon-VS scattering in the range of experimentally observed densities, and that first-principles k calculations using these densities can explain the range of measured values found in the literature. Furthermore, measured k values for bulk MoS2 are more consistent because VS defects are not as prevalent.

Abstract: The exciton states in a strained (In,Ga)As/GaAs nanocup are theoretically determined. We explore how the nanocup bottom thickness (t) affects the magnetic field dependence of the exciton energy. Strain distribution is computed by the continuum mechanical model under the approximation of isotropic elasticity. The exciton wave functions are expanded into products of the electron and hole envelope functions. For small t, the exciton ground state has zero orbital momentum and exhibits small oscillations of the second derivative when the magnetic field increases. When t approaches the value of the cup height, however, the exciton levels exhibit angular momentum transitions, whose behavior is similar to that for type-II quantum dots. Small oscillations of the oscillator strength for exciton recombination are found when the magnetic field increases. An increase in thickness of the nanocup bottom has only a small effect on those oscillations for the optically active exciton states, but the exciton ground state becomes dark when the magnetic field increases. Hence, the results of our calculations show that an increase in thickness of the nanocup bottom transforms the exciton ground energy level dependence on magnetic field from the one characteristic of type-I rings to the one characteristic of type-II dots.

Abstract: Large-scale atomistic simulations using the reactive empirical bond order force field approach is implemented to investigate thermal and mechanical properties of single-layer (SL) and multilayer (ML) molybdenum disulfide (MoS2). The amplitude of the intrinsic ripples of SL MoS2 are found to be smaller than those exhibited by graphene (GE). Furthermore, because of the van der Waals interaction between layers, the out-of-plane thermal fluctuations of ML MoS2 decreases rapidly with increasing number of layers. This trend is confirmed by the buckling transition due to uniaxial stress which occurs for a significantly larger applied tension as compared to graphene. For SL MoS2, the melting temperature is estimated to be 3700 K which occurs through dimerization followed by the formation of small molecules consisting of two to five atoms. When different types of vacancies are inserted in the SL MoS2 it results in a decrease of both the melting temperature as well as the stiffness.

Abstract: The formation of lattice defects and their effect on the electronic properties of single-layer FeCl2 are investigated by means of first-principles calculations. Among the vacancy defects, namely mono-, di-, and three-Cl vacancies and mono-Fe vacancy, the formation of mono-Cl vacancy is the most preferable. Comparison of two different antisite defects reveals that the formation of the Fe-antisite defect is energetically preferable to the Cl-antisite defect. While a single Cl vacancy leads to a 1 mu(B) decrease in the total magnetic moment of the host lattice, each Fe vacant site reduces the magnetic moment by 4 mu(B). However, adsorption of an excess Cl atom on the surface changes the electronic structure to a ferromagnetic metal or to a ferromagnetic semiconductor depending on the adsorption site without changing the ferromagnetic state of the host lattice. Both Cl-antisite and Fe-antisite defected domains change the magnetic moment of the host lattice by -1 mu(B) and +3 mu(B), respectively. The electronic ground state of defected structures reveals that (i) single-layer FeCl2 exhibits half-metallicity under the formation of vacancy and Cl-antisite defects; (ii) ferromagnetic metallicity is obtained when a single Cl atom is adsorbed on upper-Cl and Fe sites, respectively; and (iii) ferromagnetic semiconducting behavior is found when a Cl atom is adsorbed on a lower-Cl site or a Fe-antisite defect is formed. Simulated scanning electron microscope images show that atomic-scale identification of defect types is possible from their electronic charge density. Further investigation of the periodically Fe-defected structures reveals that the formation of the single-layer FeCl3 phase, which is a dynamically stable antiferromagnetic semiconductor, is possible. Our comprehensive analysis on defects in single-layer FeCl2 will complement forthcoming experimental observations.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: We calculate the Curie temperature of layered ferromagnets, chromium tri-iodide (CrI3), chromium tri-bromide (CrBr3), chromium germanium tri-telluride (CrGeTe3), and the Ned temperature of a layered antiferromagnet iron di-chloride (FeCl2), using first-principles density functional theory calculations and Monte Carlo simulations. We develop a computational method to model the magnetic interactions in layered magnetic materials and calculate their critical temperature. We provide a unified method to obtain the magnetic exchange parameters (J) for an effective Heisenberg Hamiltonian from first principles, taking into account both the magnetic ansiotropy as well as the out-of-plane interactions. We obtain the magnetic phase change behavior, in particular the critical temperature, from the susceptibility and the specific-heat, calculated using the three-dimensional Monte Carlo (METROPOLIS) algorithm. The calculated Curie temperatures for ferromagnetic materials (CrI3, CrBr3, and CrGeTe3), match well with experimental values. We show that the interlayer interaction in bulk CrI3 with R (3) over bar stacking is significantly stronger than the C2/m stacking, in line with experimental observations. We show that the strong interlayer interaction in R (3) over bar CrI3 results in a competition between the in-plane and the out-of-plane magnetic easy axes. Finally, we calculate the Ned temperature of FeCl2 to be 47 +/- 8 K and show that the magnetic phase transition in FeCl2 occurs in two steps with a high-temperature intralayer ferromagnetic phase transition and a low-temperature interlayer antiferromagnetic phase transition.

Abstract: The recent development in the field of two-dimensional magnetic materials urges reliable theoretical methodology for determination of magnetic properties. Among the available methods, ab initio four-state energy mapping based on density functional theory stands out as a powerful technique to calculate the magnetic exchange interaction in the Heisenberg spin model. Although the required formulas were explained in earlier works, the considered Hamiltonian in those studies always corresponded to the specific case that the off-diagonal part of J matrix is antisymmetric, which may be misleading in other cases. Therefore, using the most general form of the Heisenberg spin Hamiltonian, we here derive the generic formulas. With a proper choice of four different magnetic states, a single formula governs all elements of the exchange interaction matrix for any considered pair of spin sites.

Abstract: Based on the Gor'kov formalism for a clean s-wave superconductor, we develop an extended version of the single-band Ginzburg-Landau (GL) theory by means of a systematic expansion in the deviation from the critical temperature T(c), i.e., tau = 1 – T/T(c). We calculate different contributions to the order parameter and the magnetic field: the leading contributions (proportional to tau(1/2) in the order parameter and. t in the magnetic field) are controlled by the standard GL theory, while the next-to-leading terms (proportional to tau(3/2) in the gap and proportional to tau(2) in the magnetic field) constitute the extended GL (EGL) approach. We derive the free-energy functional for the extended formalism and the corresponding expression for the current density. To illustrate the usefulness of our formalism, we calculate, in a semianalytical form, the temperature-dependent correction to the GL parameter at which the surface energy becomes zero, and analytically, the temperature dependence of the thermodynamic critical field. We demonstrate that the EGL formalism is not just a mathematical extension to the theory: variations of both the gap and the thermodynamic critical field with temperature calculated within the EGL theory are found in very good agreement with the full BCS results down to low temperatures, which dramatically improves the applicability of the formalism compared to its standard predecessor.

Abstract: Graphane, hydrogenated graphene, was very recently synthesized and predicted to have great potential applications. In this work, we propose a new promising approach for hydrogenation of graphene based on density functional theory (DFT) calculations through the application of a perpendicular electric field after substitutionally doping by nitrogen atoms. These DFT calculations show that the doping by nitrogen atoms into the graphene layer and applying an electrical field normal to the graphene surface induce dissociative adsorption of hydrogen. The dissociative adsorption energy barrier of an H2 molecule on a pristine graphene layer changes from 2.7 to 2.5 eV on N-doped graphene, and to 0.88 eV on N-doped graphene under an electric field of 0.005 au. When increasing the electric field above 0.01 au, the reaction barrier disappears. Therefore, N doping and applying an electric field have catalytic effects on the hydrogenation of graphene, which can be used for hydrogen storage purposes and nanoelectronic applications.

Abstract: We study theoretically the simultaneous effect of regular and random pinning potentials on the vortex lattice structure at filling factor of 1. This structure is determined by a competition between the square symmetry of regular pinning array, by the intervortex interaction favoring a triangular symmetry, and by the randomness trying to depin vortices from their regular positions. Both analytical and molecular-dynamics approaches are used. We construct a phase diagram of the system in the plane of regular and random pinning strengths and determine typical vortex lattice defects appearing in the system due to the disorder. We find that the total disordering of the vortex lattice can occur either in one step or in two steps. For instance, in the limit of weak pinning, a square lattice of pinned vortices is destroyed in two steps. First, elastic chains of depinned vortices appear in the film; but the vortex lattice as a whole remains still pinned by the underlying square array of regular pinning sites. These chains are composed into fractal-like structures. In a second step, domains of totally depinned vortices are generated and the vortex lattice depins from regular array.

Abstract: We theoretically examine behavior of superconductivity at parallel interfaces separating the domains of another dominant collective excitation, such as charge density waves or spin density waves. Due to their competitive coupling in a two-component Ginzburg-Landau model, suppression of the dominant order parameter at the interfacial planes allows for nucleation of the (hidden) superconducting order parameter at those planes. In such a case, we demonstrate how the number of the parallel interfacial planes and the distance between them are linked to the number and the size of the emerging superconducting gaps in the system, as well as the versatility and temperature evolution of the possible superconducting phases. These findings bear relevance to a broad selection of known layered superconducting materials, as well as to further design of artificial (e.g., oxide) superlattices, where the interplay between competing order parameters paves the way towards otherwise unattainable superconducting states, some with enhanced superconducting critical temperature.