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

Optically Active MXenes in Van der Waals Heterostructures

The vertical integration of distinct 2D materials in van der Waals (vdW) heterostructures provides the opportunity for interface engineering and modulation of electronic as well as optical properties. However, scarce experimental studies reveal many challenges for vdW heterostructures, hampering the fine-tuning of their electronic and optical functionalities. Optically active MXenes, the most recent member of the 2D family, with excellent hydrophilicity, rich surface chemistry, and intriguing optical properties, are a novel 2D platform for optoelectronics applications. Coupling MXenes with various 2D materials into vdW heterostructures can open new avenues for the exploration of physical phenomena of novel quantum-confined nanostructures and devices. Therefore, the fundamental basis and recent findings in vertical vdW heterostructures composed of MXenes as a primary component and other 2D materials as secondary components are examined. Their robust designs and synthesis approaches that can push the boundaries of light-harvesting, transition, and utilization are discussed, since MXenes provide a unique playground for pursuing an extraordinary optical response or unusual light conversion features/functionalities. The recent findings are finally summarized, and a perspective for the future development of next-generation vdW multifunctional materials enriched by MXenes is provided.

A Multitechnique Study of C2H4 Adsorption on Fe3O4(001)

The adsorption/desorption of ethene (C2H4), also commonly known as ethylene, on Fe3O4(001) was studied under ultrahigh vacuum conditions using temperature-programmed desorption (TPD), scanning tunneling microscopy, X-ray photoelectron spectroscopy, and density functional theory (DFT)-based computations. To interpret the TPD data, we have employed a new analysis method based on equilibrium thermodynamics. C2H4 adsorbs intact at all coverages and interacts most strongly with surface defects such as antiphase domain boundaries and Fe adatoms. On the regular surface, C2H4 binds atop surface Fe sites up to a coverage of 2 molecules per (√2 × √2)R45° unit cell, with every second Fe occupied. A desorption energy of 0.36 eV is determined by analysis of the TPD spectra at this coverage, which is approximately 0.1-0.2 eV lower than the value calculated by DFT + U with van der Waals corrections. Additional molecules are accommodated in between the Fe rows. These are stabilized by attractive interactions with the molecules adsorbed at Fe sites. The total capacity of the surface for C2H4 adsorption is found to be close to 4 molecules per (√2 × √2)R45° unit cell.

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r2SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew-Burke-Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r2SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

Accurate GW frontier orbital energies of 134 kilo molecules

HOMO and LUMO energies are critical molecular properties that typically require high accuracy computations for practical applicability. Until now, a comprehensive dataset containing sufficiently accurate HOMO and LUMO energies has been unavailable. In this study, we introduce a new dataset of HOMO/LUMO energies for QM9 compounds, calculated using the GW method. The GW method offers adequate HOMO/LUMO prediction accuracy for diverse applications, exhibiting mean unsigned errors of 100 meV in the GW100 benchmark dataset. This database may serve as a benchmark of HOMO/LUMO prediction, delta-learning, and transfer learning, particularly for larger molecules where GW is the most accurate but still numerically feasible method. We anticipate that this dataset will enable the development of more accurate machine learning models for predicting molecular properties.

The application of the SCAN density functional to color centers in diamond

Detailed characterization of deep-level color centers requires understanding their electronic and atomic structure, which is most commonly investigated utilizing the Kohn-Sham density functional theory. Standard semilocal functionals based on the generalized gradient approximation (GGA) are inclined toward an imprecise quantitative description of defects' electronic structure. Hybrid functionals provide an improved prediction of electronic properties, albeit at a much higher computational cost. In this work, we test the newly developed Strongly Constrained and Appropriately Normed (SCAN) family of meta-GGA density functionals for selected color centers in diamond. In particular, we study nitrogen-, silicon-, germanium-, and tin-vacancy centers that have been recently investigated for their use in quantum technological applications. We show that SCAN and its derivatives, the rSCAN and r2SCAN functionals, significantly improve the calculated energies of optical transitions within the delta-self-consistent-field approach, almost reaching the accuracy of the hybrid Heyd-Scuseria-Ernzerhof (HSE) functional. In the case of the NV- center, we also show that the SCAN family of functionals improves the description of the adiabatic potential energy surfaces compared to both GGA and hybrid functionals, improving calculated luminescence lineshapes. As a result of these findings, we recommend using the SCAN family of functionals as a promising alternative for studying color centers in solids.

Characterization of the Hydrogen-Bond Network in High-Pressure Water by Deep Potential Molecular Dynamics

The hydrogen-bond (H-bond) network of high-pressure water is investigated by neural-network-based molecular dynamics (MD) simulations with first-principles accuracy. The static structure factors (SSFs) of water at three densities, i.e., 1, 1.115, and 1.24 g/cm3, are directly evaluated from 512 water MD trajectories, which are in quantitative agreement with the experiments. We propose a new method to decompose the computed SSF and identify the changes in the SSF with respect to the changes in H-bond structures. We find that a larger water density results in a higher probability for one or two non-H-bonded water molecules to be inserted into the inner shell, explaining the changes in the tetrahedrality of water under pressure. We predict that the structure of the accepting end of water molecules is more easily influenced by the pressure than by the donating end. Our work sheds new light on explaining the SSF and H-bond properties in related fields.

Implementation of the meta-GGA exchange-correlation functional in numerical atomic orbital basis: With systematic testing on SCAN, rSCAN, and r2SCAN functionals

Kohn-Sham density functional theory (DFT) is nowadays widely used for electronic structure theory simulations, and the accuracy and efficiency of DFT rely on approximations of the exchange-correlation functional. By including the kinetic energy density τ, the meta-generalized-gradient approximation (meta-GGA) family of functionals achieves better accuracy and flexibility while retaining the efficiency of semi-local functionals. For example, the strongly constrained and appropriately normed (SCAN) meta-GGA functional has been proven to yield accurate results for solid and molecular systems. We implement meta-GGA functionals with both numerical atomic orbitals and plane wave bases in the ABACUS package. Apart from the exchange-correlation potential, we also discuss the evaluation of force and stress. To validate our implementation, we perform finite-difference tests and convergence tests with the SCAN, rSCAN, and r2SCAN meta-GGA functionals. We further test water hexamers, weakly interacting molecules from the S22 dataset, as well as 13 semiconductors using the three functionals. The results show satisfactory agreement with previous calculations and available experimental values.

