Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ENERGY DEPENDENCE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗

The Role of Strong Gravity and the Nuclear Equation of State on Neutron-star Common-envelope Accretion

Common-envelope evolution is important in the formation of neutron star binaries within the isolated binary formation channel. As a neutron star inspirals within the envelope of a primary massive star, it accretes and spins up. Because neutron stars are in the strong-gravity regime, they have a substantial relativistic mass deficit, i.e., their gravitational mass is less than their baryonic mass. This effect causes some fraction of the accreted baryonic mass to convert into neutron star binding energy. The relativistic mass deficit also depends on the nuclear equation of state, since more compact neutron stars will have larger binding energies. We model the mass growth and spin-up of neutron stars inspiraling within common-envelope environments and quantify how different initial binary conditions and hadronic equations of state affect the post-common-envelope neutron star's mass and spin. From these models, we find that neutron star mass growth is suppressed by ≈15%–30%. We also find that for a given amount of accreted baryonic mass, more compact neutron stars will spin-up faster while gaining less gravitational mass, and vice versa. Here, this work demonstrates that a neutron star's strong gravity and nuclear microphysics plays a role in neutron-star-common-envelope evolution, in addition to the macroscopic astrophysics of the envelope. Strong gravity and the nuclear equation of state may thus affect both the population properties of neutron star binaries and the cosmic double neutron star merger rate.

79 ASTRONOMY AND ASTROPHYSICS↗

Energy Resilience for Mission Assurance: Case Study Scoping Document

The Energy Resilience for Mission Assurance (ERMA) project—a Department of Energy Grid Modernization Lab Consortium effort carried out via a partnership among five national laboratories— seeks to develop metrics to quantify how improvements to energy system resilience translate to improved Department of Defense (DoD) mission assurance (MA) during wide-scale, long-duration outages of the bulk power system. DoD missions are integral to national security and highly dependent on electric power. However, energy system planners—both civilian and military—lack a clear and quantifiable mapping between electric power system resilience and MA, leaving a gap in their ability to understand and consider national security outcomes within their planning efforts. The ERMA project seeks to fill this gap, providing stakeholders with new capabilities to understand the impact of electric power system resilience on MA during hazard scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Ultra-Fast High-Precision Metallic Nanoparticle Synthesis using Laser-Accelerated Protons

Laser-driven proton acceleration, as produced during the interaction of a high-intensity (I > 1 × 10 18 W/cm 2 ), short pulse (<1 ps) laser with a solid target, is a prosperous field of endeavor for manifold applications in different domains, including astrophysics, biomedicine and materials science. These emerging applications benefit from the unique features of the laser-accelerated particles such as short duration, intense flux and energy versatility, which allow obtaining unprecedented temperature and pressure conditions. In this paper, we show that laser-driven protons are perfectly suited for producing, in a single sub-ns laser pulse, metallic nanocrystals with tunable diameter ranging from tens to hundreds of nm and very high precision. Our method relies on the intense and very quick proton energy deposition, which induces in a bulk material an explosive boiling and produces nanocrystals that aggregate in a plasma plume composed by atoms detached from the proton-irradiated surface. The properties of the obtained particles depend on the deposited proton energy and on the duration of the thermodynamical process. Suitably controlling the irradiated dose allows fabricating nanocrystals of a specific size with low polydispersity that can easily be isolated in order to obtain a monodisperse nanocrystal solution. Molecular Dynamics simulations confirm our experimental results.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗

Modeling of fission gas diffusion and release for Gd 2 O 3 doped UO 2

Uranium dioxide (UO 2 ) is the primary nuclear fuel in light water reactors, and its excess neutronic reactivity can be controlled by adding burnable absorbers, such as Gd 2 O 3 . This burnable absorber has a large neutron absorption cross-section, lowering the high reactivity of the reactor's initial fuel load. However, there needs to be more understanding of how added Gd 2 O 3 influences the properties of UO 2 under irradiation. To understand the behavior of defects and fission gas in the UO 2 /Gd 2 O 3 system under irradiation, we use cluster dynamics modeling supported by density functional theory calculations. First, we calculate the formation energies of Gd point and cluster defects, and evaluate the temperature-dependent defect concentrations using the defect formation energies and entropies. We show that Gd is soluble in UO 2 , introducing a negative charge in the system. Using this information, we adapted the cluster dynamics code Centipede to model the influence of Gd on U self-diffusion and Xe diffusion in UO 2 with 10 wt% Gd 2 O 3 . Also, we analyzed the Xe diffusion as a function of Gd 2 O 3 concentration, showing that the Xe diffusivity is decreased, which means that the athermal diffusivity due to electronic stopping persists at higher temperatures. In conclusion, the decrease in Xe diffusion means that more Xe stays in the matrix, decreasing the Xe release, and lowering its influence of fission gas release on the thermomechanical properties of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Weak Dispersion of Exciton Landé Factor with Band Gap Energy in Lead Halide Perovskites: Approximate Compensation of the Electron and Hole Dependences

