Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAuxiliary-Field Quantum Monte Carlo Simulations of Neutron Matter in Chiral Effective Field Theory
149
0
0.0
(
0
)
Added by
Gabriel Wlaz{\\l}owski
Publication date
2014
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We present variational Monte Carlo calculations of the neutron matter equation of state using chiral nuclear forces. The ground-state wavefunction of neutron matter, containing non-perturbative many-body correlations, is obtained from auxiliary-field quantum Monte Carlo simulations of up to about 340 neutrons interacting on a 10^3 discretized lattice. The evolution Hamiltonian is chosen to be attractive and spin-independent in order to avoid the fermion sign problem and is constructed to best reproduce broad features of the chiral nuclear force. This is facilitated by choosing a lattice spacing of 1.5 fm, corresponding to a momentum-space cutoff of Lambda = 414 MeV/c, a resolution scale at which strongly repulsive features of nuclear two-body forces are suppressed. Differences between the evolution potential and the full chiral nuclear interaction (Entem and Machleidt Lambda = 414 MeV) are then treated perturbatively. Our results for the equation of state are compared to previous quantum Monte Carlo simulations which employed chiral two-body forces at next-to-next-to-leading order (N2LO). In addition we include the effects of three-body forces at N2LO, which provide important repulsion at densities higher than 0.02 fm^-3, as well as two-body forces at N3LO.
rate research
Read More
We show that microscopic calculations based on chiral effective field theory interactions constrain the properties of neutron-rich matter below nuclear densities to a much higher degree than is reflected in commonly used equations of state. Combined with observed neutron star masses, our results lead to a radius R = 9.7 - 13.9 km for a 1.4 M_{solar} star, where the theoretical range is due, in about equal amounts, to uncertainties in many-body forces and to the extrapolation to high densities.
The Auxiliary Field Diffusion Monte Carlo method has been applied to simulate droplets of 7 and 8 neutrons. Results for realistic nucleon-nucleon interactions, which include tensor, spin--orbit and three--body forces, plus a standard one--body confining potential, have been compared with analogous calculations obtained with Greens Function Monte Carlo methods. We have studied the dependence of the binding energy, the one--body density and the spin--orbit splittings of $^7n$ on the depth of the confining potential. The results obtained show an overall agreement between the two quantum Monte Carlo methods, although there persist differences in the evaluation of spin--orbit forces, as previously indicated by bulk neutron matter calculations. Energy density functional models, largely used in astrophysical applications, seem to provide results significantly different from those of quantum simulations. Given its scaling behavior in the number of nucleons, the Auxiliary Field Diffusion Monte Carlo method seems to be one of the best candidate to perform {sl ab initio} calculations on neutron rich nuclei.
We compute the $S$-factor of the proton-proton ($pp$) fusion reaction using chiral effective field theory ($chi$EFT) up to next-to-next-to-leading order (NNLO) and perform a rigorous uncertainty analysis of the results. We quantify the uncertainties due to (i) the computational method used to compute the $pp$ cross section in momentum space, (ii) the statistical uncertainties in the low-energy coupling constants of $chi$EFT, (iii) the systematic uncertainty due to the $chi$EFT cutoff, and (iv) systematic variations in the database used to calibrate the nucleon-nucleon interaction. We also examine the robustness of the polynomial extrapolation procedure, which is commonly used to extract the threshold $S$-factor and its energy-derivatives. By performing a statistical analysis of the polynomial fit of the energy-dependent $S$-factor at several different energy intervals, we eliminate a systematic uncertainty that can arise from the choice of the fit interval in our calculations. In addition, we explore the statistical correlations between the $S$-factor and few-nucleon observables such as the binding energies and point-proton radii of $^{2,3}$H and $^3$He as well as the $D$-state probability and quadrupole moment of $^2$H, and the $beta$-decay of $^{3}$H. We find that, with the state-of-the-art optimization of the nuclear Hamiltonian, the statistical uncertainty in the threshold $S$-factor cannot be reduced beyond 0.7%.
The density and temperature dependence of the nuclear symmetry free energy is investigated using microscopic two- and three-body nuclear potentials constructed from chiral effective field theory. The nuclear force models and many-body methods are benchmarked to properties of isospin-symmetric nuclear matter in the vicinity of the saturation density as well as the virial expansion of the neutron matter equation of state at low fugacities. The free energy per particle of isospin-asymmetric nuclear matter is calculated assuming a quadratic dependence of the interaction contributions on the isospin asymmetry. The spinodal instability at subnuclear densities is examined in detail.
We outline how auxiliary-field quantum Monte Carlo (AFQMC) can leverage graphical processing units (GPUs) to accelerate the simulation of solid state sytems. By exploiting conservation of crystal momentum in the one- and two-electron integrals we show how to efficiently formulate the algorithm to best utilize current GPU architectures. We provide a detailed description of different optimization strategies and profile our implementation relative to standard approaches, demonstrating a factor of 40 speed up over a CPU implementation. With this increase in computational power we demonstrate the ability of AFQMC to systematically converge solid state calculations with respect to basis set and system size by computing the cohesive energy of Carbon in the diamond structure to within 0.02 eV of the experimental result.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
