Multireference Perturbation Theory Research Articles

Perturbation theory in the complete degenerate active space (CDAS-PT2).

Methods based on the multireference perturbation theory (MRPT) with the one-electron zeroth-order Hamiltonian are widely used for the description of excited states, for example, due to their relatively low computational cost. However, current methods have a common drawback-use of a model space with low size. In this article, we propose the MRPT method with the model space extended to the complete active space. The one-electron zeroth-order Hamiltonian suitable for this extension is formulated. The proposed method was applied to common models, such as LiF, ethylene, and trans-butadiene. It was shown to have accuracy superior to XMCQDPT2 in most cases, especially in the case of the small active space.

Non-adiabatic electronic relaxation of tetracene from its brightest singlet excited state.

The ultrafast relaxation dynamics of tetracene following UV excitation to the bright singlet state S6 has been studied with time-resolved photoelectron spectroscopy. With the help of high-level abinitio multireference perturbation theory calculations, we assign photoelectron signals to intermediate dark electronic states S3, S4, and S5 as well as to a low-lying electronic state S2. The energetic structure of these dark states has not been determined experimentally previously. The time-dependent photoelectron yields assigned to the states S6, S5, and S4 have been analyzed and reveal the depopulation of S6 within 60fs, while S5 and S4 are populated with delays of about 50 and 80fs. The dynamics of the lower-lying states S3 and S2 seem to agree with a delayed population coinciding with the depopulation of the higher-lying states S4-S6 but could not be elucidated in full detail due to the low signal levels of the corresponding two-photon ionization probe processes.

Non-adiabatic electronic relaxation of tetracene from its brightest singlet excited state.

The ultrafast relaxation dynamics of tetracene following UV excitation to the bright singlet state S6 has been studied with time-resolved photoelectron spectroscopy. With the help of high-level abinitio multireference perturbation theory calculations, we assign photoelectron signals to intermediate dark electronic states S3, S4, and S5 as well as to a low-lying electronic state S2. The energetic structure of these dark states has not been determined experimentally previously. The time-dependent photoelectron yields assigned to the states S6, S5, and S4 have been analyzed and reveal the depopulation of S6 within 60fs, while S5 and S4 are populated with delays of about 50 and 80fs. The dynamics of the lower-lying states S3 and S2 seem to agree with a delayed population coinciding with the depopulation of the higher-lying states S4-S6 but could not be elucidated in full detail due to the low signal levels of the corresponding two-photon ionization probe processes.

Towards reliable and efficient modeling of [Cu2O2]2+-based compound electronic structures with the partially fixed reference space protocols.

This work reports a computationally efficient approach for reliable modeling of complex electronic structures based on [Cu2O2]2+ moieties. Specifically, we explore the recently developed partially fixed reference space (PFRS) protocol to minimize the active space size, taking into account the double d-shell effects. We show that the ground-state electronic structure of the core [Cu2O2]2+ model system is dominated by the d9/d10 occupations. The PFRS-crafted active spaces are further used to generate the reference wave functions for the multi-reference coupled cluster, configuration interaction, and multi-reference perturbation theory calculations. Specifically, we demonstrate that the bare [Cu2O2]2+ core can be modeled qualitatively using active spaces as small as CAS(2,2)PFRS. To obtain quantitative agreement with the reference DMRG(32,62)CI calculations, the CAS(4,4) has to be used in conjunction with the MRCCSD correction on top of it. This reliable and computationally efficient protocol is further used to model the electronic structure and properties of ammonia coordinated [Cu2O2]2+ complexes. Finally, based on the large amount of available experimental data regarding the oxo-peroxo equilibrium of [Cu2O2]2+-based systems, it is possible to formulate educated guesses regarding the effect of each experimental variable over each d-occupancy-specific state. With a large sample size of d-occupancy-specific state dependence with ligands and solvents, it should be possible to propose new ligands with specific d-occupancy and, therefore, oxidative properties based on the d-occupancy energy gaps of relatively low-cost calculations.

Toward Accurate Spin-Orbit Splittings from Relativistic Multireference Electronic Structure Theory.

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. However, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

