Search

We extend the intrinsic atomic orbital (IAO) method for the localization of molecular orbitals to calculate well-localized generalized Wannier functions in crystals in the spirit of the Pipek-Mezey method. We furthermore present a one-shot diabatic Wannierization procedure that aligns the phases of the Bloch functions, providing immediate Wannier localization, which serves as an excellent initial guess for optimization. We test our Wannier localization implementation on a number of solid-state systems, highlighting the effectiveness of the diabatic preparation, especially for localizing core bands. Partial charges of Wannier functions generated using Bloch IAOs align well with chemical intuition, which we demonstrate through the example of the adsorption of CO on a MgO surface.

First-quantized, grid-based methods for chemistry modeling are a natural and elegant fit for quantum computers. However, it is infeasible to use today's quantum prototypes to explore the power of this approach because it requires a substantial number of near-perfect qubits. Here, we use exactly emulated quantum computers with up to 36 qubits to execute deep yet resource-frugal algorithms that model 2D and 3D atoms with single and paired particles. A range of tasks is explored, from ground state preparation and energy estimation to the dynamics of scattering and ionization; we evaluate various methods within the split-operator QFT (SO-QFT) Hamiltonian simulation paradigm, including protocols previously described in theoretical papers and our own techniques. While we identify certain restrictions and caveats, generally, the grid-based method is found to perform very well; our results are consistent with the view that first-quantized paradigms will be dominant from the early fault-tolerant quantum computing era onward.

Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.

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 H 2 CS, SO 3 , HNCO, HFHF, N 4 and P 4 . 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 H 2 O. Finally, we use multi-photon statistics with a feedback control algorithm to iteratively identify quantum states that increase a particular dissociation pathway of NH 3 . These methods point to powerful new simulation tools for molecular quantum dynamics and the field of femtochemistry.

A global analytic potential energy surface for the ro-vibrational dynamics of cis- and trans-formic acid is presented, constructed using LASSO-based regression to reproduce CCSD(T)(F12*)/cc-pVTZ-F12 energies. The fit is accurate to 0.25% has an RMS deviation from the ab initio data of 9 cm -1 for the energy range 0-15000 cm -1 . Converged J = 0 vibrational eigenstates are reported, computed using vibrational configuration interaction with an internal coordinate path Hamiltonian for the torsional motion connecting the cis and trans rotamers. Methodological choices concerning the appropriate definitions of the curvilinear and diabatic bath coordinates are discussed. The zero point of the cis rotamer is 1412 cm -1 above that of the trans, which lies at 7354 cm -1 . The computed fundamentals match the bands recorded from gas-phase IR spectroscopy with an RMSD of only 3 cm -1 . A fresh assignment of the overtone spectra of both the cis and trans rotamers is presented for the energy range 0-4720 cm -1 , where 14 out of the 51 bands are reassigned on the basis of the VCI calculations.

A molecular complex of C2H2 and AuI has been generated and isolated in the gas phase through laser ablation of a gold surface in the presence of an expanding sample containing small percentages of C2H2 and CF3I in a buffer gas of argon. Rotational, B0, centrifugal distortion, ΔJ and ΔJK, and nuclear quadrupole coupling constants, χaa(Au), χbb(Au) - χcc(Au), χaa(I), and χbb(I) - χcc(I), are measured for three isotopologues of C2H2···AuI through broadband rotational spectroscopy. The complex is C2v and T-shaped with C2H2 coordinating to the gold atom via donation of electrons from the π-orbitals of ethyne. On formation of the complex, the C≡C bond of ethyne extends by 0.032(4) Å relative to r(C≡C) in isolated ethyne when the respective r0 geometries are compared. The geometry of ethyne distorts such that ∠(*-C-H) (where * indicates the midpoint of the C≡C bond) is 194.7(12)° in the r0 geometry of C2H2···AuI. Ab initio calculations at the CCSD(T)(F12*)/AVTZ level are consistent with the experimentally determined geometry and further allow calculation of the dissociation energy (De) as 136 kJ mol(-1). The χaa(Au) and χaa(I) nuclear quadrupole coupling constants of AuI and also the Au-I bond length change significantly on formation of the complex consistent with the strong interaction calculated to occur between C2H2 and AuI.

The new linear molecule Ag-C≡C-Cl has been detected and fully characterized by means of rotational spectroscopy. It was synthesized by laser ablation of a silver rod in the presence of a gaseous sample containing a low concentration of CCl4 in argon, cooled to a rotational temperature approaching ∼1-3 K through supersonic expansion, and analyzed by chirped-pulse, Fourier transform microwave spectroscopy. Six isotopologues were investigated, and for each the spectroscopic constants B0, D(J) and χ(aa)(Cl) were determined. The B0 values were interpreted to give the following bond lengths: r(Ag-C) = 2.015(14) Å and r(C-Cl) = 1.635(6) Å, with r(C≡C) = 1.2219 Å assumed from an ab initio calculation at the CCSD(T)/aug-cc-pV5Z level of theory. The Cu analogue Cu-C≡C-Cl was similarly identified and characterized.