Why Dissolving Salt in Water Decreases Its Dielectric Permittivity

The dielectric permittivity of salt water decreases on dissolving more salt. For nearly a century, this phenomenon has been explained by invoking saturation in the dielectric response of the solvent water molecules. Herein, we employ an advanced deep neural network (DNN), built using data from density functional theory, to study the dielectric permittivity of sodium chloride solutions. Notably, the decrease in the dielectric permittivity as a function of concentration, computed using the DNN approach, agrees well with experiments. Detailed analysis of the computations reveals that the dominant effect, caused by the intrusion of ionic hydration shells into the solvent hydrogen-bond network, is the disruption of dipolar correlations among water molecules. Accordingly, the observed decrease in the dielectric permittivity is mostly due to increasing suppression of the collective response of solvent waters.

Radicals in aqueous solution: assessment of density-corrected SCAN functional

We study self-interaction effects in solvated and strongly-correlated cationic molecular clusters, with a focus on the solvated hydroxyl radical. To address the self-interaction issue, we apply the DC-r2SCAN method, with the auxiliary density matrix approach. Validating our method through simulations of bulk liquid water, we demonstrate that DC-r2SCAN maintains the structural accuracy of r2SCAN while effectively addressing spin density localization issues. Extending our analysis to solvated cationic molecular clusters, we find that the hemibonded motif in the [CH3S∴CH3SH]+ cluster is disrupted in the DC-r2SCAN simulation, in contrast to r2SCAN that preserves the (three-electron-two-center)-bonded motif. Similarly, for the [SH∴SH2]+ cluster, r2SCAN restores the hemibonded motif through spin leakage, while DC-r2SCAN predicts a weaker hemibond formation influenced by solvent-solute interactions. Our findings demonstrate the potential of DC-r2SCAN combined with the auxiliary density matrix method to improve electronic structure calculations, providing insights into the properties of solvated cationic molecular clusters. This work contributes to the advancement of self-interaction corrected electronic structure theory and offers a computational framework for modeling condensed phase systems with intricate correlation effects.

Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point

The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164-168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. We explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

Accurate Simulations of Scanning Tunneling Microscope: Both Tip and Substrate States Matter

Scanning tunneling microscope (STM) provides an atomic-scale characterization tool. To this end, high-resolution measurements and accurate simulations must closely cooperate. Emerging experimental techniques, e.g., substrate spacers and tip modifications, suppress metallic couplings and improve the resolution. On the other hand, development of STM simulation methods was inactive in the past decade. Conventional simulations focus on the electronic structure of the substrate, often overlooking detailed descriptions of the tip states. Meanwhile, the overwhelming usage of periodic boundary conditions ensures effective simulations of only neutral systems. In this Perspective, we highlight the recent progress that takes the effects of both tip and substrate into account under either Tersoff-Hamann or Bardeen's approximation, which provides an accurate analysis of measured high-resolution STM results, uncovers underlying concepts, and rationally designs experimental protocols for important chemical systems. We hope this Perspective will stimulate broad interest in advanced STM simulations, highlighting the way forward for STM investigations that involve complex geometrical and electronic structures.

High-Throughput Condensed-Phase Hybrid Density Functional Theory for Large-Scale Finite-Gap Systems: The SeA Approach

High-throughput electronic structure calculations (often performed using density functional theory (DFT)) play a central role in screening existing and novel materials, sampling potential energy surfaces, and generating data for machine learning applications. By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal DFT and furnish a more accurate description of the underlying electronic structure, albeit at a computational cost that often prohibits such high-throughput applications. To address this challenge, we have constructed a robust, accurate, and computationally efficient framework for high-throughput condensed-phase hybrid DFT and implemented this approach in the PWSCF module of Quantum ESPRESSO (QE). The resulting SeA approach (SeA = SCDM + exx + ACE) combines and seamlessly integrates: (i) the selected columns of the density matrix method (SCDM, a robust noniterative orbital localization scheme that sidesteps system-dependent optimization protocols), (ii) a recently extended version of exx (a black-box linear-scaling EXX algorithm that exploits sparsity between localized orbitals in real space when evaluating the action of the standard/full-rank V^xx operator), and (iii) adaptively compressed exchange (ACE, a low-rank V^xx approximation). In doing so, SeA harnesses three levels of computational savings: pair selection and domain truncation from SCDM + exx (which only considers spatially overlapping orbitals on orbital-pair-specific and system-size-independent domains) and low-rank V^xx approximation from ACE (which reduces the number of calls to SCDM + exx during the self-consistent field (SCF) procedure). Across a diverse set of 200 nonequilibrium (H₂O)₆₄ configurations (with densities spanning 0.4-1.7 g/cm³), SeA provides a 1-2 order-of-magnitude speedup in the overall time-to-solution, i.e., ≈8-26× compared to the convolution-based PWSCF(ACE) implementation in QE and ≈78-247× compared to the conventional PWSCF(Full) approach, and yields energies, ionic forces, and other properties with high fidelity. As a proof-of-principle high-throughput application, we trained a deep neural network (DNN) potential for ambient liquid water at the hybrid DFT level using SeA via an actively learned data set with ≈8,700 (H₂O)₆₄ configurations. Using an out-of-sample set of (H₂O)₅₁₂ configurations (at nonambient conditions), we confirmed the accuracy of this SeA-trained potential and showcased the capabilities of SeA by computing the ground-truth ionic forces in this challenging system containing >1,500 atoms.

Study of polymerization of high-pressure nitrogen by ab initio molecular dynamics

