Oxygen-vacancy-induced magnetism in anti-perovskite topological Dirac semimetal Ba3SnO.

The thermodynamic, structural, magnetic and electronic properties of the pristine and intrinsic vacancy-defect-containing topological Dirac semimetal Ba3SnO are studied using first-principles density functional theory calculations. The thermodynamic stability of Ba3SnO has been evaluated with reference to its competing binary phases Ba2Sn, BaSn and BaO. Subsequently, valid limits of the atomic chemical potentials derived from the thermodynamic stability were used for assessing the formation of Ba, Sn and O vacancy defects in Ba3SnO under different synthesis environments. Based on the calculated defect-formation energies, we find that the charge-neutral oxygen vacancies are the most favourable type of vacancy defect under most chemical environments. The calculated electronic properties of pristine Ba3SnO show that inclusion of spin-orbit coupling in exchange-correlation potentials computed using generalized gradient approximation yields a semimetallic band structure exhibiting twin Dirac cones along the Γ-X path of the Brillouin zone. The effect of spin-polarization and spin-orbit coupling on the physical properties of intrinsic vacancy defects containing Ba3SnO has been examined in detail. Using Bader charges, electron localization function (ELF), electronic density of states (DOS) and spin density, we show that the isolated oxygen vacancy is a magnetic defect in anti-perovskite Ba3SnO. Our results show that the origin of magnetism in Ba3SnO is the accumulation of unpaired charges at the oxygen vacancy sites, which couple strongly with the 5d states of the Ba atom. Owing to the metastability observed in earlier theoretically predicted magnetic topological semimetals, the present study reveals the important role of intrinsic vacancy defects in giving rise to magnetism and also provides opportunities for engineering the electronic structure of a Dirac semimetal.

Superior Anchoring of Sodium Polysulfides to the Polar C2N 2D Material: A Potential Electrode Enhancer in Sodium-Sulfur Batteries.

Despite the high theoretical specific energy in rechargeable sodium-sulfur batteries, the shuttle effect severely hampers its capacity and reversibility, which could be overcome by introducing an anchoring material. We, herein, use first-principles calculations to study the low-cost, easily synthesized, environmentally friendly, and stable two-dimensional polar nitrogenated holey graphene (C2N) and nonpolar polyaniline (C3N) to investigate their performance as anchoring materials and the mechanism behind the binding to identify the best candidate to improve the performance of sodium-sulfur batteries. We gain insight into the interaction, including the lowest-energy configurations, binding energies, binding nature, charge transfer, and electronic properties. Sodium primarily contributes to binding with the nanosheets, which is in accordance with their characteristics as anchoring materials. Sodium polysulfides (NaPSs) and the S8 cluster adsorb at the pores of C2N, where there are six electron lone pairs, one for each N atom. The polar C2N binds the NaPSs much strongly than the nonpolar C3N. In contrast to C3N, the charge population substantially modifies by adsorbing NaPSs on C2N, with a substantial charge transfer from the sulfur atoms. The calculated work function of 6.04 eV for pristine C2N, comparable with the previously reported values, decreases on adsorption of the NaPSs formed from battery discharging. We suggest that the inclusion of C2N into sulfur electrodes could also improve their issue with poor conductivity.

Communication: Origin of the difference between carbon nanotube armchair and zigzag ends.

In this work, we have found that the difference between armchair and zigzag ends of carbon nanotubes (CNTs) does not pertain at close study for individual bonds and thus alternative strategies need to be developed to reach the ultimate goals in selective growth. Based on first-principles simulations, the difference between binding strengths for CNTs of different chirality was investigated using hydrogen dissociation energies at their passivated ends. When all H atoms are removed collectively we find the well-known difference: that armchair bonds are much weaker than zigzag ones, which is typically seen for both CNT ends and graphene edges. However, when individual H atoms are removed we find almost no difference in hydrogen dissociation energies, small difference in bond lengths, which by association means small difference in C-C and M-C binding energies. We show convincingly that the difference in binding energy between armchair and zigzag ends is due to a fragment stabilization effect that is only manifested when all (or several neighbouring) bonds are broken. This is because at armchair ends/edges neighbouring dangling bonds can pair-up to form C≡C triple bonds that constitute a considerable stabilization effect compared to the isolated dangling bonds at zigzag ends/edges. Consequently, in many processes, e.g., catalytic growth where bonds are normally created/broken sequentially, not collectively, the difference between armchair and zigzag ends/edges cannot be used to discriminate growth of one type over the other to achieve chiral selective growth. Strategies are discussed to realize chirality selective growth in the light of the results presented, including addition of C2-fragments to favor armchair tubes.

