Adsorption and vibrational spectroscopy of CO on the surface of MgO from periodic local coupled-cluster theory.

The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of +14.7 cm-1, in excellent agreement with the experimental shift of +14.3 cm-1. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by 5-15 kJ mol-1 and predict a negative frequency shift of about -20 cm-1, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.

Theory of acoustic polarons in the two-dimensional SSH model applied to the layered superatomic semiconductor Re6Se8Cl2.

Layered superatomic semiconductors, whose building blocks are atomically precise molecular clusters, exhibit interesting electronic and vibrational properties. In recent work [Tulyagankhodjaev et al., Science 382, 438 (2023)], transient reflection microscopy revealed quasi-ballistic exciton dynamics in Re6Se8Cl2, which was attributed to the formation of polarons due to coupling with acoustic phonons. Here, we characterize the electronic, excitonic, and phononic properties with periodic density functional theory. We further parameterize a polaron Hamiltonian with nonlocal (Su-Schrieffer-Heeger) coupling to an acoustic phonon to study the polaron ground state binding energy and dispersion relation with variational wavefunctions. We calculate a polaron binding energy of about 10 meV at room temperature, and the maximum group velocity of our polaron dispersion relation is 1.5 km/s, which is similar to the experimentally observed exciton transport velocity.

Performance of periodic EOM-CCSD for bandgaps of inorganic semiconductors and insulators.

We calculate bandgaps of 12 inorganic semiconductors and insulators composed of atoms from the first three rows of the Periodic Table using periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Our calculations are performed with atom-centered triple-zeta basis sets and up to 64 k-points in the Brillouin zone. We analyze the convergence behavior with respect to the number of orbitals and number of k-points sampled using composite corrections and extrapolations to produce our final values. When accounting for electron-phonon corrections to experimental bandgaps, we find that EOM-CCSD has a mean signed error of -0.12 eV and a mean absolute error of 0.42 eV; the largest outliers are C (error of -0.93 eV), BP (-1.00 eV), and LiH (+0.78 eV). Surprisingly, we find that the more affordable partitioned EOM-MP2 theory performs as well as EOM-CCSD.

Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics.

In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.

Interplay between Local Moment and Itinerant Magnetism in the Layered Metallic Antiferromagnet TaFe1.14Te3.

Two-dimensional antiferromagnets have garnered considerable interest for the next generation of functional spintronics. However, many bulk materials from which two-dimensional antiferromagnets are isolated are limited by their air sensitivity, low ordering temperatures, and insulating transport properties. TaFe1+yTe3 aims to address these challenges with increased air stability, metallic transport, and robust antiferromagnetism. Here, we synthesize TaFe1+yTe3 (y = 0.14), identify its structural, magnetic, and electronic properties, and elucidate the relationships between them. Axial-dependent high-field magnetization measurements on TaFe1.14Te3 reveal saturation magnetic fields ranging between 27 and 30 T with saturation magnetic moments of 2.05-2.12 µB. Magnetotransport measurements confirm that TaFe1.14Te3 is metallic with strong coupling between magnetic order and electronic transport. Angle-resolved photoemission spectroscopy measurements across the magnetic transition uncover a complex interplay between itinerant electrons and local magnetic moments that drives the magnetic transition. We demonstrate the ability to isolate few-layer sheets of TaFe1.14Te3, establishing TaFe1.14Te3 as a potential platform for two-dimensional spintronics.

Highly Accurate Electronic Structure of Metallic Solids from Coupled-Cluster Theory with Nonperturbative Triple Excitations.

Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T))-often considered the "gold standard" of main-group quantum chemistry-is inapplicable to three-dimensional metals due to an infrared divergence, preventing its application to many important problems in materials science. We study the full, nonperturbative inclusion of triple excitations (CCSDT) and propose a new, iterative method, which we call ring-CCSDT, that resums the essential triple excitations with the same N^{7} run-time scaling as CCSD(T). CCSDT and ring-CCSDT are used to calculate the correlation energy of the uniform electron gas at metallic densities and the structural properties of solid lithium. Inclusion of connected triple excitations is shown to be essential to achieving high accuracy. We also investigate semiempirical CC methods based on spin-component scaling and the distinguishable cluster approximation and find that they enhance the accuracy of their parent ab initio methods.

Accurate Quantum Chemical Reaction Energies for Lithium-Mediated Electrolyte Decomposition and Evaluation of Density Functional Approximations.

