Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional

Tunable Redox Temperature of a Co3-xMnxO4 (0 ≤ x ≤ 3) Continuous Solid Solution for Thermochemical Energy Storage

Heat-storage technologies are well suited to improve the energy efficiency of power plants and the recovery of process heat. A good option for high storage capacities, especially at high temperatures, is storing thermal energy by reversible thermochemical reactions. In particular, the Co₃O₄/CoO and Mn₂O₃/Mn₃O₄ redox-active couples are known to be very promising systems. However, cost and toxicity issues for Co oxides and the sluggish oxidation rate (leading to poor reversibility) for Mn oxide hinder the applicability of these single oxides. Considering, instead, binary Co-Mn oxide mixtures could mitigate the above-mentioned shortcomings. To examine this in detail, here, we combine first-principles atomistic calculations and experiments to provide a structural characterization and observe the thermal behavior of novel mixed-metal oxides based on cobalt/manganese metals with the spinel structure Co_3-xMn_xO₄. We show that novel Co_3-xMn_xO₄ phases indeed enhance the enthalpy of the redox reactions, facilitate reversibility, and mitigate energy losses when compared to pure metal oxide systems. Our results expand therefore the limited list of currently available thermochemical heat-storage materials and pave the way toward the implementation of tunable redox temperature materials for practical applications.

Improved description of hematite surfaces by the SCAN functional

Controversies on the surface termination of α-Fe₂O₃ (0001) focus on its surface stoichiometry dependence on the oxygen chemical potential. Density functional theory (DFT) calculations applying the commonly accepted Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional to a strongly correlated system predict the best matching surface termination, but would produce a delocalization error, resulting in an inappropriate bandgap, and thus are not applicable for comprehensive hematite system studies. Besides, the widely applied PBE+U scheme cannot provide evidence for existence of some of the successfully synthesized stoichiometric α-Fe₂O₃ (0001) surfaces. Hence, a better scheme is needed for hematite DFT studies. This work investigates whether the strongly constrained and appropriately normed (SCAN) approximation reported by Perdew et al. could provide an improved result for the as-mentioned problem, and whether SCAN can be applied to hematite systems. By comparing the results calculated with the PBE, SCAN, PBE+U, and SCAN+U schemes, we find that SCAN and SCAN+U improves the description of the electronic structure of different stoichiometric α-Fe₂O₃ (0001) surfaces with respect to the PBE results, and that they give a consistent prediction of the surface terminations. Besides, the bulk lattice constants and the bulk density of states are also improved with the SCAN functional. This study provides a general characterization of the α-Fe₂O₃ (0001) surfaces and rationalizes how the SCAN approximation improves the results of hematite surface calculations.

Composition-induced type I and direct bandgap transition metal dichalcogenides alloy vertical heterojunctions

While members of the 2D semiconducting transition metal dichalcogenide (TMD) family MX₂ (M = {Mo, W}, X = {S, Se}) are promising for device applications, stacked layer (vertical) heterojunctions exhibit features that make them inappropriate for light-emitting applications. Such vertical heterojunctions exhibit type II, rather than the preferred type I band alignment. Using density functional theory calculations, we identify the pseudo-binary and quaternary alloy composition range for which band alignment is type I. While broad regions of composition space lead to type I band alignment, most light-emitting devices require direct bandgaps. We demonstrate that by taking advantage of alloying and/or twisting between layers, a wide range of type I, direct bandgap stacked layer (vertical) heterojunctions are achievable. These results and the underlying method developed here provide new opportunities for TMD vertical heterojunction device optimization and opens the door to new classes of TMD vertical heterojunction device applications.

Adsorption processes on a Pd monolayer-modified Pt(111) electrode

