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

Dynamics & Spectroscopy with Neutrons-Recent Developments & Emerging Opportunities

This work provides an up-to-date overview of recent developments in neutron spectroscopic techniques and associated computational tools to interrogate the structural properties and dynamical behavior of complex and disordered materials, with a focus on those of a soft and polymeric nature. These have and continue to pave the way for new scientific opportunities simply thought unthinkable not so long ago, and have particularly benefited from advances in high-resolution, broadband techniques spanning energy transfers from the meV to the eV. Topical areas include the identification and robust assignment of low-energy modes underpinning functionality in soft solids and supramolecular frameworks, or the quantification in the laboratory of hitherto unexplored nuclear quantum effects dictating thermodynamic properties. In addition to novel classes of materials, we also discuss recent discoveries around water and its phase diagram, which continue to surprise us. All throughout, emphasis is placed on linking these ongoing and exciting experimental and computational developments to specific scientific questions in the context of the discovery of new materials for sustainable technologies.

Lattice dynamics in the conformational environment of multilayered hexagonal boron nitride (h-BN) results in peculiar infrared optical responses

Stacking mismatches in hexagonal boron nitride (h-BN) nanostructures affect their photonic, mechanical, and thermal properties. To access information about the stacked configuration of layered ensembles, highly sophisticated techniques like X-ray photoemission spectroscopy or electron microscopy are necessary. Here, instead, by taking advantage of the geometrical and chemical nature of h-BN, we show how simple structural models, based on shortened interplanar distances, can produce effective charge densities. Accounting these in the non-analytical part of the lattice dynamical description makes it possible to access information about the composition of differently stacked variants in experimental samples characterized by infrared spectroscopy. The results are obtained by density functional theory and confirmed by various functionals and pseudopotential approximations. Even though the method is shown using h-BN, the conclusions are more general and show how effective dielectric models can be considered as valuable theoretical pathways for the vibrational structure of any layered material.

