Incorporating solvent effects in DFT: insights from cation exchange in faujasites.

Zeolites are versatile materials renowned for their extra-framework cation exchange capabilities, with applications spanning diverse fields, including nuclear waste treatment. While detailed experimental characterization offers valuable insight, density functional theory (DFT) proves particularly adept at investigating ion exchange in zeolites, owing to its atomic and electronic resolution. However, the prevalent occurrence of zeolitic ion exchange in aqueous environments poses a challenge to conventional DFT modeling, traditionally conducted in a vacuum. This study seeks to enhance zeolite modeling by systematically evaluating predictive differences across varying degrees of aqueous solvent inclusion. Specifically focusing on monovalent cation exchange in Na-X zeolites, we explore diverse modeling approaches. These range from simple dehydrated systems (representing bare reference states in vacuum) to more sophisticated models that incorporate aqueous solvent effects through explicit water molecules and/or a dielectric medium. Through comparative analysis of DFT and semi-empirical DFT approaches, along with their validation against experimental results, our findings underscore the necessity to concurrently consider explicit and implicit solvent effects for accurate prediction of zeolitic ionic exchange.

Charge state-dependent symmetry breaking of atomic defects in transition metal dichalcogenides.

The functionality of atomic quantum emitters is intrinsically linked to their host lattice coordination. Structural distortions that spontaneously break the lattice symmetry strongly impact their optical emission properties and spin-photon interface. Here we report on the direct imaging of charge state-dependent symmetry breaking of two prototypical atomic quantum emitters in mono- and bilayer MoS2 by scanning tunneling microscopy (STM) and non-contact atomic force microscopy (nc-AFM). By changing the built-in substrate chemical potential, different charge states of sulfur vacancies (VacS) and substitutional rhenium dopants (ReMo) can be stabilized. Vac S - 1 as well as Re Mo 0 and Re Mo - 1 exhibit local lattice distortions and symmetry-broken defect orbitals attributed to a Jahn-Teller effect (JTE) and pseudo-JTE, respectively. By mapping the electronic and geometric structure of single point defects, we disentangle the effects of spatial averaging, charge multistability, configurational dynamics, and external perturbations that often mask the presence of local symmetry breaking.

Transition Temperature for Spin-Crossover Materials with the Mean Value Ensemble Hubbard-U Correction.

Calculation of transition temperatures T1/2 for thermally driven spin-crossover in condensed phases is challenging, even with sophisticated state-of-the-art density functional approximations. The first issue is the accuracy of the adiabatic crossover energy difference ΔEHL between the low- and high-spin states of the bistable metal-organic complexes. The other is the proper inclusion of entropic contributions to the Gibbs free energy from the electronic and vibrational degrees of freedom. We discuss the effects of treatments of both contributions upon the calculation of thermochemical properties for a set of 20 spin-crossover materials using a Hubbard-U correction obtained from a reference ensemble spin-state. The U values obtained from a simplest bimolecular representation may overcorrect, somewhat, the ΔEHL values, hence giving somewhat excessive reduction of the T1/2 results with respect to their U = 0 values in the crystalline phase. We discuss the origins of the discrepancies by analyzing different sources of uncertainties. By use of a first-coordination-sphere approximation and the assumption that vibrational contributions from the outermost atoms in a metal-organic complex are similar in both low- and high-spin states, we achieve T1/2 results with the low-cost, widely used PBE generalized gradient density functional approximation comparable to those from the more costly, more sophisticated r2SCAN meta-generalized gradient approximation. The procedure is promising for use in high-throughput materials screening, because it combines rather low computational effort requirements with freedom from user manipulation of parameters.

Tailoring the Angular Mismatch in MoS2 Homobilayers through Deformation Fields.

Ultrathin MoS2 has shown remarkable characteristics at the atomic scale with an immutable disorder to weak external stimuli. Ion beam modification unlocks the potential to selectively tune the size, concentration, and morphology of defects produced at the site of impact in 2D materials. Combining experiments, first-principles calculations, atomistic simulations, and transfer learning, it is shown that irradiation-induced defects can induce a rotation-dependent moiré pattern in vertically stacked homobilayers of MoS2 by deforming the atomically thin material and exciting surface acoustic waves (SAWs). Additionally, the direct correlation between stress and lattice disorder by probing the intrinsic defects and atomic environments are demonstrated. The method introduced in this paper sheds light on how engineering defects in the lattice can be used to tailor the angular mismatch in van der Waals (vdW) solids.

