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

The variational Gaussian wavepacket (VGW) method in combination with the replica-exchange Monte Carlo is applied to calculations of the heat capacities of quantum water clusters, (H(2)O)(8) and (H(2)O)(10). The VGW method is most conveniently formulated in Cartesian coordinates. These in turn require the use of a flexible (i.e., unconstrained) water potential. When the latter is fitted as a linear combination of Gaussians, all the terms involved in the numerical solution of the VGW equations of motion are analytic. When a flexible water model is used, a large difference in the timescales of the inter- and intramolecular degrees of freedom generally makes the system very difficult to simulate numerically. Yet, given this difficulty, we demonstrate that our methodology is still practical. We compare the computed heat capacities to those for the corresponding classical systems. As expected, the quantum effects shift the melting temperatures toward the lower values.     

Quantum statistical mechanics with Gaussians: equilibrium properties of van der Waals clusters

The variational Gaussian wave-packet method for computation of equilibrium density matrices of quantum many-body systems is further developed. The density matrix is expressed in terms of Gaussian resolution, in which each Gaussian is propagated independently in imaginary time beta=(k(B)T)(-1) starting at the classical limit beta=0. For an N-particle system a Gaussian exp[(r-q)(T)G(r-q)+gamma] is represented by its center qinR(3N), the width matrix GinR(3Nx3N), and the scale gammainR, all treated as dynamical variables. Evaluation of observables is done by Monte Carlo sampling of the initial Gaussian positions. As demonstrated previously at not-very-low temperatures the method is surprisingly accurate for a range of model systems including the case of double-well potential. Ideally, a single Gaussian propagation requires numerical effort comparable to the propagation of a single classical trajectory for a system with 9(N(2)+N)/2 degrees of freedom. Furthermore, an approximation based on a direct product of single-particle Gaussians, rather than a fully coupled Gaussian, reduces the number of dynamical variables to 9N. The success of the methodology depends on whether various Gaussian integrals needed for calculation of, e.g., the potential matrix elements or pair correlation functions could be evaluated efficiently. We present techniques to accomplish these goals and apply the method to compute the heat capacity and radial pair correlation function of Ne(13) Lennard-Jones cluster. Our results agree very well with the available path-integral Monte Carlo calculations.     

Spatial delocalization in para-H2 clusters

We applied the quantum path integral Monte Carlo method for the study of (para-H)N (N = 5-33) clusters at T = 2 K, exploring static and dynamic order, which originates from the effects of zero-point energy, kinetic energy, and thermal fluctuations in quantum clusters. Information on dynamic structure was inferred from the asymptotic tails of the cage correlation function calculated from the centroid Monte Carlo trajectory. The centroid cage correlation function decays to zero for large clusters (N = 15-33), manifesting the interchange of molecules between different solvation shells, with statistically diminishing back interchange. Further evidence for the floppiness of para-hydrogen clusters emerges from the Monte Carlo evolution of the centroid of a tagged molecule, which exhibits significant changes in the list of its first and second solvation shells due to the interchange of molecules between these shells.     

Thermodynamics and equilibrium structure of Ne38 cluster: quantum mechanics versus classical

The equilibrium properties of classical Lennard-Jones (LJ38) versus quantum Ne38 Lennard-Jones clusters are investigated. The quantum simulations use both the path-integral Monte Carlo (PIMC) and the recently developed variational-Gaussian wave packet Monte Carlo (VGW-MC) methods. The PIMC and the classical MC simulations are implemented in the parallel tempering framework. The classical heat capacity Cv(T) curve agrees well with that of Neirotti et al. [J. Chem. Phys. 112, 10340 (2000)], although a much larger confining sphere is used in the present work. The classical Cv(T) shows a peak at about 6 K, interpreted as a solid-liquid transition, and a shoulder at approximately 4 K, attributed to a solid-solid transition involving structures from the global octahedral (Oh) minimum and the main icosahedral (C5v) minimum. The VGW method is used to locate and characterize the low energy states of Ne38, which are then further refined by PIMC calculations. Unlike the classical case, the ground state of Ne38 is a liquidlike structure. Among the several liquidlike states with energies below the two symmetric states (Oh and C5v), the lowest two exhibit strong delocalization over basins associated with at least two classical local minima. Because the symmetric structures do not play an essential role in the thermodynamics of Ne38, the quantum heat capacity is a featureless curve indicative of the absence of any structural transformations. Good agreement between the two methods, VGW and PIMC, is obtained. The present results are also consistent with the predictions by Calvo et al. [J. Chem. Phys. 114, 7312 (2001)] based on the quantum superposition method within the harmonic approximation. However, because of its approximate nature, the latter method leads to an incorrect assignment of the Ne38 ground state as well as to a significant underestimation of the heat capacity.     

