Search

The promise of quantum computing to circumvent the exponential scaling of quantum chemistry has sparked a race to develop chemistry algorithms for quantum architecture. However, most works neglect the quantum-inherent shot noise, let alone the effect of current noisy devices. Here, we present a comprehensive study of quantum linear response (qLR) theory obtaining spectroscopic properties on simulated fault-tolerant quantum computers and present-day near-term quantum hardware. This work introduces novel metrics to analyze and predict the origins of noise in the quantum algorithm, proposes an Ansatz-based error mitigation technique, and highlights the significant impact of Pauli saving in reducing measurement costs and noise. Our hardware results using up to cc-pVTZ basis set serve as proof-of-principle for obtaining absorption spectra on quantum hardware in a general approach with the accuracy of classical multi-configurational methods. Importantly, our results exemplify that substantial improvements in hardware error rates and measurement speed are necessary to lift quantum computational chemistry from proof-of-concept to an actual impact in the field.

Calculating molecular properties using quantum devices can be done through the quantum linear response (qLR) or, equivalently, the quantum equation of motion (qEOM) formulations. Different parameterizations of qLR and qEOM are available, namely naive, projected, self-consistent, and state-transfer. In the naive and projected parameterizations, the metric is not the identity, and we show that it depends on the redundant orbital rotations. This dependency may lead to divergences in the excitation energies for certain choices of the redundant orbital rotation parameters in an idealized noise-less setting. Further, this leads to significant variance when calculations include statistical noise from finite quantum sampling.

The prediction of spectral properties via linear response (LR) theory is an important tool in quantum chemistry for understanding photo-induced processes in molecular systems. With the advances of quantum computing, we recently adapted this method for near-term quantum hardware using a truncated active space approximation with orbital rotation, named quantum linear response (qLR). In an effort to reduce the classic cost of this hybrid approach, we here derive and implement a reduced density matrix (RDM) driven approach of qLR. This allows for the calculation of spectral properties of moderately sized molecules with much larger basis sets than so far possible. We report qLR results for benzene and $R$-methyloxirane with a cc-pVTZ basis set and study the effect of shot noise on the valence and oxygen K-edge absorption spectra of H$_2$O in the cc-pVTZ basis.

We test the performance of the Polarizable Embedding Variational Quantum Eigensolver Self-Consistent-Field (PE-VQE-SCF) model for computing electric field gradients with comparisons to conventional complete active space self-consistent-field (CASSCF) calculations and experimental results. We compute quadrupole coupling constants for ice VIII and ice IX. We find that the inclusion of the environment is crucial for obtaining results that match the experimental data. The calculations for ice VIII are within the experimental uncertainty for both CASSCF and VQE-SCF for oxygen and lie close to the experimental value for ice IX as well. With the VQE-SCF, which is based on an Adaptive Derivative-Assembled Problem-Tailored (ADAPT) ansatz, we find that the inclusion of the environment and the size of the different basis sets do not directly affect the gate counts. However, by including an explicit environment, the wavefunction and, therefore, the optimization problem becomes more complicated, which usually results in the need to include more operators from the operator pool, thereby increasing the depth of the circuit., Comment: 18 pages, 5 figures

We explore Davidson methods for obtaining excitation energies and other linear response properties within quantum self-consistent linear response (q-sc-LR) theory. Davidson-type methods allow for obtaining only a few selected excitation energies without explicitly constructing the electronic Hessian since they only require the ability to perform Hessian-vector multiplications. We apply the Davidson method to calculate the excitation energies of hydrogen chains (up to H$_{10}$) and analyze aspects of statistical noise for computing excitation energies on quantum simulators. Additionally, we apply Davidson methods for computing linear response properties such as static polarizabilities for H$_2$, LiH, H$_2$O, OH$^-$, and NH$_3$, and show that unitary coupled cluster outperforms classical projected coupled cluster for molecular systems with strong correlation. Finally, we formulate the Davidson method for damped (complex) linear response, with application to the nitrogen K-edge X-ray absorption of ammonia, and the $C_6$ coefficients of H$_2$, LiH, H$_2$O, OH$^-$, and NH$_3$., Comment: 35 pages, 5 figures

