An official website of the United States government
Some reactions produce extremely hot nascent-products which nevertheless can form sufficiently long-lived van der Waals (vdW) complexes—with atoms or molecules from a bath gas—as to be observed via microwave spectroscopy. Theoretical calculations of such unbound resonance-states can be much more challenging than ordinary bound-state calculations depending on the approach employed. One encounters not only the floppy, and perhaps multi-welled potential energy surface (PES) characteristic of vdWs complexes, but in addition must contend with excitation of the intramolecular modes and its corresponding influence on the PES. Straightforward computation of the (resonance) rovibrational levels of interest, involves the added complication of the unbound nature of the wavefunction, often treated with techniques such as introducing a complex absorbing potential. Here, we have demonstrated that a simplified approach of making a series of vibrationally effective PESs for the intermolecular coordinates—one for each reaction product vibrational quantum number of interest—can produce vdW levels for the complex with spectroscopic accuracy. This requires constructing a series of appropriately weighted lower-dimensional PESs for which we use our freely available PES-fitting code AUTOSURF. The applications of this study are the Ar–CS and Ar–SiS complexes, which are isovalent to Ar–CO and Ar–SiO, the latter of which we considered in a previously reported study. Using a series of vibrationally effective PESs, rovibrational levels and predicted microwave transition frequencies for both complexes were computed variationally. A series of shifting rotational transition frequencies were also computed as a function of the diatom vibrational quantum number. The predicted transitions were used to guide and inform an experimental effort to make complementary observations. Comparisons are given for the transitions that are within the range of the spectrometer and were successfully recorded. Calculations of the rovibrational level pattern agree to within 0.2 % with experimental measurements.
more »
« less
Quantum eigenvector continuation for chemistry applications
Abstract A typical task for classical and quantum computing in chemistry is finding a potential energy surface (PES) along a reaction coordinate, which involves solving the quantum chemistry problem for many points along the reaction path. Developing algorithms to accomplish this task on quantum computers has been an active area of development, yet finding all the relevant eigenstates along the reaction coordinate remains a difficult problem, and determining PESs is thus a costly proposal. In this paper, we demonstrate the use of a eigenvector continuation—a subspace expansion that uses a few eigenstates as a basis—as a tool for rapidly exploring PESs. We apply this to determining the binding PES or torsion PES for several molecules of varying complexity. In all cases, we show that the PES can be captured using relatively few basis states; suggesting that a significant amount of (quantum) computational effort can be saved by making use of already calculated ground states in this manner.
more »
« less
- Award ID(s):
- 1752713
- PAR ID:
- 10473545
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Electronic Structure
- Volume:
- 5
- Issue:
- 4
- ISSN:
- 2516-1075
- Format(s):
- Medium: X Size: Article No. 045007
- Size(s):
- Article No. 045007
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The variational quantum eigensolver is one of the most promising approaches for performing chemistry simulations using noisy intermediate-scale quantum (NISQ) processors. The efficiency of this algorithm depends crucially on the ability to prepare multi-qubit trial states on the quantum processor that either include, or at least closely approximate, the actual energy eigenstates of the problem being simulated while avoiding states that have little overlap with them. Symmetries play a central role in determining the best trial states. Here, we present efficient state preparation circuits that respect particle number, total spin, spin projection, and time-reversal symmetries. These circuits contain the minimal number of variational parameters needed to fully span the appropriate symmetry subspace dictated by the chemistry problem while avoiding all irrelevant sectors of Hilbert space. We show how to construct these circuits for arbitrary numbers of orbitals, electrons, and spin quantum numbers, and we provide explicit decompositions and gate counts in terms of standard gate sets in each case. We test our circuits in quantum simulations of the$${H}_{2}$$ and$$LiH$$ molecules and find that they outperform standard state preparation methods in terms of both accuracy and circuit depth.more » « less
-
The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.more » « less
-
Molecular dynamics (MD) simulations with full-dimensional potential energy surfaces (PESs) obtained from high-level ab initio calculations are frequently used to model reaction dynamics of small molecules (i.e., molecules with up to 10 atoms). Construction of full-dimensional PESs for larger molecules is, however, not feasible since the number of ab initio calculations required grows rapidly with the increase of dimension. Only a small number of coordinates are often essential for describing the reactivity of even very large systems, and reduced-dimensional PESs with these coordinates can be built for reaction dynamics studies. While analytical methods based on transition-state theory framework are well established for analyzing the reduced-dimensionalPESs, MD simulation algorithms capable of generating trajectories on such surfaces are more rare. In this work, we present a new MD implementation that utilizes the relaxed reduced-dimensional PES for standard micro canonical (NVE) and canonical (NVT) MD simulations.The method is applied to the pyramidal inversion of a NH3molecule. The results from the MD simulations on a reduced, three-dimensional PES are validated against the ab initio MD simulations, as well as MD simulations on full-dimensional PES and experimental data.more » « less
-
We present the computational methodology that allows rigorous and efficient nine-dimensional (9D) quantum calculations of the intermolecular vibrational states of noncovalently bound trimers of diatomic molecules, with the monomers treated as rigid. The full 9D vibrational Hamiltonian of the trimer is partitioned into a 3D “frame” (or stretching) Hamiltonian and a 6D “bend” Hamiltonian. These two Hamiltonians are diagonalized separately, and a certain number of their lowest-energy eigenstates is included in the final 9D product contracted basis in which the full 9D intermolecular vibrational Hamiltonian is diagonalized. This methodology is applied to the 9D calculations of the intermolecular vibrational levels of (HF)3, a prototypical hydrogen-bonded trimer, on the rigid-monomer version of an ab initio calculated potential energy surface (PES). They are the first to include fully the stretch-bend coupling present in the trimer. The frequencies of all bending fundamentals considered from the present 9D calculations are about 10% lower than those from the earlier quantum 6D calculations that considered only the bending modes of the HF trimer. This means that the stretch-bend coupling is strong, and it is imperative to include it in any accurate treatment of the (HF)3 vibrations aiming to assess the accuracy of the PES employed. Moreover, the 9D results are in better agreement with the limited available spectroscopic data that those from the 6D calculations. In addition, the 9D results show sensitivity to the value of the HF bond length, equilibrium or vibrationally averaged, used in the calculations. The implication is that full-dimensional 12D quantum calculations will be required to obtain definitive vibrational excitation energies for a given PES. Our study also demonstrates that the nonadditive three-body interactions are very significant in (HF)3 and have to be included in order to obtain accurate intermolecular vibrational energy levels of the trimer.more » « less
