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.

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.

Ligand field design enables quantum manipulation of spins in Ni2+ complexes.

Creating the next generation of quantum systems requires control and tunability, which are key features of molecules. To design these systems, one must consider the ground-state and excited-state manifolds. One class of systems with promise for quantum sensing applications, which require water solubility, are d8 Ni2+ ions in octahedral symmetry. Yet, most Ni2+ complexes feature large zero-field splitting, precluding manipulation by commercial microwave sources due to the relatively large spin-orbit coupling constant of Ni2+ (630 cm-1). Since low lying excited states also influence axial zero-field splitting, D, a combination of strong field ligands and rigidly held octahedral symmetry can ameliorate these challenges. Towards these ends, we performed a theoretical and computational analysis of the electronic and magnetic structure of a molecular qubit, focusing on the impact of ligand field strength on D. Based on those results, we synthesized 1, [Ni(ttcn)2](BF4)2 (ttcn = 1,4,7-trithiacyclononane), which we computationally predict will have a small D (Dcalc = +1.15 cm-1). High-field high-frequency electron paramagnetic resonance (EPR) data yield spin Hamiltonian parameters: gx = 2.1018(15), gx = 2.1079(15), gx = 2.0964(14), D = +0.555(8) cm-1 and E = +0.072(5) cm-1, which confirm the expected weak zero-field splitting. Dilution of 1 in the diamagnetic Zn analogue, [Ni0.01Zn0.99(ttcn)2](BF4)2 (1') led to a slight increase in D to ∼0.9 cm-1. The design criteria in minimizing D in 1via combined computational and experimental methods demonstrates a path forward for EPR and optical addressability of a general class of S = 1 spins.

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.

Toward Linear Scaling Auxiliary-Field Quantum Monte Carlo with Local Natural Orbitals.

We develop a local correlation variant of auxiliary-field quantum Monte Carlo (AFQMC) based on local natural orbitals (LNO-AFQMC). In LNO-AFQMC, independent AFQMC calculations are performed for each localized occupied orbital using a truncated set of tailored orbitals. Because the size of this space does not grow with the system size for a target accuracy, the method has linear scaling. Applying LNO-AFQMC to molecular problems containing a few hundred to a thousand orbitals, we demonstrate convergence of total energies with significantly reduced costs. The savings are more significant for larger systems and larger basis sets. However, even for our smallest system studied, we find that LNO-AFQMC is cheaper than canonical AFQMC, in contrast with many other reduced-scaling methods. Perhaps most significantly, we show that energy differences converge much more quickly than total energies, making the method ideal for applications in chemistry and material science. Our work paves the way for linear scaling AFQMC calculations of strongly correlated systems, which would have a transformative effect on ab initio quantum chemistry.

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.

Can Spin-Component Scaled MP2 Achieve kJ/mol Accuracy for Cohesive Energies of Molecular Crystals?

Attaining kJ/mol accuracy in cohesive energy for molecular crystals is a persistent challenge in computational materials science. In this study, we evaluate second-order Møller-Plesset perturbation theory (MP2) and its spin-component scaled models for calculating cohesive energies for 23 molecular crystals (X23 data set). Using periodic boundary conditions and Brillouin zone sampling, we converge results to the thermodynamic and complete basis set limits, achieving an accuracy of about 2 kJ/mol (0.5 kcal/mol), which is rarely achieved in previous MP2 calculations for molecular crystals. When compared to experimental data, our results have a mean absolute error of 12.9 kJ/mol, comparable to Density Functional Theory with the PBE functional and TS dispersion correction. By separately scaling the opposite-spin and same-spin correlation energy components, using predetermined parameters, we reduce the mean absolute error to 9.5 kJ/mol. Further fine-tuning of these scaling parameters specifically for the X23 data set brings the mean absolute error down to 7.5 kJ/mol.

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.

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.

Room-temperature wavelike exciton transport in a van der Waals superatomic semiconductor.

The transport of energy and information in semiconductors is limited by scattering between electronic carriers and lattice phonons, resulting in diffusive and lossy transport that curtails all semiconductor technologies. Using Re6Se8Cl2, a van der Waals (vdW) superatomic semiconductor, we demonstrate the formation of acoustic exciton-polarons, an electronic quasiparticle shielded from phonon scattering. We directly imaged polaron transport in Re6Se8Cl2 at room temperature, revealing quasi-ballistic, wavelike propagation sustained for a nanosecond and several micrometers. Shielded polaron transport leads to electronic energy propagation lengths orders of magnitude greater than in other vdW semiconductors, exceeding even silicon over a nanosecond. We propose that, counterintuitively, quasi-flat electronic bands and strong exciton-acoustic phonon coupling are together responsible for the transport properties of Re6Se8Cl2, establishing a path to ballistic room-temperature semiconductors.

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.

Mobile Trions in Electrically Tunable 2D Hybrid Perovskites.

2D hybrid perovskites are currently in the spotlight of material research for light-harvesting and -emitting applications. It remains extremely challenging, however, to externally control their optical response due to the difficulties of introducing electrical doping. Here, an approach of interfacing ultrathin sheets of perovskites with few-layer graphene and hexagonal boron nitride into gate-tunable, hybrid heterostructures, is demonstrated. It allows for bipolar, continuous tuning of light emission and absorption in 2D perovskites by electrically injecting carriers to densities as high as 1012  cm-2 . This reveals the emergence of both negatively and positively charged excitons, or trions, with binding energies up to 46 meV, among the highest measured for 2D systems. Trions are shown to dominate light emission and propagate with mobilities reaching 200 cm2 V-1 s-1 at elevated temperatures. The findings introduce the physics of interacting mixtures of optical and electrical excitations to the broad family of 2D inorganic-organic nanostructures. The presented strategy to electrically control the optical response of 2D perovskites highlights it as a promising material platform toward electrically modulated light-emitters, externally guided charged exciton currents, and exciton transistors based on layered, hybrid semiconductors.

Data-Efficient Machine Learning Potentials from Transfer Learning of Periodic Correlated Electronic Structure Methods: Liquid Water at AFQMC, CCSD, and CCSD(T) Accuracy.

Obtaining the atomistic structure and dynamics of disordered condensed-phase systems from first-principles remains one of the forefront challenges of chemical theory. Here we exploit recent advances in periodic electronic structure and provide a data-efficient approach to obtain machine-learned condensed-phase potential energy surfaces using AFQMC, CCSD, and CCSD(T) from a very small number (≤200) of energies by leveraging a transfer learning scheme starting from lower-tier electronic structure methods. We demonstrate the effectiveness of this approach for liquid water by performing both classical and path integral molecular dynamics simulations on these machine-learned potential energy surfaces. By doing this, we uncover the interplay of dynamical electron correlation and nuclear quantum effects across the entire liquid range of water while providing a general strategy for efficiently utilizing periodic correlated electronic structure methods to explore disordered condensed-phase systems.

Electric fields drive bond homolysis.

Electric fields have been used to control and direct chemical reactions in biochemistry and enzymatic catalysis, yet directly applying external electric fields to activate reactions in bulk solution and to characterize them ex situ remains a challenge. Here we utilize the scanning tunneling microscope-based break-junction technique to investigate the electric field driven homolytic cleavage of the radical initiator 4-(methylthio)benzoic peroxyanhydride at ambient temperatures in bulk solution, without the use of co-initiators or photochemical activators. Through time-dependent ex situ quantification by high performance liquid chromatography using a UV-vis detector, we find that the electric field catalyzes the reaction. Importantly, we demonstrate that the reaction rate in a field increases linearly with the solvent dielectric constant. Using density functional theory calculations, we show that the applied electric field decreases the dissociation energy of the O-O bond and stabilizes the product relative to the reactant due to their different dipole moments.

Interfacial electric fields catalyze Ullmann coupling reactions on gold surfaces.

The electric fields created at solid-liquid interfaces are important in heterogeneous catalysis. Here we describe the Ullmann coupling of aryl iodides on rough gold surfaces, which we monitor in situ using the scanning tunneling microscope-based break junction (STM-BJ) and ex situ using mass spectrometry and fluorescence spectroscopy. We find that this Ullmann coupling reaction occurs only on rough gold surfaces in polar solvents, the latter of which implicates interfacial electric fields. These experimental observations are supported by density functional theory calculations that elucidate the roles of surface roughness and local electric fields on the reaction. More broadly, this touchstone study offers a facile method to access and probe in real time an increasingly prominent yet incompletely understood mode of catalysis.

Full Configuration Interaction Excited-State Energies in Large Active Spaces from Subspace Iteration with Repeated Random Sparsification.

We present a stable and systematically improvable quantum Monte Carlo (QMC) approach to calculating excited-state energies, which we implement using our fast randomized iteration method for the full configuration interaction problem (FCI-FRI). Unlike previous excited-state quantum Monte Carlo methods, our approach, which is based on an asymmetric variant of subspace iteration, avoids the use of dot products of random vectors and instead relies upon trial vectors to maintain orthogonality and estimate eigenvalues. By leveraging recent advances, we apply our method to calculate ground- and excited-state energies of challenging molecular systems in large active spaces, including the carbon dimer with 8 electrons in 108 orbitals (8e,108o), an oxo-Mn(salen) transition metal complex (28e,28o), ozone (18e,87o), and butadiene (22e,82o). In the majority of these test cases, our approach yields total excited-state energies that agree with those from state-of-the-art methods─including heat-bath CI, the density matrix renormalization group approach, and FCIQMC─to within sub-milliHartree accuracy. In all cases, estimated excitation energies agree to within about 0.1 eV.

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.

Integral-Direct Hartree-Fock and Møller-Plesset Perturbation Theory for Periodic Systems with Density Fitting: Application to the Benzene Crystal.

We present an algorithm and implementation of integral-direct, density-fitted Hartree-Fock (HF) and second-order Møller-Plesset perturbation theory (MP2) for periodic systems. The new code eliminates the formerly prohibitive storage requirements and allows us to study systems 1 order of magnitude larger than before at the periodic MP2 level. We demonstrate the significance of the development by studying the benzene crystal in both the thermodynamic limit and the complete basis set limit, for which we predict an MP2 cohesive energy of -72.8 kJ/mol, which is about 10-15 kJ/mol larger in magnitude than all previously reported MP2 calculations. Compared to the best theoretical estimate from literature, several modified MP2 models approach chemical accuracy in the predicted cohesive energy of the benzene crystal and hence may be promising cost-effective choices for future applications on molecular crystals.

Ground-State Properties of Metallic Solids from Ab Initio Coupled-Cluster Theory.

Metallic solids are an enormously important class of materials, but they are a challenging target for accurate wave function-based electronic structure theories and have not been studied in great detail by such methods. Here, we use coupled-cluster theory with single and double excitations (CCSD) to study the structure of solid lithium and aluminum using optimized Gaussian basis sets. We calculate the equilibrium lattice constant, bulk modulus, and cohesive energy and compare them to experimental values, finding accuracy comparable to common density functionals. Because the quantum chemical "gold standard" CCSD(T) (CCSD with perturbative triple excitations) is inapplicable to metals in the thermodynamic limit, we test two approximate improvements to CCSD, which are found to improve the results.

Simplified GW/BSE Approach for Charged and Neutral Excitation Energies of Large Molecules and Nanomaterials.

Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [Grimme, S. J. Chem. Phys. 2013, 138, 244104], we describe a simplified approach to excited-state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from ab initio GW calculations by only 0.2 eV. The performance of sBSE (based on the sGW input) is tested on the excitation energies of molecules in Thiel's set, for which it differs from ab initio GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms.