Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix 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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Rigid-Mode Limit of the Yokoya Matrix Formalism and the Burov-Lebedev Dispersion Equation

Transverse single-bunch instabilities of space-charge-dominated coasting beams with round and flat transverse geometries are studied using a unified dispersion-relation framework. The analysis combines the Burov-Lebedev formalism, which captures space-charge tune spread, Landau damping, and instability threshold behavior, with Yokoya’s projection method for representing coherent transverse mode structure and its dependence on beam aspect ratio. In the rigid-beam limit, the formulation reduces to a scalar dispersion relation of Burov-Lebedev paper. For non-rigid transverse oscillations, truncation of Yokoya’s Hermite-based expansion yields a finite-dimensional matrix eigenvalue problem in which space-charge and coupling impedance effects enter through Burov-Lebedev–type denominators. This approach provides a consistent basis for comparing rigid and non-rigid instability behavior in round and flat beams and for assessing the role of beam ellipticity in modifying coherent mode structure and stability thresholds.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Evaluation of data driven low-rank matrix factorization for accelerated solutions of the Vlasov equation

Low-rank methods have shown success in accelerating simulations of a collisionless plasma described by the Vlasov equation, but still rely on computationally costly linear algebra every time step. We propose a data-driven factorization method using artificial neural networks, specifically with convolutional layer architecture, that trains on existing simulation data. At inference time, the model outputs a low-rank decomposition of the distribution field of the charged particles, and we demonstrate that this step is faster than the standard linear algebra technique. Numerical experiments show that the method achieves comparable reconstruction accuracy for interpolation tasks, generalizing to unseen test data in a manner beyond just memorizing training data; patterns in factorization also inherently followed the same numerical trend as those within algebraic methods (e.g., truncated singular-value decomposition). However, when training on the first 70% of a time-series data and testing on the remaining 30%, the method fails to meaningfully extrapolate. Despite this limiting result, the technique may have benefits for simulations in a statistical steady-state or otherwise showing temporal stability. These results suggest that while the model offers a computationally efficient alternative for datasets with temporal stability, its current formulation is best suited for interpolation rather than for predicting future states in time-evolving systems. This study thus lays the groundwork for further refinement of neural network-based approaches to low-rank matrix factorization in high-dimensional plasma simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Dual Representations and H ∞ -Optimal Control of Partial Differential Equations

We consider H ∞ -optimal state-feedback control of the class of linear Partial Differential Equations (PDEs) which admit a Partial Integral Equation (PIE) representation. While linear matrix inequalities are commonly used for optimal control of Ordinary Differential Equations (ODEs), the absence of a universal state-space representation and suitable dual form prevents such methods from being applied to optimal control of PDEs. Specifically, for ODEs, the controller synthesis problem is defined in state-space, and duality is used to resolve the bilinearity of that synthesis problem. Recently, the PIE representation was proposed as a universal state-space representation for linear PDE systems. In this paper, we show that any PDE system represented by a PIE admits a dual PIE with identical stability and I/O properties. This result allows us to reformulate the stabilizing and optimal state-feedback control problems as convex optimization over the cone of positive Partial Integral (PI) operators. Operator inversion formulae then allow us to construct feedback gains for the original PDE system. The results are verified through application to several canonical problems in optimal control of PDEs and indicate the resulting bounds on H ∞ norm are not conservative.

42 ENGINEERING↗
DOE OSTI2026

Active learning emulators for nuclear two-body scattering in momentum space

In this work we extend the active learning emulators for two-body scattering in coordinate space with error estimation, recently developed by Maldonado et al. [Phys. Rev. C 112, 024002], to coupled-channel scattering in momentum space. Our full-order model (FOM) solver is based on the Lippmann-Schwinger integral equation for the scattering t-matrix as opposed to the radial Schrödinger equation. We use (Petrov-)Galerkin projections and high-fidelity calculations at a few snapshots across the parameter space of the interaction to construct efficient reduced-order models (ROMs), trained by a greedy algorithm for locally optimal snapshot selection. Both the FOM solver and the corresponding ROMs are implemented efficiently in Python using Google's JAX library. We present results for emulating scattering phase shifts in coupled and uncoupled channels and cross sections, and assess the accuracy of the developed ROMs and their computational speedup factors. We also develop emulator error estimation for both the t-matrix and the total cross section. The software framework for reproducing and extending our results is publicly available. Together with our recent advances in developing active-learning emulators for three-body scattering, these emulator frameworks set the stage for full Bayesian calibrations of chiral nuclear interactions and optical models against scattering data with quantified emulator errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Posterior Covariance Matrix Approximations

Here, the Davis equation of state (EOS) is commonly used to model thermodynamic relationships for high explosive (HE) reactants. Typically, the parameters in the EOS are calibrated, with uncertainty, using a Bayesian framework and Markov Chain Monte Carlo (MCMC) methods. However, MCMC methods are computationally expensive, especially for complex models with many parameters. This paper provides a comparison between MCMC and less computationally expensive Variational methods (Variational Bayesian and Hessian Variational Bayesian) for computing the posterior distribution and approximating the posterior covariance matrix based on heterogeneous experimental data. All three methods recover similar posterior distributions and posterior covariance matrices. This study demonstrates that for this EOS parameter calibration application, the assumptions made in the two Variational methods significantly reduce the computational cost but do not substantially change the results compared to MCMC.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

SuperScreen: An open-source package for simulating the magnetic response of two-dimensional superconducting devices