Specific adsorption of anions is an important aspect in surface electrochemistry for its influence on reaction kinetics in either a promoted or inhibited fashion. Perchloric acid is typically considered as an ideal electrolyte for investigating electrocatalytic reactions due to the lack of specific adsorption of the perchlorate anion on several metal electrodes. In this work, cyclic voltammetry and computational methods are combined to investigate the interfacial processes on a Pd monolayer deposited on Pt(111) single crystal electrode in perchloric acid solution. The “hydrogen region” of this Pd_MLPt(111) surface exhibits two voltammetric peaks: the first “hydrogen peak” at 0.246 V_RHE actually involves the replacement of hydrogen by hydroxyl, and the second “hydrogen peak” H_II at 0.306 V_RHE appears to be the replacement of adsorbed hydroxyl by specific perchlorate adsorption. The two peaks merge into a single peak when a more strongly adsorbed anion, such as sulfate, is involved. Our density functional theory calculations qualitatively support the peak assignment and show that anions generally bind more strongly to the Pd_MLPt(111) surface than to Pt(111).

Specific adsorption of anions is an important aspect in surface electrochemistry for its influence on reaction kinetics in either a promoted or inhibited fashion.

Competing stripe and magnetic phases in the cuprates from first principles

Realistic description of competing phases in complex quantum materials has proven extremely challenging. For example, much of the existing density-functional-theory-based first-principles framework fails in the cuprate superconductors. Various many-body approaches involve generic model Hamiltonians and do not account for the interplay between the spin, charge, and lattice degrees of freedom. Here, by deploying the recently constructed strongly constrained and appropriately normed (SCAN) density functional, we show how the landscape of competing stripe and magnetic phases can be addressed on a first-principles basis both in the parent insulator YBa2Cu3O6 and the near-optimally doped YBa2Cu3O7 as archetype cuprate compounds. In YBa2Cu3O7, we find many stripe phases that are nearly degenerate with the ground state and may give rise to the pseudogap state from which the high-temperature superconducting state emerges. We invoke no free parameters such as the Hubbard U, which has been the basis of much of the existing cuprate literature. Lattice degrees of freedom are found to be crucially important in stabilizing the various phases.

First-Principles Calculation of Water pKa Using the Newly Developed SCAN Functional

Acid/base chemistry is an intriguing topic that still constitutes a challenge for computational chemistry. While estimating the acid dissociation constant (or pK_a) could shed light on many chemistry processes, especially in the fields of biochemistry and geochemistry, evaluating the relative stability between protonated and nonprotonated species is often very difficult. Indeed, a prerequisite for calculating the pK_a of any molecule is an accurate description of the energetics of water dissociation. Here, we applied constrained molecular dynamics simulations, a noncanonical sampling technique, to investigate the water deprotonation process by selecting the OH distance as the reaction coordinate. The calculation is based on density functional theory and the newly developed SCAN functional, which has shown excellent performance in describing water structure. This first benchmark of SCAN on a chemical reaction shows that this functional accurately models the energetics of proton transfer reactions in an aqueous environment. After taking Coulomb long-range corrections and nuclear quantum effects into account, the estimated water pK_a is only 1.0 pK_a unit different from the target experimental value. Our results show that the combination of SCAN and constrained MD successfully reproduces the chemistry of water and constitutes a good framework for calculating the free energy of chemical reactions of interest.

Equation of state of water based on the SCAN meta-GGA density functional

By ab initio molecular dynamics simulations, the newly developed SCAN meta-GGA functional is proved better than the widely used PBE-GGA functional in describing the equation of state of water.

Structural diversity in conducting bilayer salts (CNB-EDT-TTF)4A

The family of recently described salts based on the electron donor CNB-EDT-TTF and different anions A, with general formula (CNB-EDT-TTF)₄A, constitutes an unprecedented type of molecular conductor based on a bilayer structure of the donors.

Similarities and differences between potassium and ammonium ions in liquid water: a first-principles study

Understanding ion solvation in liquid water is critical in optimizing materials for a wide variety of emerging technologies, including water desalination and purification.

Quantifying the hydration structure of sodium and potassium ions: taking additional steps on Jacob's Ladder