We study properties of nitrogen at high pressure and temperature (100-120 GPa, 2000-3000 K) where molecular and polymeric phases compete both in solid and liquid phase. We employ ab initio MD simulations with the SCAN functional and study the pressure-induced polymerization in liquid nitrogen for system sizes up to 288 atoms in order to reduce finite-size effects. The transition is studied upon both compression and decompression, and at 3000 K, it is found to take place between 110 and 115 GPa, coming close to experimental data. We also simulate the molecular crystalline phase close to the melting line and analyze its structure. We show that the molecular crystal in this regime is highly disordered, in particular, due to pronounced orientational and also translational disorder of the molecules. Its short-range order and vibrational density of states are very close to those of the molecular liquid revealing that the system likely represents a plastic crystal with high entropy.

Probing Confinement Effects on the Infrared Spectra of Water with Deep Potential Molecular Dynamics Simulations

The hydrogen-bond network of confined water is expected to deviate from that of the bulk liquid, yet probing these deviations remains a significant challenge. In this work, we combine large-scale molecular dynamics simulations with machine learning potential derived from first-principles calculations to examine the hydrogen bonding of water confined in carbon nanotubes (CNTs). We computed and compared the infrared spectrum (IR) of confined water to existing experiments to elucidate confinement effects. For CNTs with diameters >1.2 nm, we find that confinement imposes a monotonic effect on the hydrogen-bond network and on the IR spectrum of water. In contrast, confinement below 1.2 nm CNT diameter affects the water structure in a complex fashion, leading to a strong directional dependence of hydrogen bonding that varies nonlinearly with the CNT diameter. When integrated with existing IR measurements, our simulations provide a new interpretation for the IR spectrum of water confined in CNTs, pointing to previously unreported aspects of hydrogen bonding in this system. This work also offers a general platform for simulating water in CNTs with quantum accuracy on time and length scales beyond the reach of conventional first-principles approaches.

Versatile Interplay of Chalcogenide and Dichalcogenide Anions in the Thiovanadate Ba7S(VS3O)2(S2)3 and Its Selenide Derivatives: Elaboration and DFT Meta-GGA Study

Oxychalcogenides are emerging as promising alternative candidates for a variety of applications including for energy. Only few phases among them show the presence of Q-Q bonds (Q = chalcogenide anion) while they drastically alter the electronic structure and allow further structural flexibility. Four original oxy(poly)chalcogenide compounds in the system Ba-V-Q-O (Q = S, Se) were synthesized, characterized, and studied using density functional theory (DFT). The new structure type found for Ba₇V₂O₂S₁₃, which can be written as Ba₇S(VS₃O)₂(S₂)₃, was substituted to yield three selenide derivatives Ba₇V₂O₂S_9.304Se_3.696, Ba₇V₂O₂S_7.15Se_5.85, and Ba₇V₂O₂S_6.85Se_6.15. They represent original multiple-anion lattices and first members in the system Ba-V-Se-S-O. They exhibit in the first layer heteroleptic tetrahedra V⁵⁺S₃O and isolated Q^2- anions and in the second layer dichalcogenide pairs (Q₂)^2- with Q = S or Se. Selenide derivatives were attempted by targeting the selective substitution of isolated Q^2- or (Q₂)^2- (in distinct layers) or both by selenide, but it systematically led to concomitant and partial substitution of both sites. A DFT meta-GGA study showed that selective substitution yields local constraints due to rigid VO₃S and pairs. Experimentally, incorporation of selenide in both layers avoids geometrical mismatch and constraints. In such systems, we show that the interplay between the O/S anionic ratio around V⁵⁺, together with the presence/nature of the dichalcogenides (Q₂)^2- and isolated Q^2-, impacts in unique manners the band gap and provides a rich background to tune the band gap and the symmetry.

Statistical analysis of the performance of a variety of first-principles schemes for accurate prediction of binary semiconductor band gaps

Standard density functional theory (DFT) approximations tend to strongly underestimate band gaps, while the more accurate GW and hybrid functionals are much more computationally demanding and unsuitable for high-throughput screening. In this work, we have performed an extensive benchmark of several approximations with different computational complexity [G0W0@PBEsol, HSE06, PBEsol, modified Becke-Johnson potential (mBJ), DFT-1/2, and ACBN0] to evaluate and compare their performance in predicting the bandgap of semiconductors. The benchmark is based on 114 binary semiconductors of different compositions and crystal structures, for about half of which experimental band gaps are known. Surprisingly, we find that, compared with G0W0@PBEsol, which exhibits a noticeable underestimation of the band gaps by about 14%, the much computationally cheaper pseudohybrid ACBN0 functional shows a competitive performance in reproducing the experimental data. The mBJ functional also performs well relative to the experiment, even slightly better than G0W0@PBEsol in terms of mean absolute (percentage) error. The HSE06 and DFT-1/2 schemes perform overall worse than ACBN0 and mBJ schemes but much better than PBEsol. Comparing the calculated band gaps on the whole dataset (including the samples with no experimental bandgap), we find that HSE06 and mBJ have excellent agreement with respect to the reference G0W0@PBEsol band gaps. The linear and monotonic correlations between the selected theoretical schemes and experiment are analyzed in terms of the Pearson and Kendall rank coefficients. Our findings strongly suggest the ACBN0 and mBJ methods as very efficient replacements for the costly G0W0 scheme in high-throughput screening of the semiconductor band gaps.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.

Eco-friendly inorganic molecular novel antiperovskites for light-emitting application

The development of perovskite light-emitting diodes (PeLEDs) has progressed rapidly over the past several years, with high external quantum efficiencies exceeding 20%. However, the deployment of PeLEDs in commercial devices still faces severe challenges, such as environmental pollution, instability and low photoluminescence quantum yields (PLQYs). In this work, we perform high-throughput calculations to exhaustively search the unexplored and eco-friendly novel antiperovskite space (formula: X₃B[MN₄], with octahedron [BX₆] and tetrahedron [MN₄]). The novel antiperovskites have a unique structure whereby a tetrahedron can be embedded into an octahedral skeleton as a light-emitting center causing a space confinement effect, leading to the characteristics of a low-dimensional electronic structure, which then makes these materials potential light-emitting material candidates with a high PLQY and light-emitting stability. Under the guidance of newly derived tolerance, octahedral, and tetrahedral factors, 266 stable candidates are successfully screened out from 6320 compounds. Moreover, the antiperovskite materials Ba₃I_0.5F_0.5(SbS₄), Ca₃O(SnO₄), Ba₃F_0.5I_0.5(InSe₄), Ba₃O_0.5S_0.5(ZrS₄), Ca₃O(TiO₄), and Rb₃Cl_0.5I_0.5(ZnI₄) possess an appropriate bandgap, thermodynamic and kinetic stability, and excellent electronic and optical properties, making them promising light-emitting materials.

