Source author record

Stefano Gandolfi

Stefano Gandolfi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

Atomic nuclei from quantum Monte Carlo calculations with chiral EFT interactions

Quantum Monte Carlo methods are powerful numerical tools to accurately solve the Schrödinger equation for nuclear systems, a necessary step to describe the structure and reactions of nuclei and nucleonic matter starting from realistic interactions and currents. These ab-initio methods have been used to accurately compute properties of light nuclei -- including their spectra, moments, and transitions -- and the equation of state of neutron and nuclear matter. In this work we review selected results obtained by combining quantum Monte Carlo methods and recent Hamiltonians constructed within chiral effective field theory.

preprint2020arXiv

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass X-ray binaries (QLMXBs), photospheric radius expansion X-ray bursts (PREs), and X-ray timing observations of J0030+0451. With this data set, the form of the prior distribution still has an impact on the posterior mass-radius (MR) curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

preprint2019arXiv

Quasielastic lepton scattering and back-to-back nucleons in the short-time approximation

Understanding quasielastic electron- and neutrino-scattering from nuclei has taken on new urgency with current and planned neutrino oscillation experiments, and with electron scattering experiments measuring specific final states, such as those involving nucleon pairs in ``back-to-back'' configurations. Accurate many-body methods are available for calculating the response of light ($A \leq 12$) nuclei to electromagnetic and weak probes, but they are computationally intensive and only applicable to the inclusive response. In the present work we introduce a novel approach, based on realistic models of nuclear interactions and currents, to evaluate the short-time (high-energy) inclusive and exclusive response of nuclei. The approach accounts reliably for crucial two-nucleon dynamics, including correlations and currents, and provides information on back-to-back nucleons observed in electron and neutrino scattering experiments. We demonstrate that in the quasielastic regime and at moderate momentum transfers both initial- and final-state correlations, and two-nucleon currents are important for a quantitatively successful description of the inclusive response and final state nucleons. Finally, the approach can be extended to include relativistic---kinematical and dynamical---effects, at least approximately in the two-nucleon sector, and to describe the response in the resonance-excitation region.

preprint2016arXiv

Diffusion Monte Carlo study of strongly interacting two-dimensional Fermi gases

Ultracold atomic Fermi gases have been a popular topic of research, with attention being paid recently to two-dimensional (2D) gases. In this work, we perform T=0 ab initio diffusion Monte Carlo calculations for a strongly interacting two-component Fermi gas confined to two dimensions. We first go over finite-size systems and the connection to the thermodynamic limit. After that, we illustrate pertinent 2D scattering physics and properties of the wave function. We then show energy results for the strong-coupling crossover, in between the Bose-Einstein Condensation (BEC) and Bardeen-Cooper-Schrieffer (BCS) regimes. Our energy results for the BEC-BCS crossover are parametrized to produce an equation of state, which is used to determine Tan's contact. We carry out a detailed comparison with other microscopic results. Finally, we calculate the pairing gap for a range of interaction strengths in the strong coupling regime, following from variationally optimized many-body wave functions.

preprint2016arXiv

Vortex line in the unitary Fermi gas

We report diffusion Monte Carlo results for the ground state of unpolarized spin-1/2 fermions in a cylindrical container and properties of the system with a vortex-line excitation. The density profile of the system with a vortex line presents a non-zero density at the core. We calculate the ground-state energy per particle, the superfluid pairing gap, and the excitation energy per particle. These simulations can be extended to calculate the properties of vortex excitations in other strongly interacting systems, such as superfluid neutron matter using realistic nuclear Hamiltonians.

preprint2015arXiv

Hyperon Puzzle: Hints from Quantum Monte Carlo Calculations

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 $Λ$ 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.

preprint2015arXiv

Neutron Matter from Low to High Density

Neutron matter is an intriguing nuclear system with multiple connections to other areas of physics. Considerable progress has been made over the last two decades in exploring the properties of pure neutron fluids. Here we begin by reviewing work done to explore the behavior of very low density neutron matter, which forms a strongly paired superfluid and is thus similar to cold Fermi atoms, though at energy scales differing by many orders of magnitude. We then increase the density, discussing work that ties the study of neutron matter with the determination of the properties of neutron-rich nuclei and neutron-star crusts. After this, we review the impact neutron matter at even higher densities has on the mass-radius relation of neutron stars, thereby making contact with astrophysical observations.