The ability to reproduce the experimental structure of water around the sodium and potassium ions is a key test of the quality of interaction potentials due to the central importance of these ions in a wide range of important phenomena.

Aqueous solvation of the chloride ion revisited with density functional theory: impact of correlation and exchange approximations

Aqueous chloride is simulated using PBE-D3, PBE0-D3, and SCAN to investigate the impact of exchange and correlation approximations; we find the exact exchange fraction strongly impacts the energetics and polarizability of solvated chloride.

Vibrational mode frequency correction of liquid water in density functional theory molecular dynamics simulations with van der Waals correction

We develop a frequency correction scheme for the stretch and bending modes of liquid water, which substantially improves the prediction of the vibrational spectra.

Localization in the SCAN meta-generalized gradient approximation functional leading to broken symmetry ground states for graphene and benzene

SCAN over localizes orbitals leading to spin symmetry broken ground states in graphene and benzene.

A Novel Method for Calculation of Molecular Energies and Charge Distributions by Thermodynamic Formalization

The paper describes a new approach to the thermodynamic formalization for calculation of molecular energy and charge distribution in ground state by means of the variational equation of DFT. In order to thermodynamically formalize the molecular calculation, the pseudo chemical potential (PCP) is conceptualized, where a molecule is broken into multi-phase(atom) one-component(electron) systems and the energy of system is represented as PCP. Calculation of the molecular energy and atomic charge by PCP is put forward, thereafter the approach is proved to be valid and its efficiency (accuracy and calculation speed) is verified.

First-Principles Calculation of Transition Metal Hyperfine Coupling Constants with the Strongly Constrained and Appropriately Normed (SCAN) Density Functional and its Hybrid Variants

Density functional theory (DFT) is used extensively for the first-principles calculation of hyperfine coupling constants in both main-group and transition metal systems. As with many other properties, the performance of DFT for hyperfine coupling constants is of variable quality, particularly for transition metal complexes, because it strongly depends on the nature of the chemical system and the type of approximation to the exchange-correlation functional. Recently, a meta-generalized-gradient approximation (mGGA) functional was proposed that obeys all known exact constraints for such a method, known as the Strongly Constrained and Appropriately Normed (SCAN) functional. In view of its theoretically superior formulation a benchmark set of complexes is used to assess the performance of SCAN for the challenging case of transition metal hyperfine coupling constants. In addition, two global hybrid versions of the functional, SCANh and SCAN0, are described and tested. The values computed with the new functionals are compared with experiment and with those of other DFT approximations. Although the original SCAN and the SCAN-based hybrids may offer improved hyperfine coupling constants for specific systems, no uniform improvement is observed. On the contrary, there are specific cases where the new functionals fail badly due to a flawed description of the underlying electronic structure. Therefore, despite these methodological advances, systematically accurate and system-independent prediction of transition metal hyperfine coupling constants with DFT remains an unmet challenge.

Chemisorption Can Reverse Defect-Defect Interaction on Heterogeneous Catalyst Surfaces

Atomic-level understanding of roles of defect-defect interaction in the bonding of adsorbates on surfaces is critical for tailoring catalysts atom-by-atom and designing new catalysts. Here, from first-principles calculations, we propose a microscopic mechanism for the role of sulfur vacancy-vacancy interaction in hydrogen bonding on surfaces of MoS₂, a nonprecious two-dimensional catalyst for hydrogen evolution reaction. We find that before hydrogen adsorption the interaction of a sulfur vacancy with others is repulsive, originating from the antibonding-like coupling of occupied in-gap vacancy states. When the sulfur vacancy is adsorbed by a hydrogen atom, its interaction with other unadsorbed sulfur vacancies becomes attractive, which can be attributed to the decoupling of repulsive vacancy-vacancy interactions and the occupying of bonding-like coupling states between the in-gap vacancy states that are unoccupied before hydrogen adsorption. This repulsive-to-attractive reverse of vacancy-vacancy interaction reduces the hydrogen adsorption energy and explains why the hydrogen adsorption energy decreases with increasing sulfur vacancy concentration. The emerging picture enables a more general discussion of local defect effects on the adsorption of various adsorbates at different surfaces, providing guidance to improve catalytic performance through defect engineering.

