Search NASA⌕ Search

SEARCH · Search NASA

Results for “R MATRIX”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Parameterization of Direct and Doorway Processes in R-Matrix Formalism

R-matrix formalism is extended beyond compound nuclear (CN) resonant reactions to include parameterization of direct as well as doorway processes. Direct processes in the R-matrix exterior are parameterized by a unitary matrix that introduces mixing among wave function coefficients of the incoming and outgoing wave function components at the R-matrix channel surface. Doorway processes are parameterized by separating the Hilbert space of the interior R-matrix region into its doorway and CN subspaces, from which doorway state eigenenergies, reduced width amplitudes, and the strengths of their coupling to CN levels appear as new R-matrix parameters. Parameterization of generalized as well as the conventional Reich–Moore approximation for eliminated capture channels in the presence of direct, doorway, and CN processes is presented along with a complex-valued scattering length with contributions from direct, doorway, and CN capture processes. Derivation of Brune’s alternative R-matrix parameters is extended to include doorway states. This work suggests how R-matrix formalism could be extended further by adopting the concepts from related reaction formalisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

R-matrix School 2025: Introduction to R-matrix Theory

This technical memo serves as a lecture material for the R-matrix school 2025 and distributed among the participants. The manuscript discusses in details the introduction to the R-matrix theory including its algorithm to calculate reaction cross sections and derivation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Scattering matrix pole expansions for complex wave numbers in R -matrix theory

In this followup article to Ducru et al., we establish new results on scattering matrix pole expansions for complex wave numbers in R-matrix theory. In the past, two branches of theoretical formalisms emerged to describe the scattering matrix in nuclear physics: R-matrix theory and pole expansions. The two have been quite isolated from one another. Recently, our study of Brune's alternative parametrization of R-matrix theory has shown the need to extend the scattering matrix (and the underlying R-matrix operators) to complex wave numbers. Two competing ways of doing so have emerged from a historical ambiguity in the definitions of the shift S and penetration P functions: the legacy Lane and Thomas's “force closure” approach versus analytic continuation (which is the standard in mathematical physics). The R-matrix community has not yet come to a consensus as to which to adopt for evaluations in standard nuclear data libraries, such as ENDF. Here, in this article, we argue in favor of analytic continuation of R-matrix operators. We bridge R-matrix theory with the Humblet-Rosenfeld pole expansions, and discover new properties of the Siegert-Humblet radioactive poles and widths, including their invariance properties to changes in channel radii a c . We then show that analytic continuation of R-matrix operators preserves important physical and mathematical properties of the scattering matrix—canceling spurious poles and guaranteeing generalized unitarity—while still being able to close channels below thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Secondary γ-ray decays from the partial-wave $\textit{T}$ matrix with an $\textit{R}$-matrix application to 15 N($p, α_1γ$) 12 C

The secondary γ rays emitted following a nuclear reaction are frequently relatively straightforward to detect experimentally. Despite the large volume of such data, a practical formalism for describing these γ rays in terms of partial-wave $\textit{T}$-matrix elements has never been given. The partial-wave formalism is applicable when $\textit{R}$-matrix methods are used to describe the reaction in question. This paper supplies the needed framework, and it is demonstrated by the application to the 15 N($p, α_1γ$) 12 C reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

K-SHELL PHOTOIONIZATION OF NICKEL IONS USING R-MATRIX

We present R-matrix calculations of photoabsorption and photoionization cross sections across the K edge of the Li-like to Ca-like ion stages of Ni. Level-resolved, Breit-Pauli calculations were performed for the Li-like to Na-like stages. Term-resolved calculations, which include the mass-velocity and Darwin relativistic corrections, were performed for the Mg-like to Ca-like ion stages. This data set is extended up to Fe-like Ni using the distorted wave approximation as implemented by AUTOSTRUCTURE. The R-matrix calculations include the effects of radiative and Auger dampings by means of an optical potential. The damping processes affect the absorption resonances converging to the K thresholds causing them to display symmetric profiles of constant width that smear the otherwise sharp edge at the K-shell photoionization threshold. These data are important for the modeling of features found in photoionized plasmas.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2011