Electronic quantum coherence is a prerequisite for charge migration in molecules and emerging molecular quantum technologies. It also underlies new forms of spectroscopy and can enhance molecular function. Whether conical intersections can create electronic coherence and for how long electronic coherence can survive in aqueous environments has remained elusive and controversial, respectively. Here, we use X-ray spectroscopy to realize a breakthrough on both of these topics. We find that electronic relaxation from the $^1$B$_{2\rm{u}}$($\pi\pi^*$) state of pyrazine through conical intersections (CIs) induces previously unobserved electronic and vibrational quantum coherences between the $^1$B$_{3\rm{u}}$(n$\pi^*$) and $^1$A$_{\rm{u}}$(n$\pi^*$) states in the gas phase, which correspond to an electronic ring current. These coherences are entirely suppressed when pyrazine is dissolved in water. These observations, supported by the latest advances in multiconfigurational electronic-structure calculations and non-adiabatic dynamics simulations, confirm that CIs can create electronic coherences and that aqueous solvation can decohere them in less than 40 fs. This study opens the door to the investigation of the broad class of electronic coherences created during light-induced molecular dynamics and to quantify their susceptibility to aqueous solvation., Comment: 54 pages, 28 figures

The X-ray absorption (XA) spectra of aqueous ammonia and ammonium are computed using a combination of coupled cluster singles and doubles (CCSD) with different quantum mechanical and molecular mechanical embedding schemes. Specifically, we compare frozen Hartree--Fock (HF) density embedding, polarizable embedding (PE), and polarizable density embedding (PDE). Integrating CCSD with frozen HF density embedding is possible within the CC-in-HF framework, which circumvents the conventional system-size limitations of standard coupled cluster methods. We reveal similarities between PDE and frozen HF density descriptions, while PE spectra differ significantly. By including approximate triple excitations, we also investigate the effect of improving the electronic structure theory. The spectra computed using this approach show an improved intensity ratio compared to CCSD-in-HF. Charge transfer analysis of the excitations shows the local character of the pre-edge and main-edge, while the post-edge is formed by excitations delocalized over the first solvation shell and beyond., Comment: 35 pages, 12 figures

Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photo-induced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory using a truncated active-space version of the multi-configurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates dubbed ``proj LRSD'' and ``all-proj LRSD''.

Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH$_2$ in the series of STO-3G/6-31G/6-31G* basis sets and of H$_4$ and H$_2$O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H$_4$, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems., Comment: 18+14 pages, 4 figures, 1 table; comments welcome

We formulate and implement the Variational Quantum Eigensolver Self Consistent Field (VQE-SCF) algorithm in combination with polarizable embedding (PE), thereby extending PE to the regime of quantum computing. We test the resulting algorithm, PE-VQE-SCF, on quantum simulators and demonstrate that the computational stress on the quantum device is only slightly increased in terms of gate counts compared to regular VQE-SCF. On the other hand, no increase in shot noise was observed. We illustrate how PE-VQE-SCF may lead to the modeling of real chemical systems using a simulation of the reaction barrier of the Diels-Alder reaction between furan and ethene as an example., Comment: 19 pages, 5 figures, submitted to Journal of Chemical Physics

We present the first successful application of the coupled cluster approach to simulate the X-ray absorption (XA) spectrum of liquid water. The system size limitations of standard coupled cluster theory are overcome by employing a newly developed coupled cluster method for large molecular systems. This method combines coupled cluster singles, doubles, and perturbative triples in a multilevel framework (MLCC3-in-HF) and is able to describe the delicate nature of intermolecular interactions in liquid water. Using molecular geometries from state-of-the-art path-integral molecular dynamics, we obtain excellent agreement with experimental spectra. Additionally, we show that an accurate description of the electronic structure within the first solvation shell is sufficient to model the XA spectrum of liquid water. Furthermore, we present a rigorous charge transfer analysis with unprecedented reliability, achieved through MLCC3-in-HF. This analysis aligns with previous studies regarding the character of the prominent features of the spectrum., Comment: 32 pages, 4+12 figures

We present the development and implementation of the relativistic coupled cluster linear response theory (CC-LR) which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin), and take into account the finite lifetime of excited states via damped response theory. We showcase our implementation, which is capable to offload intensive tensor contractions onto graphical processing units (GPUs), in the calculation of: \textit{(a)} frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a emphasis on the calculation of valence absorption cross-sections for the I$_2$ molecule;\textit{(b)} indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H$_2$Se-H$_2$O dimer as a prototypical system containing hydrogen bonds; and \textit{(c)} optical rotations at the sodium D line for hydrogen peroxide analogues (H$_{2}$Y$_{2}$, Y=O, S, Se, Te). Thanks to this implementation, we are able show the similarities in performance--but often the significant discrepancies--between CC-LR and approximate methods such as density functional theory (DFT). Comparing standard CC response theory with the equation of motion formalism, we find that, for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region such as spin-spin couplings, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table it may become increasingly difficult to measure pure optical rotation at the sodium D line, due to the appearance of absorbing states.