Semimetal or Semiconductor: The Nature of High Intrinsic Electrical Conductivity in TiS2

As an intensively studied electrode material for secondary batteries, TiS₂ is known to exhibit high electrical conductivity without extrinsic doping. However, the origin of this high conductivity, either being a semimetal or a heavily self-doped semiconductor, has been debated for several decades. Here, combining quasi-particle GW calculations, density functional theory (DFT) study on intrinsic defects, and scanning tunneling microscopy/spectroscopy (STM/STS) measurements, we conclude that stoichiometric TiS₂ is a semiconductor with an indirect band gap of about 0.5 eV. The high conductivity of TiS₂ is therefore caused by heavy self-doping. Our DFT results suggest that the dominant donor defect that is responsible for the self-doping under thermal equilibrium is Ti interstitial, which is corroborated by our STM/STS measurements.

Tunable Schottky and Ohmic contacts in graphene and tellurene van der Waals heterostructures

We systematically investigate the effects of external electric field and interlayer coupling on the electronic structures and contact characteristics of hybrid graphene and tellurene (G/Te) van der Waals heterostructures (vdWHs) based on first-principles calculations. Our results show that the G/α-Te interface is formed by an n-type Schottky contact with a negligible Schottky barrier height (SBH), while the G/β-Te interface is formed by a p-type Schottky contact with a SBH of 0.51 eV. By applying external electric fields perpendicular to the G/Te interfaces or changing the interlayer distance between the graphene and tellurene monolayers, both Schottky barriers and contact types (n-type Schottky, p-type Schottky, and Ohmic) at the G/Te interfaces can be effectively modulated. The changes in charge transfer, as well as the corresponding interface dipole and potential step with the external electric field and interlayer coupling, are revealed to account for the reason for tunable Schottky and Ohmic contacts at the G/Te interfaces. Therefore, the G/Te vdWHs show tunable Schottky and Ohmic contacts with promising applications of graphene-based field-effect transistors in future experiments.

IRMOF-74(n)-Mg: a novel catalyst series for hydrogen activation and hydrogenolysis of C-O bonds

Metal-Organic Frameworks (MOFs) that catalyze hydrogenolysis reactions are rare and there is little understanding of how the MOF, hydrogen, and substrate molecules interact. In this regard, the isoreticular IRMOF-74 series, two of which are known catalysts for hydrogenolysis of aromatic C-O bonds, provides an unusual opportunity for systematic probing of these reactions. The diameter of the 1D open channels can be varied within a common topology owing to the common secondary building unit (SBU) and controllable length of the hydroxy-carboxylate struts. We show that the first four members of the IRMOF-74(Mg) series are inherently catalytic for aromatic C-O bond hydrogenolysis and that the conversion varies non-monotonically with pore size. These catalysts are recyclable and reusable, retaining their crystallinity and framework structure after the hydrogenolysis reaction. The hydrogenolysis conversion of phenylethylphenyl ether (PPE), benzylphenyl ether (BPE), and diphenyl ether (DPE) varies as PPE > BPE > DPE, consistent with the strength of the C-O bond. Counterintuitively, however, the conversion also follows the trend IRMOF-74(III) > IRMOF-74(IV) > IRMOF-74(II) > IRMOF-74(I), with little variation in the corresponding selectivity. DFT calculations suggest the unexpected behavior is due to much stronger ether and phenol binding to the Mg(ii) open metal sites (OMS) of IRMOF-74(III), resulting from a structural distortion that moves the Mg2+ ions toward the interior of the pore. Solid-state 25Mg NMR data indicate that both H2 and ether molecules interact with the Mg(ii) OMS and hydrogen-deuterium exchange reactions show that these MOFs activate dihydrogen bonds. The results suggest that both confinement and the presence of reactive metals are essential for achieving the high catalytic activity, but that subtle variations in pore structure can significantly affect the catalysis. Moreover, they challenge the notion that simply increasing MOF pore size within a constant topology will lead to higher conversions.

