Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Relativistic wave 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗
DOE OSTI2023

Spin-dependent interactions and heavy-quark transport in the quark-gluon plasma

We extend a previously constructed T-matrix approach to the quark-gluon plasma (QGP) to include the effects of spin-dependent interactions between partons. Following earlier work within the relativistic quark model, the spin-dependent interactions figure as relativistic corrections to the Cornell potential. When applied to the vacuum spectroscopy of quarkonia, in particular their mass splittings in S- and P-wave states, the issue of the Lorentz structure of the confining potential arises. We confirm that a significant admixture of a vector interaction (to the previously assumed scalar interaction) improves the description of the experimental mass splittings. The temperature corrections to the in-medium potential are constrained by results from thermal lattice-QCD for the equation of state (EoS) and heavy-quark (HQ) free energy in a self-consistent set-up for heavy- and light-parton spectral functions in the QGP. We then deploy the refined in-medium heavy-light T-matrix to compute the charm-quark transport coefficients in the QGP. The vector component of the confining potential, through its relativistic corrections, enhances the friction coefficient for charm quarks in the QGP over previous calculations by tens of percent at low momenta and temperatures, and more at higher momenta. Our results are promising for improving the current phenomenology of open heavy-flavor observables at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Investigating color screening in proton-nucleus collisions with complex potentials*

Abstract Color screening and parton inelastic scattering modify the heavy-quark antiquark potential in mediums consisting of particles from quantum chromodynamics (QCD), leading to the suppression of quarkonium production in relativistic heavy-ion collisions. Owing to the small charm/anti-charm ( ) pair production number in proton-nucleus ( pA ) collisions, the correlation between different pairs is negligible, which makes the Schrödinger equation viable for tracking the evolution of only one pair. We employ the time-dependent Schrödinger equation with an in-medium potential to study the evolution of charmonium wave functions in a hydrodynamic-like QCD medium produced in pA collisions. We explore different parametrizations of real and imaginary parts of the potential and calculate the nuclear modification factors ( ) of and in TeV energy p -Pb collisions at the Large Hadron Collider (LHC). Comparing strong and weak screening scenarios with experimental data in this approach, we arrive at the conclusion that color screening is weak at temperatures close to the deconfined phase transition. Moreover, the imaginary part of the potential is crucial in describing the experimental data, which is consistent with widely studied semi-classical approaches, where dissociation rates are essential.

Physics↗
NASA NTRS2010

Merging Black Holes

Black-hole mergers take place in regions of very strong and dynamical gravitational fields, and are among the strongest sources of gravitational radiation. Probing these mergers requires solving the full set of Einstein's equations of general relativity numerically. For more than 40 years, progress towards this goal has been very slow, as numerical relativists encountered a host of difficult problems. Recently, several breakthroughs have led to dramatic progress, enabling stable and accurate calculations of black-hole mergers. This article presents an overview of this field, including impacts on astrophysics and applications in gravitational wave data analysis.

Centrella, Joan↗
DOE OSTI2025

Finite- and infinite-volume study of DDπ scattering

We develop a comprehensive framework for extracting the pole position and properties of the doubly-charmed tetraquark $T$$^{+}_{cc}$(3875) from lattice QCD data using the relativistic three-particle formalism. This approach incorporates the effect of the one-pion exchange diagram in DDπ and DD ∗ scattering, making it applicable at energies coinciding with the left-hand cut in the partial-wave projected DD ∗ amplitude. We present an example application of this framework to existing lattice QCD data at mπ = 280 MeV. We solve the integral equations describing the DDπ reaction, use LSZ reduction to determine the corresponding DD ∗ amplitude, and find the values of the infinite-volume two- and three-body K matrices that lead to agreement with lattice DD ∗ phase shifts within their uncertainties. Using these K matrices in the three-particle quantization condition, we describe the finite- volume DD ∗ spectrum and find good agreement with the lattice QCD energies. Our results suggest that, at this pion mass, the tetraquark appears as a pair of subthreshold complex poles whose precise location strongly depends on the value of the DDπ three-particle K matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2006

Self-Consistent Ring Current Modeling with Propagating Electromagnetic Ion Cyclotron Waves in the Presence of Heavy Ions

