Search NASA⌕ Search

SEARCH · Search NASA

Results for “interstellar magnetic fields”

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

Global Characteristics of the IBEX Ribbon and Its Temporal Variation: Implications for Relationships with Solar Cycle Variability and the Interstellar Magnetic Field

The Interstellar Boundary Explorer (IBEX) has revealed the existence of the energetic neutral atom (ENA) ribbon, a narrow and nearly circular band of enhanced flux that enables probing plasma and magnetic field conditions outside the heliopause (HP). In this study, we investigate the global structure and temporal variation of the IBEX ribbon using 14 yr of IBEX-Hi ENA observations. We derive geometric parameters of the ribbon: center location, angular radius, eccentricity, and major axis orientation. We find that the ribbon exhibits systematic variations tied to the solar cycle. The ribbon center shows a solar cycle–modulated latitudinal motion at low energies. The angular radius increases with ENA energy and displays a pronounced enlargement between 2013 and 2018, coinciding with changes in high-latitude fast solar wind from polar coronal holes. The eccentricity remains consistently around 0.3 with some temporal evolution, while the major axis rotates gradually over time, becoming increasingly perpendicular to the ecliptic meridional plane. A notable hemispheric asymmetry has emerged over time: the ribbon in the northern hemisphere remains remarkably stable in position, whereas the southern ribbon exhibits temporal migration. By combining these results with solar wind observations, global heliosphere, and ribbon modeling, we suggest that this asymmetry is explained by a coupled effect of magnetic field draping and the latitude-dependent timing of fast solar wind from polar coronal holes. These results establish the most detailed long-term geometric characterization of the IBEX ribbon to date and provide observational support for models of secondary ENA production and magnetic field draping outside the HP.