Striated 2D Lattice with Sub-nm 1D Etch Channels by Controlled Thermally Induced Phase Transformations of PdSe2

2D crystals are typically uniform and periodic in-plane with stacked sheet-like structure in the out-of-plane direction. Breaking the in-plane 2D symmetry by creating unique lattice structures offers anisotropic electronic and optical responses that have potential in nanoelectronics. However, creating nanoscale-modulated anisotropic 2D lattices is challenging and is mostly done using top-down lithographic methods with ≈10 nm resolution. A phase transformation mechanism for creating 2D striated lattice systems is revealed, where controlled thermal annealing induces Se loss in few-layered PdSe₂ and leads to 1D sub-nm etched channels in Pd₂ Se₃ bilayers. These striated 2D crystals cannot be described by a typical unit cells of 1-2 Å for crystals, but rather long range nanoscale periodicity in each three directions. The 1D channels give rise to localized conduction states, which have no bulk layered counterpart or monolayer form. These results show how the known family of 2D crystals can be extended beyond those that exist as bulk layered van der Waals crystals by exploiting phase transformations by elemental depletion in binary systems.

Bidirectional heterostructures consisting of graphene and lateral MoS2/WS2 composites: a first-principles study

First-principles calculations have been performed to explore the structural and electronic properties of bidirectional heterostructures composed of graphene and (MoS₂) _X /(WS₂)_4-X (X = 1, 2, 3) lateral composites and compare them with those of heterobilayers formed by graphene and pristine MS₂ (M = Mo, W). The band gaps of the lateral heterostructures lie between those of pristine MoS₂ and WS₂. The weak coupling between the two layers can induce a tiny band-gap opening of graphene and formation of an n-type Schottky contact at the G-(MoS₂) _X /(WS₂)_4-X interface. Moreover, the combination ratio of MoS₂/WS₂ can control the electronic properties of G-(MoS₂) _X /(WS₂)_4-X . By applying external electric fields, the band gaps of (MoS₂) _X /(WS₂)_4-X (X = 0, 1, 2, 3, 4) monolayers undergo a direct-indirect transition, and semiconductor-metal transitions can be found in WS₂. External electric fields can also be used effectively to tune the binding energies, charge transfers, and band structures (the types of Schottky and Ohmic contacts) of G-(MoS₂) _X /(WS₂)_4-X heterostructures. These findings suggest that G-(MoS₂) _X /(WS₂)_4-X heterostructures can serve as high-performance nano-electronic devices.

Thermoelectric power factor of pure and doped ZnSb via DFT based defect calculations

The power factor of pure p-type ZnSb has been calculated via ab initio simulations assuming that the carrier concentrations are due to the doping effect of intrinsic zinc vacancies. With a vacancy concentration close to the experimental solubility limit we were able to perfectly reproduce the Power Factor measured in polycrystalline ZnSb samples. The methodology has then been successfully extended for predicting the effect of extrinsic doping elements on the thermoelectric properties of ZnSb. Germanium and tin seem to be promising p-type doping elements. In addition, we give, for the first time, an explanation of why it is difficult to synthesize polycrystalline n-type ZnSb samples. Indeed, compensative effects between intrinsic defects (zinc vacancies) and doping elements (Ga, or In) explain the existence of an optimal (and relatively high) dopant concentration necessary to convert ZnSb into an n-type semiconductor.

Machine Learning the Physical Nonlocal Exchange-Correlation Functional of Density-Functional Theory

