Accurate Wave Functions Research Articles

Gaussians for Electronic and Rovibrational Quantum Dynamics.

The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wave packets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to treat the electronic and nuclear degrees of freedom on equal quantum-mechanical footing by not invoking the BO approximation at all. Explicitly correlated Gaussian (ECG) basis functions have proved successful for non-BO calculations of stationary molecular states and energies, reproducing rovibrational absorption spectra with very high accuracy. In this Article, we present a proof-of-principle study of the ability of fully flexible ECGs (FFECGs) to capture the intricate electronic and rovibrational dynamics generated by short, high-intensity laser pulses. By fitting linear combinations of FFECGs to accurate wave function histories obtained on a large real-space grid for a regularized 2D model of the hydrogen atom and for the 2D Morse potential, we demonstrate that FFECGs provide a very compact description of laser-driven electronic and rovibrational dynamics.

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models.

Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.

Aufbau Suppressed Coupled Cluster Theory for Electronically Excited States.

We introduce an approach to improve single-reference coupled cluster theory in settings where the Aufbau determinant is absent from or plays only a small role in the true wave function. Using a de-excitation operator that can be efficiently hidden within a similarity transform, we create a coupled cluster wave function in which de-excitations work to suppress the Aufbau determinant and produce wave functions dominated by other determinants. Thanks to an invertible and fully exponential form, the approach is systematically improvable, size consistent, size extensive, and, interestingly, size intensive in a granular way that should make the adoption of some ground state techniques, such as local correlation, relatively straightforward. In this initial study, we apply the general formalism to create a state-specific method for orbital-relaxed, singly excited states. We find that this approach matches the accuracy of similar-cost equation-of-motion methods in valence excitations while offering improved accuracy for charge transfer states. We also find the approach to be more accurate than excited-state-specific perturbation theory in both types of states.

Accurate Quantum States for a 2D-Dipole.

Edge dislocations are crucial in understanding both mechanical and electrical transport in solid and are modeled as line distributions of dipole moments. The calculation of the electronic spectrum for the two dimensional dipole, represented by the potential energy V(r,θ)=pcosθ/r, has been the topic of several studies that show significant difficulties in obtaining accurate results. In this work, we demonstrate that the source of these difficulties is a logarithmic contribution to the behavior of the wave function at the origin that was neglected by previous authors. By taking into account this non-analytic deviation of the solution of Schrödinger's equation, superior results, with the expected rate of convergence, are obtained. This goal is accomplished by "adapting" general algorithms for solving partial derivative differential equations to include the desired asymptotic behavior. We illustrate this principle for the variational principle and finite difference methods. Accurate energies and wave functions are obtained not only for the ground state but also for the first eleven excited states and are useful for designing nanoelectronic devices. This paper demonstrates that augmentary knowledge about analytic properties of the solutions leads to the improved convergence and stability of numerical methods.

Towards a transferable fermionic neural wavefunction for molecules.

Deep neural networks have become a highly accurate and powerful wavefunction ansatz in combination with variational Monte Carlo methods for solving the electronic Schrödinger equation. However, despite their success and favorable scaling, these methods are still computationally too costly for wide adoption. A significant obstacle is the requirement to optimize the wavefunction from scratch for each new system, thus requiring long optimization. In this work, we propose a neural network ansatz, which effectively maps uncorrelated, computationally cheap Hartree-Fock orbitals, to correlated, high-accuracy neural network orbitals. This ansatz is inherently capable of learning a single wavefunction across multiple compounds and geometries, as we demonstrate by successfully transferring a wavefunction model pre-trained on smaller fragments to larger compounds. Furthermore, we provide ample experimental evidence to support the idea that extensive pre-training of such a generalized wavefunction model across different compounds and geometries could lead to a foundation wavefunction model. Such a model could yield high-accuracy ab-initio energies using only minimal computational effort for fine-tuning and evaluation of observables.

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.

Benchmarking DFT-based excited-state methods for intermolecular charge-transfer excitations.

Intermolecular charge-transfer is a highly important process in biology and energy-conversion applications where generated charges need to be transported over several moieties. However, its theoretical description is challenging since the high accuracy required to describe these excited states must be accessible for calculations on large molecular systems. In this benchmark study, we identify reliable low-scaling computational methods for this task. Our reference results were obtained from highly accurate wavefunction calculations that restrict the size of the benchmark systems. However, the density-functional theory based methods that we identify as accurate can be applied to much larger systems. Since targeting charge-transfer states requires the unambiguous classification of an excited state, we first analyze several charge-transfer descriptors for their reliability concerning intermolecular charge-transfer and single out the charge-transfer distance calculated based on the variation of electron density upon excitation (DCT) as an optimal choice for our purposes. In general, best results are obtained for orbital-optimized methods and among those, the maximum overlap method proved to be the most numerically stable variant when using the initial MOs as reference orbitals. Favorable error cancellation with optimally-tuned range-separated hybrid functionals and a rather small basis set can provide an economical yet reasonable wavefunction when using time-dependent density functional theory, which provides relevant information about the excited-state character to be used in the orbital-optimized methods. The qualitative agreement makes these fast calculations attractive for high-throughput screening applications.

