SEARCH · Engineering Papers

Results for “linear 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 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗

DOE OSTI2026

Floquet Spin Splitting and Spin Generation in Antiferromagnets

In antiferromagnetic spintronics, accessing the spin degrees of freedom is essential for generating spin currents and manipulating magnetic order, which generally requires lifting spin degeneracy. This is typically achieved through relativistic spin-orbit coupling or nonrelativistic spin splitting in altermagnets. Here, we propose an alternative approach: a dynamical spin splitting induced by an optical field in antiferromagnets. By coupling the driven system to a thermal bath, we demonstrate the emergence of steady-state pure spin currents as well as linear-response longitudinal and transverse spin currents. Crucially, thermal bath engineering enables a nonrelativistic Edelstein effect—the generation of a net spin accumulation—without relying on spin-orbit coupling. Furthermore, our results provide a broadly applicable and experimentally tunable route to control spins in antiferromagnets, offering new opportunities for spin generation and manipulation in antiferromagnetic spintronics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Stress Tensor flows, birefringence in non-linear electrodynamics and supersymmetry

We identify the unique stress tensor deformation which preserves zero-birefringence conditions in non-linear electrodynamics, which is a 4d 4 d version of the T\overline{T} T T ¯ operator. We study the flows driven by this operator in the three Lagrangian theories without birefringence - Born-Infeld, Plebanski, and reverse Born-Infeld - all of which admit ModMax-like generalizations using a root- T\overline{T} T T ¯ -like flow that we analyse in our paper. We demonstrate one way of making this root- T\overline{T} T T ¯ -like flow manifestly supersymmetric by writing the deforming operator in \mathcal{N} = 1 𝒩 = 1 superspace and exhibit two examples of superspace flows. We present scalar analogues in d = 2 d = 2 with similar properties as these theories of electrodynamics in d = 4 d = 4 . Surprisingly, the Plebanski-type theories are fixed points of the classical T\overline{T} T T ¯ -like flows, while the Born-Infeld-type examples satisfy new flow equations driven by relevant operators constructed from the stress tensor. Finally, we prove that any theory obtained from a classical stress-tensor-squared deformation of a conformal field theory gives rise to a related “subtracted” theory for which the stress-tensor-squared operator is a constant.

Physics↗

DOE OSTI2021

Excitation spectra of quantum matter without quasiparticles. II. Random t - J models

Here, we present numerical solutions of the spectral functions of large dimension t -J models with random nearest neighbor exchange and global SU(M) spin rotation symmetry. The solutions are obtained from the saddle-point equations of the large dimension limit, followed by the large M limit. The same saddle point equations also apply to the model with all-to-all and random hopping and exchange. Such a t -J model realizes a deconfined critical state proposed as a model for the optimally doped cuprates. The large M theory involves Green’s functions for fractionalized spinons and holons carrying emergent U(1) gauge charges, obeying relations similar to those of the Sachdev-Ye-Kitaev (SYK) models. The low frequency spectral functions are compared with an analytic analysis of the operator scaling dimensions with good agreement. We also compute the low frequency and temperature behavior of experimentally observable gauge-invariant observables: the electron Green’s function, the local spin susceptibility and the optical conductivity, along with the temperature dependence of the d.c. resistivity. The time reparameterization soft mode (equivalent to the boundary graviton in holographically dual models of two-dimensional quantum gravity) makes important contributions to all observables and provides a linear-in-temperature contribution to the d.c. resistivity.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Space-Time Fluid (Unabridged)

Purpose: This essay is a retelling of general relativity in a language in which space-time geometry is expressed as a fluid. This trivial and useful reformulation gives 1) a non-perturbative covariant description of cosmological inhomogeneities and 2) a simple formula describing how cosmic inhomogeneities are generated on super-horizon scales. Methods: Equating the Ricci curvature with the associated matter stress-energy gives a description of space-time geometry in terms of fluid properties. These locally measurable (covariant) non-perturbative quantities are in some ways superior to commonly used "gauge invariant" quantities. The dynamics of a quantity (kurvature) which describes cosmological inhomogeneities is described in detail. A detailed comparison is made of space-time fluid dynamics with that of a classical (Newtonian physics) fluid. Results: The fluid lexicon permits an unambiguous definition of the velocity of space-time. The evolution of the space-time fluid is in many ways identical with that of the classical fluid when expressed in Lagrangian coordinates. Kurvature is a measure of the specific binding energy of the fluid and is a most useful covariant measure of cosmological inhomogeneities. For plausible matter models kurvature will increase, even on super-horizon scales, due to non-linear hydrodynamic effects rather than gravity. This phenomena is also exhibited by classical fluids. Conclusion: The space-time fluid representation of geometrodynamics gives a simple and useful description of the evolution of cosmological inhomogeneities.