V4C3 MXene: a Type-II Nodal Line Semimetal with Potential as High-Performing Anode Material for Mg-Ion Battery.

We have used density functional theory simulations to explore the topological characteristics of a new MXene-like material, V4C3, and its oxide counterpart, assessing their potential as anode materials for Mg-ion batteries. Our research reveals that V4C3 monolayer is a topological type-II nodal line semimetal, protected by time reversal and spatial inversion symmetries. This type-II nodal line is marked by unique drumhead-like edge states that appear either inside or outside the loop circle, contingent upon the edge ending. Intriguingly, even with an increase in metallicity due to oxygen functionalization, the topological features of V4C3 remain intact. Consequently, the monolayer V4C3 has a topologically enhanced electrical conductivity that amplifies further upon oxygen functionalization. During the charging phase, a remarkable storage concentration led to a peak specific capacity of 894.73 mAh g-1 for V4C3, which only decreases to 789.33 mAh g-1 for V4C3O2. When compared to V2C, V4C3 displays a significantly lower specific capacity loss due to functionalization, demonstrating its superior electrochemical properties. Additionally, V4C3 and V4C3O2 exhibit moderate average open-circuit voltages (0.54 V for V4C3 and 0.58 V for V4C3O2) and energy barriers for intercalation migration (ranging between 0.29-0.63 eV), which are desirable for anode materials. Thus, our simulation results support V4C3 potential as an efficient anode material for Mg-ion batteries.

Realization of either physisorption or chemisorption of 2H-tetraphenylporphyrin on the Cu(111) from density functional theory.

The adsorption of organic molecules to surfaces is a central issue to achieve fully-functional molecular devices, for which porphyrins are well-studied due to their chemical stability and functional diversity. Herein, we investigate both the physical and the chemical adsorption of the free-base tetraphenylporphyrin 2H-TPP on the Cu(111) surface within the framework of density functional theory and find that the most stable physisorbed configuration is more weakly bound by -0.31 eV than the chemisorbed configuration. We use the electron localization function to investigate the difference in binding mechanisms between strong physisorption and weak chemisorption. We have computed a reaction barrier of 0.12 eV in going from physical binding to chemical bonding to the surface, and a barrier of 50 meV in going between neighboring physical binding sites. Our results support the possibility of realizing free-base porphyrins either physisorbed or chemisorbed on Cu(111) depending on the deposition procedure and experimental conditions.

Dynamics of growing carbon nanotube interfaces probed by machine learning-enabled molecular simulations.

Carbon nanotubes (CNTs), hollow cylinders of carbon, hold great promise for advanced technologies, provided their structure remains uniform throughout their length. Their growth takes place at high temperatures across a tube-catalyst interface. Structural defects formed during growth alter CNT properties. These defects are believed to form and heal at the tube-catalyst interface but an understanding of these mechanisms at the atomic-level is lacking. Here we present DeepCNT-22, a machine learning force field (MLFF) to drive molecular dynamics simulations through which we unveil the mechanisms of CNT formation, from nucleation to growth including defect formation and healing. We find the tube-catalyst interface to be highly dynamic, with large fluctuations in the chiral structure of the CNT-edge. This does not support continuous spiral growth as a general mechanism, instead, at these growth conditions, the growing tube edge exhibits significant configurational entropy. We demonstrate that defects form stochastically at the tube-catalyst interface, but under low growth rates and high temperatures, these heal before becoming incorporated in the tube wall, allowing CNTs to grow defect-free to seemingly unlimited lengths. These insights, not readily available through experiments, demonstrate the remarkable power of MLFF-driven simulations and fill long-standing gaps in our understanding of CNT growth mechanisms.

