We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus ^{16}O to identify the most influential low-energy constants in the next-to-next-to-leading order ...chiral Hamiltonian with two- and three-nucleon forces. For this purpose, we develop a subspace-projected coupled-cluster method using eigenvector continuation Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.032501. With this method, we compute the binding energy and charge radius of ^{16}O at more than 10^{6} different values of the 16 low-energy constants in one hour on a standard laptop computer. For relatively small subspace projections, the root-mean-square error is about 1% compared to full-space coupled-cluster results. We find that 58(1)% of the variance in energy can be apportioned to a single contact term in the ^{3}S_{1} wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segrè chart.
In fish, maximum O2 consumption rate (ṀO2,max) and aerobic scope can be expanded following exhaustive exercise in hyperoxia; however, the mechanisms explaining this are yet to be identified. Here, in ...exhaustively exercised rainbow trout (Oncorhynchus mykiss), we assessed the influence of hyperoxia on ṀO2,max, aerobic scope, cardiac function and blood parameters to address this knowledge gap. Relative to normoxia, ṀO2,max was 33% higher under hyperoxia, and this drove a similar increase in aerobic scope. Cardiac output was significantly elevated under hyperoxia at ṀO2,max because of increased stroke volume, indicating that hyperoxia released a constraint on cardiac contractility apparent with normoxia. Thus, hyperoxia improved maximal cardiac performance, thereby enhancing tissue O2 delivery and allowing a higher ṀO2,max. Venous blood O2 partial pressure (PvO2) was elevated in hyperoxia at ṀO2,max, suggesting a contribution of improved luminal O2 supply in enhanced cardiac contractility. Additionally, despite reduced haemoglobin and higher PvO2, hyperoxia treated fish retained a higher arterio-venous O2 content difference at ṀO2,max. This may have been possible because of hyperoxia offsetting declines in arterial oxygenation that are known to occur following exhaustive exercise in normoxia. If this occurs, increased contractility at ṀO2,max with hyperoxia may also relate to an improved O2 supply to the compact myocardium via the coronary artery. Our findings show ṀO2,max and aerobic scope may be limited in normoxia following exhaustive exercise as a result of constrained maximal cardiac performance and highlight the need to further examine whether or not exhaustive exercise protocols are suitable for eliciting ṀO2,max and estimating aerobic scope in rainbow trout.
We use Bayesian methods and Hamiltonian Monte Carlo (HMC) sampling to infer the posterior probability density function (PDF) for the low-energy constants (LECs) up to next-to-next-to-next-to-leading ...order (N3LO) in a chiral effective field theory (χEFT) description of the nucleon-nucleon interaction. In a first step, we condition the inference on neutron-proton and proton-proton scattering data and account for uncorrelated χEFT truncation errors. We demonstrate how to successfully sample the 31-dimensional space of LECs at N3LO using a revised HMC inference protocol. In a second step we extend the analysis by means of importance sampling and an empirical determination of the neutron-neutron scattering length to infer the posterior PDF for the leading charge-dependent contact LEC in the S01 neutron-neutron interaction channel. While doing so we account for the χEFT truncation error via a conjugate prior. We use the resulting posterior PDF to sample the posterior predictive distributions for the effective range parameters in the S01 wave as well as the strengths of charge-symmetry breaking and charge-independence breaking. We conclude that empirical point-estimate results of isospin breaking in the S01 channel are consistent with the PDFs obtained in our Bayesian analysis and that, when accounting for χEFT truncation errors, one must go to next-to-next-to-leading order to confidently detect isospin breaking effects.
The number of low-energy constants (LECs) in chiral effective field theory (chi EFT) grows rapidly with increasing chiral order, necessitating the use of Markov chain Monte Carlo techniques for ...sampling their posterior probability density function. For this we introduce a Hamiltonian Monte Carlo (HMC) algorithm and sample the LEC posterior up to next-to-next-to-leading order (NNLO) in the two-nucleon sector of chi EFT. We find that the sampling efficiency of HMC is three to six times higher compared to an affine-invariant sampling algorithm. We analyze the empirical coverage probability and validate that the NNLO model yields predictions for two-nucleon scattering data with largely reliable credible intervals, provided that one ignores the leading-order EFT expansion parameter when inferring the variance of the truncation error. We also find that the NNLO truncation error dominates the error budget.
This review presents some of the challenges in constructing models of atomic nuclei starting from theoretical descriptions of the strong interaction between nucleons. The focus is on statistical ...computing and methods for analyzing the link between bulk properties of atomic nuclei, such as radii and binding energies, and the underlying microscopic description of the nuclear interaction. The importance of careful model calibration and uncertainty quantification of theoretical predictions is highlighted.
Coronary perfusion and cardiac autonomic regulation may benefit myocardial oxygen delivery and thermal performance of the teleost heart, and thus influence whole animal heat tolerance. Yet, no study ...has examined how coronary perfusion affects cardiac output during warming in vivo. Moreover, while β-adrenergic stimulation could protect cardiac contractility, and cholinergic decrease in heart rate may enhance myocardial oxygen diffusion at critically high temperatures, previous studies in rainbow trout (
Oncorhynchus mykiss
) using pharmacological antagonists to block cholinergic and β-adrenergic regulation showed contradictory results with regard to cardiac performance and heat tolerance. This could reflect intra-specific differences in the extent to which altered coronary perfusion buffered potential negative effects of the pharmacological blockade. Here, we first tested how cardiac performance and the critical thermal maximum (CT
max
) were affected following a coronary ligation. We then assessed how these performances were influenced by pharmacological cholinergic or β-adrenergic blockade, hypothesising that the effects of the pharmacological treatment would be more pronounced in coronary ligated trout compared to trout with intact coronaries. Coronary blockade reduced CT
max
by 1.5 °C, constrained stroke volume and cardiac output across temperatures, led to earlier cardiac failure and was associated with reduced blood oxygen-carrying capacity. Nonetheless, CT
max
and the temperatures for cardiac failure were not affected by autonomic blockade. Collectively, our data show that coronary perfusion improves heat tolerance and cardiac performance in trout, while evidence for beneficial effects of altered cardiac autonomic tone during warming remains inconclusive.
We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate. Theoretical error estimates of this low-energy process are important for a ...reliable interpretation of forthcoming experimental results by the MuSun Collaboration. Specifically, we estimate the three dominant sources of uncertainties that impact theoretical calculations of this rate: those resulting from uncertainties in the pool of fit data used to constrain the coupling constants in the nuclear interaction, those due to the truncation of the effective field theory, and those due to uncertainties in the axial radius of the nucleon. For the capture rate into the S01 channel, we find an uncertainty of approximately 4.6s-1 due to the truncation in the effective field theory and an uncertainty of 3.9s-1 due to the uncertainty in the axial radius of the nucleon, both of which are similar in size to the targeted experimental precision of a few percent.
We present a Bayesian analysis of renormalization-group invariant nucleon-nucleon interactions at leading order in chiral effective field theory (χEFT) with momentum cutoffs in the range 400–4000 ...MeV. We use history matching to identify relevant regions in the parameter space of low-energy constants (LECs) and subsequently infer the posterior probability density of their values using Markov chain Monte Carlo. All posteriors are conditioned on experimental data for neutron-proton scattering observables and we estimate the χ EFT truncation error in an uncorrelated limit. We do not detect any significant cutoff dependence in the posterior predictive distributions for two-nucleon observables. For all cutoff values we find a multimodal LEC posterior with an insignificant mode harboring a bound 1S0 state. The 3P0 and 3P2 phase shifts emerging from the Bayesian analysis are less constrained and typically more repulsive compared to the results of a phase shift optimization. We expect that our inference will impact predictions for nuclei. This work demonstrates how to perform inference in the presence of limit-cycle-like behavior and spurious bound states, and lays the foundation for a Bayesian analysis of renormalization-group invariant χEFT interactions beyond leading order.
In this work we present a framework that allows one to solve the Faddeev equations for three-nucleon scattering using the wave-packet continuum-discretization method. We perform systematic benchmarks ...using results in the literature and study in detail the convergence of this method with respect to the number of wave packets. We compute several different elastic neutron-deuteron scattering cross-section observables for a variety of energies using chiral nucleon-nucleon interactions. For the optimized next-to-next-to-leading order interaction N2LOopt we find good agreement with data for nucleon scattering-energies ELab≤70 MeV and a slightly larger maximum of the neutron analyzing power Ay(n) at ELab=10 and 21 MeV compared with other interactions. This work represents a first step towards a systematic inclusion of three-nucleon scattering observables in the construction of next-generation nuclear interactions.