An important concern related to the performance of Li-ion batteries is the formation of a solid electrolyte interphase on the surface of the anode. This film is formed from the decomposition of electrolytes and can have important effects on the stability and performance. Here, we evaluate the decomposition pathway of ethylene carbonate and related organic electrolyte molecules using a series of density functional approximations and correlated wave function (WF) methods, including the coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)] and auxiliary-field quantum Monte Carlo (AFQMC). We find that the transition state barrier associated with ring opening varies widely across different functionals, ranging from 3.01 to 17.15 kcal/mol, which can be compared to the value of 12.84 kcal/mol predicted by CCSD(T). This large variation underscores the importance of benchmarking against accurate WF methods. A performance comparison of all of the density functionals used in this study reveals that the M06-2X-D3 (a meta-hybrid GGA), CAM-B3LYP-D3 (a range-separated hybrid), and B2GP-PLYP-D3 (a double hybrid) perform the best, with average errors of about 1.50-1.60 kcal/mol compared to CCSD(T). We also compared the performance of the WF methods that are more scalable than CCSD(T), finding that DLPNO-CCSD(T) and phaseless AFQMC with a DFT trial wave function exhibit average errors of 1.38 and 1.74 kcal/mol, respectively.

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals.

Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

Dark-Exciton Driven Energy Funneling into Dielectric Inhomogeneities in Two-Dimensional Semiconductors.

The optoelectronic and transport properties of two-dimensional transition metal dichalcogenide semiconductors (2D TMDs) are highly susceptible to external perturbation, enabling precise tailoring of material function through postsynthetic modifications. Here, we show that nanoscale inhomogeneities known as nanobubbles can be used for both strain and, less invasively, dielectric tuning of exciton transport in bilayer tungsten diselenide (WSe2). We use ultrasensitive spatiotemporally resolved optical scattering microscopy to directly image exciton transport, revealing that dielectric nanobubbles are surprisingly efficient at funneling and trapping excitons at room temperature, even though the energies of the bright excitons are negligibly affected. Our observations suggest that exciton funneling in dielectric inhomogeneities is driven by momentum-indirect (dark) excitons whose energies are more sensitive to dielectric perturbations than bright excitons. These results reveal a new pathway to control exciton transport in 2D semiconductors with exceptional spatial and energetic precision using dielectric engineering of dark state energetic landscapes.

Full-frequency dynamical Bethe-Salpeter equation without frequency and a study of double excitations.

The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is O(N6), where N is the system size. To avoid the common static screening approximation, we show that the full-frequency dynamical BSE can be exactly reformulated as a frequency-independent eigenvalue problem in an expanded space of single and double excitations. When combined with an iterative eigensolver and the density fitting approximation to the electron repulsion integrals, this reformulation yields a dynamical BSE algorithm whose computational time is O(N5), which we verify numerically. Furthermore, the reformulation provides direct access to excited states with dominant double excitation character, which are completely absent in the spectrum of the statically screened BSE. We study the 21Ag state of butadiene, hexatriene, and octatetraene and find that GW/BSE overestimates the excitation energy by about 1.5-2 eV and significantly underestimates the double excitation character.

Accurate thermochemistry of covalent and ionic solids from spin-component-scaled MP2.

We study the performance of spin-component-scaled second-order Møller-Plesset perturbation theory (SCS-MP2) for the prediction of the lattice constant, bulk modulus, and cohesive energy of 12 simple, three-dimensional covalent and ionic semiconductors and insulators. We find that SCS-MP2 and the simpler scaled opposite-spin MP2 (SOS-MP2) yield predictions that are significantly improved over the already good performance of MP2. Specifically, when compared to experimental values with zero-point vibrational corrections, SCS-MP2 (SOS-MP2) yields mean absolute errors of 0.015 (0.017) Å for the lattice constant, 3.8 (3.7) GPa for the bulk modulus, and 0.06 (0.08) eV for the cohesive energy, which are smaller than those of leading density functionals by about a factor of two or more. We consider a reparameterization of the spin-scaling parameters and find that the optimal parameters for these solids are very similar to those already in common use in molecular quantum chemistry, suggesting good transferability and reliable future applications to surface chemistry on insulators.

Polytypism, Anisotropic Transport, and Weyl Nodes in the van der Waals Metal TaFeTe4.