We train a neural network as the universal exchange-correlation functional of density-functional theory that simultaneously reproduces both the exact exchange-correlation energy and the potential. This functional is extremely nonlocal but retains the computational scaling of traditional local or semilocal approximations. It therefore holds the promise of solving some of the delocalization problems that plague density-functional theory, while maintaining the computational efficiency that characterizes the Kohn-Sham equations. Furthermore, by using automatic differentiation, a capability present in modern machine-learning frameworks, we impose the exact mathematical relation between the exchange-correlation energy and the potential, leading to a fully consistent method. We demonstrate the feasibility of our approach by looking at one-dimensional systems with two strongly correlated electrons, where density-functional methods are known to fail, and investigate the behavior and performance of our functional by varying the degree of nonlocality.

A charge optimized many-body potential for iron/iron-fluoride systems

A classical interatomic potential for iron/iron-fluoride systems is developed in the framework of the charge optimized many-body (COMB) potential. This interatomic potential takes into consideration the effects of charge transfer and many-body interactions depending on the chemical environment. The potential is fitted to a training set composed of both experimental and ab initio results of the cohesive energies of several Fe and FeF2 crystal phases, the two fluorine molecules F2 and the F2-1 dissociation energy curve, the Fe and FeF2 lattice parameters of the ground state crystalline phase, and the elastic constants of the body centered cubic Fe structure. The potential is tested in an NVT ensemble for different initial structural configurations as the crystal ground state phases, F2 molecules, iron clusters, and iron nanospheres. In particular, we model the FeF2/Fe bilayer and multilayer interfaces, as well as a system of square FeF2 nanowires immersed in an iron solid. It has been shown that there exists a reordering of the atomic positions for F and Fe atoms at the interface zone; this rearrangement leads to an increase in the charge transfer among the atoms that make the interface and put forward a possible mechanism of the exchange bias origin based on asymmetric electric charge transfer in the different spin channels.

Observation of Topological Edge States at the Step Edges on the Surface of Type-II Weyl Semimetal TaIrTe4

Topological materials harbor topologically protected boundary states. Recently, TaIrTe₄, a ternary transition-metal dichalcogenide, was identified as a type-II Weyl semimetal with the minimal nonzero number of Weyl points allowed for a time-reversal invariant Weyl semimetal. Monolayer TaIrTe₄ was proposed to host topological edge states, which, however, lacks of experimental evidence. Here, we report on the topological edge states localized at the monolayer step edges of the type-II Weyl semimetal TaIrTe₄ using scanning tunneling microscopy. One-dimensional electronic states that show substantial robustness against the edge irregularity are observed at the step edges. Theoretical calculations substantiate the topologically nontrivial nature of the edge states and their robustness against the edge termination and layer stacking. The observation of topological edge states at the step edges of TaIrTe₄ surfaces suggests that monolayer TaIrTe₄ is a two-dimensional topological insulator, providing TaIrTe₄ as a promising material for topological physics and devices.

K-Shell Core-Electron Excitations in Electronic Stopping of Protons in Water from First Principles

Understanding the role of core-electron excitation in liquid water under proton irradiation has become important due to the growing use of proton beams in radiation oncology. Using a first-principles, nonequilibrium simulation approach based on real-time, time-dependent density functional theory, we determine the electronic stopping power, the velocity-dependent energy transfer rate from irradiating ions to electrons. The electronic stopping power curve agrees quantitatively with experimental data over the velocity range available. At the same time, significant differences are observed between our first-principles result and commonly used perturbation theoretic models. Excitations of the water molecules' oxygen core electrons are a crucial factor in determining the electronic stopping power curve beyond its maximum. The core-electron contribution is responsible for as much as one third of the stopping power at the high proton velocity of 8.0 a.u. (1.6 MeV). K-shell core-electron excitations not only provide an additional channel for the energy transfer-they also significantly influence the valence electron excitations. In the excitation process, generated holes remain highly localized within a few angstroms around the irradiating proton path, whereas electrons are excited away from the path. In spite of their great contribution to the stopping power, K-shell electrons play a rather minor role in terms of the excitation density; only 1% of the hole population composes K-shell holes, even at the high proton velocity of 8.0 a.u. The excitation behavior revealed is distinctly different from that of photon-based ionizing radiation such as x or γ rays.