K-Shell Photoionization of Nickel Ions Using R-Matrix

We present R-matrix calculations of photoabsorption and photoionization cross sections across the K edge of the Li-like to Ca-like ions stages of Ni. Level-resolved, Breit-Pauli calculations were performed for the Li-like to Na-like stages. Term-resolved calculations, which include the mass-velocity and Darwin relativistic corrections, were performed for the Mg-like to Ca-like ion stages. This data set is extended up to Fe-like Ni using the distorted wave approximation as implemented by AUTOSTRUCTURE. The R-matrix calculations include the effects of radiative and Auger dampings by means of an optical potential. The damping processes affect the absorption resonances converging to the K thresholds causing them to display symmetric profiles of constant width that smear the otherwise sharp edge at the K-shell photoionization threshold. These data are important for the modeling of features found in photoionized plasmas.

Witthoeft, M. C.↗
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 OSTI2022

Phenomenological R -Matrix parameterization of direct, doorway, and compound nuclear reactions [Abstract]

Although formal expressions for scattering matrix accounting for direct, doorway, and compound nuclear (CN) resonant reactions have been derived several decades ago in both the transition ( T -)matrix formalism and the reactance ( K -)matrix formalism, the absence of corresponding expressions in phenomenological R -matrix formalism has limited the application of the latter to CN resonant reactions only. We remove this limitation by parameterizing direct, doorway, and CN resonant reactions in a phenomenological R -matrix scattering matrix, and provide a parameterization for a corresponding Reich-Moore approximation of eliminated capture channels. Direct reactions induce (previously neglected) mixing among the incoming or outgoing R -matrix channel wave functions, parameterized by real and orthonormal channel-rotation matrix, M , whereby the original scattering matrix U is transformed into M T UM . Any real and orthonormal matrix, M , can be equivalently expressed as e η , where η is a real and skew-symmetric 2 rotation-generating matrix that subsequently yields a more intuitive parameterization of eliminated direct capture reactions in Reich-Moore approximation. A phenomenological R -matrix parameterization of doorway reactions is inferred by equating the expression for reactance ( K -)matrix, given in terms of Brune’s alternative R -matrix parameterization, to a corresponding expression derived using Feshbach’s projection operator formalism. Assuming that all doorway states, just like CN states, are confined within spheres defined by R -matrix channel radii, a new R -matrix-like term induced by doorway states is gleaned, wherein each doorway state is parameterized by its energy, width, and the strength of its coupling to each CN state. Since a Reich-Moore approximation for retained-channel scattering matrix ought to approximate the effect of eliminated capture channels taking place via direct, doorway, or CN reactions, each of the three kinds of reactions contributing to the capture entails a corresponding Reich-Moore parameterization in a first-order approximation: direct contribution is parameterized by introducing finite diagonal elements of a retained-channel rotation-generating matrix, doorway contribution is parameterized by doorway capture widths, while CN contribution is parameterized by conventional Reich-Moore capture widths. We will present evidence of direct and doorway reactions observed in recent measurements of resolved resonance cross sections at the Gaerttner LINAC Center at Rensselaer Polytechnic Institute, and will outline a path for implementing this new R -matrix parameterization into the SAMMY nuclear data evaluation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Level-resolved R-matrix calculations for the electron-impact excitation of Ne{sup 3+} and Ne{sup 6+}

