Finite temperature path integral Monte Carlo simulations of structural and dynamical properties of Ar(N)-CO(2) clusters

We report finite temperature quantum mechanical simulations of structural and dynamical properties of Ar(N)-CO(2) clusters using a path integral Monte Carlo algorithm. The simulations are based on a newly developed analytical Ar-CO(2) interaction potential obtained by fitting ab initio results to an anisotropic two-dimensional Morse∕Long-range function. The calculated distributions of argon atoms around the CO(2) molecule in Ar(N)-CO(2) clusters with different sizes are consistent to the previous studies of the configurations of the clusters. A first-order perturbation theory is used to quantitatively predict the CO(2) vibrational frequency shift in different clusters. The first-solvation shell is completed at N = 17. Interestingly, our simulations for larger Ar(N)-CO(2) clusters showed several different structures of the argon shell around the doped CO(2) molecule. The observed two distinct peaks (2338.8 and 2344.5 cm(-1)) in the υ(3) band of CO(2) may be due to the different arrangements of argon atoms around the dopant molecule.     

Predicting atomic dopant solvation in helium clusters: the MgHe(n) case

We present a quantum Monte Carlo study of the solvation and spectroscopic properties of the Mg-doped helium clusters MgHe(n) with n=2-50. Three high-level [MP4, CCSD(T), and CCSDT] MgHe interaction potentials have been used to study the sensitivity of the dopant location on the shape of the pair interaction. Despite the similar MgHe well depth, the pair distribution functions obtained in the diffusion Monte Carlo simulations markedly differ for the three pair potentials, therefore indicating different solubility properties for Mg in He(n). Moreover, we found interesting size effects for the behavior of the Mg impurity. As a sensitive probe of the solvation properties, the Mg excitation spectra have been simulated for various cluster sizes and compared with the available experimental results. The interaction between the excited 1P Mg atom and the He moiety has been approximated using the diatomics-in-molecules method and the two excited 1pi and 1sigma MgHe potentials. The shape of the simulated MgHe50 spectra shows a substantial dependency on the location of the Mg impurity, and hence on the MgHe pair interaction employed. To unravel the dependency of the solvation behavior on the shape of the computed potentials, exact density-functional theory has been adapted to the case of doped He(n) and various energy distributions have been computed. The results indicate the shape of the repulsive part of the MgHe potential as an important cause of the different behaviors.     

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.     

Excited Li and Na in He(n): influence of the dimer potential energy curves

The X(2)Σ ground and the A(2)Π and B(2)Σ first two excited states of Li-He and Na-He are determined using high level complete active space self-consistent field-multireference configuration interaction ab initio method. The obtained potentials differ from the ones proposed by Pascale [Phys. Rev. A 28, 632 (1983)], more strongly for the ground than for the excited states. Quantum diffusion Monte Carlo studies of small Li(?)He(n) and Na(?)He(n) with n ≤ 5 are performed using a diatomics-in-molecule approach to model the non-pair additive interaction potential. The sensitivity of our results to the A(2)Π and B(2)Σ potentials used is assessed by an analysis of the structure and of the energetics of the clusters. For these small clusters, the physical conclusions are essentially independent of the diatomic curves employed.     

Anionic microsolvation in helium droplets: OH- (He)N structures from classical and quantum calculations

Diffusion Monte Carlo calculations are carried out for clusters of OH- (1Sigma+) with N 4He atoms, N varying up to 15, while classical configurations from a genetic algorithm optimization are obtained up to N=20. The overall interaction potential is assembled from ab initio data for the partners using the sum-of-potentials scheme. In contrast with the cationic dopants' behavior, the results indicate a very marked spatial delocalization and quantum features of the solvent adatoms surrounding the anionic impurity, thus making classical calculations of solvent's spatial locations of only limited use. In spite of the generally known repulsive interaction of negative charges in He droplets, the calculations show that this polar molecular anion is solvated by a liquidlike solvent layer, reminiscent of what happens in pure helium droplets.     

Solubility of metal atoms in helium droplets: exploring the effect of the well depth using the coinage metals Cu and Ag

We report a theoretical investigation of the solution properties of Cu and Ag atoms dissolved in He clusters. Employing our recent ab initio ground state pair potential for Me-He (Me = Ag, Cu), we simulated the species Me@He (n) (n = 2-100) by means of diffusion Monte Carlo (DMC) obtaining exact information on their energetics and the structural properties. In particular, we investigated the sensitivity of structural details on the well depth of the two interaction potentials. Whereas Ag structures the first He solvation layer similarly, to some extent, to a positive ion such as Na(+), Cu appears to require the onset of a second solvation shell for a similar dense structure to be formed despite an interaction well of 28.4 μhartree. An additional signature of the different solution behavior between Ag and Cu appears also in the dependence of the energy required to evaporate a single He atom on the size of the MeHe(n) clusters. The absorption spectrum for the (2)P ← (2)S excitation of the metals was also simulated employing the semi-classical Lax approximation to further characterize Me@He(n) (n = 2-100) using novel accurate interaction potentials between He and the lowest (2)P state of Ag and Cu in conjunction with the Diatomic-in-Molecules approach. The results indicated that Ag exciplexes should not form via a direct vertical excitation into an attractive region of the excited manifolds and that there is an interesting dependence of the shape of the Cu excitation bands on the local structure of the first solvation shell.     