Topological Nodal Surface and Quadratic Dirac Semimetal States and van Hove Singularities in ScH3 and LuH3 Superconductors.

The coexistence of non-trivial topology and superconductivity in a material may induce a novel physical phenomenon known as topological superconductivity. Topological superconductors have been the subject of intense research, yet there are severe limitations in their application due to a lack of suitable materials. Topological nodal surface semimetals with nearly flat nodal surfaces near the Fermi level can be promising materials to achieve topological superconductivity. Here, we use first-principles calculations to examine the topological electronic characteristics of two new superconductors, ScH3 and LuH3, at both ambient and high pressures. Our studies show that both ScH3 and LuH3 have van Hove singularities, which confirms their superconductivity. Interestingly, both materials host topological nodal surface states under the protection of time reversal and spatial inversion symmetries in the absence of spin-orbit coupling (SOC). These nodal surfaces are distinguished by a pair of unique drum-head-like surface states not previously observed in nodal surface semimetals. Moreover, the nodal surfaces transform into essential spin-orbit quadratic Dirac points when SOC is included. Our findings demonstrate that ScH3 and LuH3 are good candidates to investigate the exotic properties of both nodal surface semimetals (NSSMs) and quadratic Dirac semimetal states and also provide a platform to explore the coexistence of topology and superconductivity in NSSMs with promising applications in high-speed electronics and topological quantum computing.

Physical binding energies using the electron localization function in 4-hydroxyphenylboronic acid co-crystals with aza donors.

Binding energies are traditionally simulated using cluster models by computation of each synthon for each individual co-crystal former. However, our investigation of the binding strengths using the electron localization function (ELF) reveals that these can be determined directly from the crystal supercell computations. We propose a new modeling protocol for the computation of physical binding energies directly from bulk simulations using ELF analysis. In this work, we establish a correlation between ELF values and binding energies calculated for co-crystals of 4-hydroxyphenylboronic acid (4HPBA) with four different aza donors using density functional theory with varying descriptions of dispersion. Boronic acids are gaining significant interest in the field of crystal engineering, but theoretical studies on their use in materials are still very limited. Here, we present a systematic investigation of the non-covalent interactions in experimentally realized co-crystals. Prior diffraction studies on these complexes have shown the competitive nature between the boronic acid functional group and the para-substituted phenolic group forming heteromeric interactions with aza donors. We determine the stability of the co-crystals by simulating their lattice energies, and the different dispersion descriptions show similar trends in lattice energies and lattice parameters. Our study bolsters the experimental observation of the boronic acid group as a competitive co-crystal former in addition to the well-studied phenolic group. Further research on correlating ELF values for physical binding could potentially transform this approach to a viable alternative for the computation of binding energies.

Exploring the Superior Anchoring Performance of the Two-Dimensional Nanosheets B2C4P2 and B3C2P3 for Lithium-Sulfur Batteries.

Potential anchoring materials in lithium-sulfur batteries help overcome the shuttle effect and achieve long-term cycling stability and high-rate efficiency. The present study investigates the two-dimensional nanosheets B2C4P2 and B3C2P3 by employing density functional theory calculations for their promise as anchoring materials. The nanosheets B2C4P2 and B3C2P3 bind polysulfides with adsorption energies in the range from -2.22 to -0.75 and -2.43 to -0.74 eV, respectively. A significant charge transfer occurs from the polysulfides, varying from -0.74 to -0.02e and -0.55 to -0.02e for B2C4P2 and B3C2P3, respectively. Upon anchoring the polysulfides, the band gap of B3C2P3 reduces, leading to enhanced electrical conductivity of the sulfur cathode. Finally, the calculated barrier energies of B2C4P2 and B3C2P3 for Li2S indicate fast diffusion of Li when recharged. These enthralling characteristics propose that the nanosheets B2C4P2 and B3C2P3 could reduce the shuttle effect in Li-S batteries and significantly improve their cycle performance, suggesting their promise as anchoring materials.