Mean Value Ensemble Hubbard-U Correction for Spin-Crossover Molecules.

High-throughput searches for spin-crossover molecules require Hubbard-U corrections to common density functional exchange-correlation (XC) approximations. However, the Ueff values obtained from linear response or based on previous studies overcorrect the spin-crossover energies. We demonstrate that employing a linearly mixed ensemble average spin state as the reference configuration for the linear response calculation of Ueff resolves this issue. Validation on a commonly used set of spin-crossover complexes shows that these ensemble Ueff values consistently are smaller than those calculated directly on a pure spin state, irrespective of whether that be low- or high-spin. Adiabatic crossover energies using this methodology for a generalized gradient approximation XC functional are closer to the expected target energy range than with conventional Ueff values. Based on the observation that the Ueff correction is similar for different complexes that share transition metals with the same oxidation state, we devise a set of recommended averaged Ueff values for high-throughput calculations.

Photoluminescence Induced by Substitutional Nitrogen in Single-Layer Tungsten Disulfide.

The electronic and optical properties of two-dimensional materials can be strongly influenced by defects, some of which can find significant implementations, such as controllable doping, prolonged valley lifetime, and single-photon emissions. In this work, we demonstrate that defects created by remote N2 plasma exposure in single-layer WS2 can induce a distinct low-energy photoluminescence (PL) peak at 1.59 eV, which is in sharp contrast to that caused by remote Ar plasma. This PL peak has a critical requirement on the N2 plasma exposure dose, which is strongest for WS2 with about 2.0% sulfur deficiencies (including substitutions and vacancies) and vanishes at 5.6% or higher sulfur deficiencies. Both experiments and first-principles calculations suggest that this 1.59 eV PL peak is caused by defects related to the sulfur substitutions by nitrogen, even though low-temperature PL measurements also reveal that not all the sulfur vacancies are remedied by the substitutional nitrogen. The distinct low-energy PL peak suggests that the substitutional nitrogen defect in single-layer WS2 can potentially serve as an isolated artificial atom for creating single-photon emitters, and its intensity can also be used to monitor the doping concentrations of substitutional nitrogen.

Ligand Optimization of Exchange Interaction in Co(II) Dimer Single Molecule Magnet by Machine Learning.

Designing single-molecule magnets (SMMs) for potential applications in quantum computing and high-density data storage requires tuning their magnetic properties, especially the strength of the magnetic interaction. These properties can be characterized by first-principles calculations based on density functional theory (DFT). In this work, we study the experimentally synthesized Co(II) dimer (Co2(C5NH5)4(µ-PO2(CH2C6H5)2)3) SMM with the goal to control the exchange energy, ΔEJ, between the Co atoms through tuning of the capping ligands. The experimentally synthesized Co(II) dimer molecule has a very small ΔEJ < 1 meV. We assemble a DFT data set of 1081 ligand substitutions for the Co(II) dimer. The ligand exchange provides a broad range of exchange energies, ΔEJ, from +50 to -200 meV, with 80% of the ligands yielding a small ΔEJ < 10 meV. We identify descriptors for the classification and regression of ΔEJ using gradient boosting machine learning models. We compare one-hot encoded, structure-based, and chemical descriptors consisting of the HOMO/LUMO energies of the individual ligands and the maximum electronegativity difference and bond order for the ligand atom connecting to Co. We observe a similar overall performance with the chemical descriptors outperforming the other descriptors. We show that the exchange coupling, ΔEJ, is correlated to the difference in the average bridging angle between the ferromagnetic and antiferromagnetic states, similar to the Goodenough-Kanamori rules.

The 2021 room-temperature superconductivity roadmap.

Designing materials with advanced functionalities is the main focus of contemporary solid-state physics and chemistry. Research efforts worldwide are funneled into a few high-end goals, one of the oldest, and most fascinating of which is the search for an ambient temperature superconductor (A-SC). The reason is clear: superconductivity at ambient conditions implies being able to handle, measure and access a single, coherent, macroscopic quantum mechanical state without the limitations associated with cryogenics and pressurization. This would not only open exciting avenues for fundamental research, but also pave the road for a wide range of technological applications, affecting strategic areas such as energy conservation and climate change. In this roadmap we have collected contributions from many of the main actors working on superconductivity, and asked them to share their personal viewpoint on the field. The hope is that this article will serve not only as an instantaneous picture of the status of research, but also as a true roadmap defining the main long-term theoretical and experimental challenges that lie ahead. Interestingly, although the current research in superconductor design is dominated by conventional (phonon-mediated) superconductors, there seems to be a widespread consensus that achieving A-SC may require different pairing mechanisms.In memoriam, to Neil Ashcroft, who inspired us all.

