Search NASA⌕ Search

SEARCH · Search NASA

Results for “General relativity equations & solutions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Exact solution of the DeWitt-Brehme-Hobbs equation in copropagating electromagetic and gravitational waves

An accelerated charge interacts with its own electromagnetic field, a phenomenon known as electromagnetic radiation reaction. The DeWitt-Brehme-Hobbs (DWBH) equation describes the motion of a charged mass in the presence of combined electromagnetic and gravitational fields, taking into account electromagnetic radiation-reaction effects. Here, we find the first exact analytical solution of the DWBH equation in the case of a charged mass in the presence of copropagating and otherwise arbitrary electromagnetic and gravitational plane waves. As a consequence of the Penrose limit, the scenario considered here can be seen as a local limit around ultrarelativistic trajectories in a general curved spacetime. Finally, the paradigmatic example of an electromagnetic wave in the presence of a constant-amplitude gravitational wave is worked out explicitly and it is shown how the presence of the gravitational wave can qualitatively change electromagnetic radiation-reaction effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Einstein-Maxwell theory and the Weyl double copy

The Weyl double copy relates vacuum solutions in general relativity to Abelian gauge fields in Minkowski spacetime. In a previous work, we showed how the Weyl double copy can be extended to provide a treatment of external gravitational sources consistent with the classical Kerr-Schild double copy. Using this generalization, here we provide a complete double copy analysis of electrovacuum Petrov type D spacetimes. This includes the first analysis of the charged C-metric, whose single copy interpretation invokes the two-potential formalism of electrodynamics. Here, we also present the first double copy prescription for the Ricci spinor, which for non-accelerating spacetimes, takes a form similar to the original double copy relation for the Weyl spinor.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Type-II spacetimes and the double copy for fluids metrics

In our previous paper [Keeler , From Navier-Stokes to Maxwell via Einstein, .] we discussed type-D and type-N fluid-dual spacetimes and provided their associated single copies in the context of the Weyl double copy. In this work we extend our analysis to more general fluids thereby requiring the application of the double-copy picture to type-II spacetimes. By combining our previous type-D and type-N fluids via their associated stream functions we construct an example of a viable type-II double copy. We show that the gravity duals of these fluids perturbatively satisfy Einstein’s equations. We use a near-horizon expansion to identify the type-II double copy for the fluid-dual spacetimes. We show a Maxwell spinor ansatz containing a heterogeneous bispinor component is necessary to provide a viable type-II double copy at the lowest order.

Fluid-gravity correspondence↗
DOE OSTI2023

Halo formation from Yukawa forces in the very early Universe

If long-range attractive forces exist and are stronger than gravity then cosmic halo formation can begin in the radiation-dominated era. We study a simple realization of this effect in a system where dark matter fermions have Yukawa interactions mediated by scalar particles, analogous to the Higgs boson in the standard model. We develop a self-consistent description of the system including exact background dynamics of the scalar field, and precise modeling of the fermion density fluctuations. For the latter, we provide accurate approximations for the linear growth as well as quantitative modeling of the nonlinear evolution using N-body simulations. We find that halo formation occurs exponentially fast and on scales substantially larger than simple estimates predict. Here, the final fate of these halos remains uncertain, but could be annihilation, dark stars, primordial black holes, or even the existence of galaxy-sized halos at matter-radiation equality. More generally, our results demonstrate the importance of mapping scalar-mediated interactions onto structure formation outcomes and constraints for beyond the standard model theories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum field theory, worldline theory, and spin magnitude change in orbital evolution

A previous paper [Z. Bern et al., Binary dynamics through the fifth power of spin at 𝑂⁡(𝐺 2 ), Phys. Rev. Lett. 130, 201402 (2023)] identified a puzzle stemming from the amplitudes-based approach to spinning bodies in general relativity: additional Wilson coefficients appear compared to current worldline approaches to conservative dynamics of generic astrophysical objects, including neutron stars. In this paper we clarify the nature of analogous Wilson coefficients in the simpler theory of electrodynamics. We analyze the original field-theory construction, identifying definite-spin states some of which have negative norms, and relating the additional Wilson coefficients in the classical theory to transitions between different quantum spin states. We produce a new version of the theory which also has additional Wilson coefficients, but no negative-norm states. We match, through 𝒪⁡(𝛼 2 ) and 𝒪⁡(𝑆 2 ), the Compton amplitudes of these field theories with those of a modified worldline theory with extra degrees of freedom introduced by releasing the spin supplementary condition. We build an effective two-body Hamiltonian that matches the impulse and spin kick of the modified field theory and of the worldline theory, displaying additional Wilson coefficients compared to standard worldline approaches. The results are then compactly expressed in terms of an eikonal formula. Our key conclusion is that, contrary to standard approaches, while the magnitude of the spin tensor is still conserved, the magnitude of the spin vector can change under conserved Hamiltonian dynamics and this change is governed by the additional Wilson coefficients. For specific values of Wilson coefficients the results are equivalent to those from a definite spin obeying the spin supplementary condition, but for generic values they are physically inequivalent. These results warrant detailed studies of the corresponding issues in general relativity.

classical black holes↗
DOE OSTI2025

General relativistic hydrodynamic simulations of binary strange star mergers

