R-Matrix Theory from the Green-Bloch Perspective
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
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.
It is often the case that nuclear reactions that are important for societal applications or basic science are difficult to measure directly in existing facilities for accelerated beams. The difficulty might be associated with an exceedingly small cross section, as the ones obtained at beam energies well below the Coulomb barrier, or with the impossibility of devising a short-lived target, as is the case for neutron-induced reactions on unstable isotopes. The fruitful line of experimental research addressing this issue with alternative (indirect) reactions has been developed in parallel to the theory needed to make the connection between the observed data and the desired cross sections of the (direct) reaction under study. Within this context, we present here an explicit connection between the indirect cross section and the R-matrix parameters which describe the energy-dependent cross section of the desired direct reaction.
A resolved resonance region evaluation of 140,142 Ce was conducted by Oak Ridge National Laboratory. Requested by the US Nuclear Criticality Safety Program, this evaluation is based on recent high-resolution transmission and capture high-resolution measurements of nat Ce and 142 Ce conducted at JRC-Geel at the Geel Linear Accelerator facility. It is also based on recently measured thermal constants available from the EX FOR database. Starting from the resonance parameters from the ENDF/B-VIII.0 library and following a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance in formation was derived by the fit of these experimental datasets using the Reich–Moore approximation of the R-matrix theory, as implemented in the SAMMY code system. The resolved resonance region upper energy limit for 140 Ce was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. This new evaluation was found to be in good agreement not only with several integral quantities of interest to the reactor physics community, but also with the stellar Maxwellian-averaged cross section.
The INDEN for Light Elements network (INDEN-LE) held a consultants’ meeting from 20 to 23 June 2022, to review the status of the evaluations undertaken by the network as well as developments in R-matrix theory and new measurements. The summaries of the presentations and discussions can be found in this report.
Comparison is made with results of new calculations in a 25-state R-matrix theory and with results in a previous eight-state R-matrix calculation. The presence of rich resonance structure is confirmed in both experiment and theory.
Background: The preference for light nuclear systems to coagulate into -particle clusters has been well-studied. The possibility of a linear chain configuration of -particles would allow for a new way to study this phenomenon. Purpose: A rotational band of states in C has been claimed showing a linear chain structure. The mirror system, O, has been studied here to examine how this linear chain structure is affected by replacing the valence neutrons with protons. Method: A beam of C was incident into a chamber filled with He:CO gas with the tracks recorded inside the TexAT Time Projection Chamber and the recoil -particles detected by a silicon detector array to measure the cross section. Results: The experimental cross section was compared with previous studies and fit using R-Matrix theory with the previously-observed O states being transformed to the C using mirror symmetry. The measured cross section does not replicate the claimed states, with the predicted cross section exceeding that observed at several energies and angles. Conclusion: A series of possibilities are highlighted with the most likely being that the originally-seen C states did not constitute a rotational band with a potentially incorrect spin assignment due to the limitations of the angular correlation method with non-zero spin particles. The work highlights the difficulties in measuring broad resonances corresponding to a linear chain state in a high level density.
Boundary condition constraint method for generalizing R-matrix theory, permitting use of shell model as basis for nuclear-reaction calculations
This article, on electron-molecule collisions, is dedicated to the legacy of my good friend and sometime collaborator, Aaron Temkin on his retirement from the NASA-Goddard Space Flight Center after many years of work at the highest intellectual level in the theoretical treatment of electron-atom and electron-molecule scattering. Aaron's contributions to the manner in which we think about electron-molecule collisions is clear to all of us who have worked in this field. I doubt that the great progress that has occurred in the computational treatment of such complex collision problems could have happened without these contributions. For a brief historical account, see the discussion of Temkin's contribution to electron-molecule scattering in the first article of this volume by Dr. A. K. Bhatia. In this article, I will concentrate on the application of the so called, non-adiabatic R-matrix theory, to vibrational excitation and dissociative attachment, although I will also present some results applying the Linear Algebraic and Kohn-Variational methods to vibrational excitation. As a starting point for almost all computationally effective approaches to electron-molecule collisions, is the fixed nuclei approximation. That is, one recognizes, just as one does with molecular bound states, that there is a separation of electronic(fast) and nuclear(s1ow) degrees of freedom. This separation makes it possible to "freeze" the nuclei in space, calculate the collision parameters for the frozen molecule and then, somehow to add back the vibrations and rotations. The manner in which this is done, depends on the details of the collision problem. It is the work of Aaron and a number of other researchers that has provided the guidance necessary to resolve these issues.
Explore the source record for details and available documents.
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.
We present recent advances in the R-matrix formalism as well as the Bayesian evaluation framework for improved nuclear data evaluations. The advances in the R matrix formalism include: 1) direct processes, 2) doorway, as well as multistep, processes, and 3) various forms of the Reich-Moore approximation for eliminated capture channels. Furthermore, to address unreasonably small posterior uncertainties often encountered in nuclear data evaluations of large data sets using the conventional form of the Bayes’ theorem, we introduce imperfections (of the data or the model) as a formal evaluation tool for taming the evaluated uncertainties in harmony with Bayes’ theorem. These theoretical advances were motivated by the nuclear data evaluations of differential resolved resonance cross section data using the code SAMMY, as well as the integral benchmark experiments using the SCALE code system, being performed at Oak Ridge National Laboratory for the Nuclear Criticality Safety Program. Some pedagogical applications of the new formalism, as well as a snapshot of the SAMMY modernization efforts, will be presented.
The position operator 𝑟̂ appears as 𝑖∂ 𝑝 in wave mechanics, while its matrix form (e.g., under a Bloch basis) is well known diverging in diagonals, causing difficulties in basis transformation, observable yielding, etc. We aim to find a convergent r-matrix (CRM) to improve the existing divergent r-matrix (DRM), and investigate its influence at both the conceptual and the application levels. A key modification is increasing the familiar substitution of 𝑟̂ by 𝑖∂ 𝑝 to 𝑖∑ 𝑗 ∂ 𝑘 𝑗 , namely the N-th Weyl algebra. Resolving the divergence makes r-matrix rigorously defined, and we are able to show r-matrix is distinct from a spin matrix in terms of its defining principles, transformation behavior, and the observable it yields. Conceptually, the CRM fills the logical gap between the r-matrix and the Berry connection (this unremarked vagueness has caused the diagonal divergence). In application, we focus on transport, and discover that the Hermitian matrix is not identical with the associative Hermitian operator, i.e., 𝑟 𝑚,𝑛 = 𝑟$^∗_{𝑛,𝑚}$ ⇎ 𝑟̂ = 𝑟̂ † , which subtly affects the celebrated Berry curvature formula for adiabatic current. We also discuss how such a non-representation CRM can contribute to building a unified transport theory.
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.
In this work, we investigate α-clustering in 18 Ne and compare it with what is known about clustering in the mirror nucleus 18 O. The excitation function for resonant elastic scattering of 14 O + α was measured in inverse kinematics using the TexAT active-target detector system, covering excitation energies from 8 to 17 MeV. The data were analyzed using a multi-channel R-matrix approach. Detailed spectroscopic information—including excitation energies, spin and parity assignments, and partial α and total widths—was extracted from the R-matrix analysis and compared with theoretical models and previously available data. Correspondence between levels in 18 O and 18 Ne is established. An extensive shell-model analysis of the 18 O– 18 Ne mirror pair shows very good agreement between theory and experiment, especially in the distribution of clustering strength. The comparison reveals that certain states, particularly at higher excitation energies, exhibit significantly stronger α-cluster character than predicted. This suggests that these states have a collective structure aligned with the corresponding α-reaction channel.
The goal of this research grant was to calculate accurate oscillator strengths and electron collisional excitation strengths for inelastic transitions in atomic species of relevance to Planetary Atmospheres. Large scale configuration-interaction atomic structure calculations have been performed to obtain oscillator strengths and transition probabilities for transitions among the fine-structure levels and R-matrix method has been used in the calculations of electron-ion collision cross sections of C II, S I, S II, S III, and Ar II. A number of strong features due to ions of sulfur have been detected in the spectra of Jupiter satellite Io. The electron excitation cross sections for the C II and S II transitions are studied in collaboration with the experimental atomic physics group at the Jet Propulsion Laboratory. There is excellent agreement between experiment and theory which provide an accurate and broad-base test of the ability of theoretical methods used in the calculation of atomic processes. Specifically, research problems have been investigated for: electron impact excitation cross sections of C II: electron impact excitation cross sections of S III; energy levels and oscillator strengths for transitions in S III; collision strengths for electron collisional excitation of S II; electron impact excitation of inelastic transitions in Ar II; oscillator strengths of fine-structure transitions in neutral sulfur; cross sections for inelastic scattering of electrons from atomic nitrogen; and excitation of atomic ions 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.
We investigate the alpha-decay width of a near proton-threshold resonance in 11 B by the excitation functions of 7 Li(α, α) and 7 Li(α, α’) reactions. This resonance is an example of loosely bound atomic nuclei, understood as small open quantum systems, that exhibit properties significantly influenced by their coupling to the continuum. The present experiment focuses on the formation of a specific state in 11 B at excitation energy in the region near 11.4 MeV by bombarding a 7 Li target with a 4 He beam at energies ranging from 3.92 to 4.56 MeV in the laboratory frame. The study provides data complementary to previous observations of the proton emission in the β-decay of the neutron-rich halo nucleus 11 Be. An R-matrix fit to the data provides resonance energy and partial widths consistent with J π =1/2 + for a narrow near-threshold state in 11 B at 11.400(20) MeV, that is 171(20) keV above the proton emission threshold, with partial α-decay width Γ α =2.5$^{+3.5}_{–1.5}$ keV in elastic scattering and Γ α' =15.8$^{+1.9}_{–0.4}$ keV in inelastic scattering. These findings contribute to the understanding of the structure of this open quantum system resonance.