Layered van der Waals (vdW) materials belonging to the MM'Te4 structure class have recently received intense attention due to their ability to host exotic electronic transport phenomena, such as in-plane transport anisotropy, Weyl nodes, and superconductivity. Here we report two new vdW materials with strongly anisotropic in-plane structures featuring stripes of metallic TaTe2 and semiconducting FeTe2, α-TaFeTe4 and ß-TaFeTe4. We find that the structure of α-TaFeTe4 produces strongly anisotropic in-plane electronic transport (anisotropy ratio of up to 250%), outcompeting all other vdW metals, and demonstrate that it can be mechanically exfoliated to the two-dimensional (2D) limit. We also explore the possibility that broken inversion symmetry in ß-TaFeTe4 produces Weyl points in the electronic band structure. Eight Weyl nodes slightly below the Fermi energy are computationally identified for ß-TaFeTe4, indicating they may contribute to the transport behavior of this polytype. These findings identify the TaFeTe4 polytypes as an ideal platform for investigation of 2D transport anisotropy and chiral charge transport as a result of broken symmetry.

Tunable Cr4+ Molecular Color Centers.

The inherent atomistic precision of synthetic chemistry enables bottom-up structural control over quantum bits, or qubits, for quantum technologies. Tuning paramagnetic molecular qubits that feature optical-spin initialization and readout is a crucial step toward designing bespoke qubits for applications in quantum sensing, networking, and computing. Here, we demonstrate that the electronic structure that enables optical-spin initialization and readout for S = 1, Cr(aryl)4, where aryl = 2,4-dimethylphenyl (1), o-tolyl (2), and 2,3-dimethylphenyl (3), is readily translated into Cr(alkyl)4 compounds, where alkyl = 2,2,2-triphenylethyl (4), (trimethylsilyl)methyl (5), and cyclohexyl (6). The small ground state zero field splitting values (<5 GHz) for 1-6 allowed for coherent spin manipulation at X-band microwave frequency, enabling temperature-, concentration-, and orientation-dependent investigations of the spin dynamics. Electronic absorption and emission spectroscopy confirmed the desired electronic structures for 4-6, which exhibit photoluminescence from 897 to 923 nm, while theoretical calculations elucidated the varied bonding interactions of the aryl and alkyl Cr4+ compounds. The combined experimental and theoretical comparison of Cr(aryl)4 and Cr(alkyl)4 systems illustrates the impact of the ligand field on both the ground state spin structure and excited state manifold, laying the groundwork for the design of structurally precise optically addressable molecular qubits.

Simulations of Trions and Biexcitons in Layered Hybrid Organic-Inorganic Lead Halide Perovskites.

Behaving like atomically precise two-dimensional quantum wells with non-negligible dielectric contrast, the layered hybrid organic-inorganic lead halide perovskites (HOIPs) have strong electronic interactions leading to tightly bound excitons with binding energies on the order of 500 meV. These strong interactions suggest the possibility of larger excitonic complexes like trions and biexcitons, which are hard to study numerically due to the complexity of the layered HOIPs. Here, we propose and parametrize a model Hamiltonian for excitonic complexes in layered HOIPs and we study the correlated eigenfunctions of trions and biexcitons using a combination of diffusion Monte Carlo and very large variational calculations with explicitly correlated Gaussian basis functions. Binding energies and spatial structures of these complexes are presented as a function of the layer thickness. The trion and biexciton of the thinnest layered HOIP have binding energies of 35 and 44 meV, respectively, whereas a single exfoliated layer is predicted to have trions and biexcitons with equal binding energies of 48 meV. We compare our findings to available experimental data and to that of other quasi-two-dimensional materials.

Full-frequency GW without frequency.

Efficient computer implementations of the GW approximation must approximate a numerically challenging frequency integral; the integral can be performed analytically, but doing so leads to an expensive implementation whose computational cost scales as O(N6), where N is the size of the system. Here, we introduce a new formulation of the full-frequency GW approximation by exactly recasting it as an eigenvalue problem in an expanded space. This new formulation (1) avoids the use of time or frequency grids, (2) naturally obviates the need for the common "diagonal" approximation, (3) enables common iterative eigensolvers that reduce the canonical scaling to O(N5), and (4) enables a density-fitted implementation that reduces the scaling to O(N4). We numerically verify these scaling behaviors and test a variety of approximations that are motivated by this new formulation. The new formulation is found to be competitive with conventional O(N4) methods based on analytic continuation or contour deformation. In this new formulation, the relation of the GW approximation to configuration interaction, coupled-cluster theory, and the algebraic diagrammatic construction is made especially apparent, providing a new direction for improvements to the GW approximation.

Tight distance-dependent estimators for screening two-center and three-center short-range Coulomb integrals over Gaussian basis functions.