Diffusion quantum Monte Carlo for equilibrium structures and harmonic frequencies of ethane and ozone molecules

Application of the Ornstein-Uhlenbeck diffusion quantum Monte Carlo method in combination with a trial wave function constructed from the floating spherical Gaussian orbitals and spherical Gaussian geminals to studies on the equilibrium structures and harmonic frequencies of ethane and ozone is presented. These Monte Carlo computed results are compared with those of experiments as well as the coupled cluster methods with the correlation consistent basis sets for the two molecules. For ozone, we also compare the Monte Carlo results with the results from multireference calculations.     

N2O分解反应的蒙特卡罗模拟

N2O作为大气污染物之一早已受到人们的注意，许多人研究了N2O在金属氧化物上的分解反应得到了一些基本的实验事实[1]：（1）反应速率与原料气中N2O的分压近似成比例．(2)在反应初期，N2的生成速率由大变小，O2的生成速率由小变大，达到平衡时，N2生成速率是O2生成速率的二倍．(3)过渡应答法研究表明，把原料气由N2O和He的混合气体切换成纯He后，N2立即停止发生，而O2的生成仍能维持一段时间．基于上述实验事实，人们对Nzo的催化分解反应提出了一些可能的机理，主要有下面两种看法：另一方面，人们还发现几O的催化分解反应是一个振荡…     

Implementation of dynamical nucleation theory effective fragment potentials method for modeling aerosol chemistry

In this work, the dynamical nucleation theory (DNT) model using the ab initio based effective fragment potential (EFP) is implemented for evaluating the evaporation rate constant and molecular properties of molecular clusters. Predicting the nucleation rates of aerosol particles in different chemical environments is a key step toward understanding the dynamics of complex aerosol chemistry. Therefore, molecular scale models of nanoclusters are required to understand the macroscopic nucleation process. On the basis of variational transition state theory, DNT provides an efficient approach to predict nucleation kinetics. While most DNT Monte Carlo simulations use analytic potentials to model critical sized clusters, or use ab initio potentials to model very small clusters, the DNTEFP Monte Carlo method presented here can treat up to critical sized clusters using the effective fragment potential (EFP), a rigorous nonempirical intermolecular model potential based on ab initio electronic structure theory calculations, improvable in a systematic manner. The DNTEFP method is applied to study the evaporation rates, energetics, and structure factors of multicomponent clusters containing water and isoprene. The most probable topology of the transition state characterizing the evaporation of one water molecule from a water hexamer at 243 K is predicted to be a conformer that contains six hydrogen bonds, with a topology that corresponds to two water molecules stacked on top of a quadrangular (H(2)O)(4) cluster. For the water hexamer in the presence of isoprene, an increase in the cluster size and a 3-fold increase in the evaporation rate are predicted relative to the reaction in which one water molecule evaporates from a water hexamer cluster.     

Asynchronous multicanonical basin hopping method and its application to cobalt nanoclusters

