SEARCH · Engineering Papers

Results for “Radiation Scattering Equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Back-to-Back Inclusive Dijets in Deep Inelastic Scattering at Small x : Complete NLO Results and Predictions

We compute the back-to-back dijet cross section in deep inelastic scattering at small x to next-to-leading order (NLO) in the color glass condensate effective field theory. Our result can be factorized into a convolution of the Weizsäcker-Williams gluon transverse-momentum-dependent distribution function (WW gluon TMD) with a universal soft factor and an NLO coefficient function. The soft factor includes both double and single logarithms in the ratio of the relative transverse momentum P ⊥ of the dijet pair to the dijet momentum imbalance q ⊥ ; its renormalization group (RG) evolution is resummed into the Sudakov factor. Likewise, the WW TMD obeys a nonlinear RG equation in x that is kinematically constrained to satisfy both the lifetime and rapidity ordering of the projectile. Exact analytical expressions are obtained for the NLO coefficient function of transversely and longitudinally polarized photons. Our results allow for the first quantitative separation of the dynamics of Sudakov suppression from that of gluon saturation. They can be extended to other final states and provide a framework for precision tests of novel QCD many-body dynamics at the Electron-Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Nucleon off-shell structure and the free neutron valence structure from A=3 inclusive electron scattering measurements

Understanding the differences between the distribution of quarks bound in protons and neutrons is key for constraining the mechanisms of SU(6) spin-flavor symmetry breaking in Quantum Chromodynamics (QCD). While vast amounts of proton structure measurements were done, data on the structure of the neutron is much more spars as experiments typically extract the structure of neutrons from measurements of light atomic nuclei using model-dependent corrections for nuclear effects. Recently the MARATHON collaboration performed such an extraction by measuring inclusive deep-inelastic electron-scattering on helium-3 and tritium mirror nuclei where nuclear effects are expected to be similar and thus be suppressed in the helium-3 to tritium ratio. Here we evaluate the model dependence of this extraction by examining a wide range of models including the effect of using instant-form and light-cone nuclear wave functions and several different parameterizations of nucleon modification effects, including those with and without isospin dependence. We find that, while the data cannot differentiate among the different models of nuclear structure and nucleon modification, they consistently prefer a neutron-to-proton structure function ratio of at $x_B \rightarrow 1$ of $\sim 0.4$ with a typical uncertainty ($1\sigma$) of $\sim0.05$ and $\sim0.10$ for isospin-independent and isospin-dependent modification models, respectively. While strongly favoring SU(6) symmetry breaking models based on perturbative QCD and the Schwinger-Dyson equation calculation, the MARATHON data do not completely rule out the scalar di-quark models if an isospin-dependent modification exist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Quark sivers function at small x: spin-dependent odderon and the sub-eikonal evolution

We apply the formalism developed earlier for studying transverse momentum dependent parton distribution functions (TMDs) at small Bjorken x to construct the small- x asymptotics of the quark Sivers function. First, we explicitly construct the complete fundamental “polarized Wilson line” operator to sub-sub-eikonal order: this object can be used to study a variety of quark TMDs at small x . We then express the quark Sivers function in terms of dipole scattering amplitudes containing various components of the “polarized Wilson line” and show that the dominant (eikonal) term which contributes to the quark Sivers function at small x is the spin-dependent odderon, confirming the re- cent results of Dong, Zheng and Zhou. Our conclusion is also similar to the case of the gluon Sivers function derived by Boer, Echevarria, Mulders and Zhou. We also analyze the sub-eikonal corrections to the quark Sivers function using the constructed “polarized Wilson line” operator. We derive new small- x evolution equations re-summing double-logarithmic powers of α s ln 2 (1 /x ) with α s the strong coupling constant. We solve the corresponding novel evolution equations in the large- N c limit, obtaining a sub-eikonal correction to the spin-dependent odderon contribution. We conclude that the quark Sivers function at small x receives contributions from two terms and is given by ${f}_{1T}^{\perp q}\left(x,{k}_T^2\right)={C}_O\left(x,{k}_T^2\right)\frac{1}{x}+{C}_1\left({k}_T^2\right){\left(\frac{1}{x}\right)}^0+\cdots$ with the function C O ( x, ${k}_T^2$) varying slowly with x and the ellipsis denoting the subasymptotic and sub-sub-eikonal (order- x ) corrections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Electric dipole excitations near the neutron separation energies in 96 Mo

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly 96 Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Probing the core of the strong nuclear interaction