Exit wave function reconstruction from two defocus images using neural network

Wave function reconstruction from one or two defocus images is promising for live atomic resolution imaging in transmission electron microscopy. However, a robust and accurate reconstruction method we still need more attention. Here, we present a neural-network-based wave function reconstruction method, EWR-NN, that enables accurate wave function reconstruction from only two defocus images. Results from both simulated and two different experimental defocus series show that the EWR-NN method has better performance than the widely-used iterative wave function reconstruction (IWFR) method. Influence of image number, defocus deviation, residual image shifts and noise level were considered to validate the performance of EWR-NN under practical conditions. It is seen that these factors will not influence the arrangement of atom columns in the reconstructed phase images, while they can alter the absolute values of all-atom columns and degrade the contrast of the phase images.

Photoionization from the Ground and Excited Vibrational States of and Its Deuterated Isotopologues

Photoionization cross sections and rate coefficients have been calculated for all bound vibrational levels of the 1sσ g state of , HD+, and . The Born–Oppenheimer approximation is employed in our calculation of vibrationally resolved photoionization cross sections. Vibrationally resolved and local thermal equilibrium photoionization rate coefficients are presented for photon temperatures less than 50,000 K and are found to be several orders of magnitude larger than previous results in the literature. Analytic fits for the vibrationally resolved and local thermal equilibrium photoionization rate coefficients are provided. Near-threshold oscillations in the vibrationally resolved photoionization are observed. A benchmark set of photoionization cross sections are presented. Fixed-nuclei photoionization cross sections are calculated using two-center true continuum wave functions and are verified by comparison with previous calculations and are found to be in excellent agreement in all cases. Data files for our set of benchmark cross sections, rate coefficients, and fitting parameters for , HD+, and are available on Zenodo under an open-source Creative Commons Attribution license at doi:10.5281/zenodo.8304061.

The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures.

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.

Coupling properties of subjects in a three-body interaction system with an intermediary entity and symmetrical participants on both sides

This paper studies the phenomenon of light-excited entanglement between nucleon spins, which is realized by using a transient light beam to excite the spin of electrons to produce entanglement between the two nucleon spins. This standard three-body interaction system is composed of an intermedium equivalent to a spin-1 quasiparticle and two spin-1/2 nucleons. We describe the entire system, its corresponding effective Hamiltonian, and the eigenvalues and eigenstates. The Hamiltonian is analyzed under the symmetric assumption, which can significantly simplify the effective Hamiltonian and make obtaining the corresponding dynamic characteristics easier. In the calculation process, degenerate perturbation theory calculates the eigenequation, whose degeneracy is not eliminated to a first approximation. We calculate the eigenequation to the second order to obtain an accurate zero-order wave function.

Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules.

We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set, close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO2, and O2 molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol.

Quantitative Prediction of Vertical Ionization Potentials from DFT via a Graph-Network-Based Delta Machine Learning Model Incorporating Electronic Descriptors.

While accurate wave function theories like CCSD(T) are capable of modeling molecular chemical processes, the associated steep computational scaling renders them intractable for treating large systems or extensive databases. In contrast, density functional theory (DFT) is much more computationally feasible yet often fails to quantitatively describe electronic changes in chemical processes. Herein, we report an efficient delta machine learning (ΔML) model that builds on the Connectivity-Based Hierarchy (CBH) scheme─an error correction approach based on systematic molecular fragmentation protocols─and achieves coupled cluster accuracy on vertical ionization potentials by correcting for deficiencies in DFT. The present study integrates concepts from molecular fragmentation, systematic error cancellation, and machine learning. First, we show that by using an electron population difference map, ionization sites within a molecule may be readily identified, and CBH correction schemes for ionization processes may be automated. As a central feature of our work, we employ a graph-based QM/ML model, which embeds atom-centered features describing CBH fragments into a computational graph to further increase accuracy for the prediction of vertical ionization potentials. In addition, we show that the incorporation of electronic descriptors from DFT, namely electron population difference features, improves model performance well beyond chemical accuracy (1 kcal/mol) to approach benchmark accuracy. While the raw DFT results are strongly dependent on the underlying functional used, for our best models, the performance is robust and much less dependent on the functional.

Elastic scattering of electrons by water: An ab initio study

In this work we devise a theoretical and computational method to compute the elastic scattering of electrons from a non-spherical potential, such as in the case of molecules and molecular aggregates. Its main feature is represented by the ability of calculating accurate wave functions for continuum states of polycentric systems via the solution of the Lippmann-Schwinger equation, including both the correlation effects and multi-scattering interference terms, typically neglected in widely used approaches, such as the Mott theory. Within this framework, we calculate the purely elastic scattering matrix elements. As a test case, we apply our scheme to the modelling of electron-water elastic scattering. The Dirac-Hartree-Fock self-consistent field method is used to determine the non-spherical molecular potential projected on a functional space spanned by Gaussian basis set. By adding a number of multi-centric radially-arranged s-type Gaussian functions, whose exponents are system-dependent and optimized to reproduce the properties of the continuum electron wave function in different energy regions, we are able to achieve unprecedented access to the description of the low energy range of the spectrum (0.001 < E < 10 eV) up to keV, finding a good agreement with experimental data and previous theoretical results. To show the potential of our approach, we also compute the total elastic scattering cross section of electrons impinging on clusters of water molecules and zundel cation. Our method can be extended to deal with inelastic scattering events and heavy-charged particles.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions.

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