Addressing the Conflict between Mobility and Stability in Oxide Thin-film Transistors

Amorphous oxide semiconductor thin-film transistors (AOS TFTs) are ever-increasingly utilized in displays. However, to bring high mobility and excellent stability together is a daunting challenge. Here, the carrier transport/relaxation bilayer stacked AOS TFTs are investigated to solve the mobility-stability conflict. The charge transport layer (CTL) is made of amorphous In-rich InSnZnO, which favors big average effective coordination number for all cations and more edge-shared structures for better charge transport. Praseodymium-doped InSnZnO is used as the charge relaxation layer (CRL), which substantially shortens the photoelectron lifetime as revealed by femtosecond transient absorption spectroscopy. The CTL and CRL with the thickness suitable for industrial production respectively afford minute potential barrier fluctuation for charge transport and fast relaxation for photo-generated carriers, resulting in transistors with an ultrahigh mobility (75.5 cm² V^-1 s^-1 ) and small negative-bias-illumination-stress/positive-bias-temperature-stress voltage shifts (-1.64/0.76 V). The design concept provides a promising route to address the mobility-stability conflict for high-end displays.

Highly degenerate 2D ferroelectricity in pore decorated covalent/metal organic frameworks

Two-dimensional (2D) ferroelectricity, a fundamental concept in low-dimensional physics, serves as the basis of non-volatile information storage and various electronic devices. Conventional 2D ferroelectric (FE) materials are usually two-fold degenerate, meaning that they can only store two logical states. In order to break such limitation, a new concept of highly degenerate ferroelectricity with multiple FE states (more than 2) coexisting in a single 2D material is proposed. This is obtained through the asymmetrical decoration of porous covalent/metal organic frameworks (COFs/MOFs). Using first-principles calculations and Monte Carlo (MC) simulations, Li-decorated 2D Cr(pyz)₂ is systematically explored as a prototype of highly degenerate 2D FE materials. We show that 2D FE Li_0.5Cr(pyz)₂ and LiCr(pyz)₂ are four-fold and eight-fold degenerate, respectively, with sizable spontaneous electric polarization that can be switched across low transition barriers. In particular, the coupling between neighbouring electric dipoles in LiCr(pyz)₂ induces novel ferroelectricity-controlled ferroelastic transition and direction-controllable hole transport channels. Moreover, three-fold and six-fold degenerate ferroelectricity is also demonstrated in P-decorated g-C₃N₄ and Ru-decorated C₂N, respectively. Our work presents a general route to obtain highly degenerate 2D ferroelectricity, which goes beyond the two-state paradigm of traditional 2D FE materials and substantially broadens the applications of 2D FE compounds.

First-Principles Calculation of MoO2 and MoO3 Electronic and Optical Properties Compared with Experimental Data

MoO₃ and MoO₂ systems have attracted particular attention for many widespread applications thanks to their electronic and optical peculiarities; from the crystallographic point of view, MoO₃ adopts a thermodynamically stable orthorhombic phase (α-MoO₃) belonging to the space group Pbmn, while MoO₂ assumes a monoclinic arrangement characterized by space group P2₁/c. In the present paper, we investigated the electronic and optical properties of both MoO₃ and MoO₂ by using Density Functional Theory calculations, in particular, the Meta Generalized Gradient Approximation (MGGA) SCAN functional together with the PseudoDojo pseudopotential, which were used for the first time to obtain a deeper insight into the nature of different Mo-O bonds in these materials. The calculated density of states, the band gap, and the band structure were confirmed and validated by comparison with already available experimental results, while the optical properties were validated by recording optical spectra. Furthermore, the calculated band-gap energy value for the orthorhombic MoO₃ showed the best match to the experimental value reported in the literature. All these findings suggest that the newly proposed theoretical techniques reproduce the experimental evidence of both MoO₂ and MoO₃ systems with high accuracy.

Extending density functional theory with near chemical accuracy beyond pure water

Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r²SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.

Immobilization of As(III) by gibbsite and catalytic oxidation to As(V): Profound impacts of doping and unraveling of associated mechanisms

As(III) is highly toxic, and its adsorption and oxidation to As(V) by minerals represent two effective approaches to remediate As(III)-contaminated sites. Gibbsite, one of the most abundant natural minerals, shows decent adsorption for As(III), and in this study, mechanisms of As(III) immobilization and oxidation by gibbsite with different dopants (M = Fe(III), Mn(III), Mn(IV)) are addressed by periodic DFT calculations. Influences of Fe(III) content and Mn oxidation state are also inspected. Although a majority remain structurally similar to those of pristine gibbsite, new adsorption configurations emerge due to doping: Inner-sphere complexes with M - As bonds for all doping, bidentate binuclear complexes for double Fe(III) doping, and physisorption with weak O_Mn-As interactions for Mn(IV) doping. As(III) adsorption affinities are significantly altered by doping and rely on dopants, while inner-sphere complexes with M-O_As bonds are always lowest-energy except doping Mn(III) that prefers trigonal bipyramidal coordination and impedes As(III) chemisorption. Doping causes strong M-3d and O_As-2p orbital interactions that facilitate As(III) adsorption whereas disappear for pristine gibbsite. Double Fe(III)- and Mn(IV)-doped gibbsite materials are effective for As(III) oxidation to As(V), and mechanisms differ significantly although all are characterized by dual electron transfers. Activation barriers for the most favorable reaction paths amount to 1.02 and 1.26-1.31 eV, respectively. Physisorbed and outer-sphere As(III) complexes exhibit comparable reactivities as chemisorbed complexes that become focus of literature reports, and may also be involved during interfacial and environmental reactions. Results rationalize experimental observations available, and provide significantly new insights that conduce to manage As-associated pollution and design efficient As(III) scavengers and oxidation catalysts.

Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO₂ . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO₃  ⋅ 3H₂ O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO₃ )OH ⋅ 2H₂ O]. The resolution of this controversy is central to understanding this material's thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance ¹³ C{¹ H} (REDOR), ¹³ C-¹ H distances are determined with precision, and the combination of ¹³ C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO₃  ⋅ 3H₂ O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO₃ . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

Solar photodegradation of Rhodamine B dye by Cu2O/TiO2 heterostructure: experimental and computational studies of degradation and toxicity

CONTEXT AND RESULTS

In this study, the heterojunction Cu₂O/TiO₂ is used for the degradation of a cationic dye, Rhodamine B, under solar light irradiation over a wide pH range. The physical and optical properties of both semiconductors Cu₂O and TiO₂ are correlated with the photo-electrochemical characterization to establish the energy diagram of the heterojunction Cu₂O/TiO₂. X-ray diffraction, UV-visible, SEM, EDX, and BET analyses are conducted for both photocatalysts. The band gap (E_g) of 3.26 eV is obtained for TiO₂ with an indirect optical transition. In the case of Cu₂O, the transition is directly allowed at 2.05 eV. According to the BET analysis, the specific surface area of TiO₂ particles is higher (82.65 m² g^-1) than that of Cu₂O (29.81 m² g^-1). The flat band potentials, determined from the Mott-Schottky plots, are 0.3 and - 0.32 V_SCE for TiO₂ and Cu₂O, respectively. The photocatalytic activity is directly proportional to the mass ratio, and the best result is obtained for the mass ratio 1:1 of Cu₂O/TiO₂.

COMPUTATIONAL AND THEORETICAL TECHNIQUES

Furthermore, a theoretical study is conducted by using density functional theory to optimize the structure, reactivity sites of the RhB molecule, and physical parameters like the energy of the frontier molecular orbitals and electronegativity and to predict the proposed mechanism of RhB degradation as well as its intermediates. Also, molecular dynamics simulation is used to determine the adsorption behavior of RhB on TiO₂ (101) and Cu₂O (111) surfaces. The ecotoxicity evaluation showed that degradation products have significantly lower acute toxicity than RhB.

New Insights into the Volume Isotope Effect of Ice Ih from Polarizable Many-Body Potentials

The anomalous volume isotope effect (VIE) of ice Ih is calculated and analyzed based on the quasi-harmonic approximation to account for nuclear quantum effects in the Helmholtz free energy. While a lot of recently developed polarizable many-body potential functions give a normal VIE contrary to experimental results, we find that one of them, MB-pol, yields the anomalous VIE in good agreement with the most recent high-resolution neutron diffraction measurements─better than DFT calculations. The short-range three-body terms in the MB-pol function, which are fitted to CCSD(T) calculations, are found to have a surprisingly large influence. A vibrational mode group decomposition of the zero-point pressure together with a hitherto unconsidered benchmark value for the intramolecular stretching modes of H₂O ice Ih obtained from Raman spectroscopy data unveils the reason for the VIE: a delicate competition between the latter and the librations.

Prediction of boridenes as high-performance anodes for alkaline metal and alkaline Earth metal ion batteries

We conducted a comprehensive density functional theory investigation using the r²SCAN-rVV10 functional on the structural stability and electrochemical properties of boridenes for their use as anode materials in rechargeable alkaline (earth) metal-ion batteries (Li⁺, Na⁺, K⁺, Mg²⁺ and Ca²⁺). According to first-principles molecular dynamics simulations and reaction thermodynamic calculations, Mo_4/3B₂(OH)₂ and Mo_4/3B₂F₂ are unstable in the presence of alkaline (earth) metal ions due to the surface-conversion reactions between the surface terminations and adsorbates. Meanwhile, the bare Mo_4/3B₂ and Mo_4/3B₂O₂ monolayers not only can accommodate alkaline (earth) metal ions, but also form stable multi-layer adsorption structures for most of the studied metal ions (Li⁺, Na⁺, K⁺, Mg²⁺ and Ca²⁺). The predicted gravimetric capacities of the bare Mo_4/3B₂ monolayer (Mo_4/3B₂O₂) are 625.9 mA h g^-1 (357.3 mA h g^-1), 247.20 mA h g^-1 (392.1 mA h g^-1), 101.8 mA h g^-1 (206.4 mA h g^-1), 667.0 mA h g^-1, and 413.0 mA h g^-1 (485.4 mA h g^-1) for Li⁺, Na⁺, K⁺, Mg²⁺ and Ca²⁺ ions, respectively. The bare Mo_4/3B₂ exhibits lower onset charging open circuit voltages for alkaline (earth) metal ions than that of Mo_4/3B₂O₂. The diffusivities of the metal ions were revealed to be high on the boridene monolayer especially for the outer fully stable adsorption layers, where the migration energy barriers were found to be less than 0.10 eV. Similar to that of MXenes, the negative electron cloud (NEC) also plays a vital role in stabilizing the observed multi-layer adsorption structures for various metal ions on either the bare Mo_4/3B₂ or Mo_4/3B₂O₂ monolayer.

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Massive Monte Carlo simulations-guided interpretable learning of two-dimensional Curie temperature

Monte Carlo (MC) simulation of the classical Heisenberg model has become the de facto tool to estimate the Curie temperature (T _C) of two-dimensional (2D) magnets. As an alternative, here we develop data-driven models for the five most common crystal types, considering the isotropic and anisotropic exchange of up to four nearest neighbors and the single-ion anisotropy. We sample the 20-dimensional Heisenberg spin Hamiltonian and conceive a bisection-based MC technique to simulate a quarter of a million materials for training deep neural networks, which yield testing R ² scores of nearly 0.99. Since 2D magnetism has a natural tendency toward low T _C, learning-from-data is combined with data-from-learning to ensure a nearly uniform final data distribution over a wide range of T _C (10-1,000 K). Global and local analysis of the features confirms the models' interpretability. We also demonstrate that the T _C can be accurately estimated by a purely first-principles-based approach, free from any empirical corrections.