Large-scale R-matrix calculations are carried out for the electron-impact excitation of Ne{sup 3+} and Ne{sup 6+}. For Ne{sup 3+}, a 581-LSJ-level R-matrix intermediate coupling frame transformation calculation is made for excitations up to the n=4 shell. For some transitions, large effective collision strength differences are found with current 23-jKJ-level Breit-Pauli R-matrix and earlier 22-LSJ-level R-matrix jj omega (JAJOM) calculations. For Ne{sup 6+}, a 171-jKJ-level Breit-Pauli R-matrix calculation is made for excitations up to the n=5 shell. For some transitions, large effective collision strength differences are found with current 46-jKJ-level Breit-Pauli R-matrix and earlier 46-LSJ-level R-matrix JAJOM calculations. Together with existing R-matrix calculations for other ion stages, high-quality excitation data are now available for astrophysical and laboratory plasma modeling along the entire Ne isonuclear sequence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Shadow poles in the alternative parametrization of R-matrix theory

In this work, we discover new, hitherto unknown, shadow poles in Brune’s alternative parametrization of R-matrix theory [C. R. Brune, Phys. Rev. C 66, 044611 (2002)]. Where these poles are, and how many, depends on how one continues R-matrix operators to complex wavenumbers (specially the shift S and penetration P functions). This has little consequence for the exact R-matrix formalism (past the last energy threshold), as we show one can still always fully reconstruct the scattering matrix with only the previously known alternative parameters (poles and corresponding resonance widths), for which there were as many poles as the number of levels Nλ. However, we generalize the alternative parametrization to the Reich-Moore formalism, and show that the choice of continuation is now critical as it changes the alternative parameters values (poles and residue widths are now complex). In order to establish nuclear libraries with alternative parameters, the nuclear community will thus have to decide what convention to adopt. We argue in favor of analytical continuation (against the legacy Lane and Thomas approach) in a follow-up article [P. Ducru, Phys. Rev. C, submitted (2020)]. We observe the first evidence of shadow poles in the alternative parametrization of R-matrix theory in isotope xenon 134 Xe spin-parity group J π = 1=2 (–) , and show how they indeed depend on the choice of continuation to complex wavenumbers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Verification of R-matrix calculations for charged-particle reactions in the resolved resonance region for the 7 Be system

R-matrix theory is used to present nuclear reactions in the resolved resonance region. It uses information on bound states and low energy resonances to accurately parametrize cross sections on the resonances as well as the non-resonant background. Since the seminal work of Lane and Thomas (1958), the method has been widely used to analyze experimental cross-section data in a broad range of fields spanning nuclear reaction dynamics, nuclear astrophysics, ion beam analysis and their applications. Different R-matrix codes have been developed and used in these different applications with very little communication among the developers or practitioners on the capabilities, achievements or limitations of the codes. A limited comparison among three R-matrix codes on neutron-induced reactions was performed by the International Atomic Energy Agency (IAEA) International Evaluation of Neutron Cross Section Standards project (2007). Since then, significant progress has been made in their implementation of the R-matrix algorithms, and R-matrix codes have enhanced capabilities. In this paper we introduce, for the first time, the results of a comprehensive effort to verify the most widely used R-matrix codes in the various fields of nuclear science and applications: AMUR, AZURE2, CONRAD, EDA, FRESCO, GECCCOS, and SAMMY. In addition to the description of the capabilities of the codes and their specifications, we discuss the results of a joint exercise which was coordinated by the International Atomic Energy Agency. The aim of the exercise was to compare calculations of charged-particle reaction cross sections for the light composite system 7Be. The calculations were performed by the codes using identical input R-matrix parameters and other specifications and were limited to charged-particle channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The R-matrix bootstrap for the 2d O( N ) bosonic model with a boundary

The S-matrix bootstrap is extended to a 1+1d theory with O(N) symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Covariance of Neutron Cross Sections for {sup 16}O through R-matrix Analysis

Through the R-matrix analysis, neutron cross sections as well as the covariance are estimated for {sup 16}O in the resolved resonance range. Although we consider the current results are still preliminary, we present the summary of the cross section analysis and the results of data uncertainty/covariance, including those for the differential cross sections. It is found that the values obtained highlight consequences of nature in the theory as well as knowledge from measurements, which gives a realistic quantification of evaluated nuclear data covariances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