The optical properties of lead halide perovskite semiconductors in vicinity of the bandgap are controlled by excitons, so that investigation of their fundamental properties is of critical importance. The exciton Landé or g-factor g X is the key parameter, determining the exciton Zeeman spin splitting in magnetic fields. The exciton, electron, and hole carrier g-factors provide information on the band structure, including its anisotropy, and the parameters contributing to the electron and hole effective masses. Here, g X is measured by reflectivity in magnetic fields up to 60 T for lead halide perovskite crystals. The materials band gap energies at a liquid helium temperature vary widely across the visible spectral range from 1.520 up to 3.213 eV in hybrid organic–inorganic and fully inorganic perovskites with different cations and halogens: FA 0.9 Cs 0.1 PbI 2.8 Br 0.2 , MAPbI 3 , FAPbBr 3 , CsPbBr 3 , and MAPb(Br 0.05 Cl 0.95 ) 3 . The exciton g-factors are found to be nearly constant, ranging from +2.3 to +2.7. Thus, the strong dependences of the electron and hole g-factors on the bandgap roughly compensate each other when combining to the exciton g-factor. The same is true for the anisotropies of the carrier g-factors, resulting in a nearly isotropic exciton g-factor. The experimental data are compared favorably with model calculation results.

36 MATERIALS SCIENCE↗

Enhancing Site Screening for Underground Hydrogen Storage: Qualitative Site Quality Assessment - SHASTA: Subsurface Hydrogen Assessment, Storage, and Technology Acceleration Project

The global shift towards renewable energy sources to mitigate fossil fuel dependency and meet carbon emission targets by 2050 has underscored the importance of innovative solutions to address energy supply-demand imbalances. Underground Hydrogen Storage (UHS) has emerged as a promising strategy to store excess renewable energy in subsurface formations for future retrieval and utilization. This report focuses on enhancing the site screening process for UHS facilities. By drawing on insights from prior research to identify key criteria influencing site suitability, our objective is to present a comprehensive set of sixteen specific criteria essential for refining the selection of UHS sites. These criteria cover various aspects such as reservoir performance, legal access, regulatory compliance, economic viability, public acceptance, and safety and security considerations. Our proposed methodology allows users to evaluate potential sites using binary responses, transitioning from isolated factors to broader considerations that facilitate a qualitative assessment of site suitability. This method enables comparative analysis and informed decision-making, supporting stakeholders in the site selection process. Our approach is designed as a guiding framework rather than a rigid template for site screening, emphasizing the importance of customizing approaches to individual circumstances. As UHS site development project evolves, adopting a comprehensive selection methodology that integrates holistic evaluations will be essential for ensuring efficient and effective site screening processes tailored to specific project needs.

08 HYDROGEN↗

Enhancing Site Screening for Underground Hydrogen Storage: Qualitative Site Quality Assessment - SHASTA: Subsurface Hydrogen Assessment, Storage, and Technology Acceleration Project

The global shift towards renewable energy sources to mitigate fossil fuel dependency and meet carbon emission targets by 2050 has underscored the importance of innovative solutions to address energy supply-demand imbalances. Underground Hydrogen Storage (UHS) has emerged as a promising strategy to store excess renewable energy in subsurface formations for future retrieval and utilization. This report focuses on enhancing the site screening process for UHS facilities. By drawing on insights from prior research to identify key criteria influencing site suitability, our objective is to present a comprehensive set of sixteen specific criteria essential for refining the selection of UHS sites. These criteria cover various aspects such as reservoir performance, legal access, regulatory compliance, economic viability, public acceptance, and safety and security considerations. Our proposed methodology allows users to evaluate potential sites using binary responses, transitioning from isolated factors to broader considerations that facilitate a qualitative assessment of site suitability. This method enables comparative analysis and informed decision-making, supporting stakeholders in the site selection process. Our approach is designed as a guiding framework rather than a rigid template for site screening, emphasizing the importance of customizing approaches to individual circumstances. As UHS site development project evolves, adopting a comprehensive selection methodology that integrates holistic evaluations will be essential for ensuring efficient and effective site screening processes tailored to specific project needs.