Stebbins, Albert [Fermilab] (ORCID:000000023807725↗

DOE OSTI2021

Intrinsic Anomalous Hall Conductivity in a Nonuniform Electric Field

We study how the intrinsic anomalous Hall conductivity is modified in two-dimensional crystals with broken time-reversal symmetry due to weak inhomogeneity of the applied electric field. Focusing on a clean noninteracting two-band system without band crossings, we derive the general expression for the Hall conductivity at small finite wave vector q to order q 2 , which governs the Hall response to the second gradient of the electric field. Using the Kubo formula, we show that the answer can be expressed through the Berry curvature, Fubini-Study quantum metric, and the rank-3 symmetric tensor which is related to the quantum geometric connection and physically corresponds to the gauge-invariant part of the third cumulant of the position operator. We further compare our results with the predictions made within the semiclassical approach. By deriving the semiclassical equations of motion, we reproduce the result obtained from the Kubo formula in some limits. We also find, however, that the conventional semiclassical description in terms of the definite position and momentum of the electron is not fully consistent because of singular terms originating from the Heisenberg uncertainty principle. We thus present a clear example of a case when the semiclassical approach inherently suffers from the uncertainty principle, implying that it should be applied to systems in nonuniform fields with extra care.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Test of Causal Nonlinear Quantum Mechanics by Ramsey Interferometry with a Trapped Ion

Quantum mechanics requires the time evolution of the wave function to be linear. While this feature has been associated with the preservation of causality, a consistent causal nonlinear theory was recently developed. Interestingly, this theory is unavoidably sensitive to the full physical spread of the wave function, rendering existing experimental tests for nonlinearities inapplicable. Here, in this work, using well-controlled motional superpositions of a trapped ion, we set a stringent limit of 5.4×10 –12 on the magnitude of the unitless scaling factor $\tilde{ε}_γ$ for the predicted causal nonlinear perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Stochastic Error Cancellation in Analog Quantum Simulation

Analog quantum simulation is a promising path towards solving classically intractable problems in many-body physics on near-term quantum devices. However, the presence of noise limits the size of the system and the length of time that can be simulated. In our work, we consider an error model in which the actual Hamiltonian of the simulator differs from the target Hamiltonian we want to simulate by small local perturbations, which are assumed to be random and unbiased. We analyze the error accumulated in observables in this setting and show that, due to stochastic error cancellation, with high probability the error scales as the square root of the number of qubits instead of linearly. We explore the concentration phenomenon of this error as well as its implications for local observables in the thermodynamic limit. Moreover, we show that stochastic error cancellation also manifests in the fidelity between the target state at the end of time-evolution and the actual state we obtain in the presence of noise. This indicates that, to reach a certain fidelity, more noise can be tolerated than implied by the worst-case bound if the noise comes from many statistically independent sources.

Analog quantum simulation↗

DOE OSTI2024

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗

DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗

DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

The Optical Signatures of Stochastic Processes in Many-Body Exciton Scattering

We review our recent quantum stochastic model for spectroscopic lineshapes in the presence of a coevolving and nonstationary background population of excitations. Starting from a field theory description for interacting bosonic excitons, we derive a reduced model whereby optical excitons are coupled to an incoherent background via scattering as mediated by their screened Coulomb coupling. The Heisenberg equations of motion for the optical excitons are then driven by an auxiliary stochastic population variable, which we take to be the solution of an Ornstein–Uhlenbeck process. Here, we present an overview of the theoretical techniques we have developed as applied to predicting coherent nonlinear spectroscopic signals. We show how direct (Coulomb) and exchange coupling to the bath give rise to distinct spectral signatures and discuss mathematical limits on inverting spectral signatures to extract the background density of states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Univariate Prediction of Hammett Parameters and Select Relative Reaction Rates Using Loewdin Atomic Charges

Loewdin charges from density functional theory calculations were used here to obtain general, univariate linear correlations for the prediction of experimental Hammett parameters and relative reaction rates. While previous studies have established that Hirshfeld and CM5 charges perform strongly as univariate predictors, the near-ubiquitous Loewdin charges have not yet been evaluated. To this end, we assess the predictive capability of Loewdin charges for three chemical systems. First, we show that Loewdin charges outperform Hirshfeld and CM5 charges for Hammett parameter prediction. Second, we see that Loewdin charges generally perform comparably to Hirshfeld charges for predicting the relative rates of olefin cleavage by photoexcited nitroarenes. The single case of poor correlation, between relative rates and the Loewdin charges on nitrogen sites, is ameliorated when considering the net charge on the NO 2 group. Third, we show that Loewdin, Hirshfeld, and CM5 charges all perform very well for generating correlations for relative reaction rates for C–H activation of 9-(4-X-phenyl)-9H-fluorene substrates by a transition metal catalyst. The equations generated throughout the study enable the prediction of Hammett parameters and relative reaction rates. Finally, these tools can accelerate synthetic and experimental studies by enabling the in silico prediction of uncharacterized chemical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