Axial Coordinated Manganese(III) Porphyrin/Tetraazaporphyrin - 4-(10-phenylanthracen-9-yl)Pyridine Dyads: Self-Assembly, Structure and Spectral Properties in Ground and Excited States.

Self-assembly of new donor-acceptor systems based on (5,10,15,20-tetraphenylporphinato)manganese(III)/(5,10,15,20-tetra-4-tert-butylphenylporphinato)manganese(III)/(octakis(4-tert-butylphenyl)tetraazaporphinato)manganese(III) acetate ((AcO)MnTPP/(AcO)MnTBPP/(AcO)MnTAP) and 4-(10-phenylanthracen-9-yl)pyridine (PyAn) was studied using fluorescence spectroscopy and mass spectrometry. It was found that the coordination complexes of 1 : 1 composition (dyads) are formed in toluene. The spectral properties, the chemical structures and redox behavior of the dyads were described using 1H NMR, IR, ESR spectroscopy and cyclic voltammetry, respectively. The dynamic processes and the characteristics in the excited state of the dyads were obtained using the femtosecond transient absorption spectroscopy method. Density functional theory (DFT), time-dependent DFT methods were used to elucidate the dyad electronic structures and to establish the differences in their frontier molecular orbitals. The analysis of the lambda parameter and the distance of hole-pair interaction was indicated more favorable charge transfer between the macrocycle and the axial PyAn fragment in (AcO)(PyAn)MnTAP. The calculated values of the zero-field splitting parameters D and E/D, together with the g tensors of the lowest spin-orbit state for (AcO)MnTPP and (AcO)(PyAn)MnTPP were obtained using the combination of DFT and Multireference Perturbation Theory (CASSCF/NEVPT2) simulations. The data obtained develop the fundamental basis in the field of photovoltaics and show the prospects for the study of molecular systems of this class.

Ab Initio Calculation of UV-vis Absorption of Parent Mg, Fe, Co, Ni, Cu, and Zn Metalloporphyrins.