We derive distance-dependent estimators for two-center and three-center electron repulsion integrals over a short-range Coulomb potential, erfc(ωr12)/r12. These estimators are much tighter than the ones based on the Schwarz inequality and can be viewed as a complement to the distance-dependent estimators for four-center short-range Coulomb integrals and for two-center and three-center full Coulomb integrals previously reported. Because the short-range Coulomb potential is commonly used in solid-state calculations, including those with the Heyd-Scuseria-Ernzerhof functional and with our recently introduced range-separated periodic Gaussian density fitting, we test our estimators on a diverse set of periodic systems using a wide range of the range-separation parameter ω. These tests demonstrate the robust tightness of our estimators, which are then used with integral screening to calculate periodic three-center short-range Coulomb integrals with linear scaling in system size.

Improving MP2 bandgaps with low-scaling approximations to EOM-CCSD.

Despite its reasonable accuracy for ground-state properties of semiconductors and insulators, second-order Møller-Plesset perturbation theory (MP2) significantly underestimates bandgaps. In this work, we evaluate the bandgap predictions of partitioned equation-of-motion MP2 (P-EOM-MP2), which is a second-order approximation to EOM coupled-cluster theory with single and double excitations. On a test set of elemental and binary semiconductors and insulators, we find that P-EOM-MP2 overestimates bandgaps by 0.3 eV on average, which can be compared to the underestimation by 0.6 eV on average exhibited by the G0W0 approximation with a Perdew-Burke-Ernzerhof reference. We show that P-EOM-MP2, when interpreted as a Green's function-based theory, has a self-energy that includes all first- and second-order diagrams and a few third-order diagrams. We find that the GW approximation performs better for materials with small gaps and P-EOM-MP2 performs better for materials with large gaps, which we attribute to their superior treatment of screening and exchange, respectively.

Fast periodic Gaussian density fitting by range separation.

We present an efficient implementation of periodic Gaussian density fitting (GDF) using the Coulomb metric. The three-center integrals are divided into two parts by range-separating the Coulomb kernel, with the short-range part evaluated in real space and the long-range part in reciprocal space. With a few algorithmic optimizations, we show that this new method-which we call range-separated GDF (RSGDF)-scales sublinearly to linearly with the number of k-points for small to medium-sized k-point meshes that are commonly used in periodic calculations with electron correlation. Numerical results on a few three-dimensional solids show about ten-fold speedups over the previously developed GDF with little precision loss. The error introduced by RSGDF is about 10-5Eh in the converged Hartree-Fock energy with default auxiliary basis sets and can be systematically reduced by increasing the size of the auxiliary basis with little extra work.

Vibrational heat-bath configuration interaction.

We introduce vibrational heat-bath configuration interaction (VHCI) as an accurate and efficient method for calculating vibrational eigenstates of anharmonic systems. Inspired by its origin in electronic structure theory, VHCI is a selected CI approach that uses a simple criterion to identify important basis states with a pre-sorted list of anharmonic force constants. Screened second-order perturbation theory and simple extrapolation techniques provide significant improvements to variational energy estimates. We benchmark VHCI on four molecules with 12-48 degrees of freedom and use anharmonic potential energy surfaces truncated at fourth and sixth orders. When compared to other methods using the same truncated potentials, VHCI produces vibrational spectra of tens or hundreds of states with sub-wavenumber accuracy at low computational cost.

Dynamical correlation energy of metals in large basis sets from downfolding and composite approaches.

Coupled-cluster theory with single and double excitations (CCSD) is a promising ab initio method for the electronic structure of three-dimensional metals, for which second-order perturbation theory (MP2) diverges in the thermodynamic limit. However, due to the high cost and poor convergence of CCSD with respect to basis size, applying CCSD to periodic systems often leads to large basis set errors. In a common "composite" method, MP2 is used to recover the missing dynamical correlation energy through a focal-point correction, but the inadequacy of finite-order perturbation theory for metals raises questions about this approach. Here, we describe how high-energy excitations treated by MP2 can be "downfolded" into a low-energy active space to be treated by CCSD. Comparing how the composite and downfolding approaches perform for the uniform electron gas, we find that the latter converges more quickly with respect to the basis set size. Nonetheless, the composite approach is surprisingly accurate because it removes the problematic MP2 treatment of double excitations near the Fermi surface. Using this method to estimate the CCSD correlation energy in the combined complete basis set and thermodynamic limits, we find that CCSD recovers 85%-90% of the exact correlation energy at rs = 4. We also test the composite approach with the direct random-phase approximation used in place of MP2, yielding a method that is typically (but not always) more cost effective due to the smaller number of orbitals that need to be included in the more expensive CCSD calculation.