Peel-and-Stick Integration of Atomically Thin Nonlayered PtS Semiconductors for Multidimensionally Stretchable Electronic Devices.

Various near-atom-thickness two-dimensional (2D) van der Waals (vdW) crystals with unparalleled electromechanical properties have been explored for transformative devices. Currently, the availability of 2D vdW crystals is rather limited in nature as they are only obtained from certain mother crystals with intrinsically possessed layered crystallinity and anisotropic molecular bonding. Recent efforts to transform conventionally non-vdW three-dimensional (3D) crystals into ultrathin 2D-like structures have seen rapid developments to explore device building blocks of unique form factors. Herein, we explore a "peel-and-stick" approach, where a nonlayered 3D platinum sulfide (PtS) crystal, traditionally known as a cooperate mineral material, is transformed into a freestanding 2D-like membrane for electromechanical applications. The ultrathin (∼10 nm) 3D PtS films grown on large-area (>cm2) silicon dioxide/silicon (SiO2/Si) wafers are precisely "peeled" inside water retaining desired geometries via a capillary-force-driven surface wettability control. Subsequently, they are "sticked" on strain-engineered patterned substrates presenting prominent semiconducting properties, i.e., p-type transport with an optical band gap of ∼1.24 eV. A variety of mechanically deformable strain-invariant electronic devices have been demonstrated by this peel-and-stick method, including biaxially stretchable photodetectors and respiratory sensing face masks. This study offers new opportunities of 2D-like nonlayered semiconducting crystals for emerging mechanically reconfigurable and stretchable device technologies.

Nitrogen-doped carbon nanotubes: growth, mechanism and structure.

Nitrogen-doped bamboo-structured carbon nanotubes have been successfully grown using a series of cobalt/molybdenum catalysts. The morphology and structure of the nanotubes were analysed by transmission electron microscopy and Raman spectroscopy. The level of nitrogen doping, as determined by X-ray photoelectron spectroscopy, was found to range between 0.5 to 2.5 at.%. The growth of bamboo-structured nanotubes in the presence of nitrogen, in preference to single-walled and multi-walled nanotubes, was due to the greater binding energy of nitrogen for cobalt in the catalyst compared to the binding strength of carbon to cobalt, as determined by density functional theory.

Role of catalysts in dehydrogenation of MgH2 nanoclusters.

A fundamental understanding of the role of catalysts in dehydrogenation of MgH(2) nanoclusters is provided by carrying out first-principles calculations based on density functional theory. It is shown that the transition metal atoms Ti, V, Fe, and Ni not only lower desorption energies significantly but also continue to attract at least four hydrogen atoms even when the total hydrogen content of the cluster decreases. In particular, Fe is found to migrate from the surface sites to the interior sites during the dehydrogenation process, releasing more hydrogen as it diffuses. This diffusion mechanism may account for the fact that a small amount of catalysts is sufficient to improve the kinetics of MgH(2), which is essential for the use of this material for hydrogen storage in fuel-cell applications.

Spin-polarization mechanisms of the nitrogen-vacancy center in diamond.

The nitrogen-vacancy (NV) center in diamond has shown great promise for quantum information due to the ease of initializing the qubit and of reading out its state. Here we show the leading mechanism for these effects gives results opposite from experiment; instead both must rely on new physics. Furthermore, NV centers fabricated in nanometer-sized diamond clusters are stable, motivating a bottom-up qubit approach, with the possibility of quite different optical properties to bulk.

Rashba Effect and Raman Spectra of Tl2O/PtS2 Heterostructure.

