Eckart-Sayvetz conditions revisited

A numerical-tensorial "hybrid" nuclear motion Hamiltonian and dipole moment operator for spectra calculation of polyatomic nonrigid molecules

The analysis and modeling of high-resolution spectra of nonrigid molecules require a specific Hamiltonian and group-theoretical formulation that differs significantly from that of more familiar rigid systems. Within the framework of Hougen-Bunker-Johns (HBJ) theory, this paper is devoted to the construction of a nonrigid Hamiltonian based on a suitable combination of numerical calculations for the nonrigid part in conjunction with the irreducible tensor operator method for the rigid part. For the first time, a variational calculation from ab initio potential energy surfaces is performed using the HBJ kinetic energy operator built from vibrational, large-amplitude motion, and rotational tensor operators expressed in terms of curvilinear and normal coordinates. Group theory for nonrigid molecules plays a central role in the characterization of the overall tunneling splittings and is discussed in the present approach. The construction of the dipole moment operator is also examined. Validation tests consisting of a careful convergence study of the energy levels as well as a comparison of results obtained from independent computer codes are given for the nonrigid molecules CH2, CH3, NH3, and H2O2. This work paves the way for the modeling of high-resolution spectra of larger nonrigid systems.

First-Principles Modeling of Preferential Solvation in Mixed-Modifier Differential Mobility Spectrometry

Differential mobility spectrometry (DMS) separates ions based on mobility differences between high and low electric field conditions. To enhance resolution, solvents such as water and acetonitrile are often used to modify the collision environment and take advantage of differing dynamic clustering behavior between analytes that coelute in hard-sphere environments (e.g., N₂). When binary solvent mixtures are used to modify the DMS environment, one solvent can have a dominant influence over the other with respect to ion trajectories. For example, for quinoline derivatives, a 9:1 water:acetonitrile solvent mixture exhibited identical behavior to an environment containing only acetonitrile as a modifier. It was hypothesized that this effect arises due to the significantly different binding strengths of the two solvents. Here, we utilize a first-principles model of DMS to study analytes in single and binary solvent mixtures and explore the effects governing the dominance of one solvent over the other. Computed DMS dispersion curves of quinoline derivatives are in excellent agreement with those measured experimentally. For mixed-modifier environments, the predicted cluster populations show a clear preferential solvation of ions with the stronger binding solvent. The influence of ion-solvent binding energies, solvent concentration, and solvent molecule size is discussed in the context of the observed DMS behavior. This work can guide the usage of binary solvent mixtures for improving ion separations in cases where compounds coelute in pure N₂ and in single-solvent modifier environments. Moreover, our results indicate that binary solvent mixtures can be used to create a relative scale for solvent binding energies.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.

Umbrella inversion of ammonia redux

Umbrella inversion of ammonia is a prototypical example of large-amplitude vibrational motion, described with a symmetric double-well potential. The transition state of the latter corresponds to a planar (D_3h) molecular geometry, whereas the two equilibrium configurations are equivalent (C_3v) pyramidal structures, with the nitrogen atom being either 'above' or 'below' the plane of the hydrogen atoms. As commonly understood, inversion motion of ammonia corresponds to the coherent, anharmonic, vibrational motion of the molecule, which shuttles back and forth between the two potential wells; that is, oscillation of the nitrogen atom from one side of the H₃ plane to the other, via coherent tunneling. However, this intuitively appealing view of umbrella inversion results from a reduced description of the dynamics, which includes only the inversion vibrational coordinate and fully neglects all the other molecular degrees of freedom. As such, this textbook picture of inversion motion ignores the fact that the two equilibrium structures of ammonia are superimposable, and can only be distinguished by labelling the identical hydrogen nuclei. A correct description of umbrella inversion, which incorporates nuclear permutations, requires the inclusion of other molecular modes. Indeed, it is well known that the quantum symmetrization postulate engenders entanglement between ammonia's nuclear-spin, inversion, and rotation. Using the explicit expressions of the corresponding zeroth-order eigenstates, we clearly show that the inversion density of any multilevel wavepacket of ammonia, including the case of perfectly aligned molecules, is symmetrically distributed between the two potential wells, at all times. This follows from a rigorous demonstration based on the evaluation of the expectation values of the inversion coordinate or equivalent projection operators. However, provided that these wavepackets involve inversion-rotation levels with opposite parity, the inversion density may exhibit dynamical spatial localization. In the latter case, the space-fixed inversion density or, equivalently, the expectation values of the projections of the inversion coordinate on the space-fixed axes, may oscillate between opposite directions in the space-fixed frame. Nevertheless, in all cases, localization of ammonia in a single potential well is impossible, even partially or transiently. This is equivalent to saying that the nitrogen atom has the same probability (one-half) to be on either side of the H₃ plane, for any wavepacket of the molecule and at all times-a conclusion which is in perfect accord with the principle of the indistinguishability of identical particles (nuclei).

