Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPERIMENTAL DATA”

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 595 records · Page 33

Comparison of spatial dynamics and point kinetics approaches in multiphysics modeling of the molten salt reactor experiment

In this work, we present validation test results of fully coupled neutronics and thermal-hydraulics models of the Molten Salt Reactor Experiment (MSRE) against experimental data of the zero power pump transients and the natural circulation tests at low power. To capture the strong coupling between neutronics and thermal-hydraulics due to fuel circulation, and to account for the delayed neutron precursor (DNP) distribution, the porous media thermal-hydraulics solver Pronghorn was fully coupled to the spatial neutron dynamics code Griffin, which solves the neutron diffusion equation, and to the 0-D point kinetics solver Squirrel, using a 2-D homogenized representation of the MSRE. The validation test results show very good agreement with experimental data for both point kinetics and spatial dynamics simulations, capturing the strong feedback effect and DNP losses in the MSRE. The 0-D code Squirrel accurately predicted the time-dependent behavior in the MSRE given the steady-state spatial dynamics solution of Griffin.

42 - ENGINEERING↗

AQuaRef: machine learning accelerated quantum refinement of protein structures

Cryo-EM and X-ray crystallography provide crucial experimental data for obtaining atomic-detail models of biomacromolecules. Refining these models relies on library-based stereochemical data, which, in addition to being limited to known chemical entities, do not include meaningful noncovalent interactions. Quantum mechanical (QM) calculations could alleviate these issues but are too expensive for large molecules. Here we present a novel AI-enabled Quantum Refinement (AQuaRef) based on AIMNet2 machine learned interatomic potential (MLIP) mimicking QM at substantially lower computational costs. By refining 41 cryo-EM and 30 X-ray structures, we show that this approach yields atomic models with superior geometric quality compared to standard techniques, while maintaining an equal or better fit to experimental data. Notably, AQuaRef aids in determining proton positions, as illustrated in the challenging case of short hydrogen bonds in the parkinsonism-associated human protein DJ-1 and its bacterial homolog YajL.