The self-consistent treatment of the RC ion dynamics and EMIC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm. That is why the modeling of EMIC waves is critical and timely issue in magnetospheric physics. To describe the RC evolution itself this study uses the ring current-atmosphere interaction model (RAM). RAM solves the gyration and bounce-averaged Boltzmann-Landau equation inside of geosynchronous orbit. Originally developed at the University of Michigan, there are now several branches of this model currently in use as describe by Liemohn namely those at NASA Goddard Space Flight Center This study will generalize the self-consistent theoretical description of RC ions and EMIC waves that has been developed by Khazanov and include the heavy ions and propagation effects of EMIC waves in the global dynamic of self-consistent RC - EMIC waves coupling. The results of our newly developed model that will be presented at GEM meeting, focusing mainly on the dynamic of EMIC waves and comparison of these results with the previous global RC modeling studies devoted to EMIC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

Khazanov, George V.↗
DOE OSTI2021

Applying the Euclidean-signature semi-classical method to the quantum Taub models with a cosmological constant and aligned electromagnetic field

We prove the existence of a countably infinite number of “excited” states for the Lorentzian-signature Taub–Wheeler–DeWitt (WDW) equation when a cosmological constant is present using the Euclidean-signature semi-classical method. We also find a “ground” state solution when both an aligned electromagnetic field and cosmological constant are present; as a result, conjecture that the Euclidean-signature semi classical method can be used to prove the existence of a countably “infinite” number of “excited” states when the two aforementioned matter sources are present. Afterward, we prove the existence of asymptotic solutions to the vacuum Taub–WDW equation using the “no boundary” and “wormhole” solutions of the Taub Euclidean-signature Hamilton–Jacobi equation and compare their mathematical properties. We then discuss the fascinating qualitative properties of the wave functions we have computed. By utilizing the Euclidean-signature semi-classical method in the above manner, we further show its ability to prove the existence of solutions to Lorentzian-signature equations without having to invoke a Wick rotation. This feature of not needing to apply a Wick rotation makes this method potentially very useful for tackling a variety of problems in bosonic relativistic field theory and quantum gravity

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Radiation from highly relativistic geodesics

A number of recent works are reviewed concerning the generation and emission of gravitational waves. It is shown that at high frequencies, the generation of gravitational radiation is a local phenomenon. Two examples are described illustrating this generation when a high-energy particle collides against the space-time curvature. One, after Matzner and Nutku, uses a method of virtual photons; the other, after Chrzanowski and Misner, is based on the W.K.B. approximation, corresponding to geometric optics, for the inhomogeneous wave equation. This method uses a factorized integral representation of the Green function which is valid asymptotically to infinity in space.

Misner, C. W.↗
NASA NTRS1982

Electromagnetic cyclotron-loss-cone instability associated with weakly relativistic electrons

The amplification of fast extraordinary mode waves at frequencies very close to the electron cyclotron frequency, due to the presence of a population of energetic electrons with a loss-cone type distribution, is studied. Low-energy background electrons are included in the analysis. Two types of loss-cone distribution functions are considered, and it is found that the maximum growth rates for both distribution functions are of the same order of magnitude. When the thermal effects of the energetic electrons are included in the dispersion equation, the real frequencies of the waves are lower than those obtained by using the cold plasma approximation. This effect tends to enhance the growth rate. An idealized case including a parallel electric field such that the distribution function of the trapped energetic electrons is modified is also considered. It is assumed that the parallel electric field can remove the low-energy background electrons away from the source region of radiation. Both these effects increase the growth rate.

Wong, H. K.↗
NASA NTRS2013

Modeling Relativistic Electron Precipitation Bremsstrahlung X-Ray Intensities at 10-100 km Manned Vehicle Altitudes

Relativisitic electron precipitation (REP) events occur when beams or bunches of relativistic electrons of magnetospheric origin enter the Earth's atmosphere, typically at auroral latitudes. REP events are associated with a variety of space weather effects, including production of transitional and bremsstrahlung radiation, catalytic depletion of stratospheric ozone, and scintillation of transionospheric radio waves. This study examines the intensities of x‐rays produced at airliner, manned balloon, and space reuseable launch vehicles (sRLVs). The monoenergetic beam is modeled in cylindrical symetry using the paraxial ray equation. Bremsstrahlung photon production is calculated using the traditional Sauter‐Elwert cross‐section, providing x‐ray emission spectra differential in energy and angle. Attenuation is computed for a plane‐stratified standard atmosphere, and the loss processes include photoionization, Rayleigh and Compton scattering, electron‐positron pair production, and photonuclear interaction. Peak altitudes of electron energy deposition and bremsstrahlung x‐ray production were calculated for beams of energies from 1 MeV through 100 MeV. The altitude peak of bremsstrahlung deposition was consistently and significantly lower that that of the electron deposition due to the longer mean free paths of x‐rays compared to electrons within the atmosphere. For example, for a nadir‐directed monoenergetic 5 MeV beam, the peak deposition altitude was calculated to be 42 km, but the resulting bremsstrahlung deposition peaked at 25 km. This has implications for crew and passenger safety, especially with the growth of the space tourism industry. A survey of results covering the 1‐100 MeV spectrum for the three altitude ranges of interest will be presented.