Synthesis of borophane polymorphs through hydrogenation of borophene.

Synthetic two-dimensional polymorphs of boron, or borophene, have attracted attention because of their anisotropic metallicity, correlated-electron phenomena, and diverse superlattice structures. Although borophene heterostructures have been realized, ordered chemical modification of borophene has not yet been reported. Here, we synthesize "borophane" polymorphs by hydrogenating borophene with atomic hydrogen in ultrahigh vacuum. Through atomic-scale imaging, spectroscopy, and first-principles calculations, the most prevalent borophane polymorph is shown to possess a combination of two-center-two-electron boron-hydrogen and three-center-two-electron boron-hydrogen-boron bonds. Borophane polymorphs are metallic with modified local work functions and can be reversibly returned to pristine borophene through thermal desorption of hydrogen. Hydrogenation also provides chemical passivation because borophane reduces oxidation rates by more than two orders of magnitude after ambient exposure.

Scalable Substitutional Re-Doping and its Impact on the Optical and Electronic Properties of Tungsten Diselenide.

Reliable, controlled doping of 2D transition metal dichalcogenides will enable the realization of next-generation electronic, logic-memory, and magnetic devices based on these materials. However, to date, accurate control over dopant concentration and scalability of the process remains a challenge. Here, a systematic study of scalable in situ doping of fully coalesced 2D WSe2 films with Re atoms via metal-organic chemical vapor deposition is reported. Dopant concentrations are uniformly distributed over the substrate surface, with precisely controlled concentrations down to <0.001% Re achieved by tuning the precursor partial pressure. Moreover, the impact of doping on morphological, chemical, optical, and electronic properties of WSe2 is elucidated with detailed experimental and theoretical examinations, confirming that the substitutional doping of Re at the W site leads to n-type behavior of WSe2 . Transport characteristics of fabricated back-gated field-effect-transistors are directly correlated to the dopant concentration, with degrading device performances for doping concentrations exceeding 1% of Re. The study demonstrates a viable approach to introducing true dopant-level impurities with high precision, which can be scaled up to batch production for applications beyond digital electronics.

Implicit self-consistent electrolyte model in plane-wave density-functional theory.

The ab initio computational treatment of electrochemical systems requires an appropriate treatment of the solid/liquid interfaces. A fully quantum mechanical treatment of the interface is computationally demanding due to the large number of degrees of freedom involved. In this work, we develop a computationally efficient model where the electrode part of the interface is described at the density-functional theory (DFT) level, and the electrolyte part is represented through an implicit solvation model based on the Poisson-Boltzmann equation. We describe the implementation of the linearized Poisson-Boltzmann equation into the Vienna Ab initio Simulation Package, a widely used DFT code, followed by validation and benchmarking of the method. To demonstrate the utility of the implicit electrolyte model, we apply it to study the surface energy of Cu crystal facets in an aqueous electrolyte as a function of applied electric potential. We show that the applied potential enables the control of the shape of nanocrystals from an octahedral to a truncated octahedral morphology with increasing potential.

Computational Study of Low Interlayer Friction in Tin+1Cn (n = 1, 2, and 3) MXene.

The friction of adjacent Tin+1Cn (n = 1, 2, and 3) MXene layers is investigated using density functional theory (DFT) calculations and classical molecular dynamics simulations with ReaxFF potentials. The calculations reveal the sliding pathways in all three MXene systems with low energy barriers. The friction coefficients for interlayer sliding are evaluated using static calculations. Both DFT and ReaxFF methods predict friction coefficients between 0.24 and 0.27 for normal loads less than 1.2 GPa. The effect of titanium (Ti) vacancies in sublayers and terminal oxygen (O) vacancies at surfaces on the interlayer friction is further investigated using the ReaxFF potential. These defects are found to increase the friction coefficients by increasing surface roughness and creating additional attractive forces between adjacent layers. However, these defective MXenes still maintain friction coefficients below 0.31. We also consider functionalized Ti3C2 MXene terminated with -OH and -OCH3 and find that compared to the -O-terminated surface both groups further reduce the interlayer friction coefficient to 0.10-0.14.