08 HYDROGEN↗

Impact of Seasonal Snow-Cover Change on the Observed and Simulated State of the Atmospheric Boundary Layer in a High-Altitude Mountain Valley

The structure and evolution of the atmospheric boundary layer (ABL) under clear-sky fair weather conditions over mountainous terrain is dominated by the diurnal cycle of the surface energy balance and thus strongly depends on surface snow cover. We use data from three passive ground-based infrared spectrometers deployed in the East River Valley in Colorado's Rocky Mountains to investigate the response of the thermal ABL structure to changes in surface energy balance during the seasonal transition from low to high snow cover. Temperature profiles were retrieved from the infrared radiances using the optimal estimation physical retrieval Tropospheric Remotely Observed Profiling via Optimal Estimation. A nocturnal surface inversion formed in the valley during clear-sky days, which was subsequently mixed out during daytime with the development of a convective boundary layer when snow cover was low. Over high snow cover, a very shallow convective boundary layer formed, above which the inversion persisted through the daytime hours. We compare these observations to NOAA's operational High-Resolution-Rapid-Refresh model and find large warm biases on clear-sky days resulting from the model's inability to form strong nocturnal inversions and to maintain the stable stratification in the valley during daytime when there was snow on the ground. We suggest several factors contributing to the large model errors. These are (a) the inability of the model to represent well-developed thermally driven flows likely due to the too coarse horizontal grid spacing (3 km), (b) too much convective mixing during daytime, and (c) too strong vertical coupling between the valley atmosphere and the free troposphere.

54 ENVIRONMENTAL SCIENCES↗

Formation energy puzzle in intermetallic alloys: Random phase approximation fails to predict accurate formation energies

We performed density functional calculations to estimate the formation energies of intermetallic alloys. We used two semilocal approximations, the generalized gradient approximation (GGA) by Perdew-Burke- Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) meta-GGA. In addition, we utilized two nonlocal DFT functionals, the hybrid HSE06, and the state-of-the-art random phase approximation (RPA). The nonlocal functionals such as HSE06 and RPA yield accurate formation energies of binary alloys with completely-filled d-band metals, where semilocal functionals underperform. The accuracy at the nonlocal functionals is greatly reduced when a partially-filled d-band metal is present in an alloy, while PBE-GGA outperforms in these cases. We show that the accurate prediction of formation energies by any DFT method depends on its ability to predict the accurate electronic properties, e.g., valence d-band contribution to the density of states (DOS). The SCAN meta-GGA often corrects the PBE-DOS, however, it does not provide accurate formation energies compared to PBE. This is assumed to be due to the lack of proper error cancellation that should be expected due to the similar bulk nature of both alloys and their constituents, which may improve with the modification of meta-GGA ingredients. RPA yields too negative formation energies of alloys with partially-filled d-band metals. RPA results can be corrected by restoring the exchange-correlation kernel, thereby improving the short-range electron-electron correlation in metallic densities.

36 MATERIALS SCIENCE↗

Climate Change Mitigation Potential of Wind Energy

Global wind resources greatly exceed current electricity demand and the levelized cost of energy from wind turbines has shown precipitous declines. Accordingly, the installed capacity of wind turbines grew at an annualized rate of about 14% during the last two decades and wind turbines now provide ~6–7% of the global electricity supply. This renewable electricity generation source is thus already playing a role in reducing greenhouse gas emissions from the energy sector. Here we document trends within the industry, examine projections of future installed capacity increases and compute the associated climate change mitigation potential at the global and regional levels. Key countries (the USA, UK and China) and regions (e.g., EU27) have developed ambitious plans to expand wind energy penetration as core aspects of their net-zero emissions strategies. The projected climate change mitigation from wind energy by 2100 ranges from 0.3–0.8 °C depending on the precise socio-economic pathway and wind energy expansion scenario followed. The rapid expansion of annual increments to wind energy installed capacity by approximately two times current rates can greatly delay the passing of the 2 °C warming threshold relative to pre-industrial levels. To achieve the required expansion of this cost-effective, low-carbon energy source, there is a need for electrification of the energy system and for expansion of manufacturing and installation capacity.

17 WIND ENERGY↗

Performance of multifunctional piezoelectric energy harvesting microgyroscopes with material degradation

