Electron and nuclear dynamics following ionisation of modified bismethylene-adamantane

We have simulated the coupled electron and nuclear dynamics using the Ehrenfest method upon valence ionisation of modified bismethylene-adamantane (BMA) molecules where there is an electron transfer between the two π bonds. We have shown that the nuclear motion significantly affects the electron dynamics after a few fs when the electronic states involved are close in energy. We have also demonstrated how the non-stationary electronic wave packet determines the nuclear motion, more precisely the asymmetric stretching of the two π bonds, illustrating “charge-directed reactivity”. Taking into account the nuclear wave packet width results in the dephasing of electron dynamics with a half-life of 8 fs; this eventually leads to the equal delocalisation of the hole density over the two methylene groups and thus symmetric bond lengths

Molecular second-quantized Hamiltonian:Electron correlation and non-adiabatic coupling treated on an equal footing

We introduce a new theoretical and computational framework for treating molecular quantum mechanics without the Born–Oppenheimer approximation. The molecular wavefunction is represented in a tensor-product space of electronic and vibrational basis functions, with electronic basis chosen to reproduce the mean-field electronic structure at all geometries. We show how to transform the Hamiltonian to a fully second-quantized form with creation/annihilation operators for electronic and vibrational quantum particles, paving the way for polynomial-scaling approximations to the tensor-product space formalism. In addition, we make a proof-of-principle application of the new Ansatz to the vibronic spectrum of C2

Direct quantum dynamics using variational Gaussian wavepackets and Gaussian process regression

The method of direct variational quantum nuclear dynamics in a basis of Gaussian wavepackets, combined with the potential energy surfaces fitted on-the-fly using Gaussian process regression, is described together with its implementation. Enabling exact and efficient analytic evaluation of Hamiltonian matrix elements, this approach allows for black-box quantum dynamics of multidimensional anharmonic molecular systems. Example calculations of intra-molecular proton transfer on the electronic ground state of salicylaldimine are provided, and future algorithmic improvements as well as the potential for multiple-state non-adiabatic dynamics are discussed

Toward accurate QM/MM reaction barriers with large QM regions using domain based pair natural orbital coupled cluster theory

The hydroxylation reaction catalyzed by p-hydroxybenzoate hydroxylase and the Baeyer–Villiger reaction catalyzed by cyclohexanone monooxygenase are investigated by means of quantum mechanical/molecular mechanical (QM/MM) calculations at different levels of QM theory. The geometries of the stationary points along the reaction profile are obtained from QM/MM geometry optimizations, in which the QM region is treated by density functional theory (DFT). Relative energies are determined from single-point QM/MM calculations using the domain-based local pair natural orbital coupled cluster DLPNO-CCSD(T) method as QM component. The results are compared with single-point DFT/MM energies obtained using popular density functionals and with available experimental and computational data. It is found that the choice of the QM method strongly affects the computed energy profiles for these reactions. Different density functionals provide qualitatively different energy barriers (variations of the order of 10 kcal/mol in both reactions), thus limiting the confidence in DFT/MM computational predictions of energy profiles. On the other hand, the use of the DLPNO-CCSD(T) method in conjunction with large QM regions and basis sets makes it possible to achieve high accuracy. A critical discussion of all the technical aspects of the calculations is given with the aim of aiding computational chemists in the application of the DLPNO-CCSD(T) methodology in QM/MM calculations

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework

A complete description of tunnelling using direct quantum dynamics simulation: Salicylaldimine proton transfer

We demonstrate here conclusively that the variational multiconfiguration Gaussian (vMCG) method converges to the grid based full quantum dynamics multiconfiguration time-dependent Hartree result for a tunnelling problem in many dimensions, using the intramolecular proton transfer in salicylaldimine as a model system. The 13-dimensional model potential energy surface was obtained from Hartree Fock energies with the 6-31G* basis set and the expectation value of the flux operator along the transition mode was used as a benchmark characteristic. As well as showing excellent convergence of the vMCG method on the model surface using a local harmonic approximation and a moderate number of basis functions, we show that the direct dynamics version of the vMCG also performs very well, usually needs the same number of Gaussians to converge, and converges to exact results if those are obtained on an accurately fitted surface. Finally, we make an important observation that the width of the Gaussian basis functions must be chosen very carefully to obtain accurate results with the use of the frozen-width approximation