The strong nuclear interaction between nucleons (protons and neutrons) is the effective force that holds the atomic nucleus together. This force stems from fundamental interactions between quarks and gluons (the constituents of nucleons) that are described by the equations of quantum chromodynamics. However, as these equations cannot be solved directly, nuclear interactions are described using simplified models, which are well constrained at typical inter-nucleon distances but not at shorter distances. This limits our ability to describe high-density nuclear matter such as that in the cores of neutron stars. Here we use high-energy electron scattering measurements that isolate nucleon pairs in short-distance, high-momentum configurations, accessing a kinematical regime that has not been previously explored by experiments, corresponding to relative momenta between the pair above 400 mega electron volts per c (c, speed of light in vacuum). As the relative momentum between two nucleons increases and their separation thereby decreases, we observe a transition from a spin-dependent tensor force to a predominantly spin-independent scalar force. These results demonstrate the usefulness of using such measurements to study the nuclear interaction at short distances and also support the use of point-like nucleon models with two- and three-body effective interactions to describe nuclear systems up to densities several times higher than the central density of the nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as (n, γ) cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron (n), a proton (p), and the target nucleus (A) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force (n–p–A 3BF). Here, our goal is to quantify systematic uncertainties stemming from the reduction of deuteron-nucleus (d + A) dynamics to a picture of three pointlike nuclear clusters interacting via pairwise nucleon-nucleus forces, using as testing grounds d + α scattering and the 6 Li ground state. We particularly focus on quantifying uncertainties arising from the full antisymmetrization of the (A + 2)-body system with the target nucleus fixed in its ground state. We adopt the ab initio no-core shell model coupled with the resonating group method (NCSM/RGM) to compute microscopic n–α and p–α interactions, and use them in a three-body description of the d + α system by means of momentum-space Faddeev-type equations. Simultaneously, we also carry out ab initio calculations of d + α scattering and 6 Li ground state by means of six-body NCSM/RGM calculations to serve as a benchmark for the three-body model predictions given by the Faddeev calculations. By comparing the Faddeev and NCSM/RGM results, we show that the irreducible n–p–α 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach yields a 6 Li ground state that is approximately 600 keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for d + α scattering yield a 3 + resonance that is located approximately 400 keV higher in energy compared to the NCSM/RGM result. The shape of the d + α angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the 3 + resonance energy. The Faddeev three-body model predictions for d + α scattering and 6 Li using microscopic n–α and p–α potentials differ from those computed microscopically with the NCSM/RGM. These discrepancies are due to the n–p–α 3BF, which arises from two-nucleon exchange terms in the microscopic d–α interaction and are not accounted for in the three-body model Faddeev calculations. This study lays the foundation for future parametrizations of the 3BF due to Pauli exclusion principle effects in improved three-body calculations of deuteron-induced reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Investigating radiatively driven, magnetized plasmas with a university scale pulsed-power generator

We present first results from a novel experimental platform that is able to access physics relevant to topics including indirect-drive magnetized inertial confinement fusion, laser energy deposition, various topics in atomic physics, and laboratory astrophysics (for example, the penetration of B-fields into high energy density plasmas). This platform uses the x rays from a wire array Z-pinch to irradiate a silicon target, producing an outflow of ablated plasma. The ablated plasma expands into ambient, dynamically significant B-fields (∼5 T), which are supported by the current flowing through the Z-pinch. The outflows have a well-defined (quasi-1D) morphology, enabling the study of fundamental processes typically only available in more complex, integrated schemes. Experiments were fielded on the MAGPIE pulsed-power generator (1.4 MA, 240 ns rise time). On this machine, a wire array Z-pinch produces an x-ray pulse carrying a total energy of ∼15 kJ over ∼30 ns. This equates to an average brightness temperature of around 10 eV on-target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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 OSTI2024

Low-lying dipole response of Ni 64

