Quantum Monte Carlo methods have proved valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab initio calculations reproduce many low-lying states, moments, and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. The nuclear interactions and currents are reviewed along with a description of the continuum quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. A variety of results are presented, including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. Low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars are also described. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

The onset of hyperons in the core of neutron stars and the consequent softening of the equation of state have been questioned for a long time. Controversial theoretical predictions and recent astrophysical observations of neutron stars are the grounds for the so-called hyperon puzzle. We calculate the equation of state and the neutron star mass-radius relation of an infinite systems of neutrons and $\mathrm{\Lambda }$ particles by using the auxiliary field diffusion Monte Carlo algorithm. We find that the three-body hyperon-nucleon interaction plays a fundamental role in the softening of the equation of state and for the consequent reduction of the predicted maximum mass. We have considered two different models of three-body force that successfully describe the binding energy of medium mass hypernuclei. Our results indicate that they give dramatically different results on the maximum mass of neutron stars, not necessarily incompatible with the recent observation of very massive neutron stars. We conclude that stronger constraints on the hyperon-neutron force are necessary in order to properly assess the role of hyperons in neutron stars.

We calculated the equation of state of neutron matter at zero temperature by means of the auxiliary field diffusion Monte Carlo (AFDMC) method combined with a fixed-phase approximation. The calculation of the energy was carried out by simulating up to 114 neutrons in a periodic box. Special attention was given to reducing finite-size effects at the energy evaluation by adding to the interaction the effect due to the truncation of the simulation box, and by performing several simulations using different numbers of neutrons. The finite-size effects due to kinetic energy were also checked by employing the twist-averaged boundary conditions. We considered a realistic nuclear Hamiltonian containing modern two- and three-body interactions of the Argonne and Urbana family. The equation of state can be used to compare and calibrate other many-body calculations and to predict properties of neutron stars.

We present a method to investigate the kinetics of protein folding and the dynamics underlying the formation of secondary and tertiary structures during the entire reaction. By writing the solution of the Fokker-Planck equation in terms of a path integral, we derive a Hamilton-Jacobi variational principle from which we are able to compute the most probable pathway of folding. The method is applied to the folding of the Villin headpiece subdomain simulated using a Go model. An initial collapsing phase driven by the initial configuration is followed by a rearrangement phase, in which secondary structures are formed and all computed paths display strong similarities. This completely general method does not require the prior knowledge of any reaction coordinate and is an efficient tool to perform simulations of the entire folding process with available computers.

We show how nuclear effective field theory (EFT) and ab initio nuclear-structure methods can turn input from lattice quantum chromodynamics (LQCD) into predictions for the properties of nuclei. We argue that pionless EFT is the appropriate theory to describe the light nuclei obtained in LQCD simulations carried out at pion masses heavier than the physical pion mass. We solve the EFT using the effective-interaction hyperspherical harmonics and auxiliary-field diffusion Monte Carlo methods. Fitting the three leading-order EFT parameters to the deuteron, dineutron, and triton LQCD energies at $m_{\pi }\approx 800\mathrm{MeV}$, we reproduce the corresponding alpha-particle binding and predict the binding energies of mass-5 and mass-6 ground states.

International audienceWe extend the prediction range of Pionless Effective Field Theory with an analysis of the ground state of 16 O in leading order. To renormalize the theory, we use as input both experimental data and lattice QCD predictions of nuclear observables, which probe the sensitivity of nuclei to increased quark masses. The nuclear many-body Schrödinger equation is solved with the Auxiliary Field Diffusion Monte Carlo method. For the first time in a nuclear quantum Monte Carlo calculation, a linear optimization procedure, which allows us to devise an accurate trial wave function with a large number of variational parameters, is adopted. The method yields a binding energy of 4 He which is in good agreement with experiment at physical pion mass and with lattice calculations at larger pion masses. At leading order we do not find any evidence of a 16 O state which is stable against breakup into four 4 He, although higher-order terms could bind 16 O