The nuclear-spin-forbidden rovibrational transitions of water from first principles.

The water molecule occurs in two nuclear-spin isomers that differ by the value of the total nuclear spin of the hydrogen atoms, i.e., I = 0 for para-H2O and I = 1 for ortho-H2O. Spectroscopic transitions between rovibrational states of ortho and para water are extremely weak due to the tiny hyperfine nuclear-spin-rotation interaction of only ∼30 kHz and, so far, have not been observed. We report the first comprehensive theoretical investigation of the hyperfine effects and ortho-para transitions in H2 16O due to nuclear-spin-rotation and spin-spin interactions. We also present the details of our newly developed general variational approach to the simulation of hyperfine effects in polyatomic molecules. Our results for water suggest that the strongest ortho-para transitions with room-temperature intensities on the order of 10-31 cm/molecule are about an order of magnitude larger than previously predicted values and should be detectable in the mid-infrared ν2 and near-infrared 2ν1 + ν2 and ν1 + ν2 + ν3 bands by current spectroscopy experiments.

Using collocation and a hierarchical basis to solve the vibrational Schrödinger equation.

We show that it is possible to compute vibrational energy levels of polyatomic molecules with a collocation method and a basis of products of one-dimensional harmonic oscillator functions pruned so that it does not include functions for which the indices of many of the one-dimensional functions are nonzero. Functions with many nonzero indices are coupled only by terms that depend simultaneously on many coordinates, and they are typically small. The collocation equation is derived without invoking differences of interpolation operators, which simplifies implementation of the method. This, however, requires inverting a matrix whose elements are values of the pruned basis functions at the collocation points. The collocation points are the points on a Smolyak grid whose size is equal to the size of the pruned basis set. The Smolyak grid is built from symmetrized Leja points. Because both the basis and the grid are not tensor products, the inverse is not straightforward. It can be done by using so-called hierarchical 1-D basis functions. They are defined so that the matrix whose elements are the 1-D hierarchical basis functions evaluated at points is lower triangular. We test the method by applying it to compute 100 energy levels of CH2NH with an iterative eigensolver.

Field-Induced Diastereomers for Chiral Separation.

A novel approach for the state-specific enantiomeric enrichment and the spatial separation of enantiomers is presented. Our scheme utilizes techniques from strong-field laser physics-specifically an optical centrifuge in conjunction with a static electric field-to create a chiral field with defined handedness. Molecular enantiomers experience unique rotational excitation dynamics, and this can be exploited to spatially separate the enantiomers using electrostatic deflection. Notably, the rotational-state-specific enantiomeric enhancement and its handedness are fully controllable. To explain these effects, the conceptual framework of field-induced diastereomers of a chiral molecule is introduced and computationally demonstrated through robust quantum-mechanical simulations on the prototypical chiral molecule propylene oxide (C_{3}H_{6}O), for which ensembles with an enantiomeric excess of up to 30% were obtained.

Ro-vibronic transition intensities for triatomic molecules from the exact kinetic energy operator; electronic spectrum for the C̃ 1B2 ← X̃ 1A1 transition in SO2.

A procedure for calculating ro-vibronic transition intensities for triatomic molecules within the Born-Oppenheimer approximation is reported. Ro-vibrational energy levels and wavefunctions are obtained with the DVR3D suite, which solves the nuclear motion problem with an exact kinetic energy operator. Absolute transition intensities are calculated both with the Franck-Condon approximation and with a full transition dipole moment surface. The theoretical scheme is tested on C̃ 1B2 ← X̃ 1A1 ro-vibronic transitions of SO2. Ab initio potential energy and dipole moment surfaces are generated for this purpose. The calculated ro-vibronic transition intensities and cross sections are compared with the available experimental and theoretical data.

Simulating electric field interactions with polar molecules using spectroscopic databases.

Ro-vibrational Stark-associated phenomena of small polyatomic molecules are modelled using extensive spectroscopic data generated as part of the ExoMol project. The external field Hamiltonian is built from the computed ro-vibrational line list of the molecule in question. The Hamiltonian we propose is general and suitable for any polar molecule in the presence of an electric field. By exploiting precomputed data, the often prohibitively expensive computations associated with high accuracy simulations of molecule-field interactions are avoided. Applications to strong terahertz field-induced ro-vibrational dynamics of PH3 and NH3, and spontaneous emission data for optoelectrical Sisyphus cooling of H2CO and CH3Cl are discussed.

Quasiperiodic energy dependence of exciton relaxation kinetics in the sexithiophene crystal.

Femtosecond kinetics of fluorescence rise in the sexithiophene crystal is studied on a microscopic model of intraband relaxation, where exciton energy is assumed to be dissipated by phonon-accompanied scattering, with the rates calculated earlier. The temporal evolution of the exciton population is described by a set of kinetic equations, solved numerically to yield the population buildup at the band bottom. Not only the time scale but also the shape of the rise curves is found to be unusually sensitive to excitation energy, exhibiting unique quasiperiodic dependence thereon, which is rationalized in terms of the underlying model. Further simulations demonstrate that the main conclusions are robust with respect to experimental factors such as finite temperature and inherent spectral broadening of the exciting pulse, while the calculated fluorescence rise times are found to be in excellent agreement with experimental data available to date. As the rise profiles are composed of a number of exponential contributions, which varies with excitation energy, the common practice of characterizing the population buildup in the emitting state by a single value of relaxation time turns out to be an oversimplification. New experiments giving further insight into the kinetics and mechanism of intraband exciton relaxation are suggested.

A uniform approach to the description of multicenter bonding.

A novel method for investigating the multicenter bonding patterns in molecular systems by means of the so-called Electron Density of Delocalized Bonds (EDDB) is introduced and discussed. The EDDB method combines the concept of Jug's bond-order orbitals and the indirect ("through-bridge") interaction formalism and opens up new opportunities for studying the interplay between different atomic interactions as well as their impact on both local and global resonance stabilization in systems of conjugated bonds. Using several illustrative examples we demonstrate that the EDDB approach allows for a reliable quantitative description of diverse multicenter delocalization phenomena (with special regard to evaluation of the aromatic stabilization in molecular systems) within the framework of a consistent theoretical paradigm.