Simulating the vibrational quantum dynamics of molecules using photonics.

Advances in control techniques for vibrational quantum states in molecules present new challenges for modelling such systems, which could be amenable to quantum simulation methods. Here, by exploiting a natural mapping between vibrations in molecules and photons in waveguides, we demonstrate a reprogrammable photonic chip as a versatile simulation platform for a range of quantum dynamic behaviour in different molecules. We begin by simulating the time evolution of vibrational excitations in the harmonic approximation for several four-atom molecules, including H2CS, SO3, HNCO, HFHF, N4 and P4. We then simulate coherent and dephased energy transport in the simplest model of the peptide bond in proteins-N-methylacetamide-and simulate thermal relaxation and the effect of anharmonicities in H2O. Finally, we use multi-photon statistics with a feedback control algorithm to iteratively identify quantum states that increase a particular dissociation pathway of NH3. These methods point to powerful new simulation tools for molecular quantum dynamics and the field of femtochemistry.

Continuous-wave cavity ringdown for high-sensitivity polarimetry and magnetometry measurements.

We report the development of a novel variant of cavity ringdown polarimetry using a continuous-wave laser operating at 532 nm for highly precise chiroptical activity and magnetometry measurements. The key methodology of the apparatus relies upon the external modulation of the laser frequency at the frequency splitting between non-degenerate left- and right-circularly polarized cavity modes. The method is demonstrated by the evaluation of the Verdet constants of crystalline CeF3 and fused silica, in addition to the observation of gas- and solution-phase optical rotations of selected chiral molecules. Specifically, optical rotations of (i) vapors of α-pinene and R-(+)-limonene, (ii) mutarotating D-glucose in water, and (iii) acidified L-histidine solutions are determined. The detection sensitivities for the gas- and solution-phase chiral activity measurements are ∼30 and ∼120µdeg over a 30 s detection period per cavity round trip pass, respectively. Furthermore, the measured optical rotations for R-(+)-limonene are compared with computations performed using the TURBOMOLE quantum chemistry package. The experimentally observed optically rotatory dispersion of this cyclic monoterpene was thus rationalized via a consideration of its room temperature conformer distribution as determined by the aforementioned single-point energy calculations.

Molecular excited state calculations with adaptive wavefunctions on a quantum eigensolver emulation: reducing circuit depth and separating spin states.

Ab initio electronic excited state calculations are necessary for the quantitative study of photochemical reactions, but their accurate computation on classical computers is plagued by prohibitive resource scaling. The Variational Quantum Deflation (VQD) is an extension of the quantum-classical Variational Quantum Eigensolver (VQE) algorithm for calculating electronic excited state energies, and has the potential to address some of these scaling challenges using quantum computers. However, quantum computers available in the near term can only support a limited number of quantum circuit operations, so reducing the quantum computational cost in VQD methods is critical to their realisation. In this work, we investigate the use of adaptive quantum circuit growth (ADAPT-VQE) in excited state VQD calculations, a strategy that has been successful previously in reducing the resources required for ground state energy VQE calculations. We also invoke spin restrictions to separate the recovery of eigenstates with different spin symmetry to reduce the number of calculations and accumulation of errors in computing excited states. We created a quantum eigensolver emulation package - Quantum Eigensolver Building on Achievements of Both quantum computing and quantum chemistry (QEBAB) - for testing the proposed adaptive procedure against two existing VQD methods that use fixed-length quantum circuits: UCCGSD-VQD and k-UpCCGSD-VQD. For a lithium hydride test case we found that the spin-restricted adaptive growth variant of VQD uses the most compact circuits out of the tested methods by far, consistently recovers adequate electron correlation energy for different nuclear geometries and eigenstates while isolating the singlet and triplet manifold. This work is a further step towards developing techniques which improve the efficiency of hybrid quantum algorithms for excited state quantum chemistry, opening up the possibility of exploiting real quantum computers for electronic excited state calculations sooner than previously anticipated.

Hydrogen-Atom Tunneling in Metaphosphorous Acid.

Invited for the cover of this issue is X. Zeng and co-workers at Soochow University, University of Stuttgart, and Max-Planck Institute for Solid State Research. The image depicts the fast tunneling transformation of the highly elusive metaphosphorous acid (HOPO). Read the full text of the article at 10.1002/chem.202000844.

Hydrogen-Atom Tunneling in Metaphosphorous Acid.

