The global minima of urea and thiourea were characterized along with other low‐lying stationary points. Each structure was optimized with the CCSD(T) method and triple‐
- NSF-PAR ID:
- 10328082
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 22
- Issue:
- 42
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 24270 to 24281
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract ζ correlation consistent basis sets followed by harmonic vibrational frequency computations. Relative energies evaluated near the complete basis set limit with both canonical and explicitly correlated CCSD(T) techniques reveal several subtle but important details about both systems. These computations resolve a discrepancy by demonstrating that the electronic energy of the C2v second‐order saddle point of urea lies at least 1.5 kcal mol−1above the C2global minimum regardless of whether the structures were optimized with MP2, CCSD, or CCSD(T). Additionally, urea effectively has one minimum instead of two because the electronic barrier for inversion at one amino group in the Cslocal minimum vanishes at the CCSD(T) CBS limit. Characterization of both systems with the same ab initio methods and large basis sets conclusively establishes that the electronic barriers to inversion at one or both NH2groups in thiourea are appreciably smaller than in urea. CCSDT(Q)/cc‐pVTZ computations show higher‐order electron correlation effects have little impact on the relative energies and are consistently offset by core correlation effects of opposite sign and comparable magnitude. -
We present molecular dynamics (MD), polarizability driven MD (α-DMD), and pump–probe simulations of Raman spectra of the protonated nitrogen dimer N4H+, and some of its isotopologues, using the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)]-F12b/aug-cc-pVTZ based potential energy surface in permutationally invariant polynomials (PIPs) of Yu et al. [J. Phys. Chem. A 119, 11623 (2015)] and a corresponding PIP-derived CCSD(T)/aug-cc-pVTZ-tr (N:spd, H:sp) polarizability tensor surface (PTS), the latter reported here for the first time. To represent the PTS in terms of a PIP basis, we utilize a recently described formulation for computing the polarizability using a many-body expansion in the orders of dipole–dipole interactions while generating a training set using a novel approach based on linear regression for potential energy distributions. The MD/α-DMD simulations reveal (i) a strong Raman activity at 260 and 2400 cm−1, corresponding to the symmetric N–N⋯H bend and symmetric N–N stretch modes, respectively; (ii) a very broad spectral region in the 500–2000 cm−1 range, assignable to the parallel N⋯H+⋯N proton transfer overtone; and (iii) the presence of a Fermi-like resonance in the Raman spectrum near 2400 cm−1 between the Σg+ N–N stretch fundamental and the Πu overtone corresponding to perpendicular N⋯H+⋯N proton transfer.
-
ABSTRACT This paper presents rate coefficients for transitions between rotational levels of the A-type and E-type nuclear spin modifications of methanol induced by collisions with molecular hydrogen. These rate coefficients are required for an accurate determination of methanol abundance in the interstellar medium, where local thermodynamic equilibrium conditions generally do not apply. Time-independent close-coupling quantum scattering calculations have been employed to calculate cross-sections and rate coefficients for the (de-)excitation of methanol in collisions with para- and ortho-H2. These calculations utilized a potential energy surface (PES) for the interaction of methanol with H2 recently computed by the explicitly correlated CCSD(T)-F12a coupled-cluster method that employed a correlation-consistent aug-cc-pVTZ basis. Rate coefficients for temperatures ranging from 3 to 250 K were calculated for all transitions among the first 76 rotational levels of both A-type and E-type methanol, whose energies are less than or equal to 170 K. These rate coefficients are compared with those by Rabli and Flower who carried out coupled-state calculations using a PES computed by second-order many-body perturbation theory. Simple radiative transfer calculations using the present set of rate coefficients are also reported and compared with such calculations using the rate coefficients previously computed by Rabli and Flower.
-
Abstract Quantum simulation of chemical systems is one of the most promising near-term applications of quantum computers. The variational quantum eigensolver, a leading algorithm for molecular simulations on quantum hardware, has a serious limitation in that it typically relies on a pre-selected wavefunction ansatz that results in approximate wavefunctions and energies. Here we present an arbitrarily accurate variational algorithm that, instead of fixing an ansatz upfront, grows it systematically one operator at a time in a way dictated by the molecule being simulated. This generates an ansatz with a small number of parameters, leading to shallow-depth circuits. We present numerical simulations, including for a prototypical strongly correlated molecule, which show that our algorithm performs much better than a unitary coupled cluster approach, in terms of both circuit depth and chemical accuracy. Our results highlight the potential of our adaptive algorithm for exact simulations with present-day and near-term quantum hardware.
-
We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(α→) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy, thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ Hamiltonian using rotation and reflection symmetry and H2 Hamiltonian within Sz=0 subspace using S2 symmetry. We have found the variationally trained unitary gives good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers.more » « less