We present an implementation of a damped response framework for calculating resonant inelastic X-ray scattering (RIXS) at the equation-of-motion coupled cluster singles and doubles (CCSD) and second-order approximate coupled cluster singles and doubles (CC2) levels of theory in the open-source program $e^T$. This framework lays the foundation for future extension to higher excitation methods (notably, the coupled cluster singles and doubles with perturbative triples, CC3) and to multilevel approaches. Our implementation adopts a fully relaxed ground state, and different variants of the core-valence separation projection technique to address convergence issues. Illustrative results are compared with those obtained within the frozen-core core-valence separated approach, available in Q-Chem, as well as with experiment. The performance of the CC2 method is evaluated in comparison with that of CCSD. It is found that, while the CC2 method is noticeably inferior to CCSD for X-ray absorption spectra, the quality of the CC2 RIXS spectra is often comparable to that of the CCSD level of theory, when the same valence excited states are probed. Finally, we present preliminary RIXS results for a solvated molecule in aqueous solution.

We present an implementation of the B term of Magnetic Circular Dichroism within the Algebraic Diagrammatic Construction (ADC) scheme of the polarization propagator and its Intermediate State Representation. As illustrative results, the MCD spectra of the ADC variants ADC(2), ADC(2)-x and ADC(3) of the molecular systems uracil, 2-thiouracil, 4-thiouracil, purine, hypoxanthine, 1,4-naphthoquinone, 9,10-anthraquinone and 1-naphthylamine, are computed and compared with results obtained using the Coupled-Cluster Singles and Approximate Doubles (CC2) method, literature TD-DFT results as well as with available experimental data.

We present an efficient implementation of ground and excited state CCSD gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition, like density-fitting, is an inner projection method, and thus similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. The computational time spent per optimization cycle are furthermore found to be significantly lower compared to other implementations based on an inner projection method. We illustrate the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory., Comment: 12 pages, 1 figure

A time-dependent equation of motion coupled cluster singles and doubles (TD-EOM-CCSD) method is implemented, which uses a reduced basis calculated with the asymmetric band Lanczos algorithm. The approach is used to study weak-field processes in small molecules induced by ultrashort valence pump and core probe pulses. We assess the reliability of the procedure by comparing TD-EOM-CCSD absorption spectra to spectra obtained from the time-dependent coupled-cluster singles and doubles (TDCCSD) method and observe that spectral features can be reproduced for several molecules, at much lower computational times. We discuss how multiphoton absorption and symmetry can be handled in the method and general features of the core-valence separation (CVS) projection technique. We also model the transient absorption of an attosecond X-ray probe pulse by the glycine molecule., Comment: 15 pages and 5 figures

We assess the performance of different protocols for simulating excited-state X-ray absorption spectra. We consider three different protocols based on equation-of-motion coupled-cluster singles and doubles, two of them combined with the maximum overlap method. The three protocols differ in the choice of a reference configuration used to compute target states. Maximum-overlap-method time-dependent density functional theory is also considered. The performance of the different approaches is illustrated using uracil, thymine, and acetylacetone as benchmark systems. The results provide guidance for selecting an electronic structure method for modeling time-resolved X-ray absorption spectroscopy., Comment: Revised version submitted to Structural Dynamics (January 2021)

An implementation of a complex solver for the solution of the response equations required to compute the complex response functions of damped response theory is presented for the resolution-of-identity (RI) coupled-cluster singles and approximate doubles CC2 method. The implementation uses a partitioned formulation that avoids the storage of double excitation amplitudes to make it applicable to large molecules. The solver is the keystone element for the development of the damped coupled-cluster response formalism for linear and nonlinear effects in resonant frequency regions at the RI-CC2 level of theory. Illustrative results are reported for the one-photon absorption cross section of C60, the electronic circular dichroism of $n$-helicenes ($n$ = 5, 6, 7), and the $C_6$ dispersion coefficients of a set of selected organic molecules and fullerenes., Comment: Submitted to J. Chem. Phys., Dec. 2020

We report an implementation of the core-valence separation approach to the 4-component relativistic Hamiltonian based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD), for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants, and the use of different Hamiltonians based on the exact 2-component (X2C) framework, on the energies of different core ionized and excited states in halogen (CH$_3$I, HX and X$^-$, X = Cl-At) and xenon containing (Xe, XeF$_2$) species. Our results show that the X2C molecular mean-field approach [Sikkema et al., J. Chem. Phys. 2009, 131, 124116], based on 4-component Dirac-Coulomb mean-field calculations ($^2$DC$^M$) is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference 4-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets yield non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine, and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations ($^2$DCG$^M$) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals., Comment: 16 pages, 4 figures in main text, 7 figures in supplemental information