Metaphosphorous acid (HOPO), a key intermediate in phosphorus chemistry, has been generated in syn- and anti-conformations in the gas phase by high-vacuum flash pyrolysis (HVFP) of a molecular precursor ethoxyphosphinidene oxide (EtOPO→C2 H4 +HOPO) at ca. 1000 K and subsequently trapped in an N2 -matrix at 2.8 K. Unlike the two conformers of the nitrogen analogue HONO, the anti-conformer of HOPO undergoes spontaneous rotamerization at 2.8 K via hydrogen-atom tunneling (HAT) with noticeable kinetic isotope effects for H/D (>104 for DOPO) and 16 O/18 O (1.19 for H18 OPO and 1.06 for HOP18 O) in N2 -matrices.

Basis set extrapolation in pair natural orbital theories.

We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.

The rotational spectrum of H2S⋯HI and an investigation by ab initio calculations of the origins of the observed doubling of rotational transitions in both H2S⋯HI and H2S⋯F2.

The rotational spectrum of the complex H2S⋯HI observed with a pulsed-jet, Fourier-transform microwave spectrometer shows that each rotational transition is split into a closely spaced doublet, a pattern similar to that observed earlier for the halogen-bonded complex H2S⋯F2. The origin of the doubling has been investigated by means of ab initio calculations conducted at the CCSD(T)(F12*)/cc-pVDZ-F12 level. Two paths were examined by calculating the corresponding energy as a function of two angles. One path involved inversion of the configuration at S through a planar transition state of C2v symmetry via changes in the angle ϕ between the C2 axis of H2S and the line joining the H and I nuclei [the potential energy function V(ϕ)]. The other was a torsional oscillation θ about the local C2 axis of H2S that also exchanges the equivalent H nuclei [the potential energy function V(θ)]. The inversion path is slightly lower in energy and much shorter in arc length and is therefore the favored tunneling pathway. In addition, calculation of V(ϕ) for the series of hydrogen- and halogen-bonded complexes H2S⋯HX (X = F, Cl, or Br) and H2S⋯XY (XY = Cl2, Br2, ClF, BrCl, or ICl) at the same level of theory revealed that doubling is unlikely to be resolved in these, in agreement with experimental observations. The barrier heights of the V(ϕ) of all ten complexes examined were found to be almost directly proportional to the dissociation energies De.

Ab initio calculation of rovibrational states for non-degenerate double-well potentials: cis-trans isomerization of HOPO.

The rovibrational spectra of metaphosphorous acid, HOPO, and its deuterated isotopologue have been studied by vibrational configuration interaction calculations, relying on the internal coordinate path Hamiltonian and the Watson Hamiltonian. Tunneling effects for the overtones of the torsional mode, which gives rise to the cis-trans isomerization, and its rovibrational transitions have been investigated in detail. Due to strong matrix effects, comparison with experimental data is hindered, and thus, the calculations provide accurate estimates for the fundamental modes of these species.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations.

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Anharmonic excited state frequencies of para-difluorobenzene, toluene and catechol using analytic RI-CC2 second derivatives.

Analytic second nuclear derivatives for excited electronic state energies have been implemented for the resolution-of-the-identity accelerated CC2, CIS(D∞) and ADC(2) models. Our efficient implementation with O(N2) memory demands enables the treatment of medium sized molecules with large basis sets and high numerical precision and thereby paves the way for semi-numerical evaluation of the higher-order derivatives required for anharmonic corrections to excited state vibrational frequencies. We compare CC2 harmonic and anharmonic excited state frequencies with experimental values for para-difluorobenzene, toluene and catechol. Basis set problems occur for out-of-plane bending vibrations due to intramolecular basis set superposition error. For non-planar molecules and in plane modes of planar molecules, the agreement between theory and experiment is better than 30 cm-1 on average and we reassign a number of experimental bands on the basis of the ab initio predictions.

Experimental and computational studies of Criegee intermediate reactions with NH3 and CH3NH2.

