Saverio Moroni

Correlation energy and spin polarization in the 2D electron gas

The ground-state energy of the two-dimensional uniform electron gas has been calculated with a fixed-node diffusion Monte Carlo method, including backflow correlations, for a wide range of electron densities as a function of spin polarization. We give a simple analytic representation of the correlation energy which fits our simulation data and includes several known high- and low-density limits. This parametrization provides a reliable local spin density energy functional for two-dimensional systems and an estimate for the spin susceptibility. Within the proposed model for the correlation energy, a weakly first-order polarization transition occurs shortly before Wigner crystallization as the density is lowered.

Structure, rotational dynamics, and superfluidity of small OCS-doped He clusters.

The structural and dynamical properties of carbonyl sulfide (OCS) molecules solvated in helium clusters are studied using reptation quantum Monte Carlo, for cluster sizes n=3-20 He atoms. Computer simulations allow us to establish a relation between the rotational spectrum of the solvated molecule and the structure of the He solvent, and of both with the onset of superfluidity. Our results agree with a recent spectroscopic study of this system and provide a more complex and detailed microscopic picture of this system than inferred from experiments.

Condensation of Helium in Nanotube Bundles

Helium atoms are strongly attracted to the interstitial channels within a bundle of carbon nanotubes. The strong corrugation of the axial potential within a channel can produce a lattice gas system wherein the weak mutual attraction between atoms in neighboring channels induces a transition to an anisotropic condensed phase. At low temperatures, the specific heat of the adsorbate phase (with fewer than 2% of the atoms) greatly exceeds that of the host.

Quantum Monte Carlo study of helium clusters doped with nitrous oxide: quantum solvation and rotational dynamics.

Dynamical and structural properties of small (4)He(N)-N(2)O complexes have been analyzed using ground-state and finite-temperature Monte Carlo simulations. The effective rotational constants resulting from the ground-state calculations are in excellent agreement with the results of a recent spectroscopic study [Y. Xu et al., Phys. Rev. Lett. 91, 163401 (2003)]. After an initial decrease for cluster sizes up to N=8, the rotational constant increases, signaling a transition from a molecular complex to a quantum-solvated system. Such a turnaround is not present in the rotational constants extracted from the finite-temperature Monte Carlo calculations, performed for Boltzmann statistics, thus highlighting the importance of exchange effects to explain the decoupling between a solvated dopant and the helium motion.

Local-spin-density functional for multideterminant density functional theory

(Dated: February 6, 2008)Based on exact limits and quantum Monte Carlo simulations, we obtain, at any density and spinpolarization, an accurate estimate for the energy of a modified homogeneous electron gas where elec-trons repel each other only with a long-range coulombic tail. This allows us to construct an analyticlocal-spin-density exchange-correlation functional appropriate to new, multideterminantal versionsof the density functional theory, where quantum chemistry and approximate exchange-correlationfunctionals are combined to optimally describe both long- and short-range electron correlations.I. INTRODUCTION

Size-consistent variational approaches to nonlocal pseudopotentials: Standard and lattice regularized diffusion Monte Carlo methods revisited

We propose improved versions of the standard diffusion Monte Carlo (DMC) and the lattice regularized diffusion Monte Carlo (LRDMC) algorithms. For the DMC method, we refine a scheme recently devised to treat nonlocal pseudopotential in a variational way. We show that such scheme-when applied to large enough systems-maintains its effectiveness only at correspondingly small enough time-steps, and we present two simple upgrades of the method which guarantee the variational property in a size-consistent manner. For the LRDMC method, which is size-consistent and variational by construction, we enhance the computational efficiency by introducing: (i) an improved definition of the effective lattice Hamiltonian which remains size-consistent and entails a small lattice-space error with a known leading term and (ii) a new randomization method for the positions of the lattice knots which requires a single lattice-space.

Small para-hydrogen clusters doped with carbon monoxide: Quantum Monte Carlo simulations and observed infrared spectra

The structures and rotational dynamics of clusters of a single carbon monoxide molecule solvated in para-hydrogen, (paraH(2))(N)-CO, have been simulated for sizes up to N=17 using the reptation Monte Carlo technique. The calculations indicate the presence of two series of R(0) rotational transitions with J=1<--0 for cold clusters, similar to those predicted and observed in the case of He(N)-CO. Infrared spectra of these clusters have been observed in the region of the C-O stretch ( approximately 2143 cm(-1)) in a pulsed supersonic jet expansion using a tunable diode laser probe. With the help of the calculations, the observed R(0) rotational transitions have been assigned up to N=9 for the b-type series and N=14 for the a-type series. Theory and experiment agree rather well, except that theory tends to overestimate the b-type energies. The (paraH(2))(12)-CO cluster is calculated to be particularly stable and (relatively) rigid, corresponding to completion of the first solvation shell, and it is observed to have the strongest a-type transition.

Computational spectroscopy of helium-solvated molecules : Effective inertia, from small he clusters toward the nanodroplet regime

Accurate computer simulations of the rotational dynamics of linear molecules solvated in He clusters indicate that the large-size (nanodroplet) regime is attained quickly for light rotors (HCN) and slowly for heavy ones (OCS, N2O, and CO2), thus challenging previously reported results. Those results spurred the view that the different behavior of light rotors with respect to heavy ones-including a smaller reduction of inertia upon solvation of the former-would result from the lack of adiabatic following of the He density upon molecular rotation. We have performed computer experiments in which the rotational dynamics of OCS and HCN molecules was simulated using a fictitious inertia appropriate to the other molecule. These experiments indicate that the approach to the nanodroplet regime, as well as the reduction of the molecular inertia upon solvation, is determined by the anistropy of the potential, more than by the molecular weight. Our findings are in agreement with recent infrared and/or microwave experimental data which, however, are not yet totally conclusive by themselves.

Rotational dynamics of CO solvated in small He clusters: A quantum Monte Carlo study

The rotational dynamics of CO single molecules solvated in small He clusters (CO @ HeN) has been studied using reptation quantum Monte Carlo simulations for cluster sizes up to N = 30. Our results are in good agreement with the rotovibrational features of the infrared spectrum recently determined for this system and provide a deep insight into the relation between the structure of the cluster and its dynamics. Simulations for large N also provide a prediction of the effective moment of inertia of CO in the He nanodroplet regime, which has not been measured so far.

Rotational spectrum of cyanoacetylene solvated with helium atoms

The high resolution microwave spectra of He(N)-HCCCN clusters were studied in the size ranges of 1-18 and 25-31. In the absence of an accompanying infrared study, rotational excitation energies were computed by the reptation quantum Monte Carlo method and used to facilitate the search and assignment of R(0) transitions from N > 6, as well as R(1) transitions with N > 1. The assignments in the range of 25-31 are accurate to +/-2 cluster size units, with an essentially certain relative ordering. The rotational transition frequencies decrease with N = 1-6 and then show oscillatory behavior for larger cluster sizes, which is now recognized to be a manifestation of the onset and microscopic evolution of superfluidity. For cluster sizes beyond completion of the first solvation shell the rotational frequencies increase significantly above the large-droplet limit. This behavior, common to other linear molecules whose interaction with He features a strong nearly equatorial minimum, is analyzed using path integral Monte Carlo simulations. The He density in the incipient second solvation shell is shown to open a new channel for long permutation cycles, thus increasing the decoupling of the quantum solvent from the rotation of the dopant molecule.