Theoretical investigations on charge transport properties of tetrabenzo[a,d,j,m]coronene derivatives using different density functional theory functionals (B3LYP, M06-2X, and wB97XD)

Charge transport rate is one of the key parameters determining the performance of organic electronic devices. Based on density functional theory, exchange-correlation functionals which adequately account for non-covalent interactions, such as M06-2X and wB97XD, should significantly improve the accuracy of charge transport rate calculations for large systems with non-covalent interactions. In this work, the B3LYP hybrid functional, the variant hybrid functional M06-2X, and the long-range-corrected wB97XD functional were used to perform geometry optimizations and charge transport rate calculations on 11 variants of tetrabenzo[ a,d,j,m]coronene, including tetrabenzo[ a,d,j,m]coronene itself and its tetra-substituted and octa-substituted derivatives. Our results indicate that the molecular geometries of these benzocoronene semiconductors are large quasi-planar conjugated π systems, and the incorporation of different substituents significantly affects their frontier molecular orbitals. The hole carrier mobility ( µ+) and electron carrier mobility ( µ−) of the methoxy-substituted derivatives (TBC(OCH3)4and TBC(OCH3)8) were relatively low. The results of the tetrabenzo[ a,d,j,m]coronene molecules studied were consistent with using the aforementioned M06-2X, wB97XD, and B3LYP methods. We found that the octa-substituted derivatives (TBCF8, TBCCl8, TBC(CH3)8, and TBC(CN)8) could be used as p-type organic semiconductor materials.

Halogen Bond Structure and Dynamics from Molecular Simulations

Halogen bonding has emerged as an important noncovalent interaction in a myriad of applications, including drug design, supramolecular assembly, and catalysis. The current understanding of the halogen bond is informed by electronic structure calculations on isolated molecules and/or crystal structures that are not readily transferable to liquids and disordered phases. To address this issue, we present a first-principles simulation-based approach for quantifying halogen bonds in molecular systems rooted in an understanding of nuclei-nuclei and electron-nuclei spatial correlations. We then demonstrate how this approach can be used to quantify the structure and dynamics of halogen bonds in condensed phases, using solid and liquid molecular chlorine as prototypical examples with high concentrations of halogen bonds. We close with a discussion of how the knowledge generated by our first-principles approach may inform the development of classical empirical models, with a consistent representation of halogen bonding.

An ultralow-density porous ice with the largest internal cavity identified in the water phase diagram

The recent back-to-back findings of low-density porous ice XVI and XVII have rekindled the century-old field of the solid-state physics and chemistry of water. Experimentally, both ice XVI and XVII crystals can be produced by extracting guest atoms or molecules enclosed in the cavities of preformed ice clathrate hydrates. Herein, we examine more than 200 hypothetical low-density porous ices whose structures were generated according to a database of zeolite structures. Hitherto unreported porous EMT ice, named according to zeolite nomenclature, is identified to have an extremely low density of 0.5 g/cm3 and the largest internal cavity (7.88 Å in average radius). The EMT ice can be viewed as dumbbell-shaped motifs in a hexagonal close-packed structure. Our first-principles computations and molecular dynamics simulations confirm that the EMT ice is stable under negative pressures and exhibits higher thermal stability than other ultralow-density ices. If all cavities are fully occupied by hydrogen molecules, the EMT ice hydrate can easily outperform the record hydrogen storage capacity of 5.3 wt % achieved with sII hydrogen hydrate. Most importantly, in the reconstructed temperature-pressure (T-P) phase diagram of water, the EMT ice is located at deeply negative pressure regions below ice XVI and at higher temperature regions next to FAU. Last, the phonon spectra of empty-sII, FAU, EMT, and other zeolite-like ice structures are computed by using the dispersion corrected vdW-DF2 functional. Compared with those of ice XI (0.93 g/cm3), both the bending and stretching vibrational modes of the EMT ice are blue-shifted due to their weaker hydrogen bonds.