preprint2015arXiv

New insights on the hyperon puzzle from quantum Monte Carlo calculations

In the last years auxiliary field diffusion Monte Carlo has been used to assess the properties of hypernuclear systems, from light- to medium-heavy hypernuclei and hyper-neutron matter. One of the main findings is the key role played by the three-body hyperon-nucleon-nucleon interaction in the determination of the hyperon separation energy of hypernuclei and as a possible solution to the hyperon puzzle. However, there are still aspects of the employed hypernuclear potential that remain to be carefully investigated. In this paper we show that the isospin dependence of the Lambda-NN force, which is crucial in determining the NS structure, is poorly constrained by the available experimental data.

preprint2015arXiv

Using Neutron Star Observations to Determine Crust Thicknesses, Moments of Inertia, and Tidal Deformabilities

We perform a systematic assessment of models for the equation of state (EOS) of dense matter in the context of recent neutron star mass and radius measurements to obtain a broad picture of the structure of neutron stars. We demonstrate that currently available neutron star mass and radius measurements provide strong constraints on moments of inertia, tidal deformabilities, and crust thicknesses. A measurement of the moment of inertia of PSR J0737-3039A with 10% error, without any other information from observations, will constrain the EOS over a range of densities to within 50%$-$60%. We find tidal deformabilities between 0.6 and $6\times 10^{36}$ g cm$^{2}$ s$^{2}$ (to 95% confidence) for $M=1.4~\mathrm{M}_{\odot}$, and any measurement which constrains this range will provide an important constraint on dense matter. The crustal fraction of the moment of inertia can be as large as 10% for $M=1.4~\mathrm{M}_{\odot}$ permitting crusts to have a large enough moment of inertia reservoir to explain glitches in the Vela pulsar even with a large amount of superfluid entrainment. Finally, due to the uncertainty in the equation of state, there is at least a 40% variation in the thickness of the crust for a fixed mass and radius, which implies that future simulations of the cooling of a neutron star crust which has been heated by accretion will need to take this variation into account.

preprint2014arXiv

Accurate determination of the interaction between $Λ$ hyperons and nucleons from auxiliary field diffusion Monte Carlo calculations

An accurate assessment of the hyperon-nucleon interaction is of great interest in view of recent observations of very massive neutron stars. The challenge is to build a realistic interaction that can be used over a wide range of masses and in infinite matter starting from the available experimental data on the binding energy of light hypernuclei. To this end, accurate calculations of the hyperon binding energy in a hypernucleus are necessary. We present a quantum Monte Carlo study of $Λ$ and $ΛΛ$ hypernuclei up to $A=91$. We investigate the contribution of two- and three-body $Λ$-nucleon forces to the $Λ$ binding energy. Ground state energies are computed solving the Schrödinger equation for non-relativistic baryons by means of the auxiliary field diffusion Monte Carlo algorithm extended to the hypernuclear sector. We show that a simple adjustment of the parameters of the $ΛNN$ three-body force yields a very good agreement with available experimental data over a wide range of hypernuclear masses. In some cases no experiments have been performed yet, and we give new predictions. The newly fitted $ΛNN$ force properly describes the physics of medium-heavy $Λ$ hypernuclei, correctly reproducing the saturation property of the hyperon separation energy.

preprint2014arXiv

From hypernuclei to the Inner Core of Neutron Stars: A Quantum Monte Carlo Study

Auxiliary Field Diffusion Monte Carlo (AFDMC) calculations have been employed to revise the interaction between $Λ$-hyperons and nucleons in hypernuclei. The scheme used to describe the interaction, inspired by the phenomenological Argonne-Urbana forces, is the $ΛN+ΛNN$ potential firstly introduced by Bodmer, Usmani et al.. Within this framework, we performed calculations on light and medium mass hypernuclei in order to assess the extent of the repulsive contribution of the three-body part. By tuning this contribution in order to reproduce the $Λ$ separation energy in $^5_Λ$He and $^{17}_{~Λ}$O, experimental findings are reproduced over a wide range of masses. Calculations have then been extended to $Λ$-neutron matter in order to derive an analogous of the symmetry energy to be used in determining the equation of state of matter in the typical conditions found in the inner core of neutron stars.

