SEARCH · Search NASA
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.
This report summarizes our most important recent re search, building upon several years of productive research. This resulted in breakthroughs in our understanding of plasma confined by the strong field of a magnetic dipole, dissertations from several graduate students, and learning experiences from our undergraduate students.
Missions like NASA's Imaging X-ray Polarimetry Explorer (IXPE) are poised to provide an unprecedented view of the Universe in polarized X-rays. Polarization probes physical anisotropies, a fact exploited by particle physicists to look for the anisotropic a E∙B operator in the axion-like-particle (ALP) Lagrangian. Such studies have typically focused on polarization in the radio and microwaves, through local or cosmic birefringence effects. To such polarization studies we add X-rays emanating from magnetars — a class of neutron stars with near-critical strength magnetic fields — that are important targets for IXPE. ALPs produced in the neutron star core convert to X-rays in the magnetosphere; such X-rays are polarized along the direction parallel to the dipolar magnetic field at the point of conversion. Here, we develop the full theoretical formalism for ALP-induced polarization in the presence of dipolar magnetic fields. For uncorrelated photon and ALP production mechanisms, we completely disentangle the ALP contributions to the Stokes parameters in terms of the ALP intensity, the ALP-to-photon conversion probability, and the ALP-induced birefringence. In the proper limit, our results demonstrate that the inclusion of ALPs suppresses the observed degree of circular polarization compared to its pure astrophysical value. Our results can also be used to impose limits on ALP couplings with IXPE polarization data from magnetars 4U 0142+61 and 1RXS J170849.0-400910, the subject of upcoming work.
We report on Chandra X-ray Observatory (CXO) observations of the pulsar wind nebula (PWN) associated with PSR B0355+54 (eight observations with a 395 ks total exposure, performed over an eight month period). We investigated the spatial and spectral properties of the emission coincident with the pulsar, compact nebula (CN), and extended tail. We find that the CN morphology can be interpreted in a way that suggests a small angle between the pulsar spin axis and our line of sight, as inferred from the radio data. On larger scales, emission from the 7{sup ′} (≈2 pc) tail is clearly seen. We also found hints of two faint extensions nearly orthogonal to the direction of the pulsar’s proper motion. The spectrum extracted at the pulsar position can be described with an absorbed power-law + blackbody model. The nonthermal component can be attributed to magnetospheric emission, while the thermal component can be attributed to emission from either a hot spot (e.g., a polar cap) or the entire neutron star surface. Surprisingly, the spectrum of the tail shows only a slight hint of cooling with increasing distance from the pulsar. This implies either a low magnetic field with fast flow speed, or particle reacceleration within the tail. We estimate physical properties of the PWN and compare the morphologies of the CN and the extended tail with those of other bow shock PWNe observed with long CXO exposures.
At the Earth's magnetopause that separates the hot-tenuous magnetospheric plasma from the cold dense solar wind plasma, often seen is a boundary layer where plasmas of both origins coexist. Plasma diffusions of various forms have been considered as the cause of this plasma mixing. Here, we investigate the plasma transport induced by wave-particle interaction in kinetic Alfven wave (KAW) turbulence, which is one of the candidate processes. We clarify that the physical origin of the KAW-induced cross-field diffusion is the drift motions of those particles that are in Cerenkov resonance with the wave: E Multiplication-Sign B-like drift that emerges in the presence of non-zero parallel electric field component and grad-B drift due to compressional magnetic fluctuations. We find that KAW turbulence, which has a spectral breakpoint at which an MHD inertial range transits to a dissipation range, causes selective transport for particles whose parallel velocities are specified by the local Alfven velocity and the parallel phase velocity at the spectral breakpoint. This finding leads us to propose a new data analysis method for identifying whether or not a mixed plasma in the boundary layer is a consequence of KAW-induced transport across the magnetopause. The method refers to the velocity space distribution function data obtained by a spacecraft that performs in situ observations and, in principle, is applicable to currently available dataset such as that provided by the NASA's THEMIS mission.
We compute the patterns of γ-ray emission due to curvature radiation in dissipative pulsar magnetospheres. Our ultimate goal is to construct macrophysical models that are able to reproduce the observed γ-ray light curve phenomenology recently published in the Second Fermi Pulsar Catalog. We apply specific forms of Ohm's law on the open field lines using a broad range for the macroscopic conductivity values that result in solutions ranging, from near-vacuum to near-force-free. Using these solutions, we generate model γ-ray light curves by calculating realistic trajectories and Lorentz factors of radiating particles under the influence of both the accelerating electric fields and curvature radiation reaction. We further constrain our models using the observed dependence of the phase lags between the radio and γ-ray emission on the γ-ray peak separation. We perform a statistical comparison of our model radio-lag versus peak-separation diagram and the one obtained for the Fermi standard pulsars. We find that for models of uniform conductivity over the entire open magnetic field line region, agreement with observations favors higher values of this parameter. We find, however, significant improvement in fitting the data with models that employ a hybrid form of conductivity, specifically, infinite conductivity interior to the light cylinder and high but finite conductivity on the outside. In these models the γ-ray emission is produced in regions near the equatorial current sheet but modulated by the local physical properties. These models have radio lags near the observed values and statistically best reproduce the observed light curve phenomenology. Additionally, they also produce GeV photon cut-off energies.
The Magnetic Reconnection Experiment (MRX)[1] has been utilized to assess the effectiveness of minimum variance analysis on the magnetic field (MVAB) and boundary-crossing time analysis (BCTA). The neutral sheet is swept, or jogged, in a controlled manner with respect to the stationary probes by pulsed internal coil currents. Magnetic field data from measurement points resembling data from multi-spacecraft flying though a reconnecting current sheet is used to check both techniques to deduce a proper normal vector. We examine discharges with the two-dimensional (2-D) X-line structure as well as cases in which a flux rope forms within the layer. All discharges are in a two-fluid regime in which electrons are magnetized but not ions. Boundary-crossing time analysis with four sample measurement points forming a tetrahedron generates a reasonable unit normal vector and relative velocity along the normal vector for all of the tested cases. On the other hand, MVAB sometimes fails to predict a proper normal direction. This is because the X-line magnetic geometry is fundamentally 2-D or 3-D. However, the direction along the reconnecting field determined by MVAB does not deviate much from the real magnetic geometry documented by 2-D magnetic probe arrays and one additional probe at a different toroidal location. Based on these observations, we suggest a procedure for determining a local coordinate system for data from the Magnetospheric Multi-Scale (MMS) mission when spacecraft passes through a reconnecting current sheet. The distance between measurement points on the order of the ion skin depth ($c/\omega_{pi}$) is pertinent to determination of the magnetic geometry.
The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.
Supersonic flows with high Mach number are ubiquitous in astrophysics. High-powered lasers also have the ability to drive high Mach number, radiating shock waves in laboratory plasmas, and recent experiments along these lines have made it possible to recreate analogs of high Mach-number astrophysical flows under controlled conditions. Streak cameras such as the Rochester optical streak system (ROSS) are particularly helpful in diagnosing such experiments, because they acquire spatially resolved measurements of the radiating gas continuously over a large time interval, making it easy to observe how any shock waves and ablation fronts present in the system evolve with time. This paper summarizes new ROSS observations of a laboratory analog of the collision of a stellar wind with an ablating planetary atmosphere embedded within a magnetosphere. We find good agreement between the observed ROSS data and numerical models obtained with the FLASH code, but only when the effects of optical depth are properly taken into account.
The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3-D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence-induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multifluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. Herein we focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model is found to have difficulty resolving the lower hybrid drift instability, which plays a fundamental role in heating and mixing electrons in the current layer.
Here, we present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the in situ measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on 26 July 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from subelectron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.
The coefficient of reflection of a fast magnetosonic wave incident on the magnetosphere from the solar wind is studied analytically in the framework of a plane-stratified model of the medium with allowance for the transverse inhomogeneity of the magnetosphere and a jump of the plasma parameters at the magnetopause. Three factors decisively affecting the properties of reflection are taken into account: the shear flow of the solar wind plasma relative to the magnetosphere; the presence of a magnetospheric magnetohydrodynamic waveguide caused by the transverse plasma inhomogeneity; and the presence of an Alfvén resonance deep in the magnetosphere, where the oscillation energy dissipates. If the solar wind velocity exceeds the wave phase velocity along the magnetopause, then the wave energy in the solar wind is negative and such a wave experiences overreflection. In the opposite case, the wave energy is positive and the wave is reflected only partially. The wave reflection has a pronounced resonant character: the reflection coefficient has deep narrow minima or high narrow maxima at the eigenfrequencies of the magnetospheric waveguide. For other frequencies, the reflection coefficient only slightly differs from unity. The wave energy influx into the magnetosphere is positive for waves with both positive and negative energies. For waves with a negative energy, this is a consequence of their overreflection, because the flux of negative energy carried away by the reflected wave exceeds the incident flux of negative energy.
We present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the in situ measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on 26 July 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from subelectron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.
We examined the plasma flow response to meso- and microscale magnetic dipoles by performing three-dimensional full particle-in-cell simulations. We particularly focused on the formation of a magnetosphere and its dependence on the intensity of the magnetic moment. The size of a magnetic dipole immersed in a plasma flow can be characterized by a distance L from the dipole center to the position where the pressure of the local magnetic field becomes equal to the dynamic pressure of the plasma flow under the magnetohydrodynamics (MHD) approximation. In this study, we are interested in a magnetic dipole whose L is smaller than the Larmor radius of ions r{sub iL} calculated with the unperturbed dipole field at the distance L from the center. In the simulation results, we confirmed the clear formation of a magnetosphere consisting of a magnetopause and a tail region in the density profile, although the spatial scale is much smaller than the MHD scale. One of the important findings in this study is that the spatial profiles of the plasma density as well as the current flows are remarkably affected by the finite Larmor radius effect of the plasma flow, which is different from the Earth's magnetosphere. The magnetopause found in the upstream region is located at a position much closer to the dipole center than L. In the equatorial plane, we also found an asymmetric density profile with respect to the plasma flow direction, which is caused by plasma gyration in the dipole field region. The ion current layers are created in the inner region of the dipole field, and the electron current also flows in the region beyond the ion current layer because ions with a large inertia can closely approach the dipole center. Unlike the ring current structure of the Earth's magnetosphere, the current layers in the microscale dipole fields are not circularly closed around the dipole center. Since the major current is caused by the particle gyrations, the current is independently determined to be in the direction of the electron and ion gyrations, which are the same in both the upstream and downstream regions. The present analysis on the formation of a magnetosphere in the regime of a microscale magnetic dipole is significant for understanding the solar wind response to the crustal magnetic anomalies on the Moon surface, such as were recently observed by spacecraft.
A spacecraft-charging mitigation scheme necessary for the operation of a high-power electron beam in the low-density magnetosphere is analyzed. The scheme is based on a plasma contactor, i.e. a high-density charge-neutral plasma emitted prior to and during beam emission, and its ability to emit high ion currents without strong space-charge limitations. A simple theoretical model for the transient of the spacecraft potential and contactor expansion during beam emission is presented. The model focuses on the contactor ion dynamics and is valid in the limit when the ion contactor current is equal to the beam current. The model is found in very good agreement with Particle-In-Cell simulations over a large parametric study that varies the initial expansion time of the contactor, the contactor current and the ion mass. The model highlights the physics of the spacecraft-charging mitigation scheme, indicating that the most important part of the dynamics is the evolution of the outermost ion front which is pushed away by the charge accumulated in the system by the beam. The model can be also used to estimate the long-time evolution of the spacecraft potential. For a short contactor expansion (0.3 or 0.6 ms Helium plasma or 0.8 ms Argon plasma, both with 1 mA current) it yields a peak spacecraft potential of the order of 1-3 kV. This implies that a 1-mA relativistic electron beam would be easily emitted by the spacecraft.
Parallel propagating electromagnetic electron cyclotron (EMEC) waves in the extended plasma sheet (~12R{sub S}) and in the outer magnetosphere (~18R{sub S}) of Saturn have been studied. A dispersion relation for parallel propagating relativistic EMEC waves has been applied to the magnetosphere of Saturn, and comparisons have been made with the data of Voyager 1 at these radial distances. The detailed investigations for EMEC waves have been done in the presence of the perpendicular AC electric field, using the kappa distribution function. The relativistic temporal growth rate is calculated by the method of characteristic solution with the data provided by Voyager 1. The effect of the suprathermal electron density, temperature anisotropy, frequency of AC electric field, thermal energy of ions, and relativistic factor on the temporal growth rate of EMEC wave emission has been studied. The simulation results show that the growth of parallel propagating EMEC waves is significantly affected by variations in the temperature anisotropy, electron density, ion thermal energy, and relativistic factor in both the extended plasma sheet and the outer magnetosphere of Saturn. The temperature anisotropy (T⊥/T{sub ║}), ion thermal energy (K{sub B}T{sub ║i}), and electron density (n{sub 0}) have been found to be a major source of free energy for parallel propagating EMEC waves in both regions.
It is hypothesized that dark matter is composed of particles called quark nuggets, and further that these particles have a permanent magnetic dipole moment. If this hypothesis is true, calculations predict that a magnetized quark nugget (MQN) will oscillate when encountering the Earth's magnetosphere, and emit RF radiation between 30kHz and 30MHz. To support testing this hypothesis, a loop antenna sensor was designed and developed, which is described in this report. This sensor operates between 300kHz and 3MHz and achieves about -11dBfT/vHz sensitivity at 1.5MHz.
Blobs are filament-like structures within a bulk plasma. They move radially outward in curved magnetic fields, across magnetic field lines. They account for significant amount of transport at plasma edge in tokamaks and are found in tokamaks, solar coronal filaments, planetary magnetospheres, and high altitude nuclear explosions. MHD offers some analytical descriptions of blob dynamics. In a curved geometry such as our cases of interest charged particles are subject to curvature drift, grad B drift, and E x B drifts.
Hybrid Vlasov-Fourier modeling is used to investigate the nonlinear evolution of electromagnetic ion cyclotron (EMIC) waves driven by proton temperature anisotropy in plasmas with a population of He{sup +} ions and a cold proton background. In the pure proton-electron plasma, most of the free energy is converted into high-amplitude waves and currents. In the nonlinear stage, within a few hundred proton gyroperiods after the saturation, the wave spectrum shifts toward lower wave numbers and frequencies, from {omega}{approx}0.6{Omega}{sub p} to below {omega}{approx}0.25{Omega}{sub p}. In the presence of even a small population of He{sup +} ions almost all of the free energy is used in He{sup +} heating. The wave activity in the saturated state moves from the linearly unstable upper branch to the linearly stable lower one. In the presence of a background of cold protons, the waves can propagate in the frequency stop-band. Our results demonstrate that linear stability theory cannot be used to estimate the characteristics of the expected saturated wave spectra in the terrestrial magnetosphere. Significantly, our nonlinear simulations produce wave spectra which are in close agreement with the EMIC waves observed in situ by satellites as well as by ground-based magnetometers positioned at the ends of the magnetic field lines.