Tuning the magnetism of two-dimensional hematene by ferroelectric polarization

Magnetism in two-dimensional (2D) materials, that is, a 2D version of the magnetism of three-dimensional bulk materials, and the associated novel physics have recently been the focus of many spintronics researchers. Here we investigate the manipulation of 2D magnetism at the interfaces of ferromagnetic/ferroelectric hematene/BaTiO₃(001) heterostructures (HSs) fabricated via a precisely chosen sequence. By introducing four types of interfaces of 2D hematene and three-dimensional BaTiO₃ that induce different oxygen environments, the control of magnetism is directly demonstrated from first-principles. An obvious 2D electron gas originates from the Fe-3d and O-2p hybridization; the electron gas is sensitive to the interfacial atomic displacements. Robust control of both the direction and magnitude of the net magnetization has been realized for an Fe/TiO₂ terminated bilayer HS. The electron occupancies of the d_xy and d_xz orbitals and changes to the Fe-O bond play a key role in determining the magnetism of our systems. Our work not only demonstrates the technique's potential for manipulating magnetism in 2D hematene, but also sheds light on the underlying mechanism and the fundamental properties of hematene HSs.

Indirect but Efficient: Laser-Excited Electrons Can Drive Ultrafast Polarization Switching in Ferroelectric Materials

To enhance the efficiency of next-generation ferroelectric (FE) electronic devices, new techniques for controlling ferroelectric polarization switching are required. While most prior studies have attempted to induce polarization switching via the excitation of phonons, these experimental techniques required intricate and expensive terahertz sources and have not been completely successful. Here, we propose a new mechanism for rapidly and efficiently switching the FE polarization via laser-tuning of the underlying dynamical potential energy surface. Using time-dependent density functional calculations, we observe an ultrafast switching of the FE polarization in BaTiO₃ within 200 fs. A laser pulse can induce a charge density redistribution that reduces the original FE charge order. This excitation results in both desirable and highly directional ionic forces that are always opposite to the original FE displacements. Our new mechanism enables the reversible switching of the FE polarization with optical pulses that can be produced from existing 800 nm experimental laser sources.

Origin of structural stability of ScH3 molecular nanowires and their chemical-bonding behavior: Correlation effects of the Sc 3d electrons

A new stable transition-metal trihydride (ScH₃) molecular nanowire was recently reported by Li et al. [J. Am. Chem. Soc. 139, 6290-6293 (2017)]. Of the two typical structures (T-ScH₃ and O-ScH₃), T-ScH₃ is more stable than O-ScH₃. However, the reason why O-ScH₃ is less stable than T-ScH₃ was not known. Using Perdew-Burke-Ernzerhof (PBE), PBE+U, SCAN, and HSE06, as well as crystal orbital Hamilton populations (COHPs), we investigate the orbital-projected band structures and chemical bonding of T-ScH₃ and O-ScH₃. It is found that the energies calculated by PBE, SCAN, and HSE06 indeed reveal that T-ScH₃ is more stable than O-ScH₃, and there is no occupied antibonding state at the Fermi level of the COHP curves of T-ScH₃, supporting the stable Sc-H bonding of T-ScH₃. To the contrary, the Sc-H bonding of O-ScH₃ is unstable because there exist occupied antibonding states at the Fermi level of the COHP curves of O-ScH₃. We found that the results of PBE+U are consistent with those of PBE, SCAN, and HSE06 in the case of U < U_c. However, when U > U_c, the results of PBE+U are opposite to those of PBE, SCAN, and HSE06.