Abstract: Using first-principle calculations, we systematically study strain tuning of the electronic properties of two- dimensional graphitic carbon nitride nanosheets with empirical formula CnNm. We found the following: (i) the ferromagnetic ordered state in the metal-free systems (n, m) = (4,3), (10,9), and (14,12) remains stable in the presence of strain of about 6%. However, the system (9,7) loses its ferromagnetic ordering when increasing strain. This is due to the presence of topological defects in the (9,7) system, which eliminates the asymmetry between spin up and spin down of the p(z) orbitals when strain is applied. (ii) By applying uniaxial strain, a band gap opens in systems which are initially gapless. (iii) In semiconducting systems which have an initial gap of about 1 eV, the band gap is closed with applying uniaxial strain.

Abstract: The optoelectronic properties of confined water form one of the most active research areas in the past few years. Here we present the multiscale methodology to discern the out-of-plane electronic and dipolar dielectric constants (epsilon(el)(perpendicular to) and epsilon(diP)(perpendicular to)) of strongly confined water. We reveal that epsilon(perpendicular to el) and epsilon(diP)(perpendicular to) become comparable for water confined in angstrom-scale channels (with a height of less than 15 angstrom) within graphene (GE) and hexagonal boron nitride (hBN) bilayers. Channel height (h) associated with a minimum in both epsilon(e)(l)(perpendicular to) and epsilon(dip)(perpendicular to) is linked to the formation of the ordered structure of ice for h approximate to (7 -7.5) angstrom. The recently measured total dielectric constant epsilon(T)(perpendicular to) of nanoconfined water [L. Fumagalli et al., Science 360, 1339 (2018)] is corroborated by our results. Furthermore, we evaluate the contribution from the encapsulating membranes to the dielectric properties, as a function of the interlayer spacing, i.e., the height of the confining channel for water. Finally, we conduct analysis of the optical properties of both confined water and GE membranes, and show that the electron energy loss function of confined water strongly differs from that of bulk water.

Abstract: The surface termination of MXenes greatly determines the electrochemical properties and ion kinetics on their surfaces. So far, hydroxyl-, oxygen-, and fluorine-terminated MXenes have been widely studied for energy storage applications. Recently, sulfur-functionalized MXene structures, which possess low diffusion barriers, have been proposed as candidate materials to enhance battery performance. We performed first-principles calculations on the structural, stability, electrochemical, and ion dynamic properties of Li-adsorbed sulfur-functionalized groups 3B, 4B, 5B, and 6B transition-metal (M)-based MXenes (i.e., M2CS2 with M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, and W). We performed phonon calculations, which indicated that all of the above M2CS2 MXenes, except for Sc, are dynamically stable at T = 0 K. The ground-state structure of each M2CS2 monolayer depends on the type of M atom. For instance, while sulfur prefers to sit at the FCC site on Ti2CS2, it occupies the HCP site of Cr-based MXene. We determined the Li adsorption configurations at different concentrations using the cluster expansion method. The highest maximum open-circuit voltages were computed for the group 4B element (i.e., Ti, Zr, and Hf)-based M2CS2, which are larger than 2.1 V, while their average voltages are approximately 1 V. The maximum voltage for the group 6B element (i.e., Cr, Mo, W)-based M2CS2 is less than 1 V, and the average voltage is less than 0.71 V. We found that S functionalization is helpful for capacity improvements over the O-terminated MXenes. In this respect, the computed storage gravimetric capacity may reach up to 417.4 mAh/g for Ti2CS2 and 404.5 mAh/g for V2CS2. Ta-, Cr-, Mo-, and W-based M2CS2 MXenes show very low capacities, which are less than 100 mAh/g. The Li surface diffusion energy barriers for all of the considered MXenes are less than 0.22 eV, which is favorable for high charging and discharging rates. Finally, ab initio molecular dynamic simulations performed at 400 K and bond-length analysis with respect to Li concentration verify that selected promising systems are robust against thermally induced perturbations that may induce structural transformations or distortions and undesirable Li release.