We perform fully general-relativistic simulations of binary strange star mergers considering two different approaches for thermal effects. The first uses a cold equation of state (EOS) derived from a modified version of the MIT bag model which is then supplemented by a Γ-law correction. The second approach employs a microphysical description of the finite-temperature effects. We describe results obtained with the two treatments, highlighting the influence of thermal effects. We find that the postmerger dynamics differs significantly in the two cases, leading to quantitative differences in the postmerger gravitational-wave spectrum and ejecta mass. The peak frequency of the postmerger gravitational-wave emission is consistent with the established quasi-universal relations for binary neutron star mergers and as a result, our simulations cannot distinguish between mergers of neutron stars and those of strange stars. Our models with realistic treatment of finite-temperature effects produce a significant amount of ejecta ≳0.02 M ⊙ ​. Here, the resulting flux of strangelets near the Earth, computed assuming that all neutron star mergers are in fact strange-stars mergers and that the binary considered here is representative, is in tension with experimental upper limits. As such, our results tentatively disfavor a scenario in which strange-quark matter is the lowest energy state of matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nonlinearities in Black Hole Ringdowns

The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter, we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain’s (ℓ,m)=(4,4) angular harmonic, here we show the presence of a quadratic effect across a range of binary BH mass ratios that agrees with theoretical expectations. We find that the quadratic (4,4) mode’s amplitude exhibits quadratic scaling with the fundamental (2,2) mode—its parent mode. The nonlinear mode’s amplitude is comparable to or even larger than that of the linear (4,4) mode. Therefore, correctly modeling the ringdown of higher harmonics—improving mode mismatches by up to 2 orders of magnitude—requires the inclusion of nonlinear effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Constraining gravity with a new precision 𝐸 𝐺 estimator using Planck + SDSS BOSS data

The 𝐸 𝐺 statistic is a discriminating probe of gravity developed to test the prediction of general relativity (GR) for the relation between gravitational potential and clustering on the largest scales in the observable Universe. We present a novel high-precision estimator for the 𝐸 𝐺 statistic using CMB lensing and galaxy clustering correlations that carefully matches the effective redshifts across the different measurement components to minimize corrections. A suite of detailed tests is performed to characterize the estimator’s accuracy, its sensitivity to assumptions and analysis choices, and the non-Gaussianity of the estimator’s uncertainty is characterized. After finalization of the estimator, it is applied to Planck CMB lensing and SDSS CMASS and LOWZ galaxy data. We report the first harmonic space measurement of 𝐸 𝐺 using the LOWZ sample and CMB lensing and also updated constraints using the final CMASS sample and the latest Planck CMB lensing map. We find $\hat{𝐸}$$^{Planck+CMASS}_{𝐺}$ = 0.3⁢6$^{+0.06}_{−0.05}$⁢(68.27%) and $\hat{𝐸}$$^{Planck+LOWZ}_{𝐺}$ = 0.4⁢0$^{+0.11}_{−0.09}$⁢(68.27%), with additional subdominant systematic error budget estimates of 2% and 3%, respectively. Using Ω m,0 constraints from Planck and SDSS BAO observations, Λ⁢CDM-GR predicts 𝐸$^{GR}_ {𝐺}$⁡(𝑧 =0.555) = 0.401 ± 0.005 and 𝐸$^{GR}_{𝐺}$⁡(𝑧 =0.316) = 0.452 ± 0.005 at the effective redshifts of the CMASS and LOWZ based measurements. We report the measurement to be in good statistical agreement with the Λ⁢CDM-GR prediction and report that the measurement is also consistent with the more general GR prediction of scale independence for 𝐸 𝐺 . Furthermore, this work provides a carefully constructed and calibrated statistic with which 𝐸 𝐺 measurements can be confidently and accurately obtained with upcoming survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Amplitudes, supersymmetric black hole scattering at $\mathcal{O}\left({G}^5\right)$, and loop integration

We compute the potential-graviton contribution to the scattering amplitude, the radial action, and the scattering angle of two extremal black holes in $\mathcal{N}$ = 8 supergravity at the fifth post-Minkowskian order and to next-to-leading order in a large mass expansion (first self-force order). Properties of classical unitarity cuts allow us to focus on the integration-by-parts reduction of planar integrals, while nonplanar integrals at this order are obtained from the planar ones by straightforward manipulations. We present the solution to the differential equations for all master integrals necessary to evaluate the classical scattering amplitudes of massive scalar particles at this order in all gravitational theories, in particular in $\mathcal{N}$ = 8 supergravity, and in general relativity. Despite the appearance of higher-weight generalized polylogarithms and elliptic functions in the solution to the differential equation for master integrals, the final supergravity answer is remarkably simple and contains only (harmonic) polylogarithmic functions up to weight 2. The systematic analysis of elliptic integrals discussed here, as well as the particular organization of boundary integrals in $\mathcal{N}$ = 8 observables are independent of supersymmetry and may have wider applications, including to aspects of collider physics.

Black Holes↗
DOE OSTI2022

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Herein we formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results [1]-[4], indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Poisson Equation for a (General) Homogeneous d-Dimensional Ellipsoid with Applications to Beam Envelope Tracking

This note describes the solution of the free-space Poisson equation in the interior of a $d$-dimensional homogeneous ellipsoid, and the associated space charge fields. An explicit formula (\ref{Sformula}) is provided that relates the $d\times d$ matrix describing the space charge (quadratic) potential to the $d\times d$ covariance matrix of the ellipsoid. For the cases $d=2$ and $d=3$, this result is used to determine the linear map corresponding to a space charge kick, that may be used to push the beam $6\times 6$ covariance matrix during envelope tracking. The treatment of upright ellipsoids for $d=2$ and $d=3$ is well-represented in the literature. However, the approach taken here emphasizes a general ellipsoid with arbitrary correlations in any dimension. The Appendix provides a general solution of the free-space Poisson equation in dimension $d$ for a source distribution with ellipsoidal symmetry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