In this work, two complementary real-photon scattering experiments were conducted on the proton-magic 64 Ni nucleus to study the dipole response up to its neutron-separation energy of $S_n$ = 9.7 MeV. By combining both measurements, 87 $E$1 and 23 $M$1 transitions were identified above 4.3 MeV. The results of the observed $M$1 transitions were compared to shell-model calculations using two different model spaces. It was found that the inclusion of excitations across the $Z$ = 28 shell gap in the calculations has a large impact. Furthermore, average cross sections for decays to the ground state (elastic transitions) as well as to lower-lying excited states (inelastic decays) were determined. The corresponding $E$1 channel was compared to calculations within the relativistic equation of motion (REOM) framework. Whereas the calculations of highest possible complexity reproduce the fragmentation and overall behavior of the $E$1 average elastic cross section well, the predicted absolute cross sections are approximately twice as high as the experimental upper limits even though the latter also include an estimate of the inelastic-decay channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Solution Irregularity Remediation for Spatial Discretization Error Estimation for S N Transport Solutions

The discrete ordinates linear Boltzmann transport equation is typically solved in its spatially discretized form, incurring spatial discretization error. Quantification of this error for purposes such as adaptive mesh refinement or error analysis requires an a posteriori estimator, which utilizes the numerical solution to the spatially discretized equation to compute an estimate. Because the quality of the numerical solution informs the error estimate, irregularities, present in the true solution for any realistic problem configuration, tend to cause the largest deviation in the error estimate vis-a-vis the true error. In this paper, an analytical partial singular characteristic tracking (pSCT) procedure for reducing the estimator’s error is implemented within our novel residual source estimator for a zeroth-order discontinuous Galerkin scheme, at the additional cost of a single inner iteration. Here, a metric-based evaluation of the pSCT scheme versus the standard residual source estimator is performed over the parameter range of a Method of Manufactured Solutions test suite. The pSCT scheme generates near-ideal accuracy in the estimate in problems where the dominant source of the estimator’s error is the solution irregularity, namely, problems where the true solution is discontinuous and problems where the true solution’s first derivative is discontinuous and the scattering ratio is low. In problems where the scattering ratio is high and the true solution is discontinuous in the first derivative, the error in the scattering source, which is not converged by the pSCT scheme, is greater than the error incurred due to the irregularity. Ultimately, a pSCT scheme is judged to be useful for error estimation in problems where the computational cost of the scheme is justified. In the presence of many irregularities, such a scheme may be intractable for general use, but in benchmarks, as an analytical tool, or in problems that have nondissipative discontinuities, the scheme may prove invaluable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Progress Towards the Reality of Low-Temperature Critical Experiments with Low-Temperature TEX

Integral benchmarks for criticality safety do not exist below room temperature. Incredible support for below room temperature benchmarks is present across the international community. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Similarly, the IAEA, which contains 176 member states, has the same requirement for subcriticality for fissile material transportation. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University. A six-case experiment designed based off the TEX-HEU baseline experiments has been proposed utilizing a chilled vacuum chamber capable of cooling the HEU stack down to -50°C. The six proposed configurations span the entire neutron fission energy range and can be tested at intermediate temperatures between room temperature and -40°C. Of the six configurations, the fastest case utilizes nearly 135 kg of HEU. Understanding the system is required prior to critical experiment operations. In order to test the functionality of the system, testing with surrogate fuel materials has been conducted. The following discusses the design of the low temperature TEX (LT-TEX) experiment as well as the progress made via surrogate testing towards a reality in which low-temperature integral benchmarks exist in the International Criticality Safety Benchmark Evaluation Project (ICSBEP).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Benchmarking fast-neutron response of 7 Li-enriched C 7 LYC scintillators: Simulations and experiment