Quantum Mechanical/Molecular Mechanical Study on the Mechanism of the Enzymatic Baeyer–Villiger Reaction

We report a combined quantum mechanical/molecular mechanical (QM/MM) study on the mechanism of the enzymatic Baeyer–Villiger reaction catalyzed by cyclohexanone monooxygenase (CHMO). In QM/MM geometry optimizations and reaction path calculations, density functional theory (B3LYP/TZVP) is used to describe the QM region consisting of the substrate (cyclohexanone), the isoalloxazine ring of C4a-peroxyflavin, the side chain of Arg-329, and the nicotinamide ring and the adjacent ribose of NADP+, while the remainder of the enzyme is represented by the CHARMM force field. QM/MM molecular dynamics simulations and free energy calculations at the semiempirical OM3/CHARMM level employ the same QM/MM partitioning. According to the QM/MM calculations, the enzyme–reactant complex contains an anionic deprotonated C4a-peroxyflavin that is stabilized by strong hydrogen bonds with the Arg-329 residue and the NADP+ cofactor. The CHMO-catalyzed reaction proceeds via a Criegee intermediate having pronounced anionic character. The initial addition reaction has to overcome an energy barrier of about 9 kcal/mol. The formed Criegee intermediate occupies a shallow minimum on the QM/MM potential energy surface and can undergo fragmentation to the lactone product by surmounting a second energy barrier of about 7 kcal/mol. The transition state for the latter migration step is the highest point on the QM/MM energy profile. Gas-phase reoptimizations of the QM region lead to higher barriers and confirm the crucial role of the Arg-329 residue and the NADP+ cofactor for the catalytic efficiency of CHMO. QM/MM calculations for the CHMO-catalyzed oxidation of 4-methylcyclohexanone reproduce and rationalize the experimentally observed (S)-enantioselectivity for this substrate, which is governed by the conformational preferences of the corresponding Criegee intermediate and the subsequent transition state for the migration step

Quantum Mechanical/Molecular Mechanical Study on the Enantioselectivity of the Enzymatic Baeyer–Villiger Reaction of 4‑Hydroxycyclohexanone

We report a combined quantum mechanical/molecular mechanical (QM/MM) study of the effect of mutations of the Phe434 residue in the active site of cyclohexanone monooxygenase (CHMO) on its enantioselectivity toward 4-hydroxycyclohexanone. In terms of our previously established model of the enzymatic Baeyer–Villiger reaction, enantioselectivity is governed by the preference toward the equatorial ((<i>S</i>)-selectivity) or axial ((<i>R</i>)-selectivity) conformation of the substituent at the C4 carbon atom of the cyclohexanone ring in the Criegee intermediate and the subsequent rate-limiting transition state for migration (TS2). We assess the enantiopreference by locating all relevant TS2 structures at the QM/MM level. In the wild-type enzyme we find that the axial conformation is energetically slightly more stable, thus leading to a small excess of (<i>R</i>)-product. In the Phe434Ser mutant, there is a hydrogen bond between the serine side chain and the equatorial substrate hydroxyl group that is retained during the whole reaction, and hence there is pronounced reverse (<i>S</i>)-enantioselectivity. Another mutant, Phe434Ile, is shown to preserve and enhance the (<i>R</i>)-selectivity. All these findings are in accordance with experiment. The QM/MM calculations allow us to explain the effect of point mutations on CHMO enantioselectivity for the first time at the molecular level by an analysis of the specific interactions between substrate and active-site environment in the TS2 structures that satisfy the basic stereoelectronic requirement of anti-periplanarity for the migrating σ-bond