Quantum Monte Carlo calculations of symmetric nuclear matter.

We present an accurate numerical study of the equation of state of nuclear matter based on realistic nucleon-nucleon interactions by means of auxiliary field diffusion Monte Carlo (AFDMC) calculations. The AFDMC method samples the spin and isospin degrees of freedom allowing for quantum simulations of large nucleonic systems and represents an important step forward towards a quantitative understanding of problems in nuclear structure and astrophysics.

1S0 Superfluid phase transition in neutron matter with realistic nuclear potentials and modern many-body theories.

The 1S0 pairing in neutron matter is studied using realistic two- and three-nucleon interactions. The auxiliary field diffusion Monte Carlo method and correlated basis function theory are employed to get quantitative and reliable estimates of the gap. The two methods are in good agreement up to the maximum gap density and both point to a slight reduction with respect to the standard BCS value. In fact, the maximum gap is about 2.5 MeV at kF approximately 0.8 fm(-1) in BCS and 2.2-2.4 MeV at kF approximately 0.6 fm(-1)in correlated matter. In general, the computed medium polarization effects are much smaller than those previously estimated within all theories. Truncations of Argonne to simpler forms give the same gaps in BCS, provided the truncated potentials have been refitted to the same data set. The three-nucleon interaction provides an additional increase of the gap of about 0.35 MeV.

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.

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.