Relativistic restricted active space (RAS) second-order multireference perturbation theory (MRPT2) methods, incorporating spin-orbit (SO) coupling perturbatively via state interaction (SO-MRPT2/RASSCF), were used to reproduce the absorption spectra of parent metalloporphyrins containing the Mg2+, Zn2+, Co2+, Ni2+, Cu2+, or FeCl2+ ions in the 12,500-40,000 cm-1 region. Particular attention was paid to the interaction between the porphyrin ring and the metal 3d electrons in states of different multiplicities (we used metal 3d and double d-shell or 3d' orbitals). For this class of compounds, the N-electron valence state perturbation theory (NEVPT2) method is superior to the complete active space perturbation theory (CASPT2) and successfully reproduces the energies of all four characteristic transitions (Q, B, N, and L) of closed-shell metalloporphyrins. Inclusion of SO coupling was found to have very little effect on excitation energies and oscillator strengths. For FeCl2+ porphyrin, we treated ligand-to-metal charge-transfer (LMCT; π,d), metal ligand field (d,d), and metal-to-ligand charge-transfer (MLCT; d,π*) transitions within the same framework. The broad and intense spectral features associated with its B (Soret) band are attributed to multiconfigurational LMCT (d,π*) bands involving strong metal-ligand orbital mixing.

Perturbative variational quantum algorithms for material simulations

Reducing circuit depth is essential for implementing quantum simulations of electronic structure on near-term quantum devices. In this work, we propose a variational quantum eigensolver (VQE) based perturbation theory (PT) algorithm to accurately simulate electron correlation of periodic materials with shallow ansatz circuits, which are generated from adaptive derivative-assembled pseudo-trotter or qubit-excitation-based VQE calculations using a loose convergence criteria. Here, the major part of the electron correlation is described using the VQE ansatz circuit and the remaining correlation energy is described by either multireference or similarity transformation-based PT. Numerical results demonstrate that the new algorithms are able to accurately describe electron correlation of the LiH crystal with only one circuit parameter, in contrast with ~30 parameters required in the adaptive VQE to achieve the same accuracy. Meanwhile, for fixed-depth ansätze, e.g. unitary coupled cluster, we demonstrate that the VQE-based PT provides an appealing scheme to improve their accuracy.

Bistability between π-diradical open-shell and closed-shell states in indeno[1,2-a]fluorene.

Indenofluorenes are non-benzenoid conjugated hydrocarbons that have received great interest owing to their unusual electronic structure and potential applications in nonlinear optics and photovoltaics. Here we report the generation of unsubstituted indeno[1,2-a]fluorene on various surfaces by the cleavage of two C-H bonds in 7,12-dihydroindeno[1,2-a]fluorene through voltage pulses applied by the tip of a combined scanning tunnelling microscope and atomic force microscope. On bilayer NaCl on Au(111), indeno[1,2-a]fluorene is in the neutral charge state, but it exhibits charge bistability between neutral and anionic states on the lower-workfunction surfaces of bilayer NaCl on Ag(111) and Cu(111). In the neutral state, indeno[1,2-a]fluorene exhibits one of two ground states: an open-shell π-diradical state, predicted to be a triplet by density functional and multireference many-body perturbation theory calculations, or a closed-shell state with a para-quinodimethane moiety in the as-indacene core. We observe switching between open- and closed-shell states of a single molecule by changing its adsorption site on NaCl.

Towards reliable and efficient modeling of [Cu2O2]2+-based compound electronic structures with the partially fixed reference space protocols.

This work reports a computationally efficient approach for reliable modeling of complex electronic structures based on [Cu2O2]2+ moieties. Specifically, we explore the recently developed partially fixed reference space (PFRS) protocol to minimize the active space size, taking into account the double d-shell effects. We show that the ground-state electronic structure of the core [Cu2O2]2+ model system is dominated by the d9/d10 occupations. The PFRS-crafted active spaces are further used to generate the reference wave functions for the multi-reference coupled cluster, configuration interaction, and multi-reference perturbation theory calculations. Specifically, we demonstrate that the bare [Cu2O2]2+ core can be modeled qualitatively using active spaces as small as CAS(2,2)PFRS. To obtain quantitative agreement with the reference DMRG(32,62)CI calculations, the CAS(4,4) has to be used in conjunction with the MRCCSD correction on top of it. This reliable and computationally efficient protocol is further used to model the electronic structure and properties of ammonia coordinated [Cu2O2]2+ complexes. Finally, based on the large amount of available experimental data regarding the oxo-peroxo equilibrium of [Cu2O2]2+-based systems, it is possible to formulate educated guesses regarding the effect of each experimental variable over each d-occupancy-specific state. With a large sample size of d-occupancy-specific state dependence with ligands and solvents, it should be possible to propose new ligands with specific d-occupancy and, therefore, oxidative properties based on the d-occupancy energy gaps of relatively low-cost calculations.

Performance Tests of the Second-Order Approximate Internally Contracted Multireference Coupled-Cluster Singles and Doubles Method icMRCC2.

Benchmark results are presented for the second-order approximation of the internally contracted multireference coupled-cluster method with single and double excitations, icMRCC2 [Köhn, Bargholz, J. Chem. Phys. 2019, 151, 041106], which was designed as a multireference analogue of the single-reference second-order approximate coupled-cluster method CC2 [Christiansen, Koch, Jørgensen, Chem. Phys. Lett. 1995, 243, 409-418]. Vertical excitation energies of various small to medium-sized organic molecules are investigated based on established test sets from the literature. Additionally, the spectroscopic constants of ground and excited states of diatomics and the geometric parameters of excited triatomic molecules were determined and compared to the experimental data. The results show that the method clearly extends the applicability of single-reference CC2, including doubly excited states, and also artifacts of CC2 like too low Rydberg excitations and too weak multiple bonds are eliminated. The method is computationally more demanding than standard multireference second-order perturbation theories but improves significantly in accuracy, as shown by the benchmark results. In addition, it is demonstrated that small active spaces are often sufficient to obtain accurate energies with icMRCC2. Example applications like the automerization of cyclobutadiene, the deactivation pathway of ethylene, and the excited states of an iron complex with a noninnocent nitrosyl ligand demonstrate the potential of icMRCC2 in cases with strong multireference character.

The lowest singlet states of hexatriene revisited

The two lowest excited singlet states of trans-hexatriene and cis-hexatriene are studied by multireference perturbation theory approaches (NEVPT2 and CASPT2) in their quasi-degenerate version (QD-NEVPT2 and MS-CASPT2). For these states, we report spectroscopic properties such as the vertical and adiabatic excitation energies, some features of the topology of the potential energy surfaces (PES), and the emission energies. The theoretical vertical excitation energies for the 21Ag-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${2{\\hspace{0.05cm}}{^1\ extrm{A}_g^-}}$$\\end{document} and 11Bu+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${1{\\hspace{0.05cm}}{^1\ extrm{B}_u^+}}$$\\end{document} states of trans-hexatriene are found to be almost degenerate, with a value, ≃\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\simeq$$\\end{document} 5.5–5.6 eV, higher than that normally accepted in the literature, 5.2 eV and 5.1 eV, respectively. The 21A1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${2{\\hspace{0.05cm}}{^1\ extrm{A}_1}}$$\\end{document} and 11B2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${1{\\hspace{0.05cm}}{^1\ extrm{B}_2}}$$\\end{document} states of cis-hexatriene are also almost degenerate and are estimated to be at ≃\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\simeq$$\\end{document} 5.4–5.5 and ≃\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\simeq$$\\end{document} 5.5 eV, respectively, again higher than the accepted values. The adiabatic excitation energies to the 21Ag-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${2{\\hspace{0.05cm}}{^1\ extrm{A}_g^-}}$$\\end{document} and 21A1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${2{\\hspace{0.05cm}}{^1\ extrm{A}_1}}$$\\end{document} states can be observed experimentally (in particular for the cis isomer), and our results are in excellent agreement with the experimental values. On the contrary, the vertical excitation energies for these states are not directly observable in the experimental spectra and the “experimental” values are obtained by educated guesses. We show that the hypotheses underlying these guesses are not entirely grounded.

Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory.

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.

Dynamic Correlation on the Adaptive Sampling Configuration Interaction (ASCI) Reference Function: ASCI-DSRG-MRPT2.

A balanced description of static and dynamic electron correlations is at the heart of quantum chemical methods. To obtain accurate results in strongly correlated systems using wave-function-based methods, a large active space is necessary to ensure correct descriptions of static correlations. Correcting the results for dynamic correlations is also necessary. In this work, we present implementations of second-order perturbation theory for dynamic correlations based on the adaptive sampling configuration interaction self-consistent field (ASCI-SCF) method. In particular, we implemented spin-free driven similarity renormalization group second-order multireference perturbation theory (DSRG-MRPT2). The extrapolation of the ASCI + PT2 energy based on the relaxed Hamiltonian in DSRG-MRPT2 gives a reasonable approximation of DSRG-MRPT2 based on CASSCF. We demonstrate the application of the ASCI-DSRG-MRPT2 method in evaluations of the spin-state energy gaps in iron porphyrins, polyacenes, and periacenes along with the reaction energies of methane oxidation by FeO+ and electrocyclic ring formation in cethrene.

Molecular Geometry and Electronic Structure of Copper Corroles.

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree-Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S0) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T0), which is approximately 5 kcal/mol higher in energy compared to the S0 state for both the bare and substituted copper corroles. Notably, the planarity of the T0 geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S0 and T0 geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