Interface-Driven Structural Distortions and Composition Segregation in Two-Dimensional Heterostructures.

The discovery of emergent phenomena in 2D materials has sparked substantial research efforts in the materials community. A significant experimental challenge for this field is exerting atomistic control over the structure and composition of the constituent 2D layers and understanding how the interactions between layers drive both structure and properties. While no segregation for single bilayers was observed, segregation of Pb to the surface of three bilayer thick PbSe-SnSe alloy layers was discovered within [(Pbx Sn1-x Se)1+δ ]n (TiSe2 )1 heterostructures using electron microscopy. This segregation is thermodynamically favored to occur when Pbx Sn1-x Se layers are interdigitated with TiSe2 monolayers. DFT calculations indicate that the observed segregation depends on what is adjacent to the Pbx Sn1-x Se layers. The interplay between interface- and volume-free energies controls both the structure and composition of the constituent layers, which can be tuned using layer thickness.

Two-Dimensional Intrinsic Half-Metals With Large Spin Gaps.

Through a systematic search of all layered bulk compounds combined with density functional calculations employing hybrid exchange-correlation functionals, we predict a family of three magnetic two-dimensional (2D) materials with half-metallic band structures. The 2D materials, FeCl2, FeBr2, and FeI2, are all sufficiently stable to be exfoliated from bulk layered compounds. The Fe2+ ions in these materials are in a high-spin octahedral d6 configuration leading to a large magnetic moment of 4 µB. Calculations of the magnetic anisotropy show an easy-plane for the magnetic moment. A classical XY model with nearest neighbor coupling estimates critical temperatures, Tc, for the Berezinskii-Kosterlitz-Thouless transition ranging from 122 K for FeI2 to 210 K for FeBr2. The quantum confinement of these 2D materials results in unusually large spin gaps, ranging from 4.0 eV for FeI2 to 6.4 eV for FeCl2, which should defend against spin current leakage even at small device length scales. Their purely spin-polarized currents and dispersive interlayer interactions should make these materials useful for 2D spin valves and other spintronic applications.

Topology-Scaling Identification of Layered Solids and Stable Exfoliated 2D Materials.

The Materials Project crystal structure database has been searched for materials possessing layered motifs in their crystal structures using a topology-scaling algorithm. The algorithm identifies and measures the sizes of bonded atomic clusters in a structure's unit cell, and determines their scaling with cell size. The search yielded 826 stable layered materials that are considered as candidates for the formation of two-dimensional monolayers via exfoliation. Density-functional theory was used to calculate the exfoliation energy of each material and 680 monolayers emerge with exfoliation energies below those of already-existent two-dimensional materials. The crystal structures of these two-dimensional materials provide templates for future theoretical searches of stable two-dimensional materials. The optimized structures and other calculated data for all 826 monolayers are provided at our database (https://materialsweb.org).

Evaluation and comparison of classical interatomic potentials through a user-friendly interactive web-interface.