Most microelectromechanical systems (MEMS) are made from brittle materials that are prone to cracking. Damage from the formation of cracks can occur anytime during the manufacturing process or due to mechanical loading of MEMS devices. Cracks can increase the structures flexibility and lead to changes in performance or total failure. Understanding the potential influences cracks in MEMS is necessary for the longevity of these systems. The energy harvesting microgyroscope studied in this work is susceptible to crack propagation in its silicon substrate and lead zirconate titanate piezoelectric patch. To this end, an electromechanical model is developed to capture the impacts of cracks on the performance of a damaged energy harvesting microgyroscope. Numerical methods are used to approximate the damage associated with distributed crack networks. Then, the static and dynamic responses of the damaged microgyroscope are analyzed for two possible multifunctional setups. Here, the current study revealed that cracks can significantly affect the static and dynamic behaviors of the microgyroscope energy harvester. Combinations of cracking and piezoelectric property degradation cause a decrease in the static pull-in voltage and shift the resonance frequencies of the system. Additionally, it is shown that the systems sensitivity to changes in the base rotation and energy harvesting capabilities vary widely depending on the severity of material degradation.

36 MATERIALS SCIENCE↗

Site heterogeneity and broad surface-binding isotherms in modern catalysis: Building intuition beyond the Sabatier principle

Learning the science of heterogeneous catalysis and electrocatalysis always starts with the simple case of a flat, uniform surface with an ideal adsorbate. It has of course been recognized for a century that real catalysts are more complicated. For the increasingly complex catalysts of the 21st century, this Perspective argues that surface heterogeneity and non-ideal binding isotherms are central features, and their implications need to be incorporated in current thinking. A variety of systems are described herein where catalyst complexity leads to broad, non-Langmuirian surface isotherms for the binding of hydrogen atoms – and this occurs even for ideal, flat Pt(111) surfaces. Modern catalysis employs nanoscale materials whose surfaces have substantial step, edge, corner, impurity, and other defect sites, and they increasingly have both metallic and non-metallic elements M n X m , including metal oxides, chalcogenides, pnictides, carbides, doped carbons, etc. The surfaces of such catalysts are often not crystal facets of the bulk phase underneath, and they typically have a variety of potential active sites. Catalytic surfaces in operando are often non-stoichiometric, amorphous, dynamic, and impure, and often vary from one part of the surface to another. Understanding of the issues that arise at such nanoscale, multi-element catalysts is just beginning to emerge. Yet these catalysts are widely discussed using Brønsted/Bell-Evans-Polanyi (BEP) relations, volcano plots, Tafel slopes, the Butler-Volmer equation, and other linear free energy relations (LFERs), which all depend on the implicit assumption that the active sites are “similar” and that surface adsorption is close to ideal. These assumptions underly the ubiquitous intuition based on the Sabatier Principle, that the fastest catalysis will occur when key intermediates have free energies of adsorption that are not too strong nor too weak. Current catalysis research often aims to minimize the complexity of non-ideal isotherms through experimental and computational design (e.g., the use of single crystal surfaces), and these studies are the foundation of the field. In contrast, this Perspective argues that the heterogeneity of binding sites and binding energies is an inherent strength of these catalysts. Here, this diversity makes many nanoscale catalysts inherently a high-throughput screen wrapped in a tiny package. Only by making the heterogeneity part of the foundation of catalysis models, sorting the types of active sites and dissecting non-ideal binding isotherms, will modern catalysis learn to harness the inherent diversity of real catalysts. Controlling and exploiting diversity rather than avoiding it will help to optimize complex modern catalysts and catalytic conditions.

Mayer, James M.↗

Measurements and computational analysis of the natural decay of 176 Lu