Noh, Sung Jun [New Mexico Consortium, Los Alamos, ↗

Whence the Interstellar Magnetic Field Shaping the Heliosphere?

Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B magpol , aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for "downfield" stars located over 90° from Bmagpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, "Filament A," extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

79 ASTRONOMY AND ASTROPHYSICS↗

Toward Interpreting the IBEX Ribbon with Mirror Diffusion in Interstellar Turbulent Magnetic Fields

We investigate the role of the magnetohydrodynamic (MHD) turbulence measured by Voyager in the very local interstellar medium (VLISM) in modeling the Interstellar Boundary Explorer ribbon. We demonstrate that the mirroring by compressible modes of MHD turbulence dominates over that by the mean magnetic field. Based on the new mirror diffusion mechanism identified by Lazarian & Xu for particles with large pitch angles in MHD turbulence, we find that the mirror diffusion can both confine pickup ions and preserve their initial pitch angles, and thus it accounts for the enhanced intensity of energetic neutral atoms that return to the heliosphere. The ribbon width is determined by both the range of pitch angles for effective turbulent mirroring and the field line wandering induced by Alfvénic modes. It in turn provides a constraint on the amplitude of magnetic fluctuations of fast modes. The field line wandering also affects the coherence of the ribbon structure across the sky. By extrapolating the magnetic energy spectrum measured by Voyager, we find that the injection scale of the turbulence in the VLISM must be less than ~500 au for the ribbon structure to be coherent.

79 ASTRONOMY AND ASTROPHYSICS↗

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Misalignment of Interstellar Dust Filaments

We present evidence for scale-independent misalignment of interstellar dust filaments and magnetic fields. We estimate the misalignment by comparing millimeter-wave dust-polarization measurements from Planck with filamentary structures identified in neutral-hydrogen (H I ) measurements from H I 4PI. We find that the misalignment angle displays a scale independence (harmonic coherence) for features larger than the H I 4PI beamwidth (16$\acute{.}$2). We additionally find a spatial coherence on angular scales of ${ \mathcal O }(1^\circ )$. We present several misalignment estimators formed from the auto- and cross-spectra of dust-polarization and H I -based maps, and we also introduce a map-space estimator. Applied to large regions of the high-Galactic-latitude sky, we find a global misalignment angle of ~2°, which is robust to a variety of masking choices. By dividing the sky into small regions, we show that the misalignment angle correlates with the parity-violating TB cross-spectrum measured in the Planck dust maps. The misalignment paradigm also predicts a dust EB signal, which is of relevance in the search for cosmic birefringence but as yet undetected; the measurements of EB are noisier than those of TB, and our correlations of EB with misalignment angle are found to be weaker and less robust to masking choices. We also introduce an H I -based dust-polarization template constructed from the Hessian matrix of the H I intensity, which is found to correlate more strongly than previous templates with Planck dust B modes.

79 ASTRONOMY AND ASTROPHYSICS↗

A Foreshock Model for Interstellar Shocks of Solar Origin: Voyager 1 and 2 Observations

The Voyager 1 (V1) and Voyager 2 (V2) spacecraft were launched in 1977 on a mission to explore the outer planets and reach the heliopause, the boundary between the hot solar plasma and the relatively cool interstellar plasma. V1 reached the heliopause on 2012 August 25, at 121.6 au, and V2 reached the heliopause on 2018 November 5, at 119.0 au. One of their remarkable discoveries was the detection of shocks propagating into the interstellar plasma from energetic solar events. These shocks are typically preceded by electron plasma oscillations excited by electron beams streaming along interstellar magnetic field lines ahead of the shocks. The frequencies of the plasma oscillations have now provided radial electron density profiles in the outer heliosphere and in the interstellar medium to radial distances of more than 145 au. The oscillations are typically preceded by bursts of high-energy ∼5–100 MeV electrons. These electron bursts are interpreted as being due to the reflection (and acceleration) of cosmic-ray electrons by the shock at the time the shock first contacts the magnetic field line that passes through the spacecraft. Relative timing between the cosmic rays reflected by the shock and the onset of the plasma oscillations allow us, for the first time, to estimate the energy, ∼20–100 eV, of the electron beams responsible for the plasma oscillations. These observations are combined into a self-consistent model called the foreshock model that describes the interaction of shocks of solar origin with the interstellar plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Filamentary Dust Polarization and the Morphology of Neutral Hydrogen Structures

Filamentary structures in neutral hydrogen (H$\tiny{I}$) emission are well aligned with the interstellar magnetic field, so H$\tiny{I}$ emission morphology can be used to construct templates that strongly correlate with measurements of polarized thermal dust emission. We explore how the quantification of filament morphology affects this correlation. We introduce a new implementation of the Rolling Hough Transform (RHT) using spherical harmonic convolutions, which enables efficient quantification of filamentary structure on the sphere. We use this Spherical RHT algorithm along with a Hessian-based method to construct H$\tiny{I}$-based polarization templates. We discuss improvements to each algorithm relative to similar implementations in the literature and compare their outputs. By exploring the parameter space of filament morphologies with the Spherical RHT, we find that the most informative H$\tiny{I}$ structures for modeling the magnetic field structure are the thinnest resolved filaments. For this reason, we find a ~10% enhancement in the B-mode correlation with polarized dust emission with higher-resolution H$\tiny{I}$ observations. We demonstrate that certain interstellar morphologies can produce parity-violating signatures, i.e., nonzero TB and EB, even under the assumption that filaments are locally aligned with the magnetic field. Finally, we demonstrate that B modes from interstellar dust filaments are mostly affected by the topology of the filaments with respect to one another and their relative polarized intensities, whereas E modes are mostly sensitive to the shapes of individual filaments.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on the IBEX Ribbon’s Origin from Its Evolution over a Solar Cycle

Abstract In 2009, the Interstellar Boundary Explorer (IBEX) discovered a narrow “ribbon” of energetic neutral atom emissions across the sky with properties correlated with the solar wind latitudinal structure and the interstellar magnetic field draped around the heliosphere. It is widely believed that the ribbon is formed from the escape of heliospheric ENAs into the local interstellar medium and their eventual return as secondary ENAs. However, there is no consensus on the rate of pitch angle scattering of these PUIs before they become secondary ENAs. We test two opposing limits of scattering rates (“weak” versus “strong”) by solving a time-dependent model of the ribbon that evolves with the solar cycle, and we compare them to IBEX observations over 2009–2019. First, we find that both models qualitatively reproduce the evolution of IBEX fluxes for most of the data set, with a few exceptions, although the strong (or “spatial retention”) scattering model greatly underestimates the observed fluxes. Regardless, time dependence of fluxes cannot distinguish these models. Second, the ribbon’s geometric properties, i.e., its center and radius, are significantly different between the models. The spatial retention model reproduces the observed ribbon centers as a function of energy and time slightly better than the weak scattering model, and the spatial retention model reproduces the observed ribbon radius over energy and time almost perfectly, whereas the weak scattering model compares poorly. Our analysis favors the spatial retention mechanism as the source of the IBEX ribbon, but it requires modification to increase the flux of ENAs observed at 1 au.

Zirnstein, E. J. (ORCID:0000000172400618)↗

Evolution of the Heliotail Lobes over a Solar Cycle as Measured by IBEX

Energetic neutral atom (ENA) measurements by IBEX reveal that the heliotail comprises an energy-dependent multilobe structure. We examine the heliotail evolution over 11 yr of IBEX observations covering a full solar cycle (SC). We find the following: (1) The heliotail structure persists over the entire SC, comprising three ENA-enhanced and two ENA-suppressed lobes. (2) Lobe sizes and locations are generally stable but exhibit variations in ENA fluxes driven by the SC. (3) Lobe centers follow a cyclic behavior over multiple SC phases, indicating direct signatures of slow and fast solar wind (SW) interactions in the inner heliosheath (IHS). (4) The tilted plane passing through the port–starboard lobes' centers oscillates in latitude but maintains its tilt from the ecliptic plane, likely a consequence of the interstellar magnetic field draping around the heliosphere. (5) The transition of the central heliotail from a single lobe at ~1.1 keV to two lobes above ~2 keV is SC-dependent and directly reflects the IHS plasma properties, i.e., when ENA fluxes from fast SW from the polar coronal holes change over time. (6) The central lobe exhibits a substructure that is enhanced and offset from the downwind direction, possibly indicating an asymmetric ENA emission or an asymmetry in the parent plasma distribution. These results reveal the general stability of the heliotail structure over time and distinct variations in individual lobes' properties in relation to the SC phases. Furthermore, results show the effects of multiple SC phases in the tail, reflecting different ENA travel times and source histories.

79 ASTRONOMY AND ASTROPHYSICS↗

The Davis–Chandrasekhar–Fermi method revisited

ABSTRACT Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis–Chandrasekhar–Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method.

79 ASTRONOMY AND ASTROPHYSICS↗

Observation of high-energy neutrinos from the Galactic plane

The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth’s atmosphere, is unknown. Because of deflection by interstellar magnetic fields, cosmic rays produced within the Milky Way arrive at Earth from random directions. However, cosmic rays interact with matter near their sources and during propagation, which produces high-energy neutrinos. We searched for neutrino emission using machine learning techniques applied to 10 years of data from the IceCube Neutrino Observatory. By comparing diffuse emission models to a background-only hypothesis, we identified neutrino emission from the Galactic plane at the 4.5σ level of significance. The signal is consistent with diffuse emission of neutrinos from the Milky Way but could also arise from a population of unresolved point sources.

Science & Technology - Other Topics↗

Diagnostics of Magnetohydrodynamic Modes in the Interstellar Medium through Synchrotron Polarization Statistics

One of the biggest challenges in understanding magnetohydrodynamic (MHD) turbulence is identifying the plasma mode components from observational data. Previous studies on synchrotron polarization from the interstellar medium (ISM) suggest that the dominant MHD modes can be identified via statistics of Stokes parameters, which would be crucial for studying various ISM processes such as the scattering and acceleration of cosmic rays, star formation, and dynamo. In this paper, we present a numerical study of the synchrotron polarization analysis (SPA) method through systematic investigation of the statistical properties of the Stokes parameters. We derive the theoretical basis for our method from the fundamental statistics of MHD turbulence, recognizing that the projection of the MHD modes allows us to identify the modes dominating the energy fraction from synchrotron observations. Based on the discovery, we revise the SPA method using synthetic synchrotron polarization observations obtained from 3D ideal MHD simulations with a wide range of plasma parameters and driving mechanisms, and present a modified recipe for mode identification. We propose a classification criterion based on a new SPA+ fitting procedure, which allows us to distinguish between Alfvén mode and compressible/slow mode dominated turbulence. We further propose a new method to identify fast modes by analyzing the asymmetry of the SPA+ signature and establish a new asymmetry parameter to detect the presence of fast mode turbulence. Additionally, we confirm through numerical tests that the identification of the compressible and fast modes is not affected by Faraday rotation in both the emitting plasma and the foreground.

97 MATHEMATICS AND COMPUTING↗

Spectral Properties of Globally Distributed ENA Fluxes across Diverse Regions of the Heliosphere

This study analyzes energetic neutral atom (ENA) spectral properties across distinct regions of globally distributed flux (GDF) sky maps, using Interstellar Boundary Explorer data from a full solar cycle, corrected for time dispersion. By time-shifting the data to the heliosheath using GDF source distances from D. B. Reisenfeld et al., we achieve a more accurate representation of heliosheath GDF energy spectra. We quantify ENA spectral characteristics, heliosheath line-of-sight-integrated proton pressure, and heliosheath proton temperature, comparing these to solar wind properties at 1 au and interplanetary scintillation-derived solar wind data. Our findings show that the spectral index is generally anticorrelated with heliosheath proton temperature and pressure, except in the central tail, where a partial positive correlation is observed. The lowest spectral index values occur when high-latitude heliosheath regions are dominated by fast solar wind from polar coronal holes. The south pole exhibits the flattest energy spectra due to plasma heating from both fast solar wind and a late-2014 pressure pulse. The central tail shows shorter variability (5–6 yr) for spectral index and heliosheath proton temperature, while proton pressure follows the 11 yr solar cycle. Most spectral shapes exhibit a “knee” distribution, peaking during solar maximum, with an “ankle” shape observed only at the south pole during solar cycle transitions. Asymmetry in proton pressure in the lobes is driven by the draping effect of the local interstellar magnetic field. This study provides insights into the energetic properties of GDF across the heliosphere, enhancing our understanding of the heliospheric environment.

79 ASTRONOMY AND ASTROPHYSICS↗

Gradient Technique Theory: Tracing Magnetic Field and Obtaining Magnetic Field Strength

Abstract The gradient technique is a promising tool with theoretical foundations based on the fundamental properties of MHD turbulence and turbulent reconnection. Its various incarnations use spectroscopic, synchrotron, and intensity data to trace the magnetic field and measure the media magnetization in terms of Alfvén Mach number. We provide an analytical theory of gradient measurements and quantify the effects of averaging gradients along the line of sight and over the plane of the sky. We derive analytical expressions that relate the properties of gradient distribution with the Alfvén Mach number M A . We show that these measurements can be combined with measures of sonic Mach number or line broadening to obtain the magnetic field strength. The corresponding technique has advantages to the Davis–Chandrasekhar–Fermi way of obtaining the magnetic field strength.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Field Strength from Turbulence Theory. I. Using Differential Measure Approach

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

79 ASTRONOMY AND ASTROPHYSICS↗

Observation of Cosmic-Ray Anisotropy in the Southern Hemisphere with 12 yr of Data Collected by the IceCube Neutrino Observatory

We analyzed the 7.92 × 10 11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

79 ASTRONOMY AND ASTROPHYSICS↗

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. Here, in this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, δB , and large-scale, B 0 , fields, $\langle$(δB·B 0 ) 2 $\rangle$$_{V}^{1/2}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. Furthermore, we explore the relation between the turbulent mean field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetized ISM.

79 ASTRONOMY AND ASTROPHYSICS↗

Diagnosis of 3D magnetic field and mode composition in MHD turbulence with Y -parameter

ABSTRACT Magnetic fields are crucial in numerous astrophysical processes within the interstellar medium (ISM). However, the detailed determination of magnetic field geometry is notoriously challenging. Based on the modern magnetohydrodynamic (MHD) turbulence theory, we introduce a novel statistical technique, the ‘Y-parameter’, to decipher the magnetic field inclination in the ISM and identify dominant turbulence modes. The Y-parameter, calculated as the ratio of anisotropies of different Stokes parameter combinations, displays contrasting trends with the mean-field inclination angle in Alfvénic and compressible turbulence modes. A Y-parameter value around 1.5 ± 0.5 provide a statistical boundary to determine the dominant MHD turbulence modes. We have discovered specific correlations between the Y-parameter value and the inclination angle that unveil the dominant turbulence mode. This methodology, when applied to future radio polarization surveys such as LOFAR and SKA, promises to significantly enhance our knowledge of 3D magnetic field in the ISM and improve our understanding of interstellar turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