Classical empirical potentials/force-fields (FF) provide atomistic insights into material phenomena through molecular dynamics and Monte Carlo simulations. Despite their wide applicability, a systematic evaluation of materials properties using such potentials and, especially, an easy-to-use user-interface for their comparison is still lacking. To address this deficiency, we computed energetics and elastic properties of variety of materials such as metals and ceramics using a wide range of empirical potentials and compared them to density functional theory (DFT) as well as to experimental data, where available. The database currently consists of 3248 entries including energetics and elastic property calculations, and it is still increasing. We also include computational tools for convex-hull plots for DFT and FF calculations. The data covers 1471 materials and 116 force-fields. In addition, both the complete database and the software coding used in the process have been released for public use online (presently at http://www.ctcms.nist.gov/∼knc6/periodic.html) in a user-friendly way designed to enable further material design and discovery.

Structural Changes in 2D BiSe Bilayers as n Increases in (BiSe)1+δ(NbSe2)n (n = 1-4) Heterostructures.

(BiSe)1+δ(NbSe2)n heterostructures with n = 1-4 were synthesized using modulated elemental reactants. The BiSe bilayer structure changed from a rectangular basal plane with n = 1 to a square basal plane for n = 2-4. The BiSe in-plane structure was also influenced by small changes in the structure of the precursor, without significantly changing the out-of-plane diffraction pattern or value of the misfit parameter, δ. Density functional theory calculations on isolated BiSe bilayers showed that its lattice is very flexible, which may explain its readiness to adjust shape and size depending on the environment. Correlated with the changes in the BiSe basal plane structure, analysis of scanning transmission electron microscope images revealed that the occurrence of antiphase boundaries, found throughout the n = 1 compound, is dramatically reduced for the n = 2-4 compounds. X-ray photoelectron spectroscopy measurements showed that the Bi 5d3/2, 5d5/2 doublet peaks narrowed toward higher binding energies as n increased from 1 to 2, also consistent with a reduction in the number of antiphase boundaries. Temperature-dependent electrical resistivity and Hall coefficient measurements of nominally stoichiometric samples in conjunction with structural refinements and XPS data suggest a constant amount of interlayer charge transfer independent of n. Constant interlayer charge transfer is surprising given the changes in the BiSe in-plane structure. The structural flexibility of the BiSe layer may be useful in designing multiple constituent heterostructures as an interlayer between structurally dissimilar constituents.

Enhanced Li-S Batteries Using Amine-Functionalized Carbon Nanotubes in the Cathode.

The rechargeable lithium-sulfur (Li-S) battery is an attractive platform for high-energy, low-cost electrochemical energy storage. Practical Li-S cells are limited by several fundamental issues, including the low conductivity of sulfur and its reduction compounds with Li and the dissolution of long-chain lithium polysulfides (LiPS) into the electrolyte. We report on an approach that allows high-performance sulfur-carbon cathodes to be designed based on tethering polyethylenimine (PEI) polymers bearing large numbers of amine groups in every molecular unit to hydroxyl- and carboxyl-functionalized multiwall carbon nanotubes. Significantly, for the first time we show by means of direct dissolution kinetics measurements that the incorporation of CNT-PEI hybrids in a sulfur cathode stabilizes the cathode by both kinetic and thermodynamic processes. Composite sulfur cathodes based the CNT-PEI hybrids display high capacity at both low and high current rates, with capacity retention rates exceeding 90%. The attractive electrochemical performance of the materials is shown by means of DFT calculations and physical analysis to originate from three principal sources: (i) specific and strong interaction between sulfur species and amine groups in PEI; (ii) an interconnected conductive CNT substrate; and (iii) the combination of physical and thermal sequestration of LiPS provided by the CNT=PEI composite.

Ab Initio Prediction of Piezoelectricity in Two-Dimensional Materials.

Two-dimensional (2D) materials present many unique materials concepts, including material properties that sometimes differ dramatically from those of their bulk counterparts. One of these properties, piezoelectricity, is important for micro- and nanoelectromechanical systems applications. Using symmetry analysis, we determine the independent piezoelectric coefficients for four groups of predicted and synthesized 2D materials. We calculate with density-functional perturbation theory the stiffness and piezoelectric tensors of these materials. We determine the in-plane piezoelectric coefficient d11 for 37 materials within the families of 2D metal dichalcogenides, metal oxides, and III-V semiconductor materials. A majority of the structures, including CrSe2, CrTe2, CaO, CdO, ZnO, and InN, have d11 coefficients greater than 5 pm/V, a typical value for bulk piezoelectric materials. Our symmetry analysis shows that buckled 2D materials exhibit an out-of-plane coefficient d31. We find that d31 for 8 III-V semiconductors ranges from 0.02 to 0.6 pm/V. From statistical analysis, we identify correlations between the piezoelectric coefficients and the electronic and structural properties of the 2D materials that elucidate the origin of the piezoelectricity. Among the 37 2D materials, CdO, ZnO, and CrTe2 stand out for their combination of large piezoelectric coefficient and low formation energy and are recommended for experimental exploration.

Computational Screening of 2D Materials for Photocatalysis.

Two-dimensional (2D) materials exhibit a range of extraordinary electronic, optical, and mechanical properties different from their bulk counterparts with potential applications for 2D materials emerging in energy storage and conversion technologies. In this Perspective, we summarize the recent developments in the field of solar water splitting using 2D materials and review a computational screening approach to rapidly and efficiently discover more 2D materials that possess properties suitable for solar water splitting. Computational tools based on density-functional theory can predict the intrinsic properties of potential photocatalyst such as their electronic properties, optical absorbance, and solubility in aqueous solutions. Computational tools enable the exploration of possible routes to enhance the photocatalytic activity of 2D materials by use of mechanical strain, bias potential, doping, and pH. We discuss future research directions and needed method developments for the computational design and optimization of 2D materials for photocatalysis.