Ammonia and amines are emitted into the troposphere by various natural and anthropogenic sources, where they have a significant role in aerosol formation. Here, we explore the significance of their removal by reaction with Criegee intermediates, which are produced in the troposphere by ozonolysis of alkenes. Rate coefficients for the reactions of two representative Criegee intermediates, formaldehyde oxide (CH2OO) and acetone oxide ((CH3)2COO) with NH3 and CH3NH2 were measured using cavity ring-down spectroscopy. Temperature-dependent rate coefficients, k(CH2OO + NH3) = (3.1 ± 0.5) × 10-20T2 exp(1011 ± 48/T) cm3 s-1 and k(CH2OO + CH3NH2) = (5 ± 2) × 10-19T2 exp(1384 ± 96/T) cm3 s-1 were obtained in the 240 to 320 K range. Both the reactions of CH2OO were found to be independent of pressure in the 10 to 100 Torr (N2) range, and average rate coefficients k(CH2OO + NH3) = (8.4 ± 1.2) × 10-14 cm3 s-1 and k(CH2OO + CH3NH2) = (5.6 ± 0.4) × 10-12 cm3 s-1 were deduced at 293 K. An upper limit of ≤2.7 × 10-15 cm3 s-1 was estimated for the rate coefficient of the (CH3)2COO + NH3 reaction. Complementary measurements were performed with mass spectrometry using synchrotron radiation photoionization giving k(CH2OO + CH3NH2) = (4.3 ± 0.5) × 10-12 cm3 s-1 at 298 K and 4 Torr (He). Photoionization mass spectra indicated production of NH2CH2OOH and CH3N(H)CH2OOH functionalized organic hydroperoxide adducts from CH2OO + NH3 and CH2OO + CH3NH2 reactions, respectively. Ab initio calculations performed at the CCSD(T)(F12*)/cc-pVQZ-F12//CCSD(T)(F12*)/cc-pVDZ-F12 level of theory predicted pre-reactive complex formation, consistent with previous studies. Master equation simulations of the experimental data using the ab initio computed structures identified submerged barrier heights of -2.1 ± 0.1 kJ mol-1 and -22.4 ± 0.2 kJ mol-1 for the CH2OO + NH3 and CH2OO + CH3NH2 reactions, respectively. The reactions of NH3 and CH3NH2 with CH2OO are not expected to compete with its removal by reaction with (H2O)2 in the troposphere. Similarly, losses of NH3 and CH3NH2 by reaction with Criegee intermediates will be insignificant compared with reactions with OH radicals.

Efficient and accurate description of adsorption in zeolites.

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.

Ab initio instanton rate theory made efficient using Gaussian process regression.

Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into the calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but it is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but it requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

Communication: Quasi-robust local density fitting.

The requirement that the linear density fitting error in the integral exactly vanishes introduces unphysical long range contributions to the approximate density when the auxiliary basis is incomplete. A quasi-robust density fitting formulation is presented where spatial locality is recovered at the expense of permitting a linear error that is made small by the fitting procedure, which involves optimising the Coulomb potential of the approximate charge density. The method is shown to be stable and almost as accurate as standard robust density fitting without local approximations in practical calculations using standard density fitting basis sets.

Do CCSD and approximate CCSD-F12 variants converge to the same basis set limits? The case of atomization energies.

While the title question is a clear "yes" from purely theoretical arguments, the case is less clear for practical calculations with finite (one-particle) basis sets. To shed further light on this issue, the convergence to the basis set limit of CCSD (coupled cluster theory with all single and double excitations) and of different approximate implementations of CCSD-F12 (explicitly correlated CCSD) has been investigated in detail for the W4-17 thermochemical benchmark. Near the CBS ([1-particle] complete basis set) limit, CCSD and CCSD(F12*) agree to within their respective uncertainties (about ±0.04 kcal/mol) due to residual basis set incompleteness error, but a nontrivial difference remains between CCSD-F12b and CCSD(F12*), which is roughly proportional to the degree of static correlation. The observed basis set convergence behavior results from the superposition of a rapidly converging, attractive, CCSD[F12]-CCSD-F12b difference (consisting mostly of third-order terms) and a more slowly converging, repulsive, fourth-order difference between CCSD(F12*) and CCSD[F12]. For accurate thermochemistry, we recommend CCSD(F12*) over CCSD-F12b if at all possible. There are some indications that the nZaPa family of basis sets exhibits somewhat smoother convergence than the correlation consistent family.

Cooperative hydrogen bonds form a pseudocycle stabilizing an isolated complex of isocyanic acid with urea.

The shapes of macromolecules and their complexes with small molecules are often determined by extended networks of hydrogen bonds. Here, for the first time, we provide a detailed description of a cooperative pair of hydrogen bonds to an individual molecule of urea. The structure and properties of a gas phase complex formed between urea and isocyanic acid are characterised through microwave spectroscopy and ab initio calculations at the CCSD(T)(F12*)/aug-cc-pVTZ level.