Studying excited-state-specific perturbation theory on the Thiel set.

We explore the performance of a recently introduced N5-scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set's inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles normoffers a relatively low-cost way to test for doubly excited character without the need to define an active space.

Accurate Vertical Electron Detachment Energies and Multiphoton Resonant Photoelectron Spectra of Biochromophore Anions in Aqueous Solution.

We introduce a new methodology for calculating vertical electron detachment energies (VDEs) of biologically relevant chromophores in their deprotonated anionic forms in aqueous solution. It combines a large-scale mixed DFT/EFP/MD approach with the high-level multireference perturbation theory XMCQDPT2 and the Effective Fragment Potential (EFP) method. The methodology includes a multiscale flexible treatment of inner (∼1000 water molecules) and outer (∼18000 water molecules) water shells around a charged solute, capturing both the effects of specific solvation and the properties of bulk water. VDEs are calculated as a function of system size for getting a converged value at the DFT/EFP level of theory. The XMCQDPT2/EFP approach, adapted for calculating VDEs, supports the DFT/EFP results. When corrected for a solvent polarization contribution, the XMCQDPT2/EFP method yields the most accurate estimate to date of the first VDE for aqueous phenolate (7.3 ± 0.1 eV), which agrees well with liquid-jet X-ray photoelectron spectroscopy data (7.1 ± 0.1 eV). We show that the geometry of the water shell and its size are essential for accurate VDE calculations of aqueous phenolate and its biologically relevant derivatives. By simulating photoelectron spectra of aqueous phenolate upon two-photon excitation at wavelengths resonant with the S0 → S1 transition, we also provide interpretation of recent multiphoton UV liquid-microjet photoelectron spectroscopy experiments. We show that its first VDE is consistent with our estimate of 7.3 eV, when experimental two-photon binding energies are corrected for the resonant contribution.