Krause, L. Habsh↗
NASA NTRS2021

HARM3D+NUC: A New Method for Simulating the Post-merger Phase of Binary Neutron Star Mergers with GRMHD, Tabulated EOS, and Neutrino Leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short γ-ray burst. In order to understand the radioactively powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino-leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

Ariadna Murguia-Berthier↗
NASA NTRS1973

The relativistic dynamics of a point charge in the field of a plane electromagnetic wave traveling in the direction of a uniform static magnetic field

The motion of a charged particle in electromagnetic fields of various geometric configurations and arising from a variety of sources is of intrinsic interest in electromagnetic theory. The particular configuration consisting of a plane wave propagating in the presence of a static uniform magnetic field whose direction is parallel to the wave normal is examined. The analysis presented here is treated within the context of classical electromagnetic theory. A numerical solution - at least to the approximate Lorentz-Dirac equation - is obtained.

Mitchell, T. P.↗
NASA NTRS1985

Diffusion of strongly magnetized cosmic ray particles in a turbulent medium

Cosmic ray (CR) propagation in a turbulent medium is usually considered in the diffusion approximation. Here, the diffusion equation is obtained for strongly magnetized particles in the general form. The influence of a large-scale random magnetic field on CR propagation in interstellar medium is discussed. Cosmic rays are assumed to propagate in a medium with a regular field H and an ensemble of random MHD waves. The energy density of waves on scales smaller than the free path 1 of CR particles is small. The collision integral of the general form which describes interaction between relativistic particles and waves in the quasilinear approximation is used.

Ptuskin, V. S.↗
NASA NTRS1977

Sakata-Taketani spin-0 theory with external field interactions - Lagrangian formalism and causal properties

To illustrate that a relativistic field theory need not be manifestly covariant, Lorentz-invariant Lagrangian densities are constructed that yield the equation satisfied by an interacting (two-component) Sakata-Taketani spin-0 field. Six types of external field couplings are considered, two scalars, two vectors, an antisymmetric second-rank tensor, and a symmetric second-rank tensor, with the results specialized to electromagnetic interactions. For either of the two second-rank couplings, the equation is found to describe noncausal wave propagation, a property that is apparent from the dependence of the coefficients of the space derivatives on the external field; in contrast, the noncausality of the corresponding manifestly covariant Duffin-Kemmer-Petiau spin-0 equation is not so obvious. The possibilities for generalizing the results to higher spin theories involving only the essential 2(2J + 1) components for a particle with a definite spin J and mass m are discussed in considerable detail.

Guertin, R. F.↗
DOE OSTI2024

O-Mode Instability for Non-Thermal Continuum Radiation Generated Near the Plasmapause

Analysis of the ordinary mode (O-mode) instability is performed to comprehend the nonthermal continuum (NTC) radiation near the plasmapause, taking into account the relativistic wave-electron resonance effect. The energy source is the anisotropy in the velocity of the minority suprathermal electron population. Numerical solutions demonstrate that the O-mode can be unstable with multiple narrow frequency bands located close to harmonics of the electron cyclotron frequency above the local electron plasma frequency. These waves have narrow beaming angle bands of nearly 90° relative to the ambient magnetic field. Our findings indicate that NTC radiation generated by this wave-electron resonance instability near the plasmapause can propagate nearer to the magnetic equator with multiple harmonics, which is in agreement with a recent statistical study using Van Allen Probes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
NASA NTRS2005

Ring Current-Electromagnetic Ion Cyclotron Waves Coupling

The effect of Electromagnetic Ion Cyclotron (EMIC) waves, generated by ion temperature anisotropy in Earth s ring current (RC), is the best known example of wave- particle interaction in the magnetosphere. Also, there is much controversy over the importance of EMIC waves on RC depletion. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt (RB) by EMIC wave scattering during a magnetic storm. That is why the calculation of EMIC waves must be a very critical part of the space weather studies. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch-angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) the EMIC waves saturation due to their modulation instability and LHW generation are included as an additional factor that contributes to this process; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