Molecular geometries and other properties of H2O⋯AgI and H3N⋯AgI as characterised by rotational spectroscopy and ab initio calculations.

The rotational spectra of H3N⋯AgI and H2O⋯AgI have been recorded between 6.5 and 18.5 GHz by chirped-pulse Fourier-transform microwave spectroscopy. The complexes were generated through laser vaporisation of a solid target of silver or silver iodide in the presence of an argon gas pulse containing a low concentration of the Lewis base. The gaseous sample subsequently undergoes supersonic expansion which results in cooling of rotational and vibrational motions such that weakly bound complexes can form within the expanding gas jet. Spectroscopic parameters have been determined for eight isotopologues of H3N⋯AgI and six isotopologues of H2O⋯AgI. Rotational constants, B0; centrifugal distortion constants, DJ, DJK or ΔJ, ΔJK; and the nuclear quadrupole coupling constants, χaa(I) and χbb(I) - χcc(I) are reported. H3N⋯AgI is shown to adopt a geometry that has C3v symmetry. The geometry of H2O⋯AgI is Cs at equilibrium but with a low barrier to inversion such that the vibrational wavefunction for the v = 0 state has C2v symmetry. Trends in the nuclear quadrupole coupling constant of the iodine nucleus, χaa(I), of L⋯AgI complexes are examined, where L is varied across the series (L = Ar, H3N, H2O, H2S, H3P, or CO). The results of experiments are reported alongside those of ab initio calculations at the CCSD(T)(F12*)/AVXZ level (X = T, Q).

The Dynamics of the Reaction of FeO+ and H2 : A Model for Inorganic Oxidation.

Extensive density functional theory (DFT) calculations using the B3LYP functional were used to explore the sextet and quartet energy potential energy surfaces (PESs) of the title reaction, and as a basis to fit global analytical reactive PESs. Surface-hopping dynamics on these PESs reproduce the experimentally observed reactivity and confirm that hydrogen activation rather than spin-state change is rate-limiting at low reaction energy, where the main products are Fe+ and H2 O. A change in spin state is inefficient in the product region so that excited-state 4 Fe+ is the dominant product. At higher energies, spin-allowed hydrogen atom abstraction to form FeOH+ predominates. At intermediate energy, a previously unexpected rebound mechanism contributes significantly to the reactivity.

A Structurally Characterized Fluoroalkyne.

The facile synthesis of a stable and isolable compound with a fluoroalkynyl group, M-C≡CF, is reported. Reaction of [Ru(C≡CH)(η5 -C5 Me5 )(dppe)] with an electrophilic fluorinating agent (NFSI) results in the formation of the fluorovinylidene complex [Ru(=C=CHF)(η5 -C5 Me5 )(dppe)][N(SO2 Ph)2 ]. Subsequent deprotonation with LiN(SiMe3 )2 affords the fluoroalkynyl complex [Ru(C≡CF)(η5 -C5 Me5 )(dppe)]. In marked contrast to the rare and highly reactive examples of fluoroalkynes that have been reported previously, this compound can be readily isolated and structurally characterized. This has allowed the structure and bonding in the CCF motif to be explored. Further electrophilic fluorination of this species yields the difluorovinylidene complex [Ru(C=CF2 )(η5 -C5 Me5 )(dppe)][N(SO2 Ph)2 ].

Geometry of an Isolated Dimer of Imidazole Characterised by Rotational Spectroscopy and Ab Initio Calculations.

An isolated, gas-phase dimer of imidazole is generated through laser vaporisation of a solid rod containing a 1:1 mixture of imidazole and copper in the presence of an argon buffer gas undergoing supersonic expansion. The complex is characterised through broadband rotational spectroscopy and is shown to have a twisted, hydrogen-bonded geometry. Calculations at the CCSD(T)(F12*)/cc-pVDZ-F12 level of theory confirm this to be the lowest-energy conformer of the imidazole dimer. The distance between the respective centres of mass of the imidazole monomer subunits is determined to be 5.2751(1) Å, and the twist angle γ describing rotation of one monomer with respect to the other about a line connecting the centres of mass of the monomers is determined to be 87.9(4)°. Four out of six intermolecular parameters in the model geometry are precisely determined from the experimental rotational constants and are consistent with results calculated ab initio.