The multicanonical basin hopping (MUBH) method, which uses a multicanonical weight in the basin hopping (BH) Monte Carlo method, was found to be very efficient for global optimization of large-scale systems such as Lennard-Jones clusters containing more than 150 atoms. We have implemented an asynchronous parallel version of the MUBH method using the message passing interface (MPI) to take advantage of the full usage of multiprocessors in either a homogeneous or heterogeneous computational environment. Based on the intrinsic properties of the Monte Carlo method, this MPI implementation used the task parallelism to minimize interthread data communication. For a Co nanocluster consisting of N atoms, we have applied the asynchronous multicanonical basin hopping (AMUBH) method (for 181 < N < or = 200), together with BH (for 2 < or = N < 150) and MUBH (for 150 < or = N < or = 180), to search for the molecular configuration of the global energy minimum. AMUBH becomes the only practical computational scheme for locating the energy minimum within realistic computational time for a relatively large cluster.     

Modeling the first stages of Cu precipitation in α-Fe using a hybrid atomistic kinetic Monte Carlo approach

We simulate the coherent stage of Cu precipitation in α-Fe with an atomistic kinetic Monte Carlo (AKMC) model. The vacancy migration energy as a function of the local chemical environment is provided on-the-fly by a neural network, trained with high precision on values calculated with the nudged elastic band method, using a suitable interatomic potential. To speed up the simulation, however, we modify the standard AKMC algorithm by treating large Cu clusters as objects, similarly to object kinetic Monte Carlo approaches. Seamless matching between the fully atomistic and the coarse-grained approach is achieved again by using a neural network, that provides all stability and mobility parameters for large Cu clusters, after training on atomistically informed results. The resulting hybrid algorithm allows long thermal annealing experiments to be simulated, within a reasonable CPU time. The results obtained are in very good agreement with several series of experimental data available from the literature, spanning over different conditions of temperature and alloy composition. We deduce from these results and relevant parametric studies that the mobility of Cu clusters containing one vacancy plays a central role in the precipitation mechanism.     

Theoretical modeling of ionization energies of argon clusters: nuclear delocalization effects

Temperature dependence of vertical ionization energies is modeled for small argon clusters (N ≤ 13) using classical parallel-tempering Monte Carlo methods and extended interaction models based on the diatomics-in-molecules approach. Quantum effects at the zero temperature are also discussed in terms of zero-point nuclear vibrations, either at the harmonic approximation level or at the fully anharmonic level using the diffusion Monte Carlo calculations. Both approaches lead to a considerable improvement of the theoretical predictions of argon clusters ionization energies and represent a realistic way of modeling of ionization energies for weakly bound and floppy complexes in general. A thorough comparison with a recent electron-impact experiment [O. Echt et al., J. Chem. Phys. 123, 084313 (2005)] is presented and a novel interpretation of the experimental data is proposed.     

Multiple solvation configurations around phthalocyanine in helium droplets

Recent measurements of the emission spectrum of phthalocyanine solvated in superfluid helium nanodroplets exhibit a constant 10.3 cm(-1) splitting of each emission line relative to the absorption spectrum. This splitting has been attributed to two distinct helium environments near the surface of the phthalocyanine molecule. Rigid-body path-integral Monte Carlo provides a means of investigating the origin of the splitting on a detailed microscopic level. Path-integral Monte Carlo simulations of 4He(N)-phthalocyanine at 0.625 K with N ranging from 24 to 150 show two distinct helium configurations. One configuration is commensurate with the molecular substrate and the other is a triangular lattice. We investigate the energetics of these two configurations and use a method for calculating electronic spectral shifts for aromatic molecule-rare-gas clusters due to dispersive interactions to estimate the spectral splitting that would arise from the two helium configurations seen for N=150. The results are in reasonable agreement with the experimentally measured splitting, supporting the existence of two distinct local helium environments near the surface of the molecule in the nanodroplets.     

Quantum Monte Carlo simulations of selected ammonia clusters (n = 2-5): isotope effects on the ground state of typical hydrogen bonded systems

