Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RANDOM PHASE APPROXIMATION”

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 127 records · Page 7

β-decay feeding intensity distributions of 71,73 Ni

This paper presents the β-decay feeding intensity distribution and Gamow-Teller transition strength distribution of 71,73 Ni. These quantities were measured using the technique of total absorption spectroscopy at the National Superconducting Cyclotron Laboratory with the Summing NaI(Tl) detector. These measurements provide sensitive constraints to theoretical models used to predict β-decay properties far from stability for astrophysical applications. Specifically, for the astrophysical r process, the majority of the involved nuclei are not accessible by current facilities, and the nuclear input is mainly provided by theory. Furthermore, the present work reports on two neutron-rich nickel isotopes in the region where the weak r process is expected to be relevant in stellar nucleosynthesis. The experimental results are compared to two theoretical models, namely the shell model and the quasiparticle random-phase approximation, to help further refine theoretical calculations and aid in future r-process studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling quasielastic interactions of monoenergetic kaon decay-at-rest neutrinos

Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigation of the isoscalar response of 24 Mg to 6 Li scattering

24 Mg is a strongly deformed nucleus in the ground state. Deformation effects can be observed in the structure of the isoscalar giant monopole and quadrupole resonances. 24 Mg is also a nucleus that is well known to present different types of cluster-oscillation modes. Both giant resonances and cluster states are strongly populated by isoscalar transitions. To extract the E0, E1, and E2 transition strengths via 6 Li scattering. The 6 Li probe is a powerful tool for investigating the isoscalar nuclear response with a very favorable ratio of resonance-to-continuum background. Double-differential cross sections of 6 Li inelastic scattering, at the beam energy of 100 MeV/u, were measured in the excitation-energy range 10 – 40 MeV and scattering angles 0 - 3°. A multipole-decomposition analysis was performed for extracting the isoscalar E0 , E1, and E2 strength distributions. Results: The extracted multipole strengths were compared with predictions from consistent quasiparticle random phase approximation calculations. The theoretical predictions are in fair agreement with the experimental data. The E0 strength was also compared with results from antisymmetrized molecular dynamics calculations found in the literature. A few peaks in the experimental data might be associated with clustering in 24 Mg. Ground-state deformation effects were observed in the isoscalar giant monopole resonance (ISGMR) and isoscalar giant quadrupole resonance (ISGQR) distributions. The ISGMR strength is split in two peaks around 19 and 28 MeV. The ISGQR exhibits a pronounced peak at 20 MeV with a broadening at the low-energy region, similar to predictions from microscopic calculations. Signatures of excitation of cluster states were observed in the E0 response. Further studies including particle-decay measurements will be required to confirm the nature of the observed peaks.

20 ≤ A ≤ 38↗

Noniterative finite amplitude methods for E 1 and M 1 giant resonances

The finite amplitude method (FAM) is a very efficient approach for solving the fully self-consistent random-phase approximation (RPA) equations. We use FAM to rederive the RPA matrices for Skyrme functionals, calculate the electric dipole ( E 1) and the magnetic dipole ( M 1) giant resonances, and compare the results with available experimental and evaluated data. For the E 1 transitions in heavy nuclei, the calculations reproduce well the resonance energy of the photoabsorption cross sections. In the case of M 1 transitions, we show that the residual interaction does not affect the transition strength of double-magic nuclei, which suggests that the spin terms in the Skyrme force currently neglected in the present computation could improve the agreement between FAM and experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗

Extended Interacting Boson Model with Configuration Mixing Description For 98 Mo

An extended version of the interacting boson model (IBM), encompassing two-particle–two-hole (2p-2h) excitations originating from below a valence shell and one-particle–one-hole (1p-1h) excitations extending to the upper next-nearest major shell, is developed to investigate cross-shell excitation effects in the transitional nucleus 98 Mo. This extended IBM framework is utilized to reproduce experimentally observed properties such as low-lying positive-parity level energies, 𝐵⁡(𝐸⁢2) values, and the electric quadrupole moment of the 2$^{+}_{1}$ state, 𝜌 2 ⁡(𝐸⁢0, 0$^{+}_{2}$ → 0$^{+}_{g}$), as well as the strength functions of the isoscalar giant monopole resonance (ISGMR) and the isoscalar giant quadrupole resonance (ISGQR) in 98 Mo. It is shown that the extended consistent-𝑄 IBM, incorporating configuration mixing of both 2p-2h excitations from below the valence shell and 1p-1h excitations to the upper next-nearest major shell, describes the low-lying structural properties and the ISGMR and ISGQR strength distributions of 98 Mo quite well, with the ISGMR and the ISGQR strength distributions comparable to the results of the self-consistent quasiparticle random-phase approximation with Skyrme interactions.

90 ≤ A ≤ 149↗

Low-energy enhancement in the magnetic dipole γ-ray strength functions of heavy nuclei

A low-energy enhancement (LEE), which was observed experimentally in the gamma-ray strength function (γSF) describing the decay of compound nuclei, would have profound effects on r-process nucleosynthesis if it persists in heavy neutron-rich nuclei. The LEE was shown to be a feature of the magnetic dipole (M1) strength function in configuration-interaction shell-model calculations in medium-mass nuclei. However, its existence in heavy open-shell nuclei remains an open question. Here, using a combination of many-body methods, we identify a LEE in the M1 γSFs of heavy samarium nuclei. In particular, we use the static-path plus random-phase approximation (SPA+RPA), which includes static and small-amplitude quantal fluctuations beyond the mean field. Using the SPA+RPA strength as a prior, we apply the maximum-entropy method (MEM) to obtain finite-temperature M1 γSFs from exact imaginary-time response functions calculated with the shell model Monte Carlo (SMMC) method. We find that the slope of the LEE in samarium isotopes is roughly independent of the average initial energy over a wide range below the neutron separation energy. As the neutron number increases, strength transfers to a low-energy excitation, which we interpret as the scissors mode built on top of excited states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electric dipole polarizability of 58 Ni