Revisiting the generalized pseudospectral method: Energy degeneracies, radial expectation values, and Schwarz‐like inequalities of the shell‐confined atom

AbstractWe revisit the generalized pseudospectral (GPS) method to investigate the physical properties of shell‐confined atoms. Based on a general mapping function, the GPS method is developed for obtaining highly accurate bound state energies, wave functions, and radial quantities for the shell‐confined hydrogen atom. Besides the incidental and simultaneous degeneracies in the energy spectrum, we find that energy degeneracy appears very commonly in the shell‐confined system. The contour maps of the eigenenergies for some low‐lying bound states are obtained and the corresponding wave functions are analyzed. Radial expectation values with both positive and negative powers are calculated with high accuracy. By analyzing the Schwarz‐like inequalities, we show that the ground state density function for the outer‐confined hydrogen atom is second‐order monotonic, while for the shell‐confined hydrogen atom the radial densities for all bound states are zeroth‐order monotonic. In the extreme situation when the inner radius coalesces with the outer radius, the electron becomes a classical particle.

Quantum-Chemical and Quantum-Graph Models of the Dynamical Structure of CH5.

Experimental and computational results about the structure, dynamics, and rovibrational spectra of protonated methane have challenged a considerable number of traditional chemical concepts. Hereby theoretical and computational results are provided about the dynamical structure of CH5+. It is shown that the ground vibrational state investigated thus far by computations, forbidden by nuclear-spin statistics, has a structure similar to the first allowed vibrational state and, in fact, the structures of all vibrational states significantly below 200 cm-1 are highly similar. Spatial delocalization of the nuclei, determined by nuclear densities computed from accurate variational vibrational wave functions, turns out to be limited when viewed in the body-fixed frame, confirming that the effective structure of CH5+ is well described as a CH3+ tripod with a H2 unit on top of it. The interesting and unusual qualitative aspects of the sophisticated state-dependent variational results receive full explanation via simple quantum-graph models.

Variational tensor network operator

We propose a simple and generic construction of the variational tensor network operators to study the quantum spin systems by the synergy of ideas from the imaginary-time evolution and the variational optimization of trial wave functions. By applying these operators to simple initial states, accurate variational ground-state wave functions with extremely few parameters can be obtained. Furthermore, the framework can be applied to spontaneously study symmetry-breaking, symmetry-protected topological, and intrinsic topologically ordered phases, and we show that symmetries of the local tensors associated with these phases can emerge directly after the optimization without any gauge fixing. This provides a universal way to identify quantum phase transitions without prior knowledge of the system.

Electrochemical Hydrogenation of CO on Cu(100): Insights from Accurate Multiconfigurational Wavefunction Methods.

Copper (Cu) remains the most efficacious electrocatalyst for electrochemical CO2 reduction (CO2R). Its activity and selectivity are highly facet-dependent. We recently examined the commonly proposed rate-limiting CO hydrogenation step on Cu(111) via embedded correlated wavefunction (ECW) theory and demonstrated that only this higher-level theory yields predictions consistent with potential-dependent experimental kinetics. Here, to understand the differing activities of Cu(111) and Cu(100) in catalyzing CO2R, we explore CO hydrogenation on Cu(100) using ECW theory. We predict that the preferred pathway involves the reduction of adsorbed CO (*CO) to *COH via proton-coupled electron transfer (PCET) at working potentials, although *CHO also may form with a kinetically accessible but higher barrier. In contrast, our earlier work on Cu(111) concluded that *COH and *CHO formation via PCET are equally feasible. This work illustrates one possible origin of the facet dependence of CO2R mechanisms and products on Cu electrodes and sheds light on how the selectivity of CO2R electrocatalysts can be controlled by the surface morphology.

Precision Møller polarimetry for PREX-2 and CREX

The PREX-2 and CREX experiments in Hall A at Jefferson Lab are precision measurements of parity violating elastic electron scattering from complex nuclei. One requirement was that the incident electron beam polarization, typically ≈90%, be known with 1% precision. We commissioned and operated a Møller polarimeter on the beam line that exceeds this requirement, achieving a precision of 0.89% for PREX-2, and 0.85% for CREX. The uncertainty is purely systematic, accumulated from several different sources, but dominated by our knowledge of the target polarization. Our analysis also demonstrates the need for accurate atomic wave functions in order to correct for the Levchuk Effect. We describe the details of the polarimeter operation and analysis, as well as (for CREX) a comparison to results from a different polarimeter based on Compton scattering.