Background: Mainly because of its long half-life and despite its scientific relevance, spectroscopic measurements of 176 Lu forbidden β decays are very limited and lack formulation of shape factors. A direct precise measurement of its Q value is also presently unreported. In addition, the description of forbidden decays provides interesting challenges for nuclear theory. The comparison of precise experimental results with theoretical calculations for these decays can help to test underlying models and can aid the interpretation of data from other experiments. Purpose: Perform the first precision measurements of 176 Lu β-decay spectra and attempt the observation of its electron capture decays, as well as perform the first precision direct measurement of the 176 Lu β-decay Q value. Compare the shape of the precisely determined experimental β spectra to theoretical calculations, and compare the end point energy to that obtained from an independent Q value measurement. Method: Here, the 176 Lu β-decay spectra measurements and the search for electron capture decays were performed with an experimental setup that employed lutetium-containing scintillator crystals and a NaI(Tl) spectrometer for coincidence counting. The β decay Q value was determined via high-precision Penning trap mass spectrometry (PTMS) with the LEBIT facility at the National Superconducting Cyclotron Laboratory. The β-spectrum calculations were performed within the Fermi theory formalism with nuclear structure effects calculated using a shell model approach. Results: Both β transitions of 176 Lu were experimentally observed and corresponding shape factors formulated in their entire energy ranges. The search for electron capture decay branches led to an experimental upper limit of 6.3×10 –6 relative to its β decays. The 176 Lu β-decay and electron capture Q values were measured using PTMS to be 1193.0(6) and 108.9(8) keV, respectively. This enabled precise β end point energies of 596.2(6) and 195.3(6) keV to be determined for the primary and secondary β decays, respectively. The conserved vector current hypothesis was applied to calculate the relativistic vector matrix elements. The β-spectrum shape was shown to significantly depend on the Coulomb displacement energy and on the value of the axial vector coupling constant gA, which was extracted according to different assumptions. Conclusion: The implemented self-scintillation method has provided unmatched observations of 176 Lu, independently validated by the first direct measurements of its β-decay Q value by Penning trap mass spectrometry. Theoretical study of the main β transition led to the extraction of very different effective gA and log10f values, showing that a high-precision description of this transition would require a realistic nuclear structure with nucleus deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Future accelerator projects: new physics at the energy frontier

High-energy colliders provide direct access to the energy frontier, allowing to search for new physics at scales as high as the machine’s center-of-mass energy, perform precision measurements of the Standard Model (SM) parameters, including those related to the flavor sector, and determine the Higgs boson properties and their connection to electroweak symmetry breaking. Each proposed future collider option has its own specific science goals and capabilities, depending on the designed running energy (energies) amongst other parameters. In this paper, an overview of the discovery potential of future circular and linear colliders is presented. Results from searches for beyond the Standard Model (BSM) phenomena at proton–proton, proton–electron, electron–positron, and muon–antimuon colliders are summarized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dark Energy Survey Year 3 results: Constraints on extensions to Λ CDM with weak lensing and galaxy clustering

We constrain six possible extensions to the Λ cold dark matter (CDM) model using measurements from the Dark Energy Survey’s first three years of observations, alone and in combination with external cosmological probes. The DES data are the two-point correlation functions of weak gravitational lensing, galaxy clustering, and their cross-correlation. We use simulated data vectors and blind analyses of real data to validate the robustness of our results to astrophysical and modeling systematic errors. In many cases, constraining power is limited by the absence of theoretical predictions beyond the linear regime that are reliable at our required precision. The Λ CDM extensions are dark energy with a time-dependent equation of state, nonzero spatial curvature, additional relativistic degrees of freedom, sterile neutrinos with eV-scale mass, modifications of gravitational physics, and a binned σ 8 ( z ) model which serves as a phenomenological probe of structure growth. For the time-varying dark energy equation of state evaluated at the pivot redshift we find ( w p , w a ) = ( - 0.9 9 - 0.17 + 0.28 , - 0.9 ± 1.2 ) at 68% confidence with z p = 0.24 from the DES measurements alone, and ( w p , w a ) = ( - 1.0 3 - 0.03 + 0.04 , - 0. 4 - 0.3 + 0.4 ) with z p = 0.21 for the combination of all data considered. Curvature constraints of Ω k = 0.0009 ± 0.0017 and effective relativistic species N eff = 3.1 0 - 0.16 + 0.15 are dominated by external data, though adding DES information to external low-redshift probes tightens the Ω k constraints that can be made without cosmic microwave background observables by 20%. For massive sterile neutrinos, DES combined with external data improves the upper bound on the mass m eff by a factor of 3 compared to previous analyses, giving 95% limits of ( Δ N eff , m eff ) ≤ ( 0.28 , 0.20 eV ) when using priors matching a comparable Planck analysis. For modified gravity, we constrain changes to the lensing and Poisson equations controlled by functions Σ ( k , z ) = Σ 0 Ω Λ ( z ) / Ω Λ , 0 and μ ( k , z ) = μ 0 Ω Λ ( z ) / Ω Λ , 0 , respectively, to Σ 0 = 0. 6 - 0.5 + 0.4 from DES alone and ( Σ 0 , μ 0 ) = ( 0.04 ± 0.05 , 0.0 8 - 0.19 + 0.21 ) for the combination of all data, both at 68% confidence. Overall, we find no significant evidence for physics beyond Λ CDM .

79 ASTRONOMY AND ASTROPHYSICS↗

Global Sensitivity Analysis of Bulk Properties of an Atomic Nucleus

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)]. With this method, we compute the binding energy and charge radius of 16 O at more than 106 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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