Structure of small molecular hydrogen clusters

The ground state energies and structural properties of small (H₂) _N ,N=2?7, are calculated using the variational Monte Carlo method. These wavefunctions include both short- and long-range correlation effects that are important in the binding of van der Waals clusters. We have investigated these clusters using shadow wavefunctions and found that the coupling to shadow variables raises the energy in all cases, implying that the ground states of these small clusters are properly described as quantum liquids rather than solid structures.     

Theoretical studies on the bonding and thermodynamic properties of Ge n Si m (m+n=5) clusters: the precursors of germanium/silicon nanomaterials

Theoretical studies on the Ge n Si m clusters have been carried out using advanced ab initio approaches. The lowest energy isomers were determined for the clusters with compositions n+m=2-5. All possible isomers arising due to permutations of Ge and Si atoms were investigated. The L-shaped structure for the trimers, tetragonal with diagonal bond for tetramers, and a trigonal bipyramid for pentamers represent the energy optimized ground state geometries. The bonding analyses revealed that the trimers and tetramers are stabilized through multicenter pi bonding. In pentamers, this stabilizing factor is eliminated due to the further cluster growth. The ionization of clusters does not change their geometrical characteristics. The agreement of the calculated ionization and atomization energies with those obtained from the mass spectrometric studies (through estimated appearance potential) validated the reported structures of the clusters. The bonding properties of these species are discussed using their molecular orbital characteristics and analysis of natural bond orbital population data.     

HF in clusters of molecular hydrogen: II. Quantum solvation by H2 isotopomers, cluster rigidity, and comparison with CO-doped parahydrogen clusters

We present a comprehensive theoretical study of the quantum solvation of the HF molecule by small clusters of the H2 isotopomers, p-H2, HD, and o-D2, with up to 13 hydrogen solvent molecules. This complements our earlier work on the HF-doped parahydrogen clusters [H. Jiang and Z. Bacic, J. Chem. Phys. 122, 244306 (2005)]. The ground-state properties of the clusters are calculated exactly using the diffusion Monte Carlo method. Detailed information is obtained regarding the size and isotopomer dependences of the energetics, vibrationally averaged structures, and their rigidity. The rigidity of these clusters is investigated further by analyzing the distributions of their principal moments of inertia from the diffusion Monte Carlo simulations. The clusters are found to be rather rigid, especially when compared with the pure parahydrogen clusters of the same size. Extensive comparison is made with the quantum Monte Carlo results for the CO-doped parahydrogen clusters and significant differences are observed in the size evolution of certain properties, notably the chemical potential.     

Bosonic helium droplets with cationic impurities: onset of electrostriction and snowball effects from quantum calculations

Variational Monte Carlo and diffusion Monte Carlo calculations have been carried out for cations such as Li(+), Na(+), and K(+) as dopants of small helium clusters over a range of cluster sizes up to about 12 solvent atoms. The interaction has been modeled through a sum-of-potential picture that disregards higher order effects beyond atom-atom and atom-ion contributions. The latter were obtained from highly correlated ab initio calculations over a broad range of interatomic distances. This study focuses on two of the most striking features of the microsolvation in a quantum solvent of a cationic dopant: electrostriction and snowball effects. They are discussed here in detail and in relation with the nanoscopic properties of the interaction forces at play within a fully quantum picture of the cluster features.     

Connection between the virial equation of state and physical clusters in a low density vapor

We carry out Monte Carlo simulations of physical Lennard-Jones and water clusters and show that the number of physical clusters in vapor is directly related to the virial equation of state. This relation holds at temperatures clearly below the critical temperatures, in other words, as long as the cluster-cluster interactions can be neglected--a typical assumption used in theories of nucleation. Above a certain threshold cluster size depending on temperature and interaction potential, the change in cluster work of formation can be calculated analytically with the recently proposed scaling law. The breakdown of the scaling law below the threshold sizes is accurately modeled with the low order virial coefficients. Our results indicate that high order virial coefficients can be analytically calculated from the lower order coefficients when the scaling law for cluster work of formation is valid. The scaling law also allows the calculation of the surface tension and equilibrium vapor density with computationally efficient simulations of physical clusters. Our calculated values are in good agreement with those obtained with other methods. We also present our results for the curvature dependent surface tension of water clusters.     

Ground state of small mixed helium and spin-polarized tritium clusters: a quantum Monte Carlo study

