Neutron matter at zero temperature with an auxiliary field diffusion Monte Carlo method

A. Sarsa, S. Fantoni, Kevin Schmidt, F. Pederiva

Research output: Contribution to journalArticle

63 Citations (Scopus)

Abstract

The recently developed auxiliary field diffusion Monte Carlo method is applied to compute the equation of state and the compressibility of neutron matter. By combining diffusion Monte Carlo method for the spatial degrees of freedom and auxiliary field Monte Carlo method to separate the spin-isospin operators, quantum Monte Carlo can be used to simulate the ground state of many-nucleon systems (As ≤ 100). We use a path constraint to control the fermion sign problem. We have made simulations for realistic interactions, which include tensor and spin-orbit two-body potentials as well as three-nucleon forces. The Argonne v'8 and v'6 two-nucleon potentials plus the Urbana or Illinois three-nucleon potentials have been used in our calculations. We compare with fermion hypernetted chain results. We report on the results of a periodic box fermi hypernetted chain calculation, which is also used to estimate the finite size corrections to our quantum Monte Carlo simulations. Our auxiliary field diffusion Monte Carlo (AFDMC) results for V6, models of pure neutron matter are in reasonably good agreement with equivalent correlated basis function (CBF) calculations, providing energies per particle which are slightly lower than the CBF ones. However, the inclusion of the spin-orbit force leads to quite different results particularly at relatively high densities. The resulting equation of state from AFDMC calculations is harder than the one from previous Fermi hypernetted chain studies commonly used to determine the neutron star structure.

Original languageEnglish (US)
Pages (from-to)243081-2430814
Number of pages2187734
JournalPhysical Review C - Nuclear Physics
Volume68
Issue number2
StatePublished - Aug 2003

Fingerprint

Monte Carlo method
neutrons
nucleon potential
equations of state
fermions
two body problem
temperature
neutron stars
compressibility
boxes
simulation
degrees of freedom
inclusions
tensors
orbits
operators
ground state
estimates
interactions
energy

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Nuclear and High Energy Physics

Cite this

Neutron matter at zero temperature with an auxiliary field diffusion Monte Carlo method. / Sarsa, A.; Fantoni, S.; Schmidt, Kevin; Pederiva, F.

In: Physical Review C - Nuclear Physics, Vol. 68, No. 2, 08.2003, p. 243081-2430814.

Research output: Contribution to journalArticle

@article{7c273fa07bc0418da9f9a6ab3d86df1b,
title = "Neutron matter at zero temperature with an auxiliary field diffusion Monte Carlo method",
abstract = "The recently developed auxiliary field diffusion Monte Carlo method is applied to compute the equation of state and the compressibility of neutron matter. By combining diffusion Monte Carlo method for the spatial degrees of freedom and auxiliary field Monte Carlo method to separate the spin-isospin operators, quantum Monte Carlo can be used to simulate the ground state of many-nucleon systems (As ≤ 100). We use a path constraint to control the fermion sign problem. We have made simulations for realistic interactions, which include tensor and spin-orbit two-body potentials as well as three-nucleon forces. The Argonne v'8 and v'6 two-nucleon potentials plus the Urbana or Illinois three-nucleon potentials have been used in our calculations. We compare with fermion hypernetted chain results. We report on the results of a periodic box fermi hypernetted chain calculation, which is also used to estimate the finite size corrections to our quantum Monte Carlo simulations. Our auxiliary field diffusion Monte Carlo (AFDMC) results for V6, models of pure neutron matter are in reasonably good agreement with equivalent correlated basis function (CBF) calculations, providing energies per particle which are slightly lower than the CBF ones. However, the inclusion of the spin-orbit force leads to quite different results particularly at relatively high densities. The resulting equation of state from AFDMC calculations is harder than the one from previous Fermi hypernetted chain studies commonly used to determine the neutron star structure.",
author = "A. Sarsa and S. Fantoni and Kevin Schmidt and F. Pederiva",
year = "2003",
month = "8",
language = "English (US)",
volume = "68",
pages = "243081--2430814",
journal = "Physical Review C - Nuclear Physics",
issn = "0556-2813",
publisher = "American Physical Society",
number = "2",

}

TY - JOUR

T1 - Neutron matter at zero temperature with an auxiliary field diffusion Monte Carlo method

AU - Sarsa, A.

AU - Fantoni, S.

AU - Schmidt, Kevin

AU - Pederiva, F.

PY - 2003/8

Y1 - 2003/8

N2 - The recently developed auxiliary field diffusion Monte Carlo method is applied to compute the equation of state and the compressibility of neutron matter. By combining diffusion Monte Carlo method for the spatial degrees of freedom and auxiliary field Monte Carlo method to separate the spin-isospin operators, quantum Monte Carlo can be used to simulate the ground state of many-nucleon systems (As ≤ 100). We use a path constraint to control the fermion sign problem. We have made simulations for realistic interactions, which include tensor and spin-orbit two-body potentials as well as three-nucleon forces. The Argonne v'8 and v'6 two-nucleon potentials plus the Urbana or Illinois three-nucleon potentials have been used in our calculations. We compare with fermion hypernetted chain results. We report on the results of a periodic box fermi hypernetted chain calculation, which is also used to estimate the finite size corrections to our quantum Monte Carlo simulations. Our auxiliary field diffusion Monte Carlo (AFDMC) results for V6, models of pure neutron matter are in reasonably good agreement with equivalent correlated basis function (CBF) calculations, providing energies per particle which are slightly lower than the CBF ones. However, the inclusion of the spin-orbit force leads to quite different results particularly at relatively high densities. The resulting equation of state from AFDMC calculations is harder than the one from previous Fermi hypernetted chain studies commonly used to determine the neutron star structure.

AB - The recently developed auxiliary field diffusion Monte Carlo method is applied to compute the equation of state and the compressibility of neutron matter. By combining diffusion Monte Carlo method for the spatial degrees of freedom and auxiliary field Monte Carlo method to separate the spin-isospin operators, quantum Monte Carlo can be used to simulate the ground state of many-nucleon systems (As ≤ 100). We use a path constraint to control the fermion sign problem. We have made simulations for realistic interactions, which include tensor and spin-orbit two-body potentials as well as three-nucleon forces. The Argonne v'8 and v'6 two-nucleon potentials plus the Urbana or Illinois three-nucleon potentials have been used in our calculations. We compare with fermion hypernetted chain results. We report on the results of a periodic box fermi hypernetted chain calculation, which is also used to estimate the finite size corrections to our quantum Monte Carlo simulations. Our auxiliary field diffusion Monte Carlo (AFDMC) results for V6, models of pure neutron matter are in reasonably good agreement with equivalent correlated basis function (CBF) calculations, providing energies per particle which are slightly lower than the CBF ones. However, the inclusion of the spin-orbit force leads to quite different results particularly at relatively high densities. The resulting equation of state from AFDMC calculations is harder than the one from previous Fermi hypernetted chain studies commonly used to determine the neutron star structure.

UR - http://www.scopus.com/inward/record.url?scp=0141918669&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0141918669&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0141918669

VL - 68

SP - 243081

EP - 2430814

JO - Physical Review C - Nuclear Physics

JF - Physical Review C - Nuclear Physics

SN - 0556-2813

IS - 2

ER -