Equilibrium properties of quantum water clusters by the variational Gaussian wavepacket method

Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches

The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.

Thermally Induced Hydrogen-Bond Rearrangements in Small Water Clusters and the Persistent Water Tetramer

Small water clusters absorb heat and catalyze pivotal atmospheric reactions. Yet, experiments produced conflicting results on water cluster distribution under atmospheric conditions. Additionally, it is unclear which "phase transitions" such clusters exhibit, at what temperatures, and what are their underlying molecular mechanisms. We find that logarithmically small tails in the radial probability densities of (H₂O) _n clusters (n = 2 - 6) provide direct testimony for such transitions. Using the best available water potential (MB-pol), an advanced thermostating algorithm (g-BAOAB), and sufficiently long trajectories, we map the "bifurcation", "melting", and (hitherto unexplored) "vaporization" transitions, finding that both melting and vaporization proceed via a "monomer on a ring" conformer, exhibiting huge distance fluctuations at the vaporization temperatures (T _v). T _v may play a role in determining the atmospheric cluster size distribution such that the dimer and tetramer, with their exceptionally low/high T _v values, are under/over-represented in these distributions, as indeed observed in nondestructive mass spectrometric measurements.

Phase transitions in free water nanoparticles. Theoretical modeling of [H2O]48 and [H2O]118

Classical parallel-tempering Monte Carlo simulations of [H2O]48 and [H2O]118 have been performed in the isothermal-isobaric ensemble and a two-dimensional multiple-histogram method has been used to calculate the heat capacity of the two clusters. A semiempirical procedure is proposed for the inclusion of quantum effects and transformed heat capacity profiles are compared with state-of-the-art experimental data [C. Hock et al., Phys. Rev. Lett., 2009, 103, 073401]. A very good agreement is achieved. A detailed analysis of the simulation data is provided to gain an insight into the nature of the phase change which takes place in the two clusters at T ≈ 100 K.

The Nosé-Hoover looped chain thermostat for low temperature thawed Gaussian wave-packet dynamics

We have used a generalised coherent state resolution of the identity to map the quantum canonical statistical average for a general system onto a phase-space average over the centre and width parameters of a thawed Gaussian wave packet. We also propose an artificial phase-space density that has the same behaviour as the canonical phase-space density in the low-temperature limit, and have constructed a novel Nosé-Hoover looped chain thermostat that generates this density in conjunction with variational thawed Gaussian wave-packet dynamics. This forms a new platform for evaluating statistical properties of quantum condensed-phase systems that has an explicit connection to the time-dependent Schrödinger equation, whilst retaining many of the appealing features of path-integral molecular dynamics.

Comparison between different Gaussian series representations of the imaginary time propagator

A useful approximation for the thermal operator exp(-betaH[over ]) is based on its representation in terms of either frozen or thawed Gaussian states. Such approximate representations are leading-order terms in respective series representations of the thermal operator. A numerical study of the convergence properties of the frozen Gaussian series representation has been recently published. In this paper, we extend the previous study to include also the convergence properties of the more expensive thawed Gaussian series representation of the thermal operator. We consider three different formulations for the series representation and apply them to a quartic double-well potential to find that the thawed Gaussian series representation converges faster than the frozen Gaussian one. Further analysis is presented as to the convergence properties and the numerical efficiency of three different thawed Gaussian series representation. The unsymmetrized form converges most rapidly, however, the lower order approximations of the symmetrized forms are more accurate. Comparison with a standard discretized path-integral evaluation demonstrates that the Gaussian based perturbation series representation converges much faster.

Calorimetric observation of the melting of free water nanoparticles at cryogenic temperatures

We present an experimental study of the thermodynamics of free, size-selected water cluster anions consisting of 48 and 118 molecules. The measured caloric curves of the clusters are bulklike at low temperatures but show a well-defined, particle-size specific transition at 93+/-3 K for (H2O)48- and 118+/-3 K for (H2O)118-. At the transition temperature the heat capacity strongly increases, which marks the onset of melting.

Frozen Gaussian series representation of the imaginary time propagator theory and numerical tests

Thawed Gaussian wavepackets have been used in recent years to compute approximations to the thermal density matrix. From a numerical point of view, it is cheaper to employ frozen Gaussian wavepackets. In this paper, we provide the formalism for the computation of thermal densities using frozen Gaussian wavepackets. We show that the exact density may be given in terms of a series, in which the zeroth order term is the frozen Gaussian. A numerical test of the methodology is presented for deep tunneling in the quartic double well potential. In all cases, the series is observed to converge. The convergence of the diagonal density matrix element is much faster than that of the antidiagonal one, suggesting that the methodology should be especially useful for the computation of partition functions. As a by product of this study, we find that the density matrix in configuration space can have more than two saddle points at low temperatures. This has implications for the use of the quantum instanton theory.

Quantum transitions in Lennard-Jones clusters

The ground states of Lennard-Jones (LJ) clusters are estimated by finding the Gaussian wave packets that minimize the energy functional. A "phase diagram" for LJ_{n} as a function of size (n=31,...,45) and de Boer quantum delocalization length (Lambdain[0;0.3]) is constructed, showing the stability ranges for the two competing structural motifs, the Mackay and anti-Mackay icosahedra. An increase of Lambda has an effect similar to heating and as such may induce structural transformations.