Abstract: By using the full 3D generalized time-dependent Ginzbug-Landau equation, we study a long superconducting film of finite width and thickness under an applied transport current. We show that, for sufficiently large thickness, the vortices and the antivortices become curved before they annihilate each other. As they approach the center of the sample, their ends combine, producing a single closed vortex. We also determine the critical values of the thickness for which the closed vortex sets in for different values of the Ginzburg-Ladau parameter. Finally, we propose a model of how to detect a closed vortex experimentally.

Abstract: Flat bands in the energy spectrum have attracted a lot of attention in recent years because of their unique properties and promising applications. Special arrangement of impurities on monolayer graphene are proposed to generate multiflat bands in the electronic band structure. In addition to the single midgap states in the spectrum of graphene with low hydrogen density, we found closely spaced bands around the Fermi level with increasing impurity density, which are similar to discrete lines in the spectrum of quantum dots, as well as the unusual Landau-level energy spectrum of graphene in the presence of a strong magnetic field. The presence of flat bands crucially depends on whether or not there are odd or even electrons of H(F) atoms bound to graphene. Interestingly, we found that a fully hydrogenated (fluoridated) of a hexagon of graphene sheet with six hydrogen (fluorine) atoms sitting on top and bottom in consecutive order exhibits Dirac cones in the electronic band structure with a 20% smaller Fermi velocity as compared to the pristine graphene. Functionalizing graphene introduces various C-C bond lengths resulting in nonuniform strains. Such a nonuniform strain may induce a giant pseudomagnetic field in the system, resulting in quantum Hall effect.

Abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

Abstract: The superconducting properties of two recently synthesized ternary silicides with unit formula X₄CoSi (X = Nb, Ta) are investigated through ab initio calculations combined with Eliashberg theory. Interestingly, their crystal structure comprises interlocking honeycomb networks of Nb/Ta atoms. Nb₄CoSi is found to harbor better conditions for phonon-mediated superconductivity, as it possesses a higher density of states at the Fermi level, fostering stronger electron-phonon coupling. The superconducting critical temperatures (T_c) follow the same trend, with Nb₄CoSi having a twice higher value than Ta₄CoSi. Furthermore, the calculated T_c values (5.9 K vs 3.1 K) agree excellently with the experimentally obtained ones, establishing superconductivity in this new materials class as mediated by the electron-phonon coupling. Furthermore, my calculations show that the superconducting properties of these compounds do not simply correlate with the parameters of their honeycomb networks, contrary to proposals raised in the literature. Rather, their complete fermiology and phonon spectrum should be taken into account in order to explain their respective superconducting properties.

Abstract: We theoretically investigate the binding energy and electron-hole (e-h) overlap of excitonic states confined at the interface between two-dimensional materials with type-II band alignment, i.e., with lowest conduction and highest valence band edges placed in different materials, arranged in a side-by-side planar heterostructure. We propose a variational procedure within the effective mass approximation to calculate the exciton ground state and apply our model to a monolayer MoS2/WS2 heterostructure. The role of nonabrupt interfaces between the materials is accounted for in our model by assuming a WxMo1-xS2 alloy around the interfacial region. Our results demonstrate that (i) interface-bound excitons are energetically favorable only for small interface thickness and/or for systems under high dielectric screening by the materials surrounding the monolayer, and that (ii) the interface exciton binding energy and its e-h overlap are controllable by the interface width and dielectric environment.

Abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: Employing density functional theory-based methods, we investigate monolayer and bilayer structures of hexagonal SnS2, which is a recently synthesized monolayer metal dichalcogenide. Comparison of the 1H and 1T phases of monolayer SnS2 confirms the ground state to be the 1T phase. In its bilayer structure we examine different stacking configurations of the two layers. It is found that the interlayer coupling in bilayer SnS2 is weaker than that of typical transition-metal dichalcogenides so that alternative stacking orders have similar structural parameters and they are separated with low energy barriers. A possible signature of the stacking order in the SnS2 bilayer has been sought in the calculated absorbance and reflectivity spectra. We also study the effects of the external electric field, charging, and loading pressure on the characteristic properties of bilayer SnS2. It is found that (i) the electric field increases the coupling between the layers at its preferred stacking order, so the barrier height increases, (ii) the bang gap value can be tuned by the external E field and under sufficient E field, the bilayer SnS2 can become a semimetal, (iii) the most favorable stacking order can be switched by charging, and (iv) a loading pressure exceeding 3 GPa changes the stacking order. The E-field tunable band gap and easily tunable stacking sequence of SnS2 layers make this 2D crystal structure a good candidate for field effect transistor and nanoscale lubricant applications.

Abstract: Mixed magnetism (the coexistence of strong and weak magnetism in one material) is regarded as the origin of the giant magnetocaloric effect (GMCE). A good example is (Mn,Fe)(2)(Si,P), which is established as one of the best magnetocaloric materials available. Tuning the material properties are essential for optimizing its performance, and a straightforward way to do that is by doping. In this article, an ab initio electronic structure method was used to calculate the structure and magnetic properties of 3d-transition-metal-doped (Mn,Fe)(2)(Si,P) materials for magnetocaloric applications (transition metals are Cr, Co, Mn, Ni, Cu). For a steady performance, the material should be mechanically stable. A detailed analysis of the elastic constants shows that the mechanical stability of the (Mn,Fe)(2)(Si,P) system increases significantly by doping with boron without affecting the magnetic properties. Insights of the influence of doping enable future studies to understand and predict bettermagnetocaloric materials.

Abstract: We reveal the emergence of metallicKondo clouds in an atomistic carbon quantum dot, realized as a single-atom junction in a suitably patterned graphene nanoflake. Using density functional dynamical mean-field theory (DFDMFT) we show how correlation effects lead to striking features in the electronic structure of our device, and how those are enhanced by the electron-electron interactions when graphene is patterned at the atomistic scale. Our setup provides a well-controlled environment to understand the principles behind the orbital-selective Kondo physics and the interplay between orbital and spin degrees of freedom in carbon-based nanomaterials, which indicate new pathways for spintronics in atomically patterned graphene.

Abstract: Using the tight-binding approach, we investigate confined states in two different hybrid monolayer-bilayer systems: (i) a hexagonal monolayer area surrounded by bilayer graphene in the presence of a perpendicularly applied electric field and (ii) a hexagonal bilayer graphene dot surrounded by monolayer graphene. The dependence of the energy levels on dot size and external magnetic field is calculated. We find that the energy spectrum for quantum dots with zigzag edges consists of states inside the gap which range from dot-localized states, edge states, to mixed states coexisting together, whereas for dots with armchair edges, only dot-localized states are observed.

Abstract: The application of graphene to electronic and optoelectronic devices is limited by the absence of reliable semiconducting variants of this material. A promising candidate in this respect is graphene oxide, with a band gap on the order of similar to 5 eV, however, this has a finite density of states at the Fermi level. Here, we examine the electronic structure of three variants of half-fluorinated carbon on Sic(0001), i.e., the (6 root 3 x 6 root 3) R30 degrees C/SiC “buffer layer,” graphene on this (6 root 3 x 6 root 3) R30 degrees C/SiC buffer layer, and graphene decoupled from the SiC substrate by hydrogen intercalation. Using angle-resolved photoemission, core level photoemission, and x-ray absorption, we show that the electronic, chemical, and physical structure of all three variants is remarkably similar, exhibiting a large band gap and a vanishing density of states at the Fermi level. These results are explained in terms of first-principles calculations. This material thus appears very suitable for applications, even more so since it is prepared on a processing-friendly substrate. We also investigate two separate UV photon-induced modifications of the electronic structure that transform the insulating samples (6.2-eV band gap) into semiconducting (similar to 2.5-eV band gap) and metallic regions, respectively.