The possibility to achieve charge-to-spin conversion via Rashba spin-orbit effects provides stimulating opportunities toward the development of nanoscale spintronics. Here, we use first-principles calculations to study the electronic and spintronic properties of Tl2O/PtS2 heterostructure, for which we have confirmed the dynamical stability by its positive phonon frequencies. An unexpectedly high binding energy of -0.38 eV per unit cell depicts strong interlayer interactions between Tl2O and PtS2. Interestingly, we discover Rashba spin-splittings (with a large α R value) in the valence band of Tl2O stemming from interfacial spin-orbit effects caused by PtS2. The role of van der Waals binding on the orbital rearrangements has been studied using the electron localization function and atomic orbital projections, which explains in detail the electronic dispersion near the Fermi level. Moreover, we explain the distinct band structure alignment in momentum space but separation in real space of Tl2O/PtS2 heterostructure. Since two-dimensional (2D) Tl2O still awaits experimental confirmation, we calculate, for the first time, the Raman spectra of pristine Tl2O and the Tl2O/PtS2 heterostructure and discuss peak positions corresponding to vibrational modes of the atoms. These findings offer a promising avenue to explore spin physics for potential spintronics applications via 2D heterostructures.

Inversion of the shuttlecock shaped metal phthalocyanines MPc (M = Ge, Sn, Pb)--a density functional study.

Shuttlecock shaped metal-phthalocyanine (MPc) can adsorb on a substrate surface having the central metal atom either down or up and the possibility of reversible switching between these two adsorption configurations shows great promise for use in nanomechanical devices. Using density functional theory we investigate the mechanism of the internal conformational inversion of germanium, tin and lead phthalocyanine in terms of the geometry, energy barrier of the reaction, and redox properties of the central metal atom. We have found the same mechanism of inversion for GePc and SnPc but a different one for PbPc. Inversion proceeds through two transition states, separated by a planar local minimum, for GePc and SnPc, but through one transition state distorting the phthalocyanine macrocycle for PbPc. The energy barrier of inversion is 4.27 eV for PbPc and 2.12 and 3.16 eV for GePc and SnPc, respectively. Such high barriers are unlikely to be overcome at normal experimental conditions, and in many cases alternative explanations for switching between "up" and "down" conformation need to be sought, such as ionization assisted inversion or even flipping over of the molecules. Our calculations show that the inversion of GePc and SnPc is accompanied by reversible two electron oxidation (M(II)<--> M(IV)) of the metal atom, through intersystem crossing. The difference in mechanism of inversion for GePc (SnPc) and PbPc is assigned to the different nature of the central metal atom.

Distinguishing between chemical bonding and physical binding using electron localization function (ELF).

To distinguish between chemical bonding and physical binding is usually simple. They differ, in the normal case, in both interaction strength (binding energy) and interaction length (structure). However, chemical bonding can be weak (e.g. in some metallic bonding) and physical binding can be strong (e.g. due to permanent electrostatic moments, hydrogen binding, etc) making differentiation non-trivial. But since these are shared-electron or unshared-electron interactions, respectively, it is in principle possible to distinguish the type of interaction by analyzing the electron density around the interaction point(s)/interface. After all, the former should be a contact while the latter should be a tunneling barrier. Here, we investigate within the framework of density functional theory typical molecules and crystals to show the behaviour of the electron localization function (ELF) in different shared-electron interactions, such as chemical (covalent) and metallic bonding and compare to unshared-electron interactions typical for physical binding, such as ionic, hydrogen and Keesom, dispersion (van der Waals) binding and attempt to categorise them only by the ELF and the electron population in the interaction region. It is found that the ELF method is not only useful for the characterization of covalent bonds but a lot of information can be extracted also for weaker types of binding. Furthermore, the charge integration over the interaction region(s) and tracing the ELF profile can reveal the strength of the bonding/binding ranging from the triple bonds to weak dispersion.

Large-area 2D PtTe2/silicon vertical-junction devices with ultrafast and high-sensitivity photodetection and photovoltaic enhancement by integrating water droplets.

2D PtTe2 layers, a relatively new class of 2D crystals, have unique band structure and remarkably high electrical conductivity promising for emergent opto-electronics. This intrinsic superiority can be further leveraged toward practical device applications by merging them with mature 3D semiconductors, which has remained largely unexplored. Herein, we explored 2D/3D heterojunction devices by directly growing large-area (>cm2) 2D PtTe2 layers on Si wafers using a low-temperature CVD method and unveiled their superior opto-electrical characteristics. The devices exhibited excellent Schottky transport characteristics essential for high-performance photovoltaics and photodetection, i.e., well-balanced combination of high photodetectivity (>1013 Jones), small photo-responsiveness time (∼1 µs), high current rectification ratio (>105), and water super-hydrophobicity driven photovoltaic improvement (>300%). These performances were identified to be superior to those of previously explored 2D/3D or 2D layer-based devices with much smaller junction areas, and their underlying principles were confirmed by DFT calculations.