LiXO2(X = Co, Rh, Ir) and solar light photocatalytic water splitting for hydrogen generation

Alkali metal transition oxide LiCoO₂ has been successfully commercialized as a lithium-ion battery material, and some attention is paid to its homologous derivatives LiRhO₂ and LiIrO₂. However, the photocatalytic properties have not been explored yet for these compounds. Using the first-principles calculations, we carry out investigations on the electronic properties, light absorption, and mobility to understand the feasibility of LiXO₂(X = Co, Rh, Ir) for solar light photocatalytic hydrogen generation from water-splitting. The results show that the band edges of LiCoO₂ and LiRhO₂ meet the redox potential requirements of the water-splitting hydrogen evolution reaction. In addition, the enhanced absorptions of LiXO₂(X = Co, Rh, Ir) in the visible light range imply that they could well respond to solar light, while the significant difference in the mobilities of electrons or holes can strengthen spatial charge separation of the photoexcited electron-hole pairs. The solar-energy-to-hydrogen conversion efficiencies of LiCoO₂ and LiRhO₂ can reach 11.2% and 15.5%, respectively. The results support LiCoO₂ and LiRhO₂ as promising candidates for visible-light photocatalytic hydrogen production from water-splitting.

Integrating ferromagnetism and ferroelectricity in an iron chalcogenide monolayer: a first-principles study

Two-dimensional (2D) ferro-type materials have received great attention owing to the remarkable polarization effect in optoelectronics and spintronics. Using the first-principles method, the coupling between ferromagnetism and ferroelectricity is investigated in a multiferroic Janus 1T-FeSSe monolayer, which has a strong Stoner ferromagnetic ground state. The magnetic anisotropy energy (MAE) is apparently impacted by the out-of-plane asymmetry donated ferroelectricity, which is reflected by the asymmetry of the Z-MAE image. The easy magnetization axis of Janus FeSSe is the +y axis with a large MAE of 0.59 meV, rooting in unpaired d electrons of Fe atoms. The transformation of band splitting and Fermi surface can be effectively engineered by different magnetic polarization directions. The ferromagnetic (FM) coupling of the FeSSe monolayer is very robust under external strain within the range of -6% to 6%, while the strength of magnetic moment of Fe atoms and polarization are easily strain-engineered, the intrinsic mechanism of which can be elaborated by the GKA rules that depend on angles and distances. This multiferroic FeSSe monolayer provides a new platform for exploring the coupling of 2D ferromagnetism and ferroelectricity and designing low-dimensional multiferroic electronics.

Polaron formation and transport in Bi2WO6 studied by DFT+U and hybrid PBE0 functional approaches

Bi₂WO₆ (BWO) is considered as a promising material for photocatalytic water splitting. Its unique layered structure leads the charge separation, and transport is different from other materials. However, the charge transport mechanisms in BWO are not well understood. In this work, we investigated polaron formation and transport in BWO using the DFT+U and hybrid PBE0 functional approaches. We found that the electron will form 2-dimensional (2D)-shaped polarons among W sites in the ab plane of BWO with approximately 55% polaron density state on the central W site. This type of polaron is similar to the electron polarons in WO₃. For other W-based materials, the electrons may also form a 2D-shaped polaron. We found that the W 6s orbital plays an important role in these 2D-shaped electron polarons. The calculated mobility of electron polarons in BWO was consistent with experimental findings. For the hole state, it could form a small hole polaron on the O site with O 2p in character. However, it will not form a polaron on the Bi site, which is quite different from BiVO₄. This study provides insight for understanding polaron formation and transport in materials with W and Bi ions. It also provides understanding regarding charge separation and transport for materials with layered structures.

Catalytic formation of oxalic acid on the partially oxidised greigite Fe3S4(001) surface

Greigite (Fe₃S₄), with its ferredoxin-like 4Fe-4S redox centres, is a naturally occurring mineral capable of acting as a catalyst in the conversion of carbon dioxide (CO₂) into low molecular-weight organic acids (LMWOAs), which are of paramount significance in several soil and plant processes as well as in the chemical industry. In this paper, we report the reaction between CO₂ and water (H₂O) to form oxalic acid (H₂C₂O₄) on the partially oxidised greigite Fe₃S₄(001) surface by means of spin-polarised density functional theory calculations with on-site Coulomb corrections and long-range dispersion interactions (DFT+U-D2). We have calculated the bulk phase of Fe₃S₄ and the two reconstructed Tasker type 3 terminations of its (001) surface, whose properties are in good agreement with available experimental data. We have obtained the relevant phase diagram, showing that the Fe₃S₄(001) surface becomes 62.5% partially oxidised, by replacing S by O atoms, in the presence of water at the typical conditions of calcination [Mitchell et al. Faraday Discuss. 2021, 230, 30-51]. The adsorption and co-adsorption of the reactants on the partially oxidised Fe₃S₄(001) surface are exothermic processes. We have considered three mechanistic pathways to explain the formation of H₂C₂O₄, showing that the coupling of the C-C bond and second protonation are the elementary steps with the largest energy penalty. Our calculations suggest that the partially oxidised Fe₃S₄(001) surface is a mineral phase that can catalyse the formation of H₂C₂O₄ under favourable conditions, which has important implications for natural ecosystems and is a process that can be harnessed for the industrial manufacture of this organic acid.

Surface chemistry on a polarizable surface: Coupling of CO with KTaO3(001)

Polarizable materials attract attention in catalysis because they have a free parameter for tuning chemical reactivity. Their surfaces entangle the dielectric polarization with surface polarity, excess charge, and orbital hybridization. How this affects individual adsorbed molecules is shown for the incipient ferroelectric perovskite KTaO₃. This intrinsically polar material cleaves along (001) into KO- and TaO₂-terminated surface domains. At TaO₂ terraces, the polarity-compensating excess electrons form a two-dimensional electron gas and can also localize by coupling to ferroelectric distortions. TaO₂ terraces host two distinct types of CO molecules, adsorbed at equivalent lattice sites but charged differently as seen in atomic force microscopy/scanning tunneling microscopy. Temperature-programmed desorption shows substantially stronger binding of the charged CO; in density functional theory calculations, the excess charge favors a bipolaronic configuration coupled to the CO. These results pinpoint how adsorption states couple to ferroelectric polarization.