The electric dipole strength distribution in 58 Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate of the experimentally inaccessible high-energy contribution above 20 MeV, leading to an electric dipole polarizability α D ⁡ ( 58 Ni) = 3.48 (31)⁢ fm 3 . This serves as a test case for recent extensions of coupled-cluster calculations with chiral effective field theory interactions to nuclei with two nucleons on top of a closed-shell system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB + QRPA calculations with the Gogny interaction

The zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of-the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98 Zr . This region includes a benchmarking case and offers insights into nuclear deformation phenomena. Here, to investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the quasiparticle random phase approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98 Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the 𝐾 𝜋 = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the 𝐾 𝜋 = 0 − and 𝐾 𝜋 = 1 − components. Intrinsic transition densities and associated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e., we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Total absorption spectroscopy for the 𝛽 + decay strength distribution of 60 Ga

𝛽-decay properties play an important role in most astrophysical processes. In the absence of experimental data, astrophysical models rely on global theoretical calculations to provide the relevant properties. It is therefore important to provide strong experimental constraints when possible. In the case of 𝛽-decay, the most sensitive probe is the 𝛽-decay strength distribution. We report here on the first measurement of the latter quantity for the 𝛽 + decay of 60 Ga using the total absorption spectroscopy technique. The experimental results are compared to theoretical calculations often used in astrophysical models, namely the shell model and the quasiparticle random phase approximation (QRPA), as well as an extension of QRPA that includes higher-order nucleonic calculations. Both models are in reasonable agreement with the experimental results.

59 ≤ A ≤ 89↗

Statistical model of the stimulated forward Brillouin scattering driven by a randomized laser beam in plasma

The modeling of a spatially incoherent laser beam remains a central problem of the parametric instabilities in the context of inertial confinement fusion. This letter gives a simplified and comprehensive overview of the recent analytical developments regarding the modeling of these laser beams and a comparison with a dedicated experiment. Our model accounts for the first time for the statistical standard deviation of the gain and accurately captures the entanglement between wave mixing processes and the speckle correlations thus resolving the longstanding contradictions between the random phase approximation and the model of independent speckles. It is successfully compared to a recent laser beam spray experiment and the associated paraxial simulations, demonstrating that backscattering predictions require accounting for the beam spray. Furthermore, our framework thus provides a way to evaluate and guide the analysis of parametric instabilities in high laser energy experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Timescales of Quantum Equilibration, Dissipation and Fluctuation in Nuclear Collisions

Understanding the dynamics of equilibration processes in quantum systems as well as their interplay with dissipation and fluctuation is a major challenge in quantum many-body theory. The timescales of such processes are investigated in collisions of atomic nuclei using fully microscopic approaches. Findings from time-dependent Hartree-Fock and time-dependent random-phase approximation calculations are compared for 13 systems over a broad range of energies. The timescale for full mass equilibration ( ~ 2 × 10 - 20 s ) is found to be much larger than timescales for neutron-to-proton equilibration, kinetic energy, and angular momentum dissipations which are on the order of 10 - 21 s . Fluctuations of mass numbers in the fragments and correlations between their neutron and proton numbers build up within only a few 10 - 21 s . This suggests that dissipation is basically not impacted by mass equilibration, but is mostly driven by the exchange of nucleons between the fragments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coexistence and Interaction of Spinons and Magnons in an Antiferromagnet with Alternating Antiferromagnetic and Ferromagnetic Quantum Spin Chains

In conventional quasi-one-dimensional antiferromagnets with quantum spins, magnetic excitations are carried by either magnons or spinons in different energy regimes: they do not coexist independently, nor could they interact with each other. In this Letter, by combining inelastic neutron scattering, quantum Monte Carlo simulations and Random Phase Approximation calculations, we report the discovery and discuss the physics of the coexistence of magnons and spinons and their interactions in Botallackite-Cu 2 (OH) 3 Br. This is a unique quantum antiferromagnet consisting of alternating ferromagnetic and antiferromagnetic Spin-1/2 chains with weak interchain couplings. Furthermore, our study presents a new paradigm where one can study the interaction between two different types of magnetic quasiparticles: magnons and spinons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Crossing N = 28 Toward the Neutron Drip Line: First Measurement of Half-Lives at FRIB

Here, new half-lives for exotic isotopes approaching the neutron drip-line in the vicinity of N~28 for Z=12–15 were measured at the Facility for Rare Isotope Beams (FRIB) with the FRIB decay station initiator. The first experimental results are compared to the latest quasiparticle random phase approximation and shell-model calculations. Overall, the measured half-lives are consistent with the available theoretical descriptions and suggest a well-developed region of deformation below 48 Ca in the N=28 isotones. The erosion of the Z=14 subshell closure in Si is experimentally confirmed at N=28, and a reduction in the 38 Mg half-life is observed as compared with its isotopic neighbors, which does not seem to be predicted well based on the decay energy and deformation trends. This highlights the need for both additional data in this very exotic region, and for more advanced theoretical efforts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