Describing adsorption of benzene, thiophene, and xenon on coinage metals by using the Zaremba-Kohn theory-based model

Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, and rVV10) or empirical methods (DFT+D, DFT+vdW, and DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals (Cu, Ag, and Au) is a typical example of systems where vdW interaction is significant. However, it is difficult to find a general method that reasonably describes both adsorption energy and geometry of even the simple prototypes of cyclic and heterocyclic aromatic molecules such as benzene (C₆H₆) and thiophene (C₄H₄S), respectively, with reasonable accuracy. In this work, we present an alternative scheme based on Zaremba-Kohn theory, called DFT+vdW-dZK. We show that unlike other popular methods, DFT+vdW-dZK and particularly SCAN+vdW-dZK give an accurate description of the physisorption of a rare-gas atom (xenon) and two small albeit diverse prototype organic molecules on the (111) surfaces of the coinage metals.

Progress in Epitaxial Thin-Film Na3 Bi as a Topological Electronic Material

Trisodium bismuthide (Na₃ Bi) is the first experimentally verified topological Dirac semimetal, and is a 3D analogue of graphene hosting relativistic Dirac fermions. Its unconventional momentum-energy relationship is interesting from a fundamental perspective, yielding exciting physical properties such as chiral charge carriers, the chiral anomaly, and weak anti-localization. It also shows promise for realizing topological electronic devices such as topological transistors. Herein, an overview of the substantial progress achieved in the last few years on Na₃ Bi is presented, with a focus on technologically relevant large-area thin films synthesized via molecular beam epitaxy. Key theoretical aspects underpinning the unique electronic properties of Na₃ Bi are introduced. Next, the growth process on different substrates is reviewed. Spectroscopic and microscopic features are illustrated, and an analysis of semiclassical and quantum transport phenomena in different doping regimes is provided. The emergent properties arising from confinement in two dimensions, including thickness-dependent and electric-field-driven topological phase transitions, are addressed, with an outlook toward current challenges and expected future progress.

Defect-mediated ab initio thermodynamics of metastable γ-MoN(001)

Refractory transition metal nitrides exhibit a plethora of polymorphic expressions and chemical stoichiometries. To afford a better understanding of how defects may play a role in the structural and thermodynamics of these nitrides, using density-functional theory calculations, we investigate the influence of point and pair defects in bulk metastable γ-MoN and its (001) surface. We report favorable formation of Schottky defect pairs of neighboring Mo and N vacancies in bulk γ-MoN and apply this as a defect-mediated energy correction term to the surface energy of γ-MoN(001) within the ab initio atomistic thermodynamics approach. We also inspect the structural distortions in both bulk and surfaces of γ-MoN by using the partial radial distribution function, g(r), of Mo-N bond lengths. Large atomic displacements are found in both cases, leading to a broad spread of Mo-N bond length values when compared to their idealized bulk values. We propose that these structural and thermodynamic analyses may provide some insight into a better understanding of metastable materials and their surfaces.

Theoretical understanding of the properties of stepped iron surfaces with van der Waals interaction corrections

The stepped surfaces in nanoscale zero-valent iron (nZVI) play an essential role for environmental application. However, there is still currently a deficiency in the atomic understanding of stepped surface properties due to the limitation of the computational methodology. In this study, stepped Fe(210) and (211) surfaces were theoretically investigated using density functional theory (DFT) computations in terms of the flat Fe(110) surface. Our results suggest that the consideration of van der Waals (vdW) interaction correction is beneficial for the DFT study on Fe-based systems. The DF-cx method is found to be the most promising vdW correction method. The DF-cx results reveal that the stepped Fe(210) and Fe(211) surfaces experience significant surface relaxation and abnormal trends in their work function. Their electronic properties and reactivities of the surface atoms are strongly affected by the Fe coordination numbers and the strong adsorption strengths of oxygen on the surfaces are dependent on both the coordination number of the adsorbed atoms and the geometry of the adsorption sites.

Syntheses of five new layered quaternary chalcogenides SrScCuSe3, SrScCuTe3, BaScCuSe3, BaScCuTe3, and BaScAgTe3: crystal structures, thermoelectric properties, and electronic structures

Five new layered transition metal-based chalcogenides (SrScCuSe3, SrScCuTe3, BaScCuSe3, BaScCuTe3, and BaScAgTe3) were discovered by the exploratory solid-state method.

Born–Oppenheimer molecular dynamics simulations on structures of high-density and low-density water: a comparison of the SCAN meta-GGA and PBE GGA functionals

Using Born–Oppenheimer ab initio molecular dynamics (BOAIMD) simulations, the high-density water (HDW) and low-density water (LDW) structures based on SCAN meta-GGA are compared with those based on PBE GGA.

Tuplewise Material Representation Based Machine Learning for Accurate Band Gap Prediction

The open-access material databases allowed us to approach scientific questions from a completely new perspective with machine learning methods. Here, on the basis of open-access databases, we focus on the classical band gap problem for predicting accurately the band gap of a crystalline compound using a machine learning approach with newly developed tuplewise graph neural networks (TGNN), which is devised to automatically generate input representation of crystal structures in tuple types and to exploit crystal-level properties as one of the input features. Our method brings about a highly accurate prediction of the band gaps at hybrid functionals and GW approximation levels for multiple material data sets without heavy computational cost. Furthermore, to demonstrate the applicability of our prediction model, we provide a data set of GW band gaps for 45835 materials predicted by TGNN posing higher accuracy than standard density functional theory calculations.

Correlation matrix renormalization theory in multi-band lattice systems

An appropriate treatment of electronic correlation effects plays an important role in accurate descriptions of physical and chemical properties of real materials. The recently proposed correlation matrix renormalization theory with sum rule correction (CMR) [1] for studying correlated-electron materials has shown good performance in molecular systems and a periodic hydrogen chain in comparison with various quantum chemistry and quantum Monte Carlo calculations [2]. This work gives a detailed formulation and computational code implementation of CMR in multi-band periodic lattice systems. This lattice CMRab initiotheory is highly efficient, has no material specific adjustable parameters, and has no double counting issues faced by the hybrid approaches like LDA +U, DFT + DMFT and DFT + GA type theories. Benchmark studies on materials with s and p orbitals in this study show that CMR in its current implementation consistently performs well for these systems as the electron correlation increases from the bonding region to the bond breaking region.

A2AgCrCl6 (A = Li, Na, K, Rb, Cs) halide double perovskites: a transition metal-based semiconducting material series with appreciable optical characteristics

We have theoretically examined the geometries, electronic density of states and band structures of cubic and hexagonal A₂AgCrCl₆ (A = Cs, Rb, K, Na, Li) using meta-GGA SCAN-rVV10. The optimized lattice density was found to vary between 2.68 and 4.08 g cm^-3 for cubic-A₂AgCrCl₆, with the fundamental electronic bandgap (direct) in the range of 0.66-0.69 eV. The cell density of hexagonal A₂AgCrCl₆ was between 2.97 and 3.93 g cm^-3, but with an indirect bandgap of 0.93-1.02 eV. The valence band maximum and the conduction band minimum of A₂AgCrCl₆ were confirmed to be essentially of Cr(3d) character, but the contributions from the orbital states of Cl(3p) to the VBM were also appreciable. Cubic A₂AgCrCl₆ (A = Cs, Rb, K) was identified to possess genuine perovskite stoichiometry, evaluated using various geometry-based indices (viz. octahedral factor, tolerance factor, and global instability index). This was not so for A₂AgCrCl₆ (A = Na, Li), and was due to the small size of Na and Li cations that caused the critical strain of CrCl₆ octahedra and a significant decrease in the cell volume. However, all the five A₂AgCrCl₆ displayed nearly similar optical properties, including the nature of the oscillator peaks in the dielectric function, absorption coefficient, photoconductivity, reflectivity, and Tauc spectra. The zero-limit of the refractive index was calculated around 2.25 and 2.00 for cubic and hexagonal A₂AgCrCl₆, respectively, and the extinction coefficient was very small for all cases. The nature of the optical bandgap and transition peaks discussed in this study of cubic and hexagonal Cs₂AgCrCl₆ agreed well with the experiment. The examination of phonon band dispersion led to the conclusion that cubic-A₂AgCrCl₆ (A = Cs, Rb) are the only halide double perovskites of the entire series that are dynamically stable.

The Cs2AgRhCl6 Halide Double Perovskite: A Dynamically Stable Lead-Free Transition-Metal Driven Semiconducting Material for Optoelectronics

A-Site doping with alkali ions, and/or metal substitution at the B and B'-sites, are among the key strategies in the innovative development of A 2BB'X6 halide double perovskite semiconducting materials for application in energy and device technologies. To this end, we have investigated an intriguing series of five halide-based non-toxic systems, A 2AgRhCl6 (A = Li, Na, K, Rb, and Cs), using density functional theory at the SCAN-rVV10 level. The lattice stability and bonding properties emanating from this study of A 2AgRhCl6 matched well with those that have already been synthesized, characterized and discussed [viz. Cs2AgBiX6 (X = Cl, Br)]. Exploration of traditional and recently proposed tolerance factors has enabled us to identify A 2AgRhCl6 (A = K, Rb and Cs) as stable double perovskites. The band structure and density of states calculations suggested that the electronic transition from the top of the valence band [Cl(3p)+Rh(4d)] to the bottom of the conduction band [(Cl(3p)+Rh(4d)] is inherently direct at the X-point of the first Brillouin zone. The (non-spin polarized) bandgap of these materials was found in the range 0.57-0.65 eV with SCAN-rVV10, which were substantially smaller than those computed with hybrid HSE06 and PBE0, and quasi-particle GW methods. This, together with the appreciable refractive index and high absorption coefficient in the region covering the range 1.0-4.5 eV, enabled us to demonstrate that A 2AgRhCl6 (A = K, Rb, and Cs) are likely candidate materials for photoelectric applications. The results of our phonon calculations at the harmonic level suggested that the Cs2AgRhCl6 is the only system that is dynamically stable (no imaginary frequencies found around the high symmetry lines of the reciprocal lattice), although the elastic moduli properties suggested all five systems examined are mechanically stable.

Metal chalcogenide hollow polar bipyramid prisms as efficient sulfur hosts for Na-S batteries

Sodium sulfur batteries require efficient sulfur hosts that can capture soluble polysulfides and enable fast reduction kinetics. Herein, we design hollow, polar and catalytic bipyramid prisms of cobalt sulfide as efficient sulfur host for sodium sulfur batteries. Cobalt sulfide has interwoven surfaces with wide internal spaces that can accommodate sodium polysulfides and withstand volumetric expansion. Furthermore, results from in/ex-situ characterization techniques and density functional theory calculations support the significance of the polar and catalytic properties of cobalt sulfide as hosts for soluble sodium polysulfides that reduce the shuttle effect and display excellent electrochemical performance. The polar catalytic bipyramid prisms sulfur@cobalt sulfide composite exhibits a high capacity of 755 mAh g⁻¹ in the second discharge and 675 mAh g⁻¹ after 800 charge/discharge cycles, with an ultralow capacity decay rate of 0.0126 % at a high current density of 0.5 C. Additionally, at a high mass loading of 9.1 mg cm⁻², sulfur@cobalt sulfide shows high capacity of 545 mAh g⁻¹ at a current density of 0.5 C. This study demonstrates a hollow, polar, and catalytic sulfur host with a unique structure that can capture sodium polysulfides and speed up the reduction reaction of long chain sodium polysulfides to solid small chain polysulfides, which results in excellent electrochemical performance for sodium-sulfur batteries.

Sodium sulfur batteries require efficient sulfur hosts that can capture soluble polysulfides and enable fast reduction kinetics. Here, authors report hollow catalytic bipyramid prism CoS₂/C as efficient sulfur carriers, and investigate the reaction mechanism in the sodium sulfur battery.

Crystalline assembly of perylene in metal-organic framework thin film: J-aggregate or excimer? Insight into the electronic structure

The spatial orientation of chromophores defines the photophysical and optoelectronic properties of a material and serves as the main tunable parameter for tailoring functionality. Controlled assembly for achieving a predefined spatial orientation of chromophores is rather challenging. Metal-organic frameworks (MOFs) are an attractive platform for exploring the virtually unlimited chemical space of organic components and their self-assembly for device optimization. Here, we demonstrate the impact of interchromophore interactions on the photophysical properties of a surface-anchored MOF (SURMOF) based on 3,9-perylenedicarboxylicacid linkers. We predict the structural assembly of the perylene molecules in the MOF via robust periodic density functional theory calculations and discuss the impact of unit topology and π-π interaction patterns on spectroscopic and semiconducting properties of the MOF films. We explain the dual nature of excited states in the perylene MOF, where strong temperature-modulated excimer emission, enhanced by the formation of perylene J-aggregates, and low stable monomer emission are observed. We use band-like and hopping transport mechanisms to predict semiconducting properties of perylene SURMOF-2 films as a function of inter-linker interactions, demonstrating both p-type and n-type conduction mechanisms. Hole carrier mobility up to 7.34 cm²Vs^-1is predicted for the perylene SURMOF-2. The results show a promising pathway towards controlling excimer photophysics in a MOF while controlling charge carrier mobility on the basis of a predictive model.

Role of Exchange Interactions in the Magnetic Response and Intermolecular Recognition of Chiral Molecules

The physical origin of the so-called chirality-induced spin selectivity (CISS) effect has puzzled experimental and theoretical researchers over the past few years. Early experiments were interpreted in terms of unconventional spin-orbit interactions mediated by the helical geometry. However, more recent experimental studies have clearly revealed that electronic exchange interactions also play a key role in the magnetic response of chiral molecules in singlet states. In this investigation, we use spin-polarized closed-shell density functional theory calculations to address the influence of exchange contributions to the interaction between helical molecules as well as of helical molecules with magnetized substrates. We show that exchange effects result in differences in the interaction properties with magnetized surfaces, shedding light into the possible origin of two recent important experimental results: enantiomer separation and magnetic exchange force microscopy with AFM tips functionalized with helical peptides.

Photocatalytic water splitting for hydrogen generation driven by tetragonal, trigonal, hexagonal and cubic LiCoO2 and visible light

The photocatalytic properties of LiCoO₂ are not explored up to date although its cubic and trigonal structures are explored experimentally. Here, we investigate the feasibility of photocatalytic hydrogen production from water splitting driven by the tetragonal, trigonal, hexagonal and cubic LiCoO₂ with the irradiation of the visible light. The band structure, density of state, optical absorption and mobility are calculated by the first-principles density functional theory. The results show that the band edges of all the four structures of LiCoO₂ match to the conditions of the redox potentials of water splitting reaction and the enhanced optical absorption in the visible light range is observed. The obvious difference between the mobilities of the hole and electron are identified, especially for the cubic LiCoO₂. All the obtained results suggest that the considered structures of LiCoO₂ are promising candidates for the photocatalytic water splitting to produce hydrogen with the irradiation of the visible light.

Accurate and Numerically Efficient r2SCAN Meta-Generalized Gradient Approximation

The recently proposed rSCAN functional [ J. Chem. Phys. 2019 150, 161101] is a regularized form of the SCAN functional [ Phys. Rev. Lett. 2015 115, 036402] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.

Pressure-Dependent Intermediate Magnetic Phase in Thin Fe3GeTe2 Flakes

We investigated the evolution of ferromagnetism in layered Fe₃GeTe₂ flakes under different pressures and temperatures using in situ magnetic circular dichroism (MCD) spectroscopy. We found that the rectangular shape of the hysteresis loop under an out-of-plane magnetic field sweep can be sustained below 7 GPa. Above that pressure, an intermediate state appears in the low-temperature region signaled by an 8-shaped skewed hysteresis loop. Meanwhile, the coercive field and Curie temperature decrease with increasing pressures, implying the decrease of the exchange interaction and the magneto-crystalline anisotropy under pressures. The intermediate phase has a labyrinthine domain structure, which is attributed to the increase of the ratio of exchange interaction to magneto-crystalline anisotropy based on Jagla's theory. Moreover, our calculations reveal a weak structural transition around 6 GPa that corresponds to a significant change in the FeI-FeI bond length, which has strong influences on magnetic interaction. Detailed analysis on exchange interaction and magneto-crystalline anisotropy with pressure shows a consistent trend with experiments.

Electron-electron correlations and structural, spectral and polarization properties of tetragonal BaTiO3

To analyze the role of electron-electron correlation effects in structural (local-geometry), spectral and polarization properties of tetragonal BaTiO3we apply DFT +Uapproach. We demonstrate that the system properties drastically change when the value of the local Coulomb repulsionUcrosses the critical valueUc≈ 7 eV. In particular, the correlation effects cause a change of the ratio of the in-plane and inter-plane Ti-O bond lengths, which results in a flip of the order of the Tid-bands and change of the polarizability of the system. Since the consensus value ofUin BaTiO3is unknown, we discuss how the obtained results may be revealed in experimental data, especially in the optical response and ultrafast charge dynamics, where effectiveUis dynamically tuned.

Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA

Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew-Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.

High-Temperature Quantum Anomalous Hall Insulators in Lithium-Decorated Iron-Based Superconductor Materials

Quantum anomalous Hall (QAH) insulator is the key material to study emergent topological quantum effects, but its ultralow working temperature limits experiments. Here, by first-principles calculations, we find a family of stable two-dimensional (2D) structures generated by lithium decoration of layered iron-based superconductor materials Fe X(X=S,Se,Te), and predict room-temperature ferromagnetic semiconductors together with large-gap high-Chern-number QAH insulators in the 2D materials. The extremely robust ferromagnetic order is induced by the electron injection from Li to Fe and stabilized by strong ferromagnetic kinetic exchange in the 2D Fe layer. While in the absence of spin-orbit coupling (SOC), the ferromagnetism polarizes the system into a half Dirac semimetal state protected by mirror symmetry, the SOC effect results in a spontaneous breaking of mirror symmetry and introduces a Dirac mass term, which creates QAH states with sizable gaps (several tens of meV) and multiple chiral edge modes. We also find a 3D QAH insulator phase featured by a macroscopic number of chiral conduction channels in bulk LiOH-LiFe X. The findings open new opportunities to realize novel QAH physics and applications at high temperatures.

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to the leading order (-0.221Z^5/3 and -0.021Z ln Z, respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E_x(N, Z) ≈ -0.354N^2/3Z (as known before only for Z ≫ N ≫ 1) and E_c ≈ -0.02N ln N. These asymptotic estimates are most correct for atomic ions with large N and Z ≫ N, but we find that they are qualitatively and semi-quantitatively correct even for small N and N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

Highly reversible oxygen redox in layered compounds enabled by surface polyanions

Oxygen-anion redox in lithium-rich layered oxides can boost the capacity of lithium-ion battery cathodes. However, the over-oxidation of oxygen at highly charged states aggravates irreversible structure changes and deteriorates cycle performance. Here, we investigate the mechanism of surface degradation caused by oxygen oxidation and the kinetics of surface reconstruction. Considering Li₂MnO₃, we show through density functional theory calculations that a high energy orbital (lO_2p') at under-coordinated surface oxygen prefers over-oxidation over bulk oxygen, and that surface oxygen release is then kinetically favored during charging. We use a simple strategy of turning under-coordinated surface oxygen into polyanionic (SO₄)^2-, and show that these groups stabilize the surface of Li₂MnO₃ by depressing gas release and side reactions with the electrolyte. Experimental validation on Li_1.2Ni_0.2Mn_0.6O₂ shows that sulfur deposition enhances stability of the cathode with 99.0% capacity remaining (194 mA h g^-1) after 100 cycles at 1 C. Our work reveals a promising surface treatment to address the instability of highly charged layered cathode materials.

Spontaneously separated intermetallic Co3Mo from nanoporous copper as versatile electrocatalysts for highly efficient water splitting

Developing robust nonprecious electrocatalysts towards hydrogen/oxygen evolution reactions is crucial for widespread use of electrochemical water splitting in hydrogen production. Here, we report that intermetallic Co₃Mo spontaneously separated from hierarchical nanoporous copper skeleton shows genuine potential as highly efficient electrocatalysts for alkaline hydrogen/oxygen evolution reactions in virtue of in-situ hydroxylation and electro-oxidation, respectively. The hydroxylated intermetallic Co₃Mo has an optimal hydrogen-binding energy to facilitate adsorption/desorption of hydrogen intermediates for hydrogen molecules. Associated with high electron/ion transport of bicontinuous nanoporous skeleton, nanoporous copper supported Co₃Mo electrodes exhibit impressive hydrogen evolution reaction catalysis, with negligible onset overpotential and low Tafel slope (~40 mV dec⁻¹) in 1 M KOH, realizing current density of −400 mA cm⁻² at overpotential of as low as 96 mV. When coupled to its electro-oxidized derivative that mediates efficiently oxygen evolution reaction, their alkaline electrolyzer operates with a superior overall water-splitting output, outperforming the one assembled with noble-metal-based catalysts.

Electrochemical water splitting is an attractive energy conversion technology, but it usually suffers from low efficiency. Here, the authors report intermetallic Co₃Mo integrated on porous Cu as highly efficient electrocatalysts for alkaline HER/OER due to in-situ hydroxylation and electro-oxidation.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

A2AgCrBr6 (A = K, Rb, Cs) and Cs2AgCrX6(X = Cl, I) Double Perovskites: A Transition-Metal-Based Semiconducting Material Series with Remarkable Optics

With an interest to quest for transition metal-based halogenated double perovskites ABBX6 as high performance semiconducting materials for optoelectronics, this study theoretically examined the electronic structures, stability, electronic (density of states and band structures), transport (effective masses of charge carriers), and optical properties (dielectric function and absorption coefficients, etc.) of the series A2AgCrBr6 (A = K, Rb, Cs) using SCAN+rVV10. Our results showed that A2AgCrBr6 (A = Rb, Cs), but not K2AgCrBr6, has a stable perovskite structure, which was revealed using various traditionally recommended geometry-based indices. Despite this reservation, all the three systems were shown to have similar band structures, density of states, and carrier effective masses of conducting holes and electrons, as well as the nature of the real and imaginary parts of their dielectric function, absorption coefficient, refractive index, and photoconductivity spectra. The small changes observed in any specific property of the series A2AgCrBr6 were due to the changes in the lattice properties driven by alkali substitution at the A site. A comparison with the corresponding properties of Cs2AgCrX6 (X = Cl, I) suggested that halogen substitution at the X-site can not only significantly shift the position of the onset of optical absorption found of the dielectric function, absorption coefficient and refractive spectra of Cs2AgCrCl6 and Cs2AgCrI6 toward the high- and low-energy infrared regions, respectively; but that it is also responsible in modifying their stability, electronic, transport, and optical absorption preferences. The large value of the high frequency dielectric constants-together with the appreciable magnitude of absorption coefficients and refractive indices, small values of effective masses of conducting electrons and holes, and the indirect nature of the bandgap transitions, among others-suggested that cubic A2AgCrBr6 (A = Rb, Cs) and Cs2AgCrCl6 may likely be a set of optoelectronic materials for subsequent experimental characterizations.

Self-interaction error overbinds water clusters but cancels in structural energy differences

We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi-Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction-corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/[Formula: see text] from about 42 meV/[Formula: see text] for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction-corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

Highly dispersed Co deposited on Al2O3 particles via CoCp2 + H2 ALD

Highly dispersed cobalt atoms were deposited on porous alumina particles using atomic layer deposition (ALD) with a CoCp₂/H₂ chemistry at approximately 7 wt%. H₂ did not completely reduce the cyclopentadienyl organic ligands bound to deposited Co atoms at ALD reaction conditions. A sharp decline in Co deposited per cycle for two or more ALD cycles indicates that much of the Al₂O₃ surface is sterically blocked from further CoCp₂ deposition after the first CoCp₂ exposure. Temperature programmed reduction confirmed that the adsorbed precursor organic ligands persist after H₂ exposures during ALD and temperatures as high as 500 °C are required to fully reduce the organic ligands to CH₄. High resolution, element sensitive imaging showed that Co atoms were dispersed on the Al₂O₃ surface and could deposit in previously unobserved multiple growth morphologies, specifically layers that were continuous over several angstroms or discrete nanoparticles. Density functional theory calculations were used to examine Co atom adsorption, show the altered haptic binding of cracked Cp ligands, and to calculate the thermodynamics of Cp ligand decomposition. The lateral steric hindrance between organic ligands bound to deposited Co atoms, Cp ligand decomposition mechanism, and local Al₂O₃ surface termination all likely determine the observed Co growth morphology during initial ALD cycles.

Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent nontrivial effects on liquid-state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio-simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.

Direct Measurement of Electronic Band Structures at Oxide Grain Boundaries

Grain boundaries (GBs) modulate the macroscopic properties in polycrystalline materials because they have different atomic and electronic structures from the bulk. Despite the progress on the understanding of GB atomic structures, knowledge of the localized electronic band structures is still lacking. Here, we experimentally characterized the atomic structures and the band gaps of four typical GBs in α-Al₂O₃ by scanning transmission electron microscopy and valence electron energy-loss spectroscopy (EELS). It was found that the band gaps of the GBs are narrowed by 0.5-2.1 eV compared with that of 8.8 eV in the bulk. By combing core-loss EELS with first-principles calculations, we elucidated that the band gap reductions directly correlate with the decrease of the coordination numbers of Al and O ions at the GBs. These results provide in-depth understanding between the local atomic and electronic band structures for GBs and demonstrate a novel electronic-structure analysis for crystalline defects.

Hierarchically 3D Porous Ag Nanostructures Derived from Silver Benzenethiolate Nanoboxes: Enabling CO2 Reduction with a Near-Unity Selectivity and Mass-Specific Current Density over 500 A/g

Silver nanostructures with hierarchical porosities of multiple length scales have been synthesized through electrochemical reduction of silver benzenethiolate nanoboxes. The porous Ag nanostructures exhibit superior catalytic performance toward electrochemical reduction of CO₂. The Faradaic efficiency of reducing CO₂ to CO can be close to 100% at high cathodic potentials, benefiting from the readsorbed benzenethiolate ions on the Ag surface that can suppress the hydrogen evolution reaction (HER). Density functional theory calculations using the SCAN functional reveal that the disfavored H binding on the benzenethiolate-modified Ag surface is responsible for inhibiting the HER. The mass-specific activity of CO₂ reduction can be over 500 A/g because the multiple-scale porosities maximize the diffusion of reactive species to and away from the Ag surface. The unique multiscale porosities and surface modification of the as-synthesized Ag nanostructures make them a class of promising catalysts for electrochemical reduction of CO₂ in protic electrolytes to achieve maximum activity and selectivity.