Ab Initio Study of Structure and Transport Properties of Warm Dense Nitric Oxide

The structure, equation of state and transport properties of warm dense nitric oxide (NO) were investigated in wide density and temperature ranges by ab initio molecular dynamics simulations. Both the Perdew–Burke–Ernzerhof (PBE) and the strongly constrained and appropriately normed functional with revised Vydrov–van Voorhis nonlocal correlation (SCAN−rVV10) functionals were used in the simulations, and the pressures predicted by the SCAN−rVV10 functional were found to be systematically lower than those predicted using PBE and experimental data along the shock Hugoniot curve. Along the Hugoniot curve, as density increased, we found that the system transformed towards a mixture of atomic nitrogen and oxygen liquids with molecular NO that remained present up to the highest densities explored. The electrical conductivity along Hugoniot indicated that nonmetal to metal transition had taken place. We also calculated the electrical and thermal conductivities of nitric oxide in the warm dense matter regime, and used them to compute the Lorentz number. In addition, we also report the electronic density of states.

Localized Resolution of Identity Approach to the Analytical Gradients of Random-Phase Approximation Ground-State Energy: Algorithm and Benchmarks

We develop and implement a formalism which enables calculating the analytical gradients of particle-hole random-phase approximation (RPA) ground-state energy with respect to the atomic positions within the atomic orbital basis set framework. Our approach is based on a localized resolution of identity (LRI) approximation for evaluating the two-electron Coulomb integrals and their derivatives, and the density functional perturbation theory for computing the first-order derivatives of the Kohn-Sham (KS) orbitals and orbital energies. Our implementation allows one to relax molecular structures at the RPA level using both Gaussian-type orbitals (GTOs) and numerical atomic orbitals (NAOs). Benchmark calculations against previous implementations show that our approach delivers adequate numerical precision, highlighting the usefulness of LRI in the context of RPA gradient evaluations. A careful assessment of the quality of RPA geometries for small molecules reveals that post-KS RPA systematically overestimates the bond lengths. We furthermore optimized the geometries of the four low-lying water hexamers-cage, prism, cyclic, and book isomers, and determined the energy hierarchy of these four isomers using RPA. The obtained RPA energy ordering is in good agreement with that yielded by the coupled cluster method with single, double and perturbative triple excitations, despite that the dissociation energies themselves are appreciably underestimated. The underestimation of the dissociation energies by RPA is well corrected by the renormalized single excitation correction.

Effects of Ga Substitution on the Local Structure of Na2Zn2TeO6

In the work presented here, we prepared Ga-substituted NZTO (Na_2–xZn_2–xGa_xTeO₆, x = 0.00, 0.05, 0.10, 0.15, 0.20) layered materials with a soft chemical, citric acid-based synthesis method and characterized these by means of X-ray diffraction (XRD), ²³Na and ¹²⁵Te NMR, and by density functional theory (DFT) modeling. The influence of randomly distributed Ga cations on the ¹²⁵Te NMR spectra confirms the successful synthesis. With DFT-based linear response computations, we show that the local distribution of Na ions in the two neighboring interlayers influences the ¹²⁵Te chemical shift, consistent with observations. DFT modeling suggests that some of the Na sites are rarely occupied in pure NZTO but become favorable upon Ga substitution. There are clear indications that Ga substitution gives an uneven distribution of Na ions in neighboring interlayers and that the Na structure in one layer affects the adjacent layers.

The influence of randomly distributed Ga cations on the 125Te NMR spectra confirms the successful synthesis.

Comment on "Pushing the frontiers of density functionals by solving the fractional electron problem"

Kirkpatrick et al. (Reports, 9 December 2021, p. 1385) trained a neural network-based DFT functional, DM21, on fractional-charge (FC) and fractional-spin (FS) systems, and they claim that it has outstanding accuracy for chemical systems exhibiting strong correlation. Here, we show that the ability of DM21 to generalize the behavior of such systems does not follow from the published results and requires revisiting.