Absorption cross-sections and quantum yields for NO2 production (ΦNO2) are reported for gaseous methyl, ethyl, n-propyl, and isopropyl nitrate at 294 K. Absorption cross-sections in the wavelength range of 240-320 nm agree well with prior determinations. NO2 quantum yields at photoexcitation wavelengths of 290, 295, and 315 nm are unity within experimental uncertainties for all of the alkyl nitrates studied and are independent of bath gas (N2) pressure for total sample pressures in the range of 250-700 Torr. When averaged over all wavelengths and sample pressures, values of ΦNO2 are 1.03 ± 0.05 (methyl nitrate), 0.98 ± 0.09 (ethyl nitrate), 1.01 ± 0.04 (n-propyl nitrate), and 1.00 ± 0.05 (isopropyl nitrate), with uncertainties corresponding to 1 standard deviation. Absorption cross-sections for ethyl nitrate, isopropyl nitrate, and two unsaturated dinitrate compounds, but-3-ene-1,2-diyl dinitrate and (Z)-but-2-ene-1,4-diyl dinitrate in acetonitrile solution, are compared to gas-phase values, and over the wavelength range of 260-315 nm, the gas-phase values are well-reproduced by dividing the liquid-phase cross-sections by 2.0, 1.6, 1.7, and 2.2, respectively. Reasonable estimates of the gas-phase absorption cross-sections for low-volatility organic nitrates can therefore be obtained by halving the values for acetonitrile solutions. The quantum yield for NO2 formation from photoexcitation of but-3-ene-1,2-diyl dinitrate at 290 nm is significantly lower than those for the alkyl (mono) nitrates: a best estimate of ΦNO2 ≤ 0.25 is obtained from the experimental measurements.

H (Rydberg) atom photofragment translation spectroscopy and high-level ab initio electronic structure calculations are used to explore the photodissociation dynamics of three para-substituted thiophenols (p-YPhSH; Y = CH(3), F, and MeO). UV excitation in the wavelength range 305 > λ(phot) > 240 nm results in S-H bond fission and formation of p-YPhS radicals in their ground (X̃(2)B(1)) and first excited (Ã(2)B(2)) electronic states; the X̃/Ã state product branching ratio, Γ, varies with para-Y substituent and excitation wavelength. Excitation at λ(phot) < 265 nm results in direct population of the dissociative 1(1)πσ* potential energy surface (PES). Γ falls across the series p-CH(3)PhSH > p-FPhSH > p-MeOPhSH. Branching is ultimately determined at the conical intersection (CI) formed by the 1(1)πσ* and ground (S(0)) PESs at extended R(S-H) bond length but is sensitively dependent on the orientation of the S-H bond (relative to the ring plane) in the S(0) molecules prior to photoexcitation. Excitation at λ(phot) > 265 nm populates quasi-bound levels of the respective 1(1)ππ* states, which predissociate rapidly by tunneling under the lower diabats of the 1(1)ππ*/1(1)πσ* CI at short R(S-H). Less extreme X̃/Ã product branching ratios are measured, implicating intramolecular vibrational redistribution within the photoexcited 1(1)ππ* molecules prior to their sampling the region of the 1(1)πσ*/S(0) CI.

We have computed barrier heights of 71.8 +/- 2.0 and 216.4 +/- 2.0 kJ mol(-1) for the reactions CH4 + CH3* --> CH3* + CH4 and CH4 + CH3* --> H* + C2H6, respectively, using explicitly correlated coupled cluster theory with singles and doubles combined with standard coupled cluster theory with up to connected quadruple excitations. Transition-state theory has been used to compute the respective reaction rate constants in the temperature interval of 250-1500 K. The computed rates for the reaction to ethane are orders of magnitude slower than those used in the mechanism of Norinaga and Deutschmann (Ind. Eng. Chem. Res. 2007, 46, 3547.) for the modeling of the chemical vapor deposition of pyrolytic carbon.

Benchmark, frozen-core CCSD(T) equilibrium harmonic vibrational frequencies of 12 closed-shell and five open-shell molecules are computed to within 1 cm-1 of the basis set limit using the explicitly correlated CCSD(T)-R12 method. The convergence of the standard CCSD(T) method with the one-particle basis sets of Dunning and co-workers is examined and found to be slow, with mean and maximum absolute errors of 1.3 and 3.5 cm-1 remaining at the cc-pV6Z level. Finite basis set effects do not appear to introduce systematic errors in equilibrium harmonic frequencies, and mean absolute errors reduce by a factor of 2 for each basis set cardinal number increment. The convergence of individual equilibrium harmonic frequencies is not guaranteed to be monotonic due to the associated shift in the equilibrium structure. The inclusion of computed scalar relativistic effects and previously available corrections for core-valence correlation and higher-order excitations in the cluster operator results in an agreement with experimentally derived harmonic frequencies of 0.1, 0.3, and -0.4 cm-1 for HF, N2, and CO, respectively. F2 continues to present a challenge to computational chemistry with an error of 3.2 cm-1, primarily resulting from the high basis set dependence of the quadruples contribution.