We report results for the ground-state energy and structural properties of small (4)He-T↓ clusters consisting of up to four T↓ and eight (4)He atoms. These results have been obtained using very well-known (4)He-(4)He and T↓- T↓ interaction potentials and several models for the (4)He- T↓ interatomic potential. All the calculations have been performed with variational and diffusion Monte Carlo methods. It takes at least three atoms to form a mixed bound state. In particular, for small clusters the binding energies are significantly affected by the precise form of the (4)He- T↓ interatomic potential but the stability limits remain unchanged. The only exception is the (4)He(2)T↓ trimer whose stability in the case of the weakest (4)He- T↓ interaction potential is uncertain while it seems stable for other potentials. The mixed trimer (4)He(T↓)(2), a candidate for the Borromean state, is not bound. All other studied clusters are stable. Some of the weakest bound clusters can be classified as quantum halo as a consequence of having high probability of being in a classically forbidden region.     

Energetics driving the short-range order in CuxPd1-x/Ru(0001) monolayer surface alloys

The energetics determining the distinct short-range order in two-dimensional (2D) monolayer Cu(x)Pd(1-x) surface alloys on a Ru(0001) substrate were investigated by Monte Carlo simulations and density functional theory calculations. Using a 2D lattice gas Hamiltonian based on effective pair interaction (EPI) parameters, the EPIs were derived for different Cu concentrations with Monte Carlo (MC) simulations by comparing with the atomic distributions obtained from atomic resolution STM images and the related Warren-Cowley short-range order parameters (Hoster et al., Phys. Rev. B, 2006, 73 165413). The ground state structures and mixing energies at 0 K derived from these EPIs agree well with mixing energies determined from DFT calculations of different ordered surface alloys. Additional MC simulations yield rather low transition temperatures which explain the absence of ordered 2D phases in the experiments. The consequences of our findings for the use of alloy surfaces and surface alloys as model systems for adsorption and catalytic reaction studies are discussed.     

Molecular impurities in helium clusters

We review recently developed theories of molecular interactions with helium clusters. Ground state energies and structures are obtained from Monte Carlo and density functional calculations for a variety of impurities and a range of cluster sizes. Of particular interest are chemical potentials, location of the impurity, and shape of the helium density surrounding it. Included in the work summarized here are results for light molecules (H₂ and D₂) which interact weakly with He and results for heavier molecules (Cl₂ and SF₆) which interact more strongly with He. Finally, theoretical and experimental results are compared for SF₆ in helium 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.     

Quantum Monte Carlo study of small pure and mixed spin-polarized tritium clusters

We have investigated the stability limits of small spin-polarized clusters consisting of up to ten spin-polarized tritium T downward arrow atoms and the mixtures of T downward arrow with spin-polarized deuterium D downward arrow and hydrogen H downward arrow atoms. All of our calculations have been performed using the variational and diffusion Monte Carlo methods. For clusters with D downward arrow atoms, the released node procedure is used in cases where the wave function has nodes. In addition to the energy, we have also calculated the structure of small clusters using unbiased estimators. Results obtained for pure T downward arrow clusters are in good accordance with previous calculations, confirming that the trimer is the smallest spin-polarized tritium cluster. Our results show that mixed T downward arrow-H downward arrow clusters having up to ten atoms are unstable and that it takes at least three tritium atoms to bind one, two, or three D downward arrow atoms. Among all the considered clusters, we have found no other Borromean states except the ground state of the T downward arrow trimer.     

Importance sampling for quantum Monte Carlo in manifolds: addressing the time scale problem in simulations of molecular aggregates

Several importance sampling strategies are developed and tested for stereographic projection diffusion Monte Carlo in manifolds. We test a family of one parameter trial wavefunctions for variational Monte Carlo in stereographically projected manifolds which can be used to produce importance sampling. We use the double well potential in one dimensional Euclidean space to study systematically sampling issues for diffusion Monte Carlo. We find that diffusion Monte Carlo with importance sampling in manifolds is orders of magnitude more efficient compared to unguided diffusion Monte Carlo. Additionally, diffusion Monte Carlo with importance sampling in manifolds can overcome problems with nonconfining potentials and can suppress quasiergodicity effectively. We obtain the ground state energy and the wavefunction for the Stokmayer trimer.     

Thermodynamic properties of van der Waals fluids from Monte Carlo simulations and perturbative Monte Carlo theory

Computer simulations have been performed for fluids with van der Waals potential, that is, hard spheres with attractive inverse power tails, to determine the equation of state and the excess energy. On the other hand, the first- and second-order perturbative contributions to the energy and the zero- and first-order perturbative contributions to the compressibility factor have been determined too from Monte Carlo simulations performed on the reference hard-sphere system. The aim was to test the reliability of this "exact" perturbation theory. It has been found that the results obtained from the Monte Carlo perturbation theory for these two thermodynamic properties agree well with the direct Monte Carlo simulations. Moreover, it has been found that results from the Barker-Henderson [J. Chem. Phys. 47, 2856 (1967)] perturbation theory are in good agreement with those from the exact perturbation theory.