Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userUncertainty quantification for proton-proton fusion in chiral effective field theory
263
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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%.
rate research
Read More
The astrophysical S-factor for proton-proton fusion, S_11(E), is obtained with the nuclear matrix element analytically calculated in pionless effective field theory. To the third order, the zero-energy result S_11(0) and the first energy derivative S_11(0) are found to be (3.99 pm 0.14)* 10^-25 MeV b and S_11(0)*(11.3 pm 0.1) MeV^-1, respectively; both consistent with the current adopted values. The second energy derivative is also calculated for the first time, and the result S_11(0) = S_11(0)*(170 pm 2) MeV^-2 contributes at the level of 0.5% to the fusion rate at the solar center, which is smaller than 1% as previously estimated.
The proton-proton fusion reaction, $ppto de^+ u$, is studied in pionless effective field theory (EFT) with di-baryon fields up to next-to leading order. With the aid of the di-baryon fields, the effective range corrections are naturally resummed up to the infinite order and thus the calculation is greatly simplified. Furthermore, the low-energy constant which appears in the axial-current-di-baryon-di-baryon contact vertex is fixed through the ratio of two- and one-body matrix elements which reproduces the tritium lifetime very precisely. As a result we can perform a parameter free calculation for the process. We compare our numerical result with those from the accurate potential model and previous pionless EFT calculations, and find a good agreement within the accuracy better than 1%.
We present a systematic calculation of the cross section for the lepton-proton bremsstrahlung process l + p --> l + p + gamma in chiral perturbation theory at next-to-leading order. This process corresponds to an undetected background signal for the proposed MUSE experiment at PSI. MUSE is designed to measure elastic scattering of low-energy electrons and muons off a proton target in order to extract a precise value of the protons r.m.s. radius. We show that the commonly used peaking approximation, which is used to evaluate the radiative tail for the elastic cross section, is not applicable for muon-proton scattering at the low-energy MUSE kinematics. Furthermore, we point out a certain pathology with the standard chiral power counting scheme associated with electron scattering, whereby the next-to-next-to-leading order contribution from the pion loop diagrams is kinematically enhanced and numerically of the same magnitude as the next-to-leading order corrections. We correct a misprint in a commonly cited review article.
We employ an effective field theory (EFT) that exploits the separation of scales in the p-wave halo nucleus $^8mathrm{B}$ to describe the process $^7mathrm{Be}(p,gamma)^8mathrm{B}$ up to a center-of-mass energy of 500 keV. The calculation, for which we develop the lagrangian and power counting, is carried out up to next-to-leading order (NLO) in the EFT expansion. The power counting we adopt implies that Coulomb interactions must be included to all orders in $alpha_{rm em}$. We do this via EFT Feynman diagrams computed in time-ordered perturbation theory, and so recover existing quantum-mechanical technology such as the two-potential formalism for the treatment of the Coulomb-nuclear interference. Meanwhile the strong interactions and the E1 operator are dealt with via EFT expansions in powers of momenta, with a breakdown scale set by the size of the ${}^7$Be core, $Lambda approx 70$ MeV. Up to NLO the relevant physics in the different channels that enter the radiative capture reaction is encoded in ten different EFT couplings. The result is a model-independent parametrization for the reaction amplitude in the energy regime of interest. To show the connection to previous results we fix the EFT couplings using results from a number of potential model and microscopic calculations in the literature. Each of these models corresponds to a particular point in the space of EFTs. The EFT structure therefore provides a very general way to quantify the model uncertainty in calculations of $^7mathrm{Be}(p,gamma)^8mathrm{B}$. We also demonstrate that the only N$^2$LO corrections in $^7mathrm{Be}(p,gamma)^8mathrm{B}$ come from an inelasticity that is practically of N$^3$LO size in the energy range of interest, and so the truncation error in our calculation is effectively N$^3$LO. We also discuss the relation of our extrapolated $S(0)$ to the previous standard evaluation.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
