Anharmonic Assignment of the Water Octamer Spectrum in the OH Stretch Region.

We interface the quasi-classical trajectory approach with an ab initio potential energy surface for water to assign the vibrational spectroscopical features of the OH stretch region of the water octamer cluster, which is considered to be a precursor of ice. An attempt by Li et al. to assign their recent reference experiment involved lower-level calculations based on an ad hoc scaled harmonic approach. Differently from the conclusions of this previous assignment, which invoked the contribution of 5 conformers and a solvated form of the water heptamer in the spectrum, we find out that the spectroscopic features can be related to the 4 conformers of the octamer lying lower in energy.

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.

Quantum Critical Behavior of One-Dimensional Soft Bosons in the Continuum.

We consider a zero-temperature one-dimensional system of bosons interacting via the soft-shoulder potential in the continuum, typical of dressed Rydberg gases. We employ quantum Monte Carlo simulations, which allow for the exact calculation of imaginary-time correlations, and a stochastic analytic continuation method, to extract the dynamical structure factor. At finite densities, in the weakly interacting homogeneous regime, a rotonic spectrum marks the tendency to clustering. With strong interactions, we indeed observe cluster liquid phases emerging, characterized by the spectrum of a composite harmonic chain. Luttinger theory has to be adapted by changing the reference lattice density field. In both the liquid and cluster liquid phases, we find convincing evidence of a secondary mode, which becomes gapless only at the transition. In that region, we also measure the central charge and observe its increase towards c=3/2, as recently evaluated in a related extended Bose-Hubbard model, and we note a fast reduction of the Luttinger parameter. For two-particle clusters, we then interpret such observations in terms of the compresence of a Luttinger liquid and a critical transverse Ising model, related to the instability of the reference lattice density field towards coalescence of sites, typical of potentials which are flat at short distances. Even in the absence of a true lattice, we are able to evaluate the spatial correlation function of a suitable pseudospin operator, which manifests ferromagnetic order in the cluster liquid phase, exponential decay in the liquid phase, and algebraic order at criticality.