Zubatyuk, Roman [Carnegie Mellon University, Pitts↗

Universal correlation for falling film evaporation on a horizontal plain tube

Falling film evaporators are widely used in many industrial sectors, such as desalination, refrigeration & air-conditioning, and chemical processes. The heat transfer correlations suitable for wider operating conditions and various refrigerants are useful for designing falling film evaporators. Although many researchers developed the correlation for their own experimental data, there is a limited universal correlation for falling film evaporators. In addition, the dry-out effects were not considered in the available universal correlations. Here, therefore, this study proposes a new correlation for falling film evaporation on a single horizontal plain tube, which includes nucleate boiling, film evaporation, and dry-out effects. The proposed correlation covers a wider range of operating conditions (Heat flux is of 2.5–168kWm -2 , and film Reynolds number is up to 3500), various refrigerants (R134a, R245fa, R1234ze(E), R123, R32, R410A, R290, and R600a), and the large number of experimental data (1101 data in total). The proposed correlation is validated against the available standard correlation with regard to mean absolute deviation (MAD) and R 2 values. The MAD of the new correlation is about 10%, respectively, which is the lowest among other correlations. Similarly, the R 2 showed a value of 0.92, which is the highest compared to the existing correlations. A simplified correlation is also proposed.

42 ENGINEERING↗

Developing new analysis tools for near surface radio-based neutrino detectors

The ARIANNA experiment is an Askaryan radio detector designed to measure high-energy neutrino induced cascades within the Antarctic ice. Ultra-high-energy neutrinos above 10 16 eV have an extremely low flux, so experimental data captured at trigger level need to be classified correctly to retain as much neutrino signal as possible. We first describe two new physics-based neutrino selection methods, or "cuts", (the updown and dipole cut) that extend the previously published analysis to a specialized ARIANNA station with 8 antenna channels, which is double the number used in the prior analysis. For a standard trigger with a threshold signal to noise ratio at 4.4, the new cuts produce a neutrino efficiency of > 95% per station-year of operation, while rejecting 99.93% of the background (corresponding to 53 remaining experimental background events). When the new cuts are combined with a previously developed cut using neutrino waveform templates, all background is removed at no change of efficiency. In addition, the neutrino efficiency is extrapolated to 1,000 station-years of operation, obtaining 91%. This work then introduces a new selection method (the deep learning cut) to augment the identification of neutrino events by using deep learning methods and compares the efficiency to the physics-based analysis. The deep learning cut gives 99% signal efficiency per station-year of operation while rejecting 99.997% of the background (corresponding to 2 remaining experimental background events), which are subsequently removed by the waveform template cut at no significant change in efficiency. Here the results of the deep learning cut were verified using measured cosmic rays which shows that the simulations do not introduce artifacts with respect to experimental data. The paper demonstrates that the background rejection and signal efficiency of near surface antennas meets the requirements of a large scale future array, as considered in baseline design of the radio component of IceCube-Gen2.

97 MATHEMATICS AND COMPUTING↗

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY↗

Samarium: from a distorted-fcc phase to melting under dynamic compression using in-situ x-ray diffraction

Abstract Lattice and electronic structure interactions for f -electrons are fundamental challenges for lanthanide equation of state development. Difficulties in first-principles calculations, such as density functional theory (DFT), emphasize the need for well-characterized experimental data. Here, we measure in-situ x-ray diffraction of shocked samarium (Sm) and temperature along the Hugoniot for the first time, providing direct evidence for phase transitions. We report direct evidence of a distorted fcc (dfcc) phase at 23 GPa. Shocked samarium melts from the dfcc phase starting at 33 GPa (1333 K), with complete melt at 40 GPa (1468 K). Previous work indicated shock melt at 27 GPa (1200 K), underscoring the significance of x-ray measurements for detecting phase transitions. Interestingly, our observed melting is in sharp contrast with the melting reported by a diamond anvil cell study. These experimental data can tightly constrain first principles calculations and serve as key touchstones for equation of state modeling.

36 MATERIALS SCIENCE↗

Experimental constraints on the properties of $P_c$ states

Using new experimental data on photoproduction and $\Lambda _b^0$ decays, we derive constraints on the properties of the LHCb $P_c$ states. We conclude that $P_c(4312)$, $P_c(4380)$ and $P_c(4440)$ can be described as $\Sigma _c\bar{D}$, $\Sigma _c^*\bar{D}$ and $\Sigma _c\bar{D}^*$ molecules, but that $P_c(4457)$ does not fit into the same picture. Based on the apparent absence of additional partner states, and the striking disparity between $\Lambda _c^+\bar{D}^0$ and $\Lambda _c^+\bar{D}^{*0}$ decays, we conclude that $P_c(4440)$ has $3/2^-$ quantum numbers. Using heavy-quark symmetry we predict large experimental signals for $P_c$ states in $\eta _c p K^-$, $\Lambda _c^+\bar{D}^{*0}K^-$, and $\Sigma _c^{(*)}\bar{D}K^-$. We also argue that current experimental data on photoproduction is almost at the level of sensitivity required to observe $P_c$ states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydrodynamic description of D meson production in high-energy heavy-ion collisions *

The large values and constituent-quark-number scaling of the elliptic flow of low-p T D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-p T charm quarks, we compare the D 0 meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. D 0 mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D 0 mesons at P T < 3 GeV/c as measured in Au+Au collisions at $sqrts_{NN} = 200$ GeV. Though the experimental data on spectra are consistent with the hydrodynamic result at small p T ~ 1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum, GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Finally, our comparative study indicates that charm quarks only approach local thermal equilibrium at small p T , even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate p T .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Experimental study on kinetic oxidation of graphite IG-110 by steam

Graphite is proposed for use in High-temperature Gas-cooled Reactors (HTGRs) as the fuel matrix, neutron moderator/reflector, and core structural material. One important property of nuclear grade graphite is their resistance to oxidation in high-temperature environment. Extensive investigation has been performed in the literature for graphite oxidation by air. However, available experimental data are still limited for graphite oxidation by steam under conditions comparable to a postulated steam ingress accident in HTGRs. In this study, the oxidation rate of graphite IG-110 by steam was measured at temperatures from 850 to 1100 °C with the steam partial pressure varying from 0.5 to 20.0 kPa and the hydrogen partial pressure varying from 0 to 2.0 kPa. Further analysis confirms the oxidation process in this present study is dominated by the chemical kinetics, which lends credit to the data for being used to develop numerical models. It was observed that the increase of the kinetic oxidation rate with the steam partial pressure tends to become less apparent if the steam partial pressure keeps increasing. In addition, it was found that the partitioning of hydrogen inhibits the graphite-steam reaction process even with the steam partial pressure up to 20.0 kPa. However, this inhibiting effect starts to become saturated when the hydrogen partial pressure exceeds 1.0 kPa. The oxidation rates were fitted to the conventional Langmuir-Hinshelwood (LH) and Boltzmann-enhanced Langmuir-Hinshelwood (BLH) models by a multivariable optimization algorithm. The BLH model exhibits a better accuracy than the LH model within the specified experimental conditions. The predicted oxidation rate using the BLH model shows a mean relative difference of about 24% with the maximum difference of about 55% when compared with our experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Optimizing training trajectories in variational autoencoders via latent Bayesian optimization approach *

Unsupervised and semi-supervised ML methods such as variational autoencoders (VAE) have become widely adopted across multiple areas of physics, chemistry, and materials sciences due to their capability in disentangling representations and ability to find latent manifolds for classification and/or regression of complex experimental data. Like other ML problems, VAEs require hyperparameter tuning, e.g. balancing the Kullback–Leibler and reconstruction terms. However, the training process and resulting manifold topology and connectivity depend not only on hyperparameters, but also their evolution during training. Because of the inefficiency of exhaustive search in a high-dimensional hyperparameter space for the expensive-to-train models, here we have explored a latent Bayesian optimization (zBO) approach for the hyperparameter trajectory optimization for the unsupervised and semi-supervised ML and demonstrated for joint-VAE with rotational invariances. We have demonstrated an application of this method for finding joint discrete and continuous rotationally invariant representations for modified national institute of standards and technology database (MNIST) and experimental data of a plasmonic nanoparticles material system. The performance of the proposed approach has been discussed extensively, where it allows for any high dimensional hyperparameter trajectory optimization of other ML models.

42 ENGINEERING↗

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗

Simulation studies of the interaction of laser radiation with additively manufactured foams

The interaction of laser radiation with foams of various porosities and low densities has been the subject of several numerical and experimental studies. In all cases, the modeling of low-Z under-dense foams as uniform gases of equivalent average density using standard radiation-hydrodynamics codes has resulted in heat-front velocities that are considerably faster than those observed experimentally. It has been theoretically conjectured that this difference may be attributed to the breakdown of the foam’s morphology, leading to a dynamics of filament expansion where the ion and electron energy partitions are significantly different from those calculated using the uniform gas model. We found that 3D computer simulations employing a disconnected representation of the foam’s microstructure which allowed for the dynamics of foam element heating, expansion, and stagnation largely supported the theoretical picture. Simulations using this model for laser experiments on under-dense 2 mg cc -1 SiO 2 aerogel foams reproduced the experimental data fairly well. We used the validated model in simulations of low-density structured foam-like materials (produced via additive manufacturing) with a variety of morphologies. We found that the log-pile configurations were consistent with the analytical propagation model of Gus’kov et al. Further validation of the model was obtained by simulating experiments performed at the Jupiter Laser Facility using the log-pile and octet-truss foam morphologies. Simulations of the foam-laser interaction using a wave propagation code showed that the microstructure was able to enhance stimulated Brillouin scattering (SBS) by concentrating the light energy into density holes. In turn, this promotes laser filamentation, reducing SBS and bringing the predicted values closer to the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization and Scaling of Forced Convective Swirl in Sinusoidal Wavy-Plate-Fin Cores of Compact Heat Exchangers

The characterization and scaling of the thermal-hydraulic performance in wavy plate-fin compact heat exchanger cores, based on the understanding of the physical phenomena and heat transfer enhancement mechanism is delineated. Experimental data are presented for forced convection in air (Pr = 0.71) with flow rates in the range 50 ≤ Re ≤ 4000. A variety of wavy-fin cores that span viable applications, with geometrical attributes described by the cross section aspect ratio α (ratio of fin spacing to height), fin corrugation aspect ratio γ (ratio of 2× corrugation amplitude to wave pitch), and fin spacing ratio ζ (ratio of fin spacing to wave pitch), are considered. To characterize and correlate the vortex-flow mixing in interfin spaces, a Swirl number Sw is introduced from the balance of viscous, inertial and centrifugal forces. It is shown that the laminar, transitional and turbulent flow regimes can be explicitly identified by this Swirl number. Based on the experimental results and extended analysis, new correlations for Fanning friction factor f and Colburn factor j are developed with Sw, α, γ, and ζ as scaling parameters. Requisite expressions are devised by a superposition of both enhancement components due to the corrugated surface area enlargement and induced swirl flow field, and they are combined to cover the laminar, transitional and turbulent regimes by the method of asymptotic matching. The resulting generalized correlations are further shown to predict all available experimental data for f and j factors to within ±20% and ±15%, respectively.

42 ENGINEERING↗

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Energy dependent calculations of fission product, prompt, and delayed neutron yields for neutron induced fission on 235 U, 238 U, and 239 Pu

In this article, we perform energy-dependent calculations of independent and cumulative fission product yields for 235 U, 238 U, and 239 Pu in the first chance fission region. Starting with the primary fission fragment distributions taken from available experimental data and analytical functions based on assumptions for the excitation energy and spin-parity distributions, the Hauser-Feshbach statistical decay treatment for fission fragment de-excitation is applied to more than 1,000 fission fragments for the incident neutron energies up to 5 MeV. The calculated independent yields are then used as an input of β-decay calculations to produce the cumulative yield, and summation calculations are performed. Model parameters in these procedures are adjusted by applying the Bayesian technique at the thermal energy for 235 U and 239 Pu and in the fast energy range for 238 U. The calculated fission observable quantities, such as the energy-dependent cumulative yields, and prompt and delayed neutron yields, are compared with available experimental data. We also study the possible impact of the second chance fission opening on the energy dependence of the delayed neutron yield by on the opening of second chance fission on the energy dependence calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

40 Ca transverse response function from coupled-cluster theory

Here, we present calculations of the 40 Ca transverse response function obtained from coupled-cluster theory used in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-next-to leading order in chiral effective field theory with and without Δ degrees of freedom. We first benchmark this approach on the 4 He nucleus and compare both the transverse sum rule and the response function to earlier calculations based on different methods. As expected from the power counting of the chiral expansion of electromagnetic currents and from previous studies, our results retaining only one-body term underestimate the experimental data for 4 He by about 20%. However, when the method is applied to 40 Ca at the same order of the expansion, response functions do not lack strength and agree well with the world electron scattering data. We discuss various sources of theoretical uncertainties and comment on the comparison of our results with the available experiments.

39 ≤ A ≤ 58↗