Corresponding Active Orbital Spaces along Chemical Reaction Paths.

The accuracy of reaction energy profiles calculated with multiconfigurational electronic structure methods and corrected by multireference perturbation theory depends crucially on consistent active orbital spaces selected along the reaction path. However, it has been challenging to choose molecular orbitals that can be considered corresponding in different molecular structures. Here, we demonstrate how active orbital spaces can be selected consistently along reaction coordinates in a fully automatized way. The approach requires no structure interpolation between reactants and products. Instead, it emerges from a synergy of the Direct Orbital Selection orbital mapping ansatz combined with our fully automated active space selection algorithm autoCAS. We demonstrate our algorithm for the potential energy profile of the homolytic carbon-carbon bond dissociation and rotation around the double bond of 1-pentene in the electronic ground state. However, our algorithm also applies to electronically excited Born-Oppenheimer surfaces.

CO Inversion on a NaCl(100) Surface: A Multireference Quantum Embedding Study.

We develop a multireference quantum embedding model to investigate a recent experimental observation of the isomerization of vibrationally excited CO molecules on a NaCl(100) surface [Science 2020, 367, 175-178]. To explore this mechanism, we built a reduced potential energy surface of CO interacting with NaCl(100) using a second-order multireference perturbation theory, modeling the adsorbate-surface interaction with our previously developed active space embedding theory (ASET). We considered an isolated CO molecule on NaCl(100) and a high-coverage CO monolayer (1/1), and for both we generated potential energy surfaces parametrized by the CO stretching, adsorption, and inversion coordinates. These surfaces are used to determine stationary points and adsorption energies and to perform a vibrational analysis of the states relevant to the inversion mechanism. We found that for near-equilibrium bond lengths, CO adsorbed in the C-down configuration is lower in energy than in the O-down configuration. Stretching of the C-O bond reverses the energetic order of these configurations, supporting the accepted isomerization mechanism. The vibrational constants obtained from these potential energy surfaces show a small (< 10 cm-1) blue- and red-shift for the C-down and O-down configurations, respectively, in agreement with experimental assignments and previous theoretical studies. Our vibrational analysis of the monolayer case suggests that the O-down configuration is energetically more stable than the C-down one beyond the 16th vibrational excited state of CO, a value slightly smaller than the one from quasi-classical trajectory simulations (22nd) and consistent with the experiment. Our analysis suggests that CO-CO interactions in the monolayer play an important role in stabilizing highly vibrationally excited states in the O-down configuration and reducing the barrier between the C-down and O-down geometries, therefore playing a crucial role in the inversion mechanism.

Statistical errors in reduced density matrices sampled from quantum circuit simulation and the impact on multireference perturbation theory

In this work, we present a detailed analysis of statistical errors in reduced density matrices (RDMs) of active space wavefunctions sampled from quantum circuit simulation and the impact on results obtained by the multireference theories.