preprint2013arXiv

Computational Nuclear Quantum Many-Body Problem: The UNEDF Project

The UNEDF project was a large-scale collaborative effort that applied high-performance computing to the nuclear quantum many-body problem. UNEDF demonstrated that close associations among nuclear physicists, mathematicians, and computer scientists can lead to novel physics outcomes built on algorithmic innovations and computational developments. This review showcases a wide range of UNEDF science results to illustrate this interplay.

preprint2013arXiv

From nuclear structure to neutron stars

Recent progress in quantum Monte Carlo with modern nucleon-nucleon interactions have enabled the successful description of properties of light nuclei and neutron-rich matter. As a demonstration, we show that the agreement between theoretical calculations of the charge form factor of 12C and the experimental data is excellent. Applying similar methods to isospin-asymmetric systems allows one to describe neutrons confined in an external potential and homogeneous neutron-rich matter. Of particular interest is the nuclear symmetry energy, the energy cost of creating an isospin asymmetry. Combining these advances with recent observations of neutron star masses and radii gives insight into the equation of state of neutron-rich matter near and above the saturation density. In particular, neutron star radius measurements constrain the derivative of the symmetry energy.

preprint2013arXiv

Neutral current interactions of low-energy neutrinos in dense neutron matter

We report the results of a calculation of the response of cold neutron matter to neutral-current interactions with low energy neutrinos, carried out using an effective interaction and effective operators consistently derived within the formalism of Correlated Basis Functions. The neutrino mean free path obtained from the calculated responses turns out to be strongly affected by both short and long range correlations, leading to a sizable increase with respect to the prediction of the Fermi gas model. The consistency between the proposed approach and Landau theory of normal Fermi liquids also has been investigated, using a set of Landau parameters obtained from the matrix elements of the effective interaction.

preprint2013arXiv

Precise determination of the structure factor and contact in a unitary Fermi gas

We present a high-precision determination of the universal contact parameter in a strongly interacting Fermi gas. In a trapped gas at unitarity we find the contact to be $3.06 \pm 0.08$ at a temperature of 0.08 of the Fermi temperature in a harmonic trap. The contact governs the high-momentum (short-range) properties of these systems and this low temperature measurement provides a new benchmark for the zero temperature homogeneous contact. The experimental measurement utilises Bragg spectroscopy to obtain the dynamic and static structure factors of ultracold Fermi gases at high momentum in the unitarity and molecular Bose-Einstein condensate (BEC) regimes. We have also performed quantum Monte Carlo calculations of the static properties, extending from the weakly coupled Bardeen-Cooper-Schrieffer (BCS) regime to the strongly coupled BEC case, which show agreement with experiment at the level of a few percent.

preprint2013arXiv

Quantum Monte Carlo study of strongly interacting Fermi gases

In recent years Quantum Monte Carlo techniques provided to be a valuable tool to study strongly interacting Fermi gases at zero temperature. We have used QMC methods to investigate several properties of the two-components Fermi gas at unitarity and in the BCS-BEC crossover, both with equal and unequal masses corresponding to the $Li-K$ Fermi mixture. In this paper we present several recent QMC results, including the energy at zero and finite effective range, the contact parameter and the static structure factor, which, at low momentum, depends strongly on the phonons in the unitary Fermi gas.

preprint2012arXiv

Effective-Range Dependence of Resonantly Interacting Fermions

We extract the leading effective range corrections to the equation of state of the unitary Fermi gas from ab initio fixed-node quantum Monte Carlo (FNQMC) calculations in a periodic box using a density functional theory (DFT), and show them to be universal by considering several two-body interactions. Furthermore, we find that the DFT is consistent with the best available unbiased QMC calculations, analytic results, and experimental measurements of the equation of state. We also discuss the asymptotic effective-range corrections for trapped systems and present the first QMC results with the correct asymptotic scaling.

preprint2012arXiv

Quantum Monte Carlo approaches to nuclear and atomic physics