New insights into the fast neutron response of 7 Li-enriched C 7 LYC scintillators have been obtained through detailed Geant4 simulations of 35 Cl(n,p) and 35 Cl(n, α) reactions that are primarily responsible for the fast neutron response of C 7 LYC. The light output of a 1” x1” C 7 LYC detector was simulated using the Birks equation for light quenching in scintillators. The simulations were benchmarked against neutron scattering data from a 12 C target obtained with a white neutron beam at the Los Alamos LANSCE facility and fourteen 1”x1” C 7 LYC detectors. Incident neutron energies extracted from time-of-flight measurements, and appropriate recoil corrections for the different detector angles were used to map the response of C 7 LYC to mono-energetic neutrons. Resonances in the 35 Cl(n,p) reaction cross-section below 3 MeV incident neutron energy were directly observed in C 7 LYC. The experimental data and simulations are compared with ENDF databases, as well as recent direct cross-section measurements of neutrons scattered from a 35 Cl target. In conclusion, the results suggest a scrutiny of the 35 Cl(n,p) and 35 Cl(n, α) cross sections at low energy in the ENDF library.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

T-odd leading-twist quark TMDs at small x

We study the small-x asymptotics of the flavor non-singlet T-odd leading-twist quark transverse momentum dependent parton distributions (TMDs), the Sivers and Boer-Mulders functions. While the leading eikonal small-x asymptotics of the quark Sivers function is given by the spin-dependent odderon [1, 2], we are interested in revisiting the sub-eikonal correction considered by us earlier in [3]. We first simplify the expressions for both TMDs at small Bjorken x and then construct small-x evolution equations for the resulting operators in the large-N c limit, with N c the number of quark colors. For both TMDs, the evolution equations resum all powers of the double-logarithmic parameter α s ln 2 (1/x), where α s is the strong coupling constant, which is assumed to be small. Solving these evolution equations numerically (for the Sivers function) and analytically (for the Boer-Mulders function) we arrive at the following leading small-x asymptotics of these TMDs at large N c :$$ {\displaystyle \begin{array}{l}{f}_{1T}^{\perp NS}\left(x\ll 1,{k}_T^2\right)={C}_O\left(x,{k}_T^2\right)\frac{1}{x}+{C}_1\left(x,{k}_T^2\right){\left(\frac{1}{x}\right)}^{3.4\sqrt{\frac{\alpha_s{N}_c}{4\pi }}}\\ {}{h}_1^{\perp \textrm{NS}}\left(x\ll 1,{k}_T^2\right)=C\left(x,{k}_T^2\right){\left(\frac{1}{x}\right)}^{-1}.\end{array}} $$ The functions C O (x,$ {k}_T^2 $), C 1 (x,$ {k}_T^2 $), and C(x,$ {k}_T^2 $) can be readily obtained in our formalism: they are mildly x-dependent and do not strongly affect the power-of-x asymptotics shown above. The function C O , along with the 1/x factor, arises from the odderon exchange. For the sub-eikonal contribution to the quark Sivers function (the term with C 1 ), our result shown above supersedes the one obtained in [3] due to the new contributions identified recently in [4].

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Back-to-back dijet production in DIS at arbitrary Bjorken x: TMD gluon distributions to twist-3 accuracy

We derive the gluon transverse-momentum-dependent (TMD) operator structure of back-to-back\\\\r\\\\nquark–antiquark dijet production in deep inelastic scattering at arbitrary Bjorken-x to twist-3 ac\\\\r\\\\ncuracy. Working at leading order in the strong coupling and in the kinematic regime where the\\\\r\\\\ntransverse momentum imbalance of the jets is much smaller than their individual transverse mo\\\\r\\\\nmenta, we perform a systematic gradient expansion of the quark propagator in a background gluon\\\\r\\\\nfield. This expansion organizes multiple interactions with the target in terms of longitudinal Wilson\\\\r\\\\nlines and gauge-invariant field-strength insertions, yielding a TMD description valid beyond the\\\\r\\\\nstrict high-energy eikonal (x → 0) approximation. We obtain explicit cross sections for longitudi\\\\r\\\\nnally and transversely polarized virtual photons, identifying all contributing gluon TMD operators\\\\r\\\\nup to twist-3, including structures involving F+−, Fij, and three-gluon correlators. The full lon\\\\r\\\\ngitudinal phase eixP+z− associated with Bjorken-x is retained throughout. In the small-x limit,\\\\r\\\\nour results reproduce the known sub-eikonal expressions obtained in the Color Glass Condensate\\\\r\\\\nframework, establishing a direct connection between the general-x TMD expansion and high-energy\\\\r\\\\nfactorization. We further reduce the operator basis using equations of motion, minimizing the num\\\\r\\\\nber of independent nonperturbative matrix elements entering the cross section. This work provides\\\\r\\\\na systematic foundation for extending TMD analyses of dijet production beyond leading twist, es\\\\r\\\\ntablishing a unified operator framework valid at arbitrary Bjorken-x that smoothly interpolates\\\\r\\\\nbetween moderate- and small-x descriptions of gluon TMDs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Quantal diffusion approach for multinucleon transfer processes in the 58,64 Ni + 208 Pb reactions: Toward the production of unknown neutron-rich nuclei