Quantitative understanding of the spatial distribution of magnetic fields and Meissner screening currents in two-dimensional (2D) superconductors and mesoscopic thin film superconducting devices is critical to interpreting the results of magnetic measurements of such systems. Here, we introduce SuperScreen, an open-source Python package for simulating the response of 2D superconductors to trapped flux and applied time-independent or quasi-DC magnetic fields for any value of the effective magnetic penetration depth, Λ. Given an applied magnetic field, SuperScreen solves the 2D London equation using an efficient matrix inversion method to obtain the Meissner currents and magnetic fields in and around structures composed of one or more superconducting thin films of arbitrary geometry. Further, SuperScreen can be used to model screening effects and calculate self- and mutual-inductance in thin film superconducting devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Greedy emulators for nuclear two-body scattering

Applications of reduced basis method emulators are increasing in low-energy nuclear physics because they enable fast and accurate sampling of high-fidelity calculations, enabling robust uncertainty quantification. Here, in this paper, we develop, implement, and test two model-driven emulators based on the (Petrov-)Galerkin projection using the prototypical test case of two-body scattering with the Minnesota potential and a more realistic local chiral potential. The high-fidelity scattering equations are solved with the matrix Numerov method, a reformulation of the popular Numerov recurrence relation for solving special second-order differential equations as a linear system of coupled equations. A novel error estimator based on reduced-space residuals is applied to an active learning approach (a greedy algorithm) to choosing training samples (“snapshots”) for the emulator and contrasted with a proper orthogonal decomposition (POD) approach. Both approaches allow for computationally efficient offline-online decompositions, but the greedy approach requires many fewer snapshot calculations. These developments set the groundwork for emulating scattering observables based on chiral nucleon-nucleon and three-nucleon interactions and optical models, where computational speed-ups are necessary for Bayesian uncertainty quantification. Our emulators and error estimators are widely applicable to linear systems.

Bayesian methods↗
DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2023

Scalar Breit interaction for molecular calculations

Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Au n (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. In conclusion, the scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multidimensional coherent spectroscopy of light-driven states and their collective modes in multiband superconductors

We present a comprehensive theory of light-controlled multiband superconductivity, and apply it to predict distinctive signatures of light-driven superconducting (SC) states in terahertz multidimensional coherent spectroscopy (THz-MDCS) experiments. We first derive gauge-invariant Maxwell-Bloch equations for multiband BCS superconductors with spatial fluctuations. We consider driving electromagnetic fields determined self-consistently by Maxwell's equations. By calculating the THz-MDCS spectra measured experimentally in the clean SC limit, we identify unique signatures of finite-momentum Cooper-pairing states that live longer than the laser pulse. They are controlled by a pair of THz laser pulses with well-defined relative phase (pulse pair). The pseudospin oscillators that describe the properties of these SC states are parametrically driven by both finite-momentum Cooper pairing and by time oscillations of the order-parameter relative phase. We show that such strong parametric driving leads to drastic changes in the THz-MDCS spectral shape from the predictions of third-order nonlinear susceptibility calculations. These spectral changes strongly depend on the interband-to-intraband interaction ratio and on the collective modes of the light-driven state. For negligible interband interaction, the spectra show a transition with increasing field, from traditional pump-probe, four-wave-mixing, and third-harmonic generation peaks determined by the laser frequency to sidebands determined by the excitations of the driven system. These sidebands emerge from difference-frequency Raman processes in the nonequilibrium SC state. For interband couplings weaker than the intraband pairing, we show that the Leggett phase collective mode leads to harmonic sidebands around the traditional pump-probe peaks. Additional Higgs collective mode peaks result from light-induced inversion-symmetry breaking in a thin-film geometry. For strong interband coupling, we find a transition from a nonequilibrium finite Cooper-pair momentum state characterized by hybrid-Higgs amplitude mode peaks in THz-MDCS spectra to a driven state identified experimentally by the emergence of Floquet-type sidebands at bi-Higgs frequencies. Those dominant bi-Higgs-frequency satellites are manifestations of a new order parameter relative phase collective mode that characterizes the nonequilibrium SC state. The predicted interaction- and field-dependent transitions in the spectral profile allow us to propose THz-MDCS experiments for quantum tomography of light-driven superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Field Theory of the Fermi Function

The Fermi function F ( Z , E ) accounts for QED corrections to beta decays that are enhanced at either small electron velocity β or large nuclear charge Z . For precision applications, the Fermi function must be combined with other radiative corrections and with scale- and scheme-dependent hadronic matrix elements. We formulate the Fermi function as a field theory object and present a new factorization formula for QED radiative corrections to beta decays. We provide new results for the anomalous dimension of the corresponding effective operator complete through three loops, and resum perturbative logarithms and π enhancements with renormalization-group methods. Our results are important for tests of fundamental physics with precision beta decay and related processes. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A supersymmetric SYK model with a curious low energy behavior

We consider N = 2,4 supersymmetric SYK models that have a peculiar low energy behavior, with the entropy going like S = S 0 + (constant)T a , where a ≠ 1. The large N equations for these models are a generalization of equations that have been previously studied as an unjustified truncation of the planar diagrams describing the BFSS matrix quantum mechanics or other related matrix models. Here we reanalyze these equations in order to better understand the low energy physics of these models. We find that the scalar fields develop large expectation values which explore the low energy valleys in the potential. The low energy physics is dominated by quadratic fluctuations around these values. These models were previously conjectured to have a spin glass phase. We did not find any evidence for this phase by using the usual diagnostics, such as searching for replica symmetry breaking solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fluctuations in Hill’s equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton’s energy into Standard Model particles. This conversion often proceeds through nonperturbative mechanisms such as parametric resonance, which is described by Hill’s equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill’s equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. Furthermore, these results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.

Cosmology↗
DOE OSTI2026

A limit on the total lepton number in the Universe from BBN and the CMB

At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