Search NASA⌕ Search

SEARCH · Search NASA

Results for “R-matrix Analysis”

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.

49 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

How Bayesian methods can improve R -matrix analyses of data: The example of the d t reaction

The 3 H(d, n) 4 He reaction is of significant interest in nuclear astrophysics and nuclear applications. It is an important, early step in big-bang nucleosynthesis and a key process in nuclear fusion reactors. We use one- and two-level R-matrix approximations to analyze data on the cross section for this reaction at center-of-mass energies below 215 keV. We critically examine the data sets using a Bayesian statistical model that allows for both common-mode and additional point-to-point un- certainties. We use Markov Chain Monte Carlo sampling to evaluate this R-matrix-plus-statistical model and find two-level R-matrix results that are stable with respect to variations in the channel radii. The S factor at 40 keV evaluates to 25.36(19) MeV b (68% credibility interval). We discuss our Bayesian analysis in detail and provide guidance for future applications of Bayesian methods to R-matrix analyses. We also discuss possible paths to further reduction of the S-factor uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Near-Threshold Rovibrational Excitation of HF by Electron Impact

The pronounced threshold peaks found in cross sections for electron-impact excitation of the vibrational levels of HF have been studied using an R-matrix method, in which electron exchange, correlation effects and polarization are included. Vibration is modelled inside the R-matrix sphere using the energy-modified adiabatic phase matrix method (EMAP). Rotation of the target molecule is treated outside the R-matrix sphere, using the frame transformation method. The external rovibronic close-coupling equations are solved by analytic methods. Analysis of the multichannel threshold structures indicates that the peaks found for rovibrational excitation of nu = 1 can be interpreted as virtual states while those for nu = 2 correspond to nuclear-excited Feshbach resonances. These two different results correspond simply to a change of sign of a small parameter that is determined by the threshold analysis. Whereas the overall shapes of the excitation peaks agree with earlier theoretical results, we have found detailed rotational fine structures for scattering energies on the order of the rotational constant above and below the opening of the vibrational channels. The partial and total cross sections for rovibrational state-to-state transitions in the energy region of the peaks are analysed. Characteristic cusp structures and excitation peaks are found at the rovibrational thresholds for transitions with small angular momentum transfer (Delta)j =0, 1 and 2. Differential cross sections are similar to recent measurements. These experimental data do not show the relatively sharp rise at the vibrational thresholds indicated by theory. Theory and experiment are in good agreement for the shape and magnitude of the angular dependence of cross sections for rovibrational state-to-state transitions, especially at larger energies.

Thummel, Helmar T.↗
DOE OSTI2021

Two-photon above-threshold ionization of helium

Multiphoton ionization provides a clear window into the nature of electron correlations in the helium atom. In the present study, the final-state energy range extends up to the region near the N = 2 and N = 3 ionization thresholds, where two-photon ionization proceeds via continuum intermediate states above the lowest threshold. Our calculations are performed using multichannel quantum defect theory (MQDT) and the streamlined R-matrix method. The sum and integration over all intermediate states in the two-photon ionization amplitude is evaluated using the inhomogeneous R-matrix method developed by Robicheaux and Gao. The seamless connection of that method with MQDT allows us to present high-resolution spectra of the final-state Rydberg resonances. Our analysis classifies the resonances above the N = 2 threshold in terms of their group theory quantum numbers. Furthermore, their dominant decay channels are found to obey the previously conjectured propensity rule far more weakly for these even-parity states than was observed for the odd-parity states relevant to single-photon ionization.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Mapping Rydberg States of H 2 with the Halfium R-Matrix Method

In this article, we use the Halfium R-matrix method to investigate the Rydberg states of the H 2 molecule up to n = 20, filling the gap above the low-lying bound states already calculated with configuration interaction packages. Moreover, we show that the use of Quantum Defect Theory scaling laws, allows for a comprehensive analysis of the regular patterns resulting from the coupling between Rydberg series and doubly excited states. The results should open the door for more efficient quasi-diabatization of the potential energy curves which is required for calculating cross sections and rate coefficients of the (e + H 2 + ) collisional processes, involved in the plasma modeling for fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

R -Matrix Analysis and Statistical Properties of Dysprosium Isotopes in the Neutron Energy Ranges Up To A Few Kev

In support of the Nuclear Criticality Safety Program, a set of evaluated resonance parameters was generated for seven dysprosium isotopes in the neutron energy range from thermal up to a few keV. The evaluation methodology used the Reich-Moore approximation to fit, with the R-matrix code SAMMY, the high-resolution capture and transmission measurements on natural and enriched samples recently performed at the Rensselaer Polytechnic Institute Gaerttner LINear ACcelerator facility. Additional transmission data measured on enriched samples by Liou at the Columbia University Nevis synchrocyclotron in the mid-seventies were used to gauge the neutron widths above 15 eV. Thermal constants such as absorption and (in)coherent scattering cross sections and corresponding scattering lengths were calibrated to the National Institute of Standards and Technology’s compilation except for 161,164 Dy isotopes.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Probabilistic neural networks for improved analyses with phenomenological R -matrix