We present the first investigation of excited state dynamics by resonant Auger-Meitner spectroscopy (also known as resonant Auger spectroscopy) using the nucleobase thymine as an example. Thymine is photoexcited in the UV and probed with X-ray photon energies at and below the oxygen K-edge. After initial photoexcitation to a {\pi}{\pi}* excited state, thymine is known to undergo internal conversion to an n{\pi}* excited state with a strong resonance at the oxygen K-edge, red-shifted from the ground state {\pi}* resonances of thymine (see our previous study Wolf et al., Nat. Commun., 2017, 8, 29). We resolve and compare the Auger-Meitner electron spectra associated both with the excited state and ground state resonances, and distinguish participator and spectator decay contributions. Furthermore, we observe simultaneously with the decay of the n{\pi}* state signatures the appearance of additional resonant Auger-Meitner contributions at photon energies between the n{\pi}* state and the ground state resonances. We assign these contributions to population transfer from the n{\pi}* state to a {\pi}{\pi}* triplet state via intersystem crossing on the picosecond timescale based on simulations of the X-ray absorption spectra in the vibrationally hot triplet state. Moreover, we identify signatures from the initially excited {\pi}{\pi}* singlet state which we have not observed in our previous study.

Vibrationally resolved near-edge x-ray absorption spectra at the K-edge for a number of small molecules have been computed from anharmonic vibrational configuration interaction calculations of the Franck-Condon factors. The potential energy surfaces for ground and core-excited states were obtained at the core-valence separated CC2, CCSD, CCSDR(3), and CC3 levels of theory, employing the Adaptive Density-Guided Approach (ADGA) scheme to select the single points at which to perform the energy calculations. We put forward an initial attempt to include pair-mode coupling terms to describe the potential of polyatomic molecules

A computational expression for the Faraday A term of magnetic circular dichroism (MCD) is derived within coupled cluster response theory and alternative computational expressions for the B term are discussed. Moreover, an approach to compute the (temperature-independent) MCD ellipticity in the context of coupled cluster damped response is presented, and its equivalence with the stick-spectrum approach in the limit of infinite lifetimes is demonstrated. The damped response approach has advantages for molecular systems or spectral ranges with a high density of states. Illustrative results are reported at the coupled cluster singles and doubles level and compared to time-dependent density functional theory results., Comment: Submitted to J. Chem. Phys. on May 10, 2020

The eT program is an open source electronic structure package with emphasis on coupled cluster and multilevel methods. It includes efficient spin adapted implementations of ground and excited singlet states, as well as equation of motion oscillator strengths, for CCS, CC2, CCSD, and CC3. Furthermore, eT provides unique capabilities such as multilevel Hartree-Fock and multilevel CC2, real-time propagation for CCS and CCSD, and efficient CC3 oscillator strengths. With a coupled cluster code based on an efficient Cholesky decomposition algorithm for the electronic repulsion integrals, eT has similar advantages as codes using density fitting, but with strict error control. Here we present the main features of the program and demonstrate its performance through example calculations. Because of its availability, performance, and unique capabilities, we expect eT to become a valuable resource to the electronic structure community., Comment: 31 pages and 10 figures - supplementary information included in the uploaded files

An extensive analysis has been carried out of the performance of standard families of basis sets with the hierarchy of coupled cluster methods CC2, CCSD, CC3 and CCSDT in computing selected Oxygen, Carbon and Nitrogen K-edge (vertical) core excitation and ionization energies within a core-valence separated scheme in the molecules water, ammonia, and carbon monoxide. Complete basis set limits for the excitation energies have been estimated via different basis set extrapolation schemes. The importance of scalar relativistic effects has been established within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e)., Comment: Manuscript submitted to Adv. Quantum Chemistry

The circular and linear magnetic birefringences corresponding to the Faraday and the Cotton-Mouton effects, respectively, have been measured in xenon at $\lambda = 1064$ nm. The experimental setup is based on time dependent magnetic fields and a high finesse Fabry-Perot cavity. Our value of the Faraday effect is the first measurement at this wavelength. It is compared to theoretical predictions. Our uncertainty of a few percent yields an agreement at better than 1$\sigma$ with the computational estimate when relativistic effects are taken into account. Concerning the Cotton-Mouton effect, our measurement, the second ever published at $\lambda = 1064$ nm, agrees at better than $1\sigma$ with theoretical predictions. We also compare our error budget with those established for other experimental published values.