Competing electronic states emerging on polar surfaces

Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO2 termination of KTaO3(001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions.

How Good Is the Density-Corrected SCAN Functional for Neutral and Ionic Aqueous Systems, and What Is So Right about the Hartree-Fock Density?

Density functional theory (DFT) is the most widely used electronic structure method, due to its simplicity and cost effectiveness. The accuracy of a DFT calculation depends not only on the choice of the density functional approximation (DFA) adopted but also on the electron density produced by the DFA. SCAN is a modern functional that satisfies all known constraints for meta-GGA functionals. The density-driven errors, defined as energy errors arising from errors of the self-consistent DFA electron density, can hinder SCAN from achieving chemical accuracy in some systems, including water. Density-corrected DFT (DC-DFT) can alleviate this shortcoming by adopting a more accurate electron density which, in most applications, is the electron density obtained at the Hartree-Fock level of theory due to its relatively low computational cost. In this work, we present extensive calculations aimed at determining the accuracy of the DC-SCAN functional for various aqueous systems. DC-SCAN (SCAN@HF) shows remarkable consistency in reproducing reference data obtained at the coupled cluster level of theory, with minimal loss of accuracy. Density-driven errors in the description of ionic aqueous clusters are thoroughly investigated. By comparison with the orbital-optimized CCD density in the water dimer, we find that the self-consistent SCAN density transfers a spurious fraction of an electron across the hydrogen bond to the hydrogen atom (H^*, covalently bound to the donor oxygen atom) from the acceptor (O_A) and donor (O_D) oxygen atoms, while HF makes a much smaller spurious transfer in the opposite direction, consistent with DC-SCAN (SCAN@HF) reduction of SCAN overbinding due to delocalization error. While LDA seems to be the conventional extreme of density delocalization error, and HF the conventional extreme of (usually much smaller) density localization error, these two densities do not quite yield the conventional range of density-driven error in energy differences. Finally, comparisons of the DC-SCAN results with those obtained with the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method show that DC-SCAN represents a more accurate approach to reducing density-driven errors in SCAN calculations of ionic aqueous clusters. While the HF density is superior to that of SCAN for noncompact water clusters, the opposite is true for the compact water molecule with exactly 10 electrons.

Dynamics in the O(2 × 1) adlayer on Ru(0001): bridging timescales from milliseconds to minutes by scanning tunneling microscopy

The dynamics within an O(2 × 1) adlayer on Ru(0001) is studied by density functional theory and high-speed scanning tunneling microscopy. Transition state theory proposes dynamic oxygen species in the reduced O(2 × 1) layer at room temperature. Collective diffusion processes can result in structural reorientations of characteristic stripe patterns. Spiral high-speed scanning tunneling microscopy measurements reveal this reorientation as a function of time in real space. Measurements, ranging over several minutes with constantly high frame rates of 20 Hz resolved the gradual reorientation. Moreover, reversible fast flipping events of stripe patterns are observed. These measurements relate the observations of long-term atomic rearrangements and their underlying fast processes captured within several tens of milliseconds.

The dynamics within an O(2 × 1) adlayer on Ru(0001) is studied by density functional theory and high-speed scanning tunneling microscopy.

Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water has been a challenging task from ab initio molecular dynamics (AIMD). The difficulty mainly comes from a lack of accurate and efficient exchange-correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly-constrained and appropriately normed (SCAN) meta-GGA functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange-correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice in describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. Our AIMD simulations indicate hydroxide diffuses slower than hydronium in water, which is consistent with the experiments.

Influence of Magnetic Moment on Single Atom Catalytic Activation Energy Barriers

Design of the molecular environment of single atom catalysts (SAC) is promising for achieving high catalytic activity without expensive and scarce platinum-group metals (PGM). We utilize a first principles approach to examine how the spin state of the SAC and reactants can affect catalytic energy barriers of V, Fe, Mo, and Ta on two different graphene defects with differing magnetic moments. Spin polarized projected density of states and climbing image nudged elastic band calculations demonstrate relatively lower activation energy barriers for systems with higher spin state asymmetry near the Fermi energy; CO oxidation on Ta and V SAC have decreases in activation barrier energies of 27% and 44%, respectively.

Graphic Abstract

Li5Sn, the Most Lithium-Rich Binary Stannide: A Combined Experimental and Computational Study

From reaction of excess lithium with tin, we isolate well-crystallized Li₅Sn and solve the crystal structure from single-crystal X-ray diffraction data. The orthorhombic structure (space group Cmcm) features the same coordination polyhedra around tin and lithium as previously predicted by electronic structure calculations for this composition, however differently arranged. An extensive ab initio analysis, including thermodynamic integration using Langevin dynamics in combination with a machine-learning potential (moment tensor potential), is conducted to understand the thermodynamic stability of this Cmcm Li₅Sn structure observed in our experiments. Among the 108 Li₅Sn structures systematically derived using the structure enumeration algorithm, including the experimental Cmcm structure and those obtained in previous ab initio studies, another new structure with the space group Immm is found to be energetically most stable at 0 K. This computationally discovered Immm structure is also found to be thermodynamically more stable than the Cmcm structure at finite temperatures, indicating that the Cmcm Li₅Sn structure observed in our experiments is favored likely due to kinetic reasons rather than thermodynamics.

Metal-Metal Bonding as an Electrode Design Principle in the Low-Strain Cluster Compound LiScMo3O8

Electrode materials for Li⁺-ion batteries require optimization along several disparate axes related to cost, performance, and sustainability. One of the important performance axes is the ability to retain structural integrity though cycles of charge/discharge. Metal-metal bonding is a distinct feature of some refractory metal oxides that has been largely underutilized in electrochemical energy storage, but that could potentially impact structural integrity. Here LiScMo₃O₈, a compound containing triangular clusters of metal-metal bonded Mo atoms, is studied as a potential anode material in Li⁺-ion batteries. Electrons inserted though lithiation are localized across rigid Mo₃ triangles (rather than on individual metal ions), resulting in minimal structural change as suggested by operando diffraction. The unusual chemical bonding allows this compound to be cycled with Mo atoms below a formally +4 valence state, resulting in an acceptable voltage regime that is appropriate for an anode material. Several characterization methods including potentiometric entropy measurements indicate two-phase regions, which are attributed through extensive first-principles modeling to Li⁺ ordering. This study of LiScMo₃O₈ provides valuable insights for design principles for structural motifs that stably and reversibly permit Li⁺ (de)insertion.

Computational analysis of the optical response of ZnSe with d-orbital defects

The doping of wide band-gap semiconducting ZnSe by transition metal (TM) atoms finds applications from mid-infrared lasing, sensing, photoelectrochemical cells, to nonlinear optics. Yet understanding the response of these materials at the atomic and electronic level is lacking, particularly in comparing a range of TM dopants, which were studied primarily by phenomenological crystal-field theory. In this work, to investigate bulk ZnSe singly doped with first-row TM atoms, specifically Ti through Cu, we applied a first-principles approach and crystal-field theory to explain the origin of the infrared absorption. We show that the use of an appropriate exchange-correlation functional and a HubbardUcorrection to account for electron correlation improved the determination of the electronic transitions in these systems. We outline an approach for the calculation of the crystal-field splitting from first-principles and find it useful in providing a measure of dopant effects, also in qualitative comparison to our experimental characterization for ZnSe doped with Fe, Cr, and Ni. Our calculated absorption spectra indicate absorption signatures in the mid-infrared range, while the absorption in the visible portion of the spectrum is attributed to the ZnSe host. Our calculations will potentially motivate further experimental exploration of TM-doped ZnSe. Finally, the methods used here provide a route towards computational high-throughput screening of TM dopants in III-V materials through a combination of the electronic band structure and crystal-field theory.