Variational Monte Carlo, diffusion Monte Carlo, and stereographic projection path integral simulations are performed on eight selected species from the (NH(3))(n), (ND(3))(n), (NH(2)D)(n), and (NH(3))(n-1)(ND(3)) clusters. Each monomer is treated as a rigid body with the rotation spaces mapped by the stereographic projection coordinates. We compare the energy obtained from path integral simulations at several low temperatures with those obtained by diffusion Monte Carlo, for two dimers, and we find that at 4 K, the fully deuterated dimer energy is in excellent agreement with the ground state energy of the same. The ground state wavefunction for the (NH(3))(2-5) clusters is predominantly localized in the global minimum of the potential energy. In all simulations of mixed isotopic substitutions, we find that the heavier isotope is almost exclusively the participant in the hydrogen bond.     

Temperature measurement from scattering spectra of clusters: theoretical treatment

Scattering spectra from X-ray, electron or neutron diffraction experiments are sufficient to describe the phase behaviour of noble gas clusters and to determine their temperature. Using classical Monte Carlo simulations combined with optimized data analysis and Path Integral Monte Carlo calculations as “idealized experiments” we obtain scattering spectra of Ar- and Ne-clusters. Starting from the classical and quantum mechanical hypervirial theorems we devise a method to estimate the temperature and the caloric curves (which describe the phase behaviour of the noble gas clusters) directly from these scattering spectra using an interatomic potential function as input. As applications we studied for Ar-clusters the effect of different model potentials on the temperature estimate thus contributing to the intricate question of what experimentally is the temperature of an isolated cluster. For Ne-clusters we investigate the differences between classical and quantum mechanical treatment.     

Cluster cross sections from pickup measurements: are the established methods consistent?

Pickup of several molecules, H(2)O, HBr, and CH(3)OH, and Ar atoms on free Ar(N) clusters has been investigated in a molecular beam experiment. The pickup cross sections of the clusters with known mean sizes, ?≈ 150 and 260 were measured by two independent methods: (i) the cluster beam velocity decrease due to the momentum transfer of the picked up molecules to the clusters, and (ii) Poisson distribution of a selected cluster fragment ion as a function of the pickup pressure. In addition, the pickup cross sections were calculated using molecular dynamics and Monte Carlo simulations. The simulations support the results of the velocity measurements. On the other hand, the Poisson distributions yield significantly smaller cross sections, inconsistent with the known Ar(N) cluster sizes. These results are discussed in terms of: (i) an incomplete coagulation of guest molecules on the argon clusters when two or more molecules are picked up; and (ii) the fragmentation pattern of the embedded molecules and their clusters upon ionization on the Ar cluster. We conclude that the Poisson distribution method has to be cautiously examined, if conclusions should be drawn about the cluster cross section, or the mean cluster size ?, and the number of picked up molecules.     

Excitation levels and magic numbers of small parahydrogen clusters (N<or=40)

The excitation energies of parahydrogen clusters have been systematically calculated by the diffusion Monte Carlo technique in steps of 1 molecule from 3 to 40 molecules. These clusters possess a very rich spectra, with angular momentum excitations arriving up to L=13 for the heavier ones. No regular pattern can be guessed in terms of the angular momenta and the size of the cluster. Clusters with N=13 and 36 are characterized by a peak in the chemical potential and a large energy gap of the first excited level, which indicate the magical character of these clusters. From the calculated excitation energies, the partition function has been obtained, thus allowing for an estimate of thermal effects. An enhanced production is predicted for cluster sizes of N=13, 31, and 36, which is in agreement with the experiment.     

Methanol-water mixtures: a microsolvation study using the effective fragment potential method

The formation process of methanol-water mixtures, (MeOH/H2O)n, n = 2, 3,..., 8 is studied at the molecular level using the general effective fragment potential (EFP2) method and second-order perturbation theory (MP2). Extensive Monte Carlo/simulated annealing global optimizations were used to locate global minimum structures for each n, for both homo and hetero clusters. Mixing at the microscopic level was investigated, and some general conclusions about the microsolvation behavior of these mixtures are presented. For all of these clusters, incomplete mixing is observed at the molecular level. Low-energy (MeOH/H2O)n clusters retain much of their initial structure in the global minima of the mixed clusters.     