Here we present a method for measurement analyses based on probabilistic deep neural networks that provide several advantages over conventional analyses with phenomenological models. These include predicting physical quantities directly from data, the rapid generation of statistically robust uncertainties, and the ability to bypass some parameters that may induce ambiguities and complications in data analysis. As deep learning methods make predictions through “black boxes,” the uncertainty quantification is typically challenging. We use a probabilistic framework that provides thorough uncertainty quantification and is straightforward to follow in practice. With the network architecture based on the Transformer, we demonstrate the current method for predicting nuclear resonance parameters from scattering data using the phenomenological R-matrix model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Lowest-energy broad 𝛼-cluster resonances in 19 F

There is strong astrophysical interest in the structure of 19 F states near the α-decay threshold, as their properties are important for understanding the development of α clustering in the 20 Ne region. The emergence of clustered states, and more generally of states that couple strongly to reaction channels near their decay thresholds, is also a topic of current interest in theoretical nuclear physics. In this work, we determine the parameters of broad low-spin states in 19 F near the α-decay threshold and present a theoretical study of these states. The analysis is restricted to ℓ = 0 and 1 resonances in the α + 15 N system close to the α threshold in 19 F. Excitation functions for 15 N(α,α) elastic scattering were measured using the Thick Target Inverse Kinematics (TTIK) method, and these new data, together with older high–energy-resolution measurements, were analyzed within the R-matrix framework. The nuclear structure of 19 F was calculated using configuration-interaction methods with a recently developed effective interaction Hamiltonian. As a result, we identify a sequence of α-cluster resonances in 19 F and map the distribution of clustering strength, which is relevant for astrophysical reaction modeling and for the theoretical understanding of many-body dynamics and the emergence of clustering in loosely bound or unstable nuclei. Furthermore, the work advances theoretical insight into the origins of clustering and highlights open questions for future theoretical and experimental studies.

19F↗
DOE OSTI2026

Connecting ground-state properties of 6 Li to each other and to scattering data

We examine the relationship between the asymptotic normalization coefficient (ANC) of 6 Li and other low-energy observables in the α–deuteron system. Our analysis uses a set of calculations carried out within the ab initio no core shell model with continuum (NCSMC) using a variety of inter-nucleon interactions and basis sizes, and yielding 6 Li deuteron separation energies between 1.3 and 1.8 MeV (Hebborn et al 2022 Phys. Rev. Lett. 129 042503). These NCSMC calculations show that the square of the ANC is strongly correlated with the separation energy over this range. In this work, we investigate the origin of this correlation using the phenomenological R-matrix, a single-channel potential and a perturbative approach. We show that this correlation occurs because the depth of the α–deuteron central potential changes by only a small relative amount as the separation energy varies. We then investigate if the ANC can be accurately extracted from α–deuteron phase shifts in an ideal case in which low-energy data are available and there are no experimental errors. We find that both R-matrix and Coulomb-modified effective-range theory (CM-ERE) yield extracted ANCs close to, although not exactly equal to, the NCSMC value, provided the extrapolation is constrained by the known position of the bound-state pole and at least three terms are included in the fit function. The R-matrix approach converges faster than the CM-ERE as the number of parameters increases and is also more robust against the inclusion of low-energy and high-energy phase shift data. Finally, our study also shows that a naive quantification of uncertainties by comparing different truncations used in both theories is not accurate, and suggests the accuracy of ANCs extracted from phase shift data needs further investigation.

R-matrix↗
DOE OSTI2022

Modernization efforts for the R -Matrix code SAMMY [Abstract]