Derivation of ρ-dependent coordinate transformations for nonrigid molecules in the Hougen-Bunker-Johns formalism

In this paper, we report a series of transformations for the construction of a Hamiltonian model for nonrigid polyatomic molecules in the framework of the Hougen-Bunker-Johns formalism (HBJ). This model is expressed in normal mode coordinates for small vibrations and in a specific coordinate ρ to describe the large amplitude motion. For the first time, a general procedure linking the "true" curvilinear coordinates to ρ is proposed, allowing the expression of the potential energy part in the same coordinate representation as the kinetic energy operator, whatever the number of atoms. A Lie group-based method is also proposed for the derivation of the reference configuration in the internal axis system. This work opens new perspectives for future high-resolution spectroscopy studies of nonrigid, medium-sized molecules using HBJ-type Hamiltonians. Illustrative examples and computation of vibrational energy levels on semirigid and nonrigid molecules are given to validate this method.

Capturing Deuteration Effects in a Molecular Mechanics Force Field: Deuterated THF and the THF-Water Miscibility Gap

Deuteration is a common chemical modification used in conjunction with experiments such as neutron scattering, NMR, and Fourier-transform infrared for the study of molecular systems. Under the Born-Oppenheimer (BO) approximation, while the underlying potential energy surface remains unchanged by isotopic substitutions, isotopic substitution still alters intramolecular vibrations, which in turn may alter intermolecular interactions. Molecular mechanics (MM) force fields used in classical molecular dynamics simulations are assumed to represent local approximations of the BO potential energy surfaces, and hence, MD simulations using simple isotopic mass substitutions should capture BO-compatible isotope effects. However, standard MM force-field parameterizations do not directly fit to the local harmonic quantum mechanical (QM) Hessian that describes the BO surface, but rather to QM normal-modes and/or mass-dependent internal-coordinate derived distortion energies. Here, using tetrahydrofuran (THF)-water mixtures as our model system, we show that not only does a simple mass-substitution approach fail to capture an experimentally characterized deuteration effect (the loss of the closed-loop miscibility gap associated with the complete deuteration of THF) but also it is necessary to generate new MM force-field parameters that correctly describe isotopic dependent vibrations to capture the experimental deuteration effect. We show that the origin of this failure is a result of using mass-dependent features to fit the THF MM force field, which unintentionally biases the bonded terms of the force field to represent only the isotopologue used during the original force-field parameterization. In addition, we make use of our isotopologue-corrected force field for D₈THF to examine the molecular origins of the isotope-dependent loss of the THF-water miscibility gap.

Internal- and rho-axis systems of molecules with one large amplitude internal motion: The geometry of rho

The internal-axis system (IAS) of molecules with a large amplitude internal motion (LAM) is determined by integrating the kinematic equation of the IAS by Lie-group and Lie-algebraic methods. Numerical examples on hydrogen peroxide, nitrous acid, and acetaldehyde demonstrate the methods. By exploiting the special product structure of the solution matrix, simple methods are devised for calculating the transformation to the rho-axis system (RAS) along with the value of the parameter ρ characterizing a RAS rotational-LAM kinetic energy operator. The parameter ρ so calculated agrees exactly with that one obtained by the Floquet method as shown in the example of acetaldehyde. Geometrical interpretation of ρ is given. The advantageous property of the RAS over the IAS in retaining simple periodic boundary conditions is numerically demonstrated.

A general formulation of the quasiclassical trajectory method for reduced-dimensionality reaction dynamics calculations

Dimension reduction by freezing the unimportant coordinates is widely used in intramolecular and reaction dynamics calculations when the solution of the accurate full-dimensional nuclear Schrödinger equation is not feasible. In this paper we report on a novel form of the exact classical internal-coordinate Hamiltonian for full and reduced-dimensional vibrational motion of polyatomic molecules with the purpose of using it in quasiclassical trajectory (QCT) calculations. The derivation is based on the internal to body-fixed frame transformation, as in the t-vector formalism, however it does not require the introduction of rotational variables to allow cancellation of non-physical rotations within the body-fixed frame. The formulas needed for QCT calculations: normal mode analysis and state sampling as well as for following the dynamics and normal-mode quantum number assignment at instantaneous states are presented. The procedure is demonstrated on the CH4, CD4, CH3D and CHD3 isotopologs of methane using three reduced-dimensional models, which were previously used in quantum reactive scattering studies of the CH4 + X → CH3 + HX type reactions. The reduced-dimensional QCT methodology formulated this way combined with full-dimensional QCT calculations makes possible the classical validation of reduced-dimensional models that are used in the quantum mechanical description of the nuclear dynamics in reactive systems [A. Vikár et al., J. Phys. Chem. A, 2016, 120, 5083-5093].