Wafer-Scale Growth of 2D PtTe2 with Layer Orientation Tunable High Electrical Conductivity and Superior Hydrophobicity.

Platinum ditelluride (PtTe2) is an emerging semimetallic two-dimensional (2D) transition-metal dichalcogenide (TMDC) crystal with intriguing band structures and unusual topological properties. Despite much devoted efforts, scalable and controllable synthesis of large-area 2D PtTe2 with well-defined layer orientation has not been established, leaving its projected structure-property relationship largely unclarified. Herein, we report a scalable low-temperature growth of 2D PtTe2 layers on an area greater than a few square centimeters by reacting Pt thin films of controlled thickness with vaporized tellurium at 400 °C. We systematically investigated their thickness-dependent 2D layer orientation as well as its correlated electrical conductivity and surface property. We unveil that 2D PtTe2 layers undergo three distinct growth mode transitions, i.e., horizontally aligned holey layers, continuous layer-by-layer lateral growth, and horizontal-to-vertical layer transition. This growth transition is a consequence of competing thermodynamic and kinetic factors dictated by accumulating internal strain, analogous to the transition of Frank-van der Merwe (FM) to Stranski-Krastanov (SK) growth in epitaxial thin-film models. The exclusive role of the strain on dictating 2D layer orientation has been quantitatively verified by the transmission electron microscopy (TEM) strain mapping analysis. These centimeter-scale 2D PtTe2 layers exhibit layer orientation tunable metallic transports yielding the highest value of ∼1.7 × 106 S/m at a certain critical thickness, supported by a combined verification of density functional theory (DFT) and electrical measurements. Moreover, they show intrinsically high hydrophobicity manifested by the water contact angle (WCA) value up to ∼117°, which is the highest among all reported 2D TMDCs of comparable dimensions and geometries. Accordingly, this study confirms the high material quality of these emerging large-area 2D PtTe2 layers, projecting vast opportunities employing their tunable layer morphology and semimetallic properties from investigations of novel quantum phenomena to applications in electrocatalysis.

The importance of strong carbon-metal adhesion for catalytic nucleation of single-walled carbon nanotubes.

Density functional theory is used to show that the adhesion between single-walled carbon nanotubes (SWNTs) and the catalyst particles from which they grow needs to be strong to support nanotube growth. It is found that Fe, Co, and Ni, commonly used to catalyze SWNT growth, have larger adhesion strengths to SWNTs than Cu, Pd, and Au and are therefore likely to be more efficient for supporting growth. The calculations also show that to maintain an open end of the SWNT it is necessary that the SWNT adhesion strength to the metal particle is comparable to the cap formation energy of the SWNT end. This implies that the difference between continued and discontinued SWNT growth to a large extent depends on the carbon-metal binding strength, which we demonstrate by molecular dynamics (MD) simulations. The results highlight that first principles computations are vital for the understanding of the binding strength's role in the SWNT growth mechanism and are needed to get accurate force field parameters for MD.

Bonding of atomic phosphorus to polycyclic hydrocarbons and curved graphitic surfaces.

We present a theoretical study of the bonding of atomic phosphorus to planar hydrocarbons and to curved graphite-like surfaces. We find that bonding of phosphorus to planar polycyclic hydrocarbons induces curvature away from the phosphorus atom, as defined by the pyramidalization angle. Similarly, bonding of atomic phosphorus to the [5,5] fulvalene-circulene semifullerene and buckminsterfullerene is only possible on the convex side of the carbon surface. On the other hand, we find the interaction of atomic phosphorus with the concave side of fullerene-like surfaces to be nonbonding for both quartet and doublet spin states. We find the prerequisite for stable epoxy-type bonds within these systems is the ability of the carbon atoms to maintain or induce curvature away from the P.C=C bond.