The R-Matrix code SAMMY is a widely used nuclear data evaluation code focused on the resolved range, which includes corrections for experimental effects. The code is still mostly written in Fortran 77, and uses a memory management system suitable for the time of its initial writing (1984). A modernization effort is under way to bring the code in-line with modern software development practices. A continuous-integration testing framework was added, automating the large existing set of test cases. It is run on every commit. The memory management was updated to current standard practices suitable for modern software analysis tools. The code can be obtained from https://code.ornl.gov/RNSD/SAMMY. The resonance parameters and covariance information are now stored in C++ objects shared by SAMMY and AMPX, the processing code that generates nuclear data libraries for SCALE. This allows for easier maintenance and access to the resonance parameters inside and outside of SAMMY. This feature is already used by accessing and changing parameters in memory in the Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments project, Further plans include the switch to the ENDF reading and writing routines in AMPX, as these routines are more robust, easier to maintain, and support more features. Of note here is support for the new GNDS format. Previously it wasn’t easy to share the full covariance matrix for evaluations containing more than one isotope due to limitations on the ENDF format; this is now supported in GNDS. The data are currently available in a binary SAMMY format and can be exported to GNDS to make them more widely available and sharable. The next step will be to use the same resonance processing code at 0K in AMPX and SAMMY as one of the available Reich-Moore R-Matrix formalism. The first step toward this goal is to isolate the reconstruction into a module that takes resonance parameters as its input and does not depend on SAMMY global parameters. This goal has been achieved and it should now be possible to more easily change the resonance formalism and add enhancements as the Phenomenological R-Matrix parameterization of direct, doorway, and compound nuclear reactions discussed elsewhere on this conference. This concerted modernization and enhancement effort provides multiple advantages to the nuclear data community. It will allow parameter optimization using enhanced formalisms, including experimental effects, that better match complex experimental data. Then those evaluated parameters can immediately be passed off to AMPX to be reconstructed with the exact same cross section model and be put into a data library for subsequent testing using SCALE and the Valid Benchmark suite or other suitable benchmark suites.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Consistency of $^{16}O(n,α)$ cross sections

The evaluated $^{16}O(n,α)$ cross sections in the ENDF/B-VIII.0 nuclear data library remain uncertain because of systematic discrepancies in the measured data. In the energy region below the first excited state, R-matrix analyses rely heavily on the inverse reaction, and in particular, they rely on the measured $^{13}C(α,n)^{16}O$ cross sections reported by Bair et al. in 1973 and Harissopulos et al. in 2005. The Harissopulos cross section values are systematically lower than those previously reported by Bair. Here, drawing on the available experimental information, this paper briefly describes and demonstrates that the two sets of measured cross sections are consistent.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Trojan horse method as an indirect approach to study resonant reactions in nuclear astrophysics

The primary goal of the Trojan horse method (THM) is to analyze resonant rearrangement reactions when the density of the resonance levels is low and statistical models cannot be applied. The main difficulty of the analysis is related with the facts that in the final state the THM reaction involves three particles and that the intermediate particle, which is transferred from the Trojan horse particle to the target nucleus to form a resonance state, is virtual. Another difficulty is associated with the Coulomb interaction between the particles, especially, taking into account that the goal of the THM is to study resonant rearrangement reactions at very low energies important for nuclear astrophysics. The exact theory of such reactions with three charged particles is very complicated and is not available. This is why different approximations are used to analyze THM reactions. In this review paper we describe a new approach based on a few-body formalism that provides a solid basis for deriving the THM reaction amplitude taking into account rescattering of the particles in the initial, intermediate and final states of the THM reaction. Since the THM uses a two-step reaction in which the first step is the transfer reaction populating a resonance state, we address the theory of the transfer reactions. The theory is based on the surface-integral approach and R-matrix formalism. We also discuss application of the THM to resonant reactions populating both resonances located on the second energy sheet and subthreshold resonances, which are subthreshold bound states located at negative energies close to thresholds. We consider the application of the THM to determine the astrophysical factors of resonant radiative-capture reactions at energies so low that direct measurements can hardly be performed due to the negligibly small penetrability factor in the entry channel of the reaction. Here, we elucidated the main ideas of the THM and outline necessary conditions to perform the THM experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Investigation of secondary γ -ray angular distributions using the N 15 ( p , α 1 γ ) C * 12 reaction

The observation of secondary γ-rays provides an alternative method of measuring cross sections that populate excited final states in nuclear reactions. The angular distributions of these γ-rays also provide information on the underlying reaction mechanism. Despite a large amount of data of this type in the literature, publicly available R-matrix codes do not have the ability to calculate these types of angular distributions. In this report, the mathematical formalism derived in Brune and deBoer is implemented in the R-matrix code AZURE2 and calculations are compared with previous data from the literature for the 15 N( p, α1γ ) 12 C* reaction. In addition, new measurements, made at the University of Notre Dame Nuclear Science Laboratory using the Hybrid Array of Gamma Ray Detectors (HAGRiD), are reported that span an energy range from E p = 0.88 to 4.0 MeV. Excellent agreement between data and the phenomenological fit is obtained up to the limit of the previous fit at E p = 2.0 MeV and the R-matrix fit is extended from E x ≈ 13.5 up to E x ≈ 15.3 MeV, where 15 N+p and 12 C+α reactions are fit simultaneously for the first time. An excellent reproduction of the 15 N( p, α1γ ) 12 C* and 12 C(α, α) 12 C data is achieved, but inconsistencies and difficulty in fitting other data is encountered and discussed.

6 ≤ A ≤ 19↗