Rotational excitations of N2O in small helium clusters and the role of Bose permutation symmetry

We present a detailed study of the energetics, structures, and Bose properties of small clusters of (4)He containing a single nitrous oxide (N(2)O) molecule, from N=1 (4)He up to sizes corresponding to completion of the first solvation shell around N(2)O (N=16 (4)He). Ground state properties are calculated using the importance-sampled rigid-body diffusion Monte Carlo method, rotational excited state calculations are made with the projection operator imaginary time spectral evolution method, and Bose permutation exchange and associated superfluid properties are calculated with the finite temperature path integral method. For N< or =5 the helium atoms are seen to form an equatorial ring around the molecular axis, at N=6 helium density starts to occupy the second (local) minimum of the N(2)O-He interaction at the oxygen side of the molecule, and N=9 is the critical size at which there is onset of helium solvation all along the molecular axis. For N> or =8 six (4)He atoms are distributed in a symmetric, quasirigid ring around N(2)O. Path integral calculations show essentially complete superfluid response to rotation about the molecular axis for N> or =5, and a rise of the perpendicular superfluid response from zero to appreciable values for N> or =8. Rotational excited states are computed for three values of the total angular momentum, J=1-3, and the energy levels fitted to obtain effective spectroscopic constants that show excellent agreement with the experimentally observed N dependence of the effective rotational constant B(eff). The non-monotonic behavior of the rotational constant is seen to be due to the onset of long (4)He permutation exchanges and associated perpendicular superfluid response of the clusters for N> or =8. We provide a detailed analysis of the role of the helium solvation structure and superfluid properties in determining the effective rotational constants.     

A diffusion quantum Monte Carlo study on the lowest singlet and triplet electronic states of BN molecule

Ab initio calculation of both the lowest singlet and triplet electronic states of BN has been performed by the fixed-node Ornstein-Uhlenbeck diffusion quantum Monte Carlo method with the floating spherical Gaussian orbitals and spherical Gaussian geminals. The Monte Carlo calculation gives equilibrium bond lengths and equilibrium harmonic frequencies of 1.3317(7) A and 1529(7) cm(-1), respectively, for the lowest triplet state and 1.2751(7) A and 1709(8) cm(-1), respectively, for the lowest singlet state. Also, the Monte Carlo calculation reports an energy separation of 178(83) cm(-1) between the two electronic states and recommends the ground state is the lowest triplet state.     

Interaction of molecular nitrogen and oxygen with extraframework cations in zeolites with double six-membered rings of oxygen-bridged silicon and aluminum atoms: a DFT study

The interaction of N(2) and O(2) with extraframework cations of zeolite frameworks was studied by DFT, using the B3LYP method. The extraframework cation sites located in the vicinity of the double six-member rings (D6R) of FAU zeolites (SI, SI', SIII') were considered and clusters with composition (M(n)(+))(2/)(n)()H(12)Si(10)Al(2)O(18), M = Li(+), Na(+), K(+), Ca(2+), were selected to represent the adsorption centers. The cation sites SII in the center of single six-membered rings (S6R) were modeled by [M(I)H(12)Si(4)Al(2)O(6)](-) and M(II)H(12)Si(4)Al(2)O(6) clusters. The adsorption energy of N(2) and O(2) is the highest for Li(+) cations at the SIII' cation sites, while for the SI' and SII sites the adsorption energies decrease in the order Ca(2+) > Na(+) > Li(+). The calculated small N(2) adsorption energy for Li(+) cations at SII sites suggests that these sites do not take part in the sorption process in agreement with results of NMR studies and Monte Carlo simulations. The N(2) adsorption complexes with the extraframework cations are linear, while those of O(2) are bent regardless of the extraframework cation location. The SIII' cation sites are the most favorable ones with respect to N(2) adsorption capacity and N(2)/O(2) selectivity; the SII sites are less selective and the SI sites are not accessible.