On the Probability Density of the Nuclei in a Vibrationally Excited Molecule

Laser-Induced Electron Symmetry Restoration in Oriented Molecules Made Simple

Electron symmetry determines many important properties of molecules, from selection rules for photoelectron spectroscopy to symmetry selection rules for chemical reactions. The original electron symmetry is broken if a laser pulse changes the initial state, typically the ground state g, to a superposition of g and an excited state e with different irreducible representations (IRREPs). Quantum dynamics simulations for two examples, the oriented benzene and LiCN molecules, show that the original electron symmetry can be restored by means of a reoptimized π-laser pulse which transfers the component in the excited state e to another state e', or to several others with the same IRREP as the ground state. This method lends itself to much easier experimental applications than all previous ones because it allows the healing of electron symmetry immediately, without any attosecond constraint on the timing of the second pulse.

Orientational decoherence within molecules and emergence of the molecular shape

The question of classicality is addressed in relation with the shape of the nuclear skeleton of molecular systems. As the most natural environment, the electrons of the molecule are considered as continuously monitoring agents for the nuclei. For this picture, an elementary formalism of decoherence theory is developed and numerical results are presented for few-particle systems. The numerical examples suggest that the electron-nucleus Coulomb interaction is sufficient for inducing a blurred shape with strong quantum coherences in compounds of the lightest elements, H₂, D₂, T₂, and HeH⁺.

Representing molecular ground and excited vibrational eigenstates with nuclear densities obtained from semiclassical initial value representation molecular dynamics

We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in the protonated glycine molecule [Aieta et al., Nat Commun 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in the correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, also to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond lengths, angles, and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching, the typical harmonic nodal pattern is absent in the anharmonic distribution.

Anharmonic quantum nuclear densities from full dimensional vibrational eigenfunctions with application to protonated glycine

The interpretation of molecular vibrational spectroscopic signals in terms of atomic motion is essential to understand molecular mechanisms and for chemical characterization. The signals are usually assigned after harmonic normal mode analysis, even if molecular vibrations are known to be anharmonic. Here we obtain the quantum anharmonic vibrational eigenfunctions of the 11-atom protonated glycine molecule and we calculate the density distribution of its nuclei and its geometry parameters, for both the ground and the O-H stretch excited states, using our semiclassical method based on ab initio molecular dynamics trajectories. Our quantum mechanical results describe a molecule elongated and more flexible with respect to what previously thought. More importantly, our method is able to assign each spectral peak in vibrational spectroscopy by showing quantitatively how normal modes involving different functional groups cooperate to originate that spectroscopic signal. The method will possibly allow for a better rationalization of experimental spectroscopy.

H3 + as a five-body problem described with explicitly correlated Gaussian basis sets

Various explicitly correlated Gaussian (ECG) basis sets are considered for the solution of the molecular Schrödinger equation with particular attention to the simplest polyatomic system, H₃ ⁺. Shortcomings and advantages are discussed for plain ECGs, ECGs with the global vector representation, floating ECGs and their numerical projection, and ECGs with complex parameters. The discussion is accompanied with particle density plots to visualize the observations. In order to be able to use large complex ECG basis sets in molecular calculations, a numerically stable algorithm is developed, the efficiency of which is demonstrated for the lowest rotationally and vibrationally excited states of H₂ and H₃ ⁺.