In recent years, substantial efforts have been made for the study of multinucleon transfer reactions at energies around the Coulomb barrier both experimentally and theoretically, aiming at the production of unknown neutron-rich heavy nuclei. It is crucial to provide reliable theoretical predictions based on microscopic theories with sufficient predictive power. Purpose: This article aims to clarify the applicability of the quantal diffusion approach based on the stochastic mean-field (SMF) theory for multinucleon transfer processes. Isotope production cross sections are evaluated for the reactions of 64 Ni + 208 Pb at E c.m. = 268 MeV and 58 Ni + 208 Pb at E c.m. = 270 MeV and are compared with available experimental data. Methods: Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are carried out for a range of initial orbital angular momenta with Skyrme SLy4d functional. Quantal diffusion equations, derived based on the SMF theory, for variances and covariance of neutron and proton numbers of reaction products are solved, with microscopic drift and diffusion coefficients obtained from time evolution of occupied single-particle orbitals in TDHF. Secondary de-excitation processes, both particle evaporation and fission, are simulated by a statistical compound-nucleus de-excitation model, gemini++. Dynamics of a fast isospin equilibration process followed by a slow drift toward the mass symmetry are commonly observed, as expected. Various reaction outcomes are evaluated, including average mass and charge numbers of reaction products, total kinetic energy loss (TKEL), scattering angle, contact time, and production cross sections for primary and secondary products. By comparing with the experimental data, we find that SMF and TDHF quantitatively reproduce experimental data for few-nucleon-transfer channels around the average values. In contrast, for many-nucleon-transfer channels, we find that the SMF approach provides much better description of the experimentally measured isotopic distributions. The results underline the importance of beyond-mean-field effects, especially one-body (mean-field) fluctuations and correlations, in describing multinucleon transfer processes. Moreover, through a combined analysis of SMF with a statistical model, gemini++, we find a significant contribution of transfer-induced fission, which is consistent with the experimental observation. In some cases, the SMF approach overestimates the isotopic width, requiring further improvements of the theoretical description. Possible ways to improve the description are discussed. The SMF approach is designed to describe the quantum many-body problem according to an ensemble of mean-field trajectories, taking into account part of many-body correlations in the description. As it requires feasible computational costs comparable to the ordinary TDHF approach, together with further model improvements, it will be a promising tool in the search for optimal reaction conditions to produce yet-unknown neutron-rich heavy nuclei through the multinucleon transfer reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Gluon splitting at small x: a unified derivation for the JIMWLK, DGLAP and CSS equations

We revisit the calculation of the next-to-leading order (NLO) corrections to dijet production in electron-ion collisions at small x. We focus on the back-to-back configuration where the relative transverse momentum P ⊥ of the measured jets is much larger than both their momentum imbalance K ⊥ and the target saturation momentum Q s (x, A). In this regime, we present for the first time a complete calculation of the real NLO corrections, associated with gluon emissions outside the jet cones, at leading power in 1/P ⊥ . Our result exhibits TMD factorisation, with the same hard factor as at tree-level and a NLO correction to the Weiszäcker-Williams (WW) gluon transverse momentum dependent (TMD) distribution which involves four Wilson-line operators. By studying different kinematical regimes for K ⊥ and for the radiated gluon, we recover all the quantum evolutions that were previously identified for this process at NLO: the B-JIMWLK high-energy evolution and the CSS evolution of the gluon WW TMD, and the DGLAP evolution of the gluon PDF. When both K ⊥ and the transverse momentum transferred by the target are large compared to Q s , all the Wilson-line operators boil down to the unintegrated gluon distribution and our NLO result for the gluon TMD can be used to isolate the transverse-momentum dependent gluon splitting function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