Quantum Monte Carlo methods have proven to be valuable in the study of strongly correlated quantum systems, particularly nuclear physics and cold atomic gases. Historically, such ab initio simulations have been used to study properties of light nuclei, including spectra and form factors, low-energy scattering, and high-momentum properties including inclusive scattering and one- and two-body momentum distributions. More recently they have been used to study the properties of homogeneous and inhomogeneous neutron matter and cold atomic gases. There are close analogies between these seemingly diverse systems, including the equation of state, superfluid pairing, and linear response to external probes. In this paper, we compare and contrast results found in nuclear and cold atom physics. We show updated lattice results for the energy of the homogeneous unitary Fermi gas and comparisons with neutron matter, as well as for the dependence of the cold atom energy on the mass ratio between paired particles, which yields insights on the structure of the ground state. We also provide new lattice and continuum results for the harmonically trapped unitary gas, again comparing neutron matter and cold atoms.

preprint2012arXiv

Quantum Monte Carlo study of inhomogeneous neutron matter

We present an ab-initio study of neutron drops. We use Quantum Monte Carlo techniques to calculate the energy up to 54 neutrons in different external potentials, and we compare the results with Skyrme forces. We also calculate the rms radii and radial densities, and we find that a re-adjustment of the gradient term in Skyrme is needed in order to reproduce the properties of these systems given by the ab-initio calculation. By using the ab-initio results for neutron drops for close- and open-shell configurations, we suggest how to improve Skyrme forces when dealing with systems with large isospin-asymmetries like neutron-rich nuclei.

preprint2012arXiv

The equation of state of neutron star matter and the symmetry energy

We present an overview of microscopical calculations of the Equation of State (EOS) of neutron matter performed using Quantum Monte Carlo techniques. We focus to the role of the model of the three-neutron force in the high-density part of the EOS up to a few times the saturation density. We also discuss the interplay between the symmetry energy and the neutron star mass-radius relation. The combination of theoretical models of the EOS with recent neutron stars observations permits us to constrain the value of the symmetry energy and its slope. We show that astrophysical observations are starting to provide important insights into the properties of neutron star matter.

preprint2011arXiv

Auxiliary Field quantum Monte Carlo for Strongly Paired Fermions

We demonstrate that the inclusion of a BCS importance function dramatically increases the efficiency of the auxiliary field method for strong pairing. We calculate the ground-state energy of an unpolarized fermi gas at unitarity with up to 66 particles and lattices of up to $27^3$ sites. The method has no fermion sign problem, and an accurate result is obtained for the universal parameter $ξ$. Several different forms of the kinetic energy adjusted to the unitary limit but with different effective ranges extrapolate to the same continuum limit within error bars. The finite effective range results for different interactions are consistent with a linear term proportional to the Fermi momentum times the effective range. The new method described herein will have many applications in superfluid cold atom systems and in both electronic and nuclear structures when pairing is important.

preprint2011arXiv

Recent progress on the accurate determination of the equation of state of neutron and nuclear matter

The problem of accurately determining the equation of state of nuclear and neutron matter at density near and beyond saturation is still an open challenge. In this paper we will review the most recent progress made by means of Quantum Monte Carlo calculations, which are at present the only ab-inito method capable to treat a sufficiently large number of particles to give meaningful estimates depending only on the choice of the nucleon-nucleon interaction. In particular, we will discuss the introduction of density-dependent interactions, the study of the temperature dependence of the equation of state, and the possibility of accurately studying the effect of the onset of hyperons by developing an accurate hyperon-nucleon and hyperon-nucleon-nucleon interaction.

preprint2011arXiv

Resonantly Interacting Fermions In a Box

We use two fundamental theoretical frameworks to study the finite-size (shell) properties of the unitary gas in a periodic box: 1) an ab initio Quantum Monte Carlo (QMC) calculation for boxes containing 4 to 130 particles provides a precise and complete characterization of the finite-size behavior, and 2) a new Density Functional Theory (DFT) fully encapsulates these effects. The DFT predicts vanishing shell structure for systems comprising more than 50 particles, and allows us to extrapolate the QMC results to the thermodynamic limit, providing the tightest bound to date on the ground-state energy of the unitary gas: ξ_S <= 0.383(1). We also apply the new functional to few-particle harmonically trapped systems, comparing with previous calculations.