Search

We analyze the count rates of 40–139 keV ions that were measured in situ by the Low Energy Charged Particle instrument on Voyager 1 in order to identify the suprathermal ion anisotropies beyond the termination shock and in the very local interstellar (IS) medium (VLISM). The analysis results in a region of ∼9–10 au before the heliopause (HP) where the radial anisotropy of ions is negative, while the azimuthal ion anisotropy inside the heliosheath lies in the − T direction. In agreement to our previous analyses, we identify a positive radial anisotropy of ions up to at least ∼30 au beyond the HP, which becomes nearly zero from 2021 up to 2023 November (for ∼10 au). Notably, the anisotropy in the azimuthal direction is statistically zero throughout the upstream region, i.e., for ∼40 au past the HP, showing that the 40–139 keV ion anisotropy in the VLISM is only in the radial direction and has no azimuthal component. The presence of suprathermal ions of solar origin over such long spatial scales in the VLISM, along with the inflow of ions from IS space into the heliosheath, are important constraints for characterizing the interaction of the heliosheath with the VLISM. Our observations provide indications that V1 has entered a new regime in the VLISM since (at least) the year 2021, progressively developing characteristics akin to the pristine IS medium. Alternatively, this drop to nearly zero radial anisotropies beyond 2021 could be a manifestation of a prolonged compression/shock of solar origin. [ABSTRACT FROM AUTHOR]

We present deep 1.4–4.8 μ m JWST-NIRCam imaging of the Serpens Main star-forming region and identify 20 candidate protostellar outflows, most with bipolar structure and identified driving sources. The outflow position angles (PAs) are strongly correlated, and they are aligned within ±24° of the major axis of the Serpens filament. These orientations are further aligned with the angular momentum vectors of the two disk shadows in this region. We estimate that the probability of this number of young stars being coaligned if sampled from a uniform PA distribution is 10−4. This in turn suggests that the aligned protostars, which seem to be at similar evolutionary stages based on their outflow dynamics, formed at similar times with a similar spin inherited from a local cloud filament. Further, there is tentative evidence for a systematic change in average PA between the northwestern and southeastern cluster, as well as increased scatter in the PAs of the southeastern protostars. SOFIA-HAWC+ archival dust polarization observations of Serpens Main at 154 and 214 μ m are perpendicular to the dominant jet orientation in the northwestern region in particular. We measure and locate shock knots and edges for all of the outflows and provide an identifying catalog. We suggest that Serpens main is a cluster that formed from an isolated filament and due to its youth retains its primordial outflow alignment. [ABSTRACT FROM AUTHOR]

We report the detection of the ordered alignment between the magnetic field and kpc-scale bubbles in the nearby spiral galaxy, NGC 628. Applying the Velocity Gradient Technique on CO spectroscopic data from the ALMA-PHANGS, the magnetic field of NGC 628 is measured at the scale of 191 pc (∼4″). The large-scale magnetic field is oriented parallel to the spiral arms and curves around the galactic bubble structures in the mid-infrared emission observed by the James Webb Space Telescope. A total of 21 bubble structures have been identified at the edges of spiral arms with scales over 300 pc, which includes two kpc-scale structures. These bubbles are caused by supernova remnants and prolonged star formation and are similar to the outflow chimneys found in neutral hydrogen in galactic disks. At the edge of the bubbles, the shocks traced by the O iii emission present a curved magnetic field that parallels the bubble's shell. The magnetic field follows the bubble expansion and binds the gas in the shell to trigger further star formation. By analyzing the larger sample of 1694 bubbles, we found a distinct radial-size distribution of bubbles in NGC 628 indicating the star formation history in the galaxy. [ABSTRACT FROM AUTHOR]

Magnetic fields are now widely recognized as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modeled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organizing magnetic fields at the scale of the disk or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at a sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr versus 1–2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due to α -quenching near the midplane as the growing mean-field produces a magnetic α that opposes the kinetic α. The magnetic energy in our models of the LSD shows a slightly sublinear response to increasing resolution, indicating that we are converging toward the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution. [ABSTRACT FROM AUTHOR]

We present a cross-correlation analysis between 1 ′ resolution total intensity and polarization observations from the Atacama Cosmology Telescope (ACT) at 150 and 220 GHz and 15″ mid-infrared photometry from the Wide-field Infrared Survey Explorer (WISE) over 107 12.°5 × 12.°5 patches of sky. We detect a spatially isotropic signal in the WISE×ACT TT cross-power spectrum at 30 σ significance that we interpret as the correlation between the cosmic infrared background at ACT frequencies and polycyclic aromatic hydrocarbon (PAH) emission from galaxies in WISE, i.e., the cosmic PAH background. Within the Milky Way, the Galactic dust TT spectra are generally well described by power laws in ℓ over the range 103 < ℓ < 104, but there is evidence both for variability in the power-law index and for non-power-law behavior in some regions. We measure a positive correlation between WISE total intensity and ACT E -mode polarization at 1000 < ℓ ≲ 6000 at >3 σ in each of 35 distinct ∼100 deg2 regions of the sky, suggesting that alignment between Galactic density structures and the local magnetic field persists to subparsec physical scales in these regions. The distribution of TE amplitudes in this ℓ range across all 107 regions is biased to positive values, while there is no evidence for such a bias in the TB spectra. This work constitutes the highest- ℓ measurements of the Galactic dust TE spectrum to date and indicates that cross-correlation with high-resolution mid-infrared measurements of dust emission is a promising tool for constraining the spatial statistics of dust emission at millimeter wavelengths. [ABSTRACT FROM AUTHOR]

Knowledge of the initial conditions of high-mass star formation is critical for theoretical models, but are not well observed. Built on our previous characterization of a Galaxy-wide sample of 463 candidate high-mass starless clumps (HMSCs), here we investigate the dynamical state of a representative subsample of 44 HMSCs (radii 0.13–1.12 pc) using Green Bank Telescope NH _3 (1,1) and (2,2) data from the Radio Ammonia Mid-Plane Survey pilot data release. By fitting the two NH _3 lines simultaneously, we obtain velocity dispersion, gas kinetic temperature, NH _3 column density and abundance, Mach number, and virial parameter. Thermodynamic analysis reveals that most HMSCs have Mach number

Determining the magnitude and direction of the interstellar magnetic field ( B _ISM ) is a long-standing problem. To date, some methods to infer the direction and magnitude have utilized best-fit models to the positions of the termination shock and heliopause measured by Voyager 1 and 2. Other models use the circularity of the Interstellar Boundary Explorer (IBEX) ribbon assuming a secondary energetic neutral atom (ENA) mechanism. Previous studies have revealed that the B _ISM organizes the orientation of the heliotail with respect to the solar meridian. Here we propose a new way to infer the direction of the B _ISM based on ENA observations of the heliotail. IBEX observations of the heliotail have revealed high-latitude lobes of enhanced ENA flux at energies >2 keV. Analyses showed that the high-latitude lobes are nearly aligned with the solar meridian, while also exhibiting a rotation with solar cycle. We show, using steady-state solar wind conditions, that the inclination of the lobes reproduced with commonly used values for the angle ( α _BV ) between B _ISM and the interstellar flow in the hydrogen deflection plane (40° < α _BV < 60°) is inconsistent with the IBEX ENA observations. We report that 0° < α _BV < 20° best replicates the heliotail lobe inclinations observed by IBEX. Additionally, our model results indicate that the variation of the solar magnetic field magnitude with solar cycle causes the longitudinal rotation of the lobes observed by IBEX by affecting the inclination of the lobes.

We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8−0.3 and G39.2−0.3 in L -band (1–2 GHz) radio continuum in The H i /OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16″ beam (≲0.7 pc) from Faraday dispersion (σ ϕ ). From THOR, we find median σ ϕ of 15.9 ± 3.2 rad m−2 for G46.8−0.3 and 17.6 ± 1.6 rad m−2 for G39.2−0.3. When comparing to polarization at 6 cm, we find evidence for σ ϕ ≳ 30 rad m−2 in localized regions where we detect no polarization in THOR. We combine Faraday depth dispersion with the rotation measure (RM) structure function (SF) and find evidence for a break in the SF on scales less than the THOR beam. We estimate the RM SF of the foreground interstellar medium using the SF of extragalactic radio sources and pulsars to find that the RM fluctuations we measure originate within the SNRs for all but the largest angular scales. [ABSTRACT FROM AUTHOR]

Voyager 1 has been moving through the very local interstellar medium (VLISM) for approximately one solar cycle, from 122.58 au on 2012/DOY 238 (August 25) to 158.5 au on 2023.0. Previously, an abrupt increase ("jump") in the magnetic field strength B and proton density N by a factor of 1.35 and 1.36, respectively, was observed during an interval of ∼8 days in 2020.40. After the jump, B continued to increase to a maximum value ∼0.56 nT at ∼2021.4 and then declined until B returned to the postjump value of 0.5 nT on 2021.85, 1.45 yr after the jump. The magnetic field strength declined briefly from 0.5 nT on 2021.85 to 0.47 nT on 2021.95 and then increased sporadically to 0.52 nT at 2023.0. Thus, the magnetic field strength remained strong for at least 2.6 yr. The magnetic hump and the density hump were a compression wave propagating through the VLISM. The compression wave was generated by a region with large dynamic pressure in the solar wind that propagated through the inner heliosheath and collided with the heliopause. The magnetic field strength continued to remain strong, with slow variations, until the end of our observations at 2023.0. It is suggested that the magnetic hump evolved from the large dynamic pressure, high speeds, and density observed at 1 au between ∼2015 and ∼2017. [ABSTRACT FROM AUTHOR]

We have obtained the polarization data cube of the supernova remnant (SNR) VRO 42.05.01 at 1240 MHz using the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). Three-dimensional Faraday synthesis is applied to the FAST data to derive the Faraday depth spectrum. The peak Faraday depth map shows a large area of enhanced foreground rotation measure (RM) of ∼60 rad m−2 extending along the remnant's "wing" section, which coincides with a large-scale H i shell at −20 km s−1. The two depolarization patches within the "wing" region with RMs of 97 rad m−2 and 55 rad m−2 coincide with two H i structures in the H i shell. Faraday screen model fitting on the Canadian Galactic Plane Survey 1420 MHz full-scale polarization data reveals a distance of (0.7–0.8) d SNR in front of the SNR with enhanced regular magnetic field there. The highly piled-up magnetic field indicates that the H i shell at −20 km s−1 could originate from an old evolved SNR. [ABSTRACT FROM AUTHOR]

We present observations of polarized dust emission at 850 μ m from the L43 molecular cloud, which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ( N H 2 ∼ 10 22 –1023 cm−2) complex molecular cloud with a submillimeter-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to ∼160 ± 30 μ G in the main starless core and up to ∼90 ± 40 μ G in the more diffuse, extended region. These field strengths give magnetically super- and subcritical values, respectively, and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. [ABSTRACT FROM AUTHOR]

A "Ribbon" of enhanced energetic neutral atom (ENA) emissions was discovered by the Interstellar Boundary Explorer in 2009, redefining our understanding of the heliosphere boundaries and the physical processes occurring at the interstellar interface. The Ribbon signal is intertwined with that of a globally distributed flux (GDF) that spans the entire sky. To a certain extent, Ribbon separation methods enabled examining its evolution independent of the underlying GDF. Observations over a full solar cycle revealed the Ribbon's evolving nature, with intensity variations closely tracking those of the solar wind (SW) structure after a few years delay, accounting for the SW–ENA recycling process. In this work, we examine the Ribbon structure, namely its ENA fluxes, angular extent, width, and circularity properties for two years, 2009 and 2019, representative of the declining phases of two adjacent solar cycles. We find that, (i) the Ribbon ENA fluxes have recovered in the nose direction and south of it down to ∼25° (for energies below 1.7 keV) and not at mid and high ecliptic latitudes; (ii) the Ribbon width exhibits significant variability as a function of azimuthal angle; (iii) circularity analysis suggests that the 2019 Ribbon exhibits a statistically consistent radius with that in 2009. The Ribbon's partial recovery is aligned with the consensus of a heliosphere with its closest point being southward of the nose region. The large variability of the Ribbon width as a function of azimuth in 2019 compared to 2009 is likely indicative of small-scale processes within the Ribbon. [ABSTRACT FROM AUTHOR]

The linear polarization of thermal dust emission provides a powerful tool to probe interstellar and circumstellar magnetic fields, because aspherical grains tend to align themselves with magnetic field lines. While the Radiative Alignment Torque (RAT) mechanism provides a theoretical framework for this phenomenon, some aspects of this alignment mechanism still need to be quantitatively tested. One such aspect is the possibility that the reference alignment direction changes from the magnetic field (" B -RAT") to the radiation field k-vector (" k -RAT") in areas of strong radiation fields. We investigate this transition toward the Orion Bar PDR, using multiwavelength SOFIA HAWC+ dust polarization observations. The polarization angle maps show that the radiation field direction is on average not the preferred grain alignment axis. We constrain the grain sizes for which the transition from B -RAT to k -RAT occurs in the Orion Bar (grains ≥ 0.1 μ m toward the most irradiated locations), and explore the radiatively driven rotational disruption that may take place in the high-radiation environment of the Bar for large grains. While the grains susceptible to rotational disruption should be in suprathermal rotation and aligned with the magnetic field, k -RAT aligned grains would rotate at thermal velocities. We find that the grain size at which the alignment shifts from B -RAT to k -RAT corresponds to grains too large to survive the rotational disruption. Therefore, we expect a large fraction of grains to be aligned at suprathermal rotation with the magnetic field, and to potentially be subject to rotational disruption, depending on their tensile strength. [ABSTRACT FROM AUTHOR]

In order to study the initial conditions of massive star formation, we have previously built a sample of 463 high-mass starless clumps (HMSCs) across the inner Galactic plane covered by multiple continuum surveys. Here, we use 13CO(2–1) line data from the SEDIGISM survey, which covers 78° in longitude (−60° < l < 18°, ∣ b ∣ < 0.°5) with 30″ resolution, to investigate the global dynamical state of these parsec-scale HMSCs (207 sources with good-quality data, mass 102–105 M ⊙, size 0.1–3.6 pc). We find that most HMSCs are highly turbulent with a median Mach number   ∼ 8.2 , and 44%–55% of them are gravitationally bound (with virial parameter α vir ≲ 2) if no magnetic fields are present. A median magnetic field strength of 0.33–0.37 mG would be needed to support these bound clumps against collapse, in agreement with previous observations of magnetic fields in regions of massive star formation. Luminosity-to-mass ratio, an important tracer of evolutionary stage, is strongly correlated with dust temperature. Magnetic field strength is also correlated with density. The Larson line width–size scaling does not hold in HMSCs. This study advances our understanding of the global properties of HMSCs, and our high-resolution observations with the Atacama Large Millimeter/submillimeter Array are in progress to study the resolved properties. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Filaments are ubiquitous in the interstellar medium, yet their formation and evolution remain the topic of intense debate. In order to obtain a more comprehensive view of the 3D morphology and evolution of the Musca filament, we model the C18O(2-1) emission along the filament crest with several large-scale velocity field structures. This indicates that Musca is well described by a 3D curved cylindrical filament with longitudinal mass inflow to its center unless the filament is a transient structure with a lifetime ≲0.1 Myr. Gravitational longitudinal collapse models of filaments appear unable to explain the observed velocity field. To better understand these kinematics, we further analyze a map of the C18O(2-1) velocity field at the location of SOFIA HAWC+ dust polarization observations that trace the magnetic field in the filament. This unveils an organized magnetic field that is oriented roughly perpendicular to the filament crest. Although the velocity field is also organized, it progressively changes its orientation by more than 90° when laterally crossing the filament crest and thus appears disconnected from the magnetic field in the filament. This strong lateral change of the velocity field over the filament remains unexplained and might be associated with important longitudinal motion that can be associated to the large-scale kinematics along the filament. [ABSTRACT FROM AUTHOR]

We investigate the connection between the regulation of star formation and the cycling of baryons both within and in and out of galaxies. We use idealized numerical simulations of Milky Way–mass galaxies, in which we vary the galaxy morphology and stellar feedback strength. By following individual gas parcels through the disk, spiral arms, and massive star-forming clumps, we quantify how gas moves through the different phases of the interstellar medium (ISM) and forms stars. We show that the residence time of gas in the dense ISM phase (τ SF), the nature of spiral arms, and the clump properties depend on both the galaxy morphology and stellar feedback. We quantify signatures of the baryon cycle within galaxies using the temporal and spatial power spectrum density (PSD) of the star formation rate (SFR). Stronger stellar feedback leads to more bursty star formation while the correlation timescale of the SFH is longer, because stronger feedback dissolves the dense ISM phase, leading to a more homogeneous ISM and a decrease in τ SF. The bulge strength has a similar effect: the deep gravitational potential in a bulge-dominant galaxy imposes a strong shear force that breaks apart gas clumps in the ISM; this subsequently inhibits the fragmentation of gas and therefore the star formation in the disk, leading to a decrease in the spatial power on scales of ∼1 kpc. We conclude that measurements of the temporal and spatial PSD of the SFR can provide constraints on the baryon cycle and the star formation process. [ABSTRACT FROM AUTHOR]

We present BISTRO Survey 850 μ m dust emission polarization observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the James Clerk Maxwell Telescope (JCMT). We observe a roughly triangular network of dense filaments. We detect nine of the dense starless cores embedded within these filaments in polarization, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field dominated to being matter dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 yr. [ABSTRACT FROM AUTHOR]

The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R136 has been shown to be the main source of energy creating multiple parsec-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B -fields) of 30 Doradus inferred from the far-infrared polarimetric observations by SOFIA/HAWC+ at 89, 154, and 214 μ m. The B -field morphology is complex, showing bending structures around R136. In addition, we use high spectral and angular resolution [C ii ] observations from SOFIA/GREAT and CO(2-1) from APEX. The kinematic structure of the region correlates with the B -field morphology and shows evidence of multiple expanding-shells. Our B -field strength maps, estimated using the Davis–Chandrasekhar–Fermi method and structure-function, show variations across the cloud within a maximum of 600, 450, and 350 μ G at 89, 154, and 214 μ m, respectively. We estimated that the majority of the 30 Doradus clouds are subcritical and sub-Alfvénic. The probability distribution function of the gas density shows that the turbulence is mainly compressively driven, while the plasma beta parameter indicates supersonic turbulence. We show that the B -field is sufficient to hold the cloud structure integrity under feedback from R136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique using [C ii ] and CO(2-1) is likely to confirm these suggestions. [ABSTRACT FROM AUTHOR]

We present the first Chandra X-ray Observatory (CXO) catalog of “pulsar X-ray filaments,” or “misaligned outflows.” These are linear, synchrotron-radiating features powered by ultrarelativistic electrons and positrons that escape from bow shock pulsars. The filaments are misaligned with the (large) pulsar velocity, distinguishing them from the pulsar wind nebula (PWN) trail, which is also often visible in CXO ACIS images. Spectral fits and morphological properties are extracted for five secure filaments and three candidates using a uniform method. We present a search of archival CXO data for linear diffuse features; the known examples are recovered and a few additional weak candidates are identified. We also report on a snapshot CXO ACIS survey of pulsars with properties similar to the filament producers, finding no new filaments but some diffuse emission, including one PWN trail. Finally, we provide an updated model for the pulsar properties required to create filaments in light of these new observations.

We use polarization data from SOFIA HAWC+ to investigate the interplay between magnetic fields and stellar feedback in altering gas dynamics within the high-mass star-forming region RCW 36, located in Vela C. This region is of particular interest as it has a bipolar H ii region powered by a massive star cluster, which may be impacting the surrounding magnetic field. To determine if this is the case, we apply the histogram of relative orientations (HRO) method to quantify the relative alignment between the inferred magnetic field and elongated structures observed in several data sets such as dust emission, column density, temperature, and spectral line intensity maps. The HRO results indicate a bimodal alignment trend, where structures observed with dense gas tracers show a statistically significant preference for perpendicular alignment relative to the magnetic field, while structures probed by the photodissociation region (PDR) tracers tend to align preferentially parallel relative to the magnetic field. Moreover, the dense gas and PDR associated structures are found to be kinematically distinct such that a bimodal alignment trend is also observed as a function of line-of-sight velocity. This suggests that the magnetic field may have been dynamically important and set a preferred direction of gas flow at the time that RCW 36 formed, resulting in a dense ridge developing perpendicular to the magnetic field. However, on filament scales near the PDR region, feedback may be energetically dominating the magnetic field, warping its geometry and the associated flux-frozen gas structures, causing the observed preference for parallel relative alignment.

Synchrotron observation serves as a tool for studying magnetic fields in the interstellar medium and intracluster medium, yet its ability to unveil three-dimensional (3D) magnetic fields, meaning probing the field’s plane-of-the-sky (POS) orientation, inclination angle relative to the line of sight, and magnetization from one observational data, remains largely underexplored. Inspired by the latest insights into anisotropic magnetohydrodynamic (MHD) turbulence, we found that synchrotron emission’s intensity structures inherently reflect this anisotropy, providing crucial information to aid in 3D magnetic field studies: (i) the structure’s elongation gives the magnetic field’s POS orientation and (ii) the structure’s anisotropy degree and topology reveal the inclination angle and magnetization. Capitalizing on this foundation, we integrate a machine learning approach—convolutional neural network (CNN)—to extract this latent information, thereby facilitating the exploration of 3D magnetic fields. The model is trained on synthetic synchrotron emission maps, derived from 3D MHD turbulence simulations encompassing a range of sub-Alfvénic to super-Alfvénic conditions. We show that the CNN is physically interpretable and the CNN is capable of obtaining the POS orientation, inclination angle, and magnetization. Additionally, we test the CNN against the noise effect and the missing low-spatial frequency. We show that this CNN-based approach maintains a high degree of robustness even when only high-spatial frequencies are maintained. This renders the method particularly suitable for application to interferometric data lacking single-dish measurements. We applied this trained CNN to the synchrotron observations of a diffuse region. The CNN-predicted POS magnetic field orientation shows a statistical agreement with that derived from synchrotron polarization.

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.

We analyze the count rates of 40–139 keV ions that were measured in situ by the Low Energy Charged Particle instrument on Voyager 1 in order to identify the suprathermal ion anisotropies beyond the termination shock and in the very local interstellar (IS) medium (VLISM). The analysis results in a region of ∼9–10 au before the heliopause (HP) where the radial anisotropy of ions is negative, while the azimuthal ion anisotropy inside the heliosheath lies in the − T direction. In agreement to our previous analyses, we identify a positive radial anisotropy of ions up to at least ∼30 au beyond the HP, which becomes nearly zero from 2021 up to 2023 November (for ∼10 au). Notably, the anisotropy in the azimuthal direction is statistically zero throughout the upstream region, i.e., for ∼40 au past the HP, showing that the 40–139 keV ion anisotropy in the VLISM is only in the radial direction and has no azimuthal component. The presence of suprathermal ions of solar origin over such long spatial scales in the VLISM, along with the inflow of ions from IS space into the heliosheath, are important constraints for characterizing the interaction of the heliosheath with the VLISM. Our observations provide indications that V1 has entered a new regime in the VLISM since (at least) the year 2021, progressively developing characteristics akin to the pristine IS medium. Alternatively, this drop to nearly zero radial anisotropies beyond 2021 could be a manifestation of a prolonged compression/shock of solar origin.

Using 3D dust maps and Planck polarized dust emission data, we investigate the influence of the 3D geometry of the nearby interstellar medium (ISM) on the statistics of the dust polarization on large ( $80^{\prime} $ ) scales. We test recent models that assume that the magnetic field probed by the polarized dust emission is preferentially tangential to the Local Bubble wall, but we do not find an imprint of the Local Bubble geometry on the dust polarization fraction. We also test the hypothesis that the complexity of the 3D dust distribution drives some of the measured variation of the dust polarization fraction. We compare sight lines with similar total column densities and find that, on average, the dust polarization fraction decreases when the dust column is substantially distributed among multiple components at different distances. Conversely, the dust polarization fraction is higher for sight lines where the dust is more concentrated in 3D space. This finding is statistically significant for the dust within 1.25 kpc, but the effect disappears if we only consider dust within 270 pc. In conclusion, we find that the extended 3D dust distribution, rather than solely the dust associated with the Local Bubble, plays a role in determining the observed dust polarization fraction at 80′. This conclusion is consistent with a simple analytical prediction and remains robust under various modifications to the analysis. These results illuminate the relationship between the 3D geometry of the ISM and tracers of the interstellar magnetic field. We discuss implications for our understanding of the polarized dust foreground to the cosmic microwave background.

We report on Atacama Large Millimeter/submillimeter Array observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from northwest to southeast with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS ( J = 5 → 4) emission yields an outflow energy of approximately 3.5 × 10 ^45 erg, aligning with previous interferometric studies. Utilizing the Davis–Chandrasekhar–Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 ± 1 mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary HII region is also larger than the magnetic energy, suggesting that protostellar feedback may be the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions.

The polarization spectrum, or wavelength dependence of the polarization fraction, of interstellar dust emission provides important insights into the grain alignment mechanism of interstellar dust grains. We investigate the far-infrared polarization spectrum of a realistic simulated high-mass star-forming cloud under various models of grain alignment and emission. We find that neither a homogeneous grain alignment model nor a grain alignment model that includes collisional dealignment is able to produce the falling spectrum seen in observations. On the other hand, we find that a grain alignment model with grain alignment efficiency dependent on local temperature is capable of producing a falling spectrum that is in qualitative agreement with observations of OMC-1. For the model most in agreement with OMC-1, we find no correlation between the temperature and the slope of the polarization spectrum. However, we do find a positive correlation between the column density and the slope of the polarization spectrum. We suggest this latter correlation to be the result of wavelength-dependent polarization by absorption.

We study the anisotropy of centroid and integrated intensity maps with synthetic observations. We perform postprocess radiative transfer including the optically thick regime that was not covered in Hernández-Padilla et al. We consider the emission in various CO molecular lines that range from optically thin to optically thick ( ^12 CO, ^13 CO, C ^18 O, and C ^17 O). The results for the velocity centroids are similar to those in the optically thin case. For instance, the anisotropy observed can be attributed to the Alfvén mode, which dominates over the slow and fast modes when the line of sight is at a high inclination with respect to the mean magnetic field. A few differences arise in the models with higher opacity, where some dependence on the sonic Mach number becomes evident. In contrast to the optically thin case, maps of integrated intensity become more anisotropic in optically thick lines. In this situation the scales probed are restricted, due to absorption, to smaller scales, which are known to be more anisotropic. We discuss how the sonic Mach number can affect the latter results, with highly supersonic cases exhibiting a lower degree of anisotropy.

By using surface brightness maps of Tycho’s supernova remnant (SNR) in radio and X-rays, along with the properties of thermal and synchrotron emission, we have derived the postshock density and magnetic field (MF) strength distributions over the projection of this remnant. Our analysis reveals a density gradient oriented toward the northwest, while the MF strength gradient aligns with the Galactic plane, pointing eastward. Additionally, utilizing this MF map, we have derived the spatial distributions of the cutoff frequency and maximum energy of electrons in Tycho’s SNR. We further comment on the implications of these findings for interpreting the gamma-ray emission from Tycho’s SNR.

Utilizing Planck polarized dust emission maps at 353 GHz and large-area maps of the neutral hydrogen (H i ) cold neutral medium (CNM) fraction ( f _CNM ), we investigate the relationship between dust polarization fraction ( p _353 ) and f _CNM in the diffuse high latitude ( $\left|b\right|\gt 30^\circ $ ) sky. We find that the correlation between p _353 and f _CNM is qualitatively distinct from the p _353 –H i column density ( N _H _i ) relationship. At low column densities ( N _H _i < 4 × 10 ^20 cm ^−2 ) where p _353 and N _H _i are uncorrelated, there is a strong positive p _353 – f _CNM correlation. We fit the p _353 – f _CNM correlation with data-driven models to constrain the degree of magnetic field disorder between phases along the line of sight. We argue that an increased magnetic field disorder in the warm neutral medium (WNM) relative to the CNM best explains the positive p _353 – f _CNM correlation in diffuse regions. Modeling the CNM-associated dust column as being maximally polarized, with a polarization fraction p _CNM ∼ 0.2, we find that the best-fit mean polarization fraction in the WNM-associated dust column is 0.22 p _CNM . The model further suggests that a significant f _CNM -correlated fraction of the non-CNM column (an additional 18.4% of the H i mass on average) is also more magnetically ordered, and we speculate that the additional column is associated with the unstable medium. Our results constitute a new large-area constraint on the average relative disorder of magnetic fields between the neutral phases of the interstellar medium, and are consistent with the physical picture of a more magnetically aligned CNM column forming out of a disordered WNM.

Probing magnetic fields in astrophysical environments is both important and challenging. The Gradient Technique (GT) is a new tool for tracing magnetic fields, rooted in the properties of magnetohydrodynamic (MHD) turbulence and turbulent magnetic reconnection. In this work, we examine the performance of GT when applied to synthetic synchrotron emission and spectroscopic data obtained from sub-Alfvénic and trans-Alfvénic MHD simulations. We demonstrate the improved accuracy of GT in tracing magnetic fields in the absence of low spatial frequencies. Additionally, we apply a low-spatial-frequency filter to a diffuse neutral hydrogen region selected from the GALFA-H i survey. Our results show an increased alignment between the magnetic fields inferred from GT and the Planck 353 GHz polarization measurements.

We present deep 1.4–4.8 μ m JWST-NIRCam imaging of the Serpens Main star-forming region and identify 20 candidate protostellar outflows, most with bipolar structure and identified driving sources. The outflow position angles (PAs) are strongly correlated, and they are aligned within ±24° of the major axis of the Serpens filament. These orientations are further aligned with the angular momentum vectors of the two disk shadows in this region. We estimate that the probability of this number of young stars being coaligned if sampled from a uniform PA distribution is 10 ^−4 . This in turn suggests that the aligned protostars, which seem to be at similar evolutionary stages based on their outflow dynamics, formed at similar times with a similar spin inherited from a local cloud filament. Further, there is tentative evidence for a systematic change in average PA between the northwestern and southeastern cluster, as well as increased scatter in the PAs of the southeastern protostars. SOFIA-HAWC+ archival dust polarization observations of Serpens Main at 154 and 214 μ m are perpendicular to the dominant jet orientation in the northwestern region in particular. We measure and locate shock knots and edges for all of the outflows and provide an identifying catalog. We suggest that Serpens main is a cluster that formed from an isolated filament and due to its youth retains its primordial outflow alignment.

A direct consequence of Faraday rotation is that the polarized radio sky does not resemble the total intensity sky at long wavelengths. We analyze G137+7, which is undetectable in total intensity but appears as a depolarization feature. We use the first polarization maps from the Canadian Hydrogen Intensity Mapping Experiment. Our 400–729 MHz bandwidth and angular resolution, $17^{\prime} $ – $30^{\prime} $ , allow us to use Faraday synthesis to analyze the polarization structure. In polarized intensity and polarization angle maps, we find a tail extending 10° from the head and designate the combined object, the tadpole . Similar polarization angles, distinct from the background, indicate that the head and tail are physically associated. The head appears as a depolarized ring in single channels, but wideband observations show that it is a Faraday rotation feature. Our investigations of H i and H α find no connections to the tadpole. The tail suggests motion of either the gas or an ionizing star through the interstellar medium; the B2(e) star HD 20336 is a candidate. While the head features a coherent, ∼ −8 rad m ^−2 Faraday depth, Faraday synthesis also identifies multiple components in both the head and tail. We verify the locations of the components in the spectra using QU fitting. Our results show that approximately octave-bandwidth Faraday rotation observations at ∼600 MHz are sensitive to low-density ionized or partially ionized gas, which is undetectable in other tracers.

We revisit the relation between the variance of three-dimensional (3D) density ( ${\sigma }_{\rho }^{2}$ ) and that of the projected two-dimensional (2D) column density ( ${\sigma }_{{\rm{\Sigma }}}^{2}$ ) in turbulent media, which is of great importance in obtaining turbulence properties from observations. Earlier studies showed that ${\sigma }_{{\rm{\Sigma }}/{{\rm{\Sigma }}}_{0}}^{2}/{\sigma }_{\rho /{\rho }_{0}}^{2}={ \mathcal R }$ , where Σ/Σ _0 and ρ / ρ _0 are 2D column and 3D volume densities normalized by their mean values, respectively. The factor ${ \mathcal R }$ depends only on the density spectrum for isotropic turbulence in a cloud that has similar dimensions along and perpendicular to the line of sight. Our major findings in this paper are as follows. First, we show that the factor ${ \mathcal R }$ can be expressed in terms of N , the number of independent eddies along the line of sight. To be specific, ${\sigma }_{{\rm{\Sigma }}/{{\rm{\Sigma }}}_{0}}^{2}/{\sigma }_{\rho /{\rho }_{0}}^{2}$ is proportional to ∼1/ N due to the averaging effect arising from independent eddies along the line of sight. Second, we show that the factor ${ \mathcal R }$ needs to be modified if the dimension of the cloud in the line-of-sight direction is different from that in the perpendicular direction. However, if we express ${\sigma }_{{\rm{\Sigma }}/{{\rm{\Sigma }}}_{0}}^{2}/{\sigma }_{\rho /{\rho }_{0}}^{2}$ in terms of N , the expression remains same even in the case the cloud has different dimensions along and perpendicular to the line of sight. Third, when we plot $N{\sigma }_{{\rm{\Sigma }}/{{\rm{\Sigma }}}_{0}}^{2}$ against ${\sigma }_{\rho /{\rho }_{0}}^{2}$ , two quantities roughly lie on a single curve regardless of the sonic Mach number, which implies that we can directly obtain the latter from the former. We discuss observational implications of our findings.

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.

We present 850 μ m linear polarization and C ^18 O (3 − 2) and ^13 CO (3 − 2) molecular line observations toward the filaments (F13 and F13S) in the Cocoon Nebula (IC 5146) using the JCMT POL-2 and Heterodyne Array Receiver Program instruments. F13 and F13S are found to be thermally supercritical with identified dense cores along their crests. Our findings include that the polarization fraction decreases in denser regions, indicating reduced dust grain alignment efficiency. The magnetic field vectors at core scales tend to be parallel to the filaments, but disturbed at the high density regions. Magnetic field strengths measured using the Davis–Chandrasekhar–Fermi method are 58 ± 31 and 40 ± 9 μ G for F13 and F13S, respectively, and it reveals subcritical and sub-Alfvénic filaments, emphasizing the importance of magnetic fields in the Cocoon region. Sinusoidal C ^18 O (3 − 2) velocity and density distributions are observed along the filaments’ skeletons, and their variations are mostly displaced by ∼1/4 × the wavelength of the sinusoid, indicating core formation occurred through the fragmentation of a gravitationally unstable filament, but with shorter core spacings than predicted. Large-scale velocity fields of F13 and F13S, studied using ^13 CO (3 − 2) data, present a V-shape transverse velocity structure. We propose a scenario for the formation and evolution of F13 and F13S, along with the dense cores within them. A radiation shock front generated by a B-type star collided with a sheet-like cloud about 1.4 Myr ago. The filaments became thermally critical due to mass infall through self-gravity ∼1 Myr ago, and subsequently, dense cores formed through gravitational fragmentation, accompanied by the disturbance of the magnetic field.

The radio transient ASKAP J173608.2–321635, at the position ( ℓ , b ) = (356.°0872, −0.°0390), was serendipitously observed by The H i /OH/Recombination line survey of the Galactic center at three epochs in 2020 March, 2020 April, and 2021 February. The source was detected only on 2020 April 11 with a flux density of 20.6 ± 1.1 mJy at 1.23 GHz and in-band spectral index of α = −3.1 ± 0.2. The commensal Very Large Array Low-band Ionsophere and Transient Experiment simultaneously detected the source at 339 MHz with a flux density of 122.6 ± 20.4 mJy, indicating a spectral break below 1 GHz. The rotation measure (RM) in 2020 April was 63.9 ± 0.3 rad m ^−2 , which almost triples the range of the variable RM observed by Wang et al. to ∼130 rad m ^−2 . The polarization angle, corrected for Faraday rotation, was 97° ± 6°. The 1.23 GHz linear polarization was 76.7% ± 3.9% with wavelength-dependent depolarization, indicating a Faraday depth dispersion of ${\sigma }_{\phi }={4.8}_{-0.7}^{+0.5}\ \mathrm{rad}\ {{\rm{m}}}^{-2}$ . We find an upper limit to the circular polarization of ∣ V ∣/ I < 10.1%. Interpretation of the data in terms of diffractive scattering of radio waves by a plasma near the source indicates an electron density and a line-of-sight magnetic field strength within a factor of 3 of n _e ∼ 2 cm ^−3 and B _∥ ∼ 2 × 10 ^5 μ G . Combined with causality limits to the size of the source, these parameters are consistent with the low-frequency spectral break resulting from synchrotron self-absorption, not free–free absorption. A possible interpretation of the source is a highly supersonic neutron star interacting with a changing environment.

The magnetic field of a molecular cloud core may play a role in the formation of circumstellar disks in the core. We present magnetic field morphologies in protostellar cores of 16 targets in the Atacama Large Millimeter/submillimeter Array large program “Early Planet Formation in Embedded Disks (eDisk),” which resolved their disks with 7 au resolutions. The 0.1 pc scale magnetic field morphologies were inferred from the James Clerk Maxwell Telescope POL-2 observations. The mean orientations and angular dispersions of the magnetic fields in the dense cores are measured and compared with the radii of the 1.3 mm continuum disks and the dynamically determined protostellar masses from the eDisk program. We observe a significant correlation between the disk radii and the stellar masses. We do not find any statistically significant dependence of the disk radii on the projected misalignment angles between the rotational axes of the disks and the magnetic fields in the dense cores, nor on the angular dispersions of the magnetic fields within these cores. However, when considering the projection effect, we cannot rule out a positive correlation between disk radii and misalignment angles in three-dimensional space. Our results suggest that the morphologies of magnetic fields in dense cores do not play a dominant role in the disk formation process. Instead, the sizes of protostellar disks may be more strongly affected by the amount of mass that has been accreted onto star+disk systems, and possibly other parameters, for example, magnetic field strength, core rotation, and magnetic diffusivity.

Angular dispersion functions are typically used to estimate the fluctuations in polarization angle around the mean magnetic field orientation in dense regions, such as molecular clouds. The technique provides accurate turbulent-to-regular magnetic field ratios, $\langle {B}_{t}^{2}{\rangle }^{1/2}/{B}_{\mathrm{pos}}$ , which are often underestimated by the classic Davis–Chandrasekhar–Fermi method. We assess the technique's suitability to characterize the turbulent and regular plane-of-sky magnetic field in low-density structures of the nearby interstellar medium (ISM), particularly when the turbulence outer scale, δ , is smaller than the smallest scale observed, ℓ _min . We use optical polarization maps of three intermediate-latitude fields (∣ b ∣ ≳ 7.°5) with dimensions of 0.°3 × 0.°3, sourced from the Interstellar Polarization Survey–General ISM catalog. We decomposed the H i emission detected by the Galactic All-Sky Survey within our fields to estimate the multiphase ISM properties associated with the structure coupled to the magnetic field. We produced maps of the plane-of-sky magnetic field strength ( B _pos ), mass density ( ρ ), and turbulent velocity dispersion ( σ _v _,turb ). In the regions with well-defined structures at d < 400 pc, the average B _– ranges from ∼3 μ G to ∼9 μ G , depending on the method and physical properties. In the region where structures extend up to 1000 pc, B _pos varies from ∼1 to ∼3 μ G . The results agree with previous estimations in the local, diffuse ISM. Finally, optical starlight polarization and thermal dust polarization at 353 GHz consistently reveal a highly regular plane-of-sky magnetic field orientation unfazed by diffuse dust structures observed at 12 μ m.

Magnetic fields are now widely recognized as critical at many scales to galactic dynamics and structure, including multiphase pressure balance, dust processing, and star formation. Using imposed magnetic fields cannot reliably model the interstellar medium's (ISM) dynamical structure nor phase interactions. Dynamos must be modeled. ISM models exist of turbulent magnetic fields using small-scale dynamo (SSD). Others model the large-scale dynamo (LSD) organizing magnetic fields at the scale of the disk or spiral arms. Separately, neither can fully describe the galactic magnetic field dynamics nor topology. We model the LSD and SSD together at a sufficient resolution to use the low explicit Lagrangian resistivity required. The galactic SSD saturates within 20 Myr. We show that the SSD is quite insensitive to the presence of an LSD and is even stronger in the presence of a large-scale shear flow. The LSD grows more slowly in the presence of SSD, saturating after 5 Gyr versus 1–2 Gyr in studies where the SSD is weak or absent. The LSD primarily grows in warm gas in the galactic midplane. Saturation of the LSD occurs due to α -quenching near the midplane as the growing mean-field produces a magnetic α that opposes the kinetic α . The magnetic energy in our models of the LSD shows a slightly sublinear response to increasing resolution, indicating that we are converging toward the physical solution at 1 pc resolution. Clustering supernovae in OB associations increases the growth rates for both the SSD and the LSD, compared to a horizontally uniform supernova distribution.

We report the detection of the ordered alignment between the magnetic field and kpc-scale bubbles in the nearby spiral galaxy, NGC 628. Applying the Velocity Gradient Technique on CO spectroscopic data from the ALMA-PHANGS, the magnetic field of NGC 628 is measured at the scale of 191 pc (∼4″). The large-scale magnetic field is oriented parallel to the spiral arms and curves around the galactic bubble structures in the mid-infrared emission observed by the James Webb Space Telescope. A total of 21 bubble structures have been identified at the edges of spiral arms with scales over 300 pc, which includes two kpc-scale structures. These bubbles are caused by supernova remnants and prolonged star formation and are similar to the outflow chimneys found in neutral hydrogen in galactic disks. At the edge of the bubbles, the shocks traced by the O iii emission present a curved magnetic field that parallels the bubble’s shell. The magnetic field follows the bubble expansion and binds the gas in the shell to trigger further star formation. By analyzing the larger sample of 1694 bubbles, we found a distinct radial-size distribution of bubbles in NGC 628 indicating the star formation history in the galaxy.

We investigate the role of magnetic field on the gas dynamics in a galactic bulge region by three-dimensional simulations with radiative cooling and heating. While a high-temperature corona with T > 10 ^6 K is formed in the halo regions, the temperature near the midplane is ≲10 ^4 K following the thermal equilibrium curve determined by the radiative cooling and heating. Although the thermal energy of the interstellar gas is lost by radiative cooling, the saturation level of the magnetic field strength does not significantly depend on the radiative cooling and heating. The magnetic field strength is amplified to 10 μ G on average and reaches several hundred microgauss locally. We find the formation of magnetically dominated regions at midlatitudes in the case with the radiative cooling and heating, which is not seen in the case without radiative effect. The vertical thickness of the midlatitude regions is 50–150 pc at the radial location of 0.4–0.8 kpc from the Galactic center, which is comparable to the observed vertical distribution of neutral atomic gas. When we take the average of different components of energy density integrated over the galactic bulge region, the magnetic energy is comparable to the thermal energy. We conclude that the magnetic field plays a substantial role in controlling the dynamical and thermal properties of the galactic bulge region.

Galactic dynamo models have generally relied on input parameters that are very challenging to constrain. We address this problem by developing a model that uses observable quantities as input: the galaxy rotation curve, the surface densities of the gas, stars and star formation rate, and the gas temperature. The model can be used to estimate parameters of the random and mean components of the magnetic field, as well as the gas scale height, root-mean-square velocity and the correlation length and time of the interstellar turbulence, in terms of the observables. We use our model to derive theoretical scaling relations for the quantities of interest, finding reasonable agreement with empirical scaling relations inferred from observation. We assess the dependence of the results on different assumptions about turbulence driving, finding that agreement with observations is improved by explicitly modeling the expansion and energetics of supernova remnants. The model is flexible enough to include alternative prescriptions for the physical processes involved, and we provide links to two open-source python programs that implement it.

Magnetic fields are ubiquitous in the Universe and are thought to play an important role in various astrophysical processes. Polarization of thermal emission from dust grains aligned with the magnetic field is widely used to measure the 2D magnetic field projected onto the plane of the sky, but its component along the line of sight is not yet constrained. Here, we introduce a new method to infer 3D magnetic fields using thermal dust polarization and grain alignment physics. We first develop a physical model of thermal dust polarization using the modern grain alignment theory based on the magnetically enhanced radiative torque alignment theory. We then test this model with synthetic observations of magnetohydrodynamic simulations of a filamentary cloud with our updated POLARIS code. Combining the tested physical polarization model with synthetic polarization, we show that the B-field inclination angles can be accurately constrained by the polarization degree from synthetic observations. Compared to the true 3D magnetic fields, our method based on grain alignment physics is more accurate than the previous methods that assume uniform grain alignment. This new technique paves the way for tracing 3D B-fields using thermal dust polarization and grain alignment theory and for constraining dust properties and grain alignment physics.

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter “ Y _turb ” in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

Starlight polarimetry, when combined with accurate distance measurements, allows for exploration of the three-dimensional structure of local magnetic fields in great detail. We present optical polarimetric observations of stars in and close to the Southern Coalsack, taken from the Interstellar Polarization Survey. Located in five fields of view approximately 0.°3 × 0.°3 in size, these data represent the highest density of optical polarimetric observations in the Southern Coalsack to date. Using these data, combined with accurate distances and extinctions based on Gaia data, we are able to characterize the magnetic field of the Coalsack and disentangle contributions to the polarization caused by the Southern Coalsack and a background structure. For the Southern Coalsack, we find an average magnetic field orientation of θ ∼ 75° with respect to the Galactic north pole and an average plane-of-sky magnetic field strength of approximately B _POS = 10 μ G, using the Davis–Chandrasekhar–Fermi method. These values are in agreement with some earlier estimates of the Coalsack’s magnetic field. In order to study the distant structure, we introduce a simple method to separate and isolate the polarization of distant stars from foreground contribution. For the distant structure, which we estimate to be located at a distance of approximately 1.3–1.5 kpc, we find an average magnetic field orientation of θ ∼ 100° and estimate a field strength of B _POS ∼ 10 μ G, although this will remain highly uncertain until the precise nature of the distant structure can be uncovered.

We present CO( J = 1–0) multiline observations toward the L914 dark cloud in the vicinity of the Cygnus X region, using the 13.7 m millimeter telescope of the Purple Mountain Observatory. The CO observations reveal in the L914 cloud a long filament with an angular length of ∼3.°6, corresponding to approximately 50 pc at the measured distance of ∼ 760 pc. Furthermore, a group of hair-like striations are discovered in two subregions of the L914 cloud, which are connected with the dense ridge of the filament. These striations display quasiperiodic characteristics in both the CO intensity images and position–velocity diagrams. Two of the striations also show increasing velocity gradients and dispersions toward the dense ridge, which could be fitted by accretion flows under gravity. Based on the Planck 353 GHz dust polarization data, we find that the striations are well aligned with the magnetic fields. Moreover, both the striations and magnetic fields are perpendicular to the dense ridge, constructing a bimodal configuration. Using the classic method, we estimate the strength of the magnetic field and further evaluate the relative importance of gravity, turbulence, and the magnetic field, finding that the L914 cloud is strongly magnetized. Our results suggest that magnetic fields play an important role in the formation of filamentary structures by channeling the material along the striations toward the dense ridge. The comparison between the observations and simulations suggests that striations could be a product of the magnetohydrodynamic process.

Direct-sampling observations of interstellar neutral gas, including hydrogen and deuterium, have been performed for more than one cycle of solar activity by IBEX. The IBEX viewing is restricted to directions perpendicular to the spacecraft–Sun line, which limits the observations to several months each year. This restriction is removed in the forthcoming mission called Interstellar Mapping and Acceleration Probe. The IMAP-Lo instrument will have the capability of adjusting the angle of its boresight with the spacecraft rotation axis. We continue a series of studies of the resulting science opportunities. We adopt the schedule of adjusting the boresight angle suggested by Kubiak et al. and focus on interstellar hydrogen and deuterium during solar maximum. Based on an extensive set of simulations, we identify the times during the calendar year and the elongation angles of the boresight needed to measure the abundance of D/H at the termination shock and to unambiguously observe interstellar H without a contribution from interstellar He. Furthermore, IMAP-Lo will be able to resolve the primary and secondary populations, in particular, to view the secondary population with little contribution from the primary. We show that the expected signal is sensitive to details of radiation pressure, particularly its dependence on the radial speed of the atoms, and to details of the behavior of the distribution function of the primary and secondary populations at the heliopause. Therefore, IMAP-Lo will be able to provide the observations needed to address compelling questions in heliospheric physics and even in general astrophysics.

We report 850 μ m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B -fields In STar forming Regions Observations large program on the James Clerk Maxwell Telescope. These data reveal a well-structured nonuniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30° from north to east. Field strength estimates and a virial analysis of the major clumps indicate that NGC 2264C is globally dominated by gravity, while in 2264D, magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and longitudinal collapsing, driven by the region's global gravity.

The Davis–Chandrasekhar–Fermi (DCF) method is widely used to evaluate magnetic fields in star-forming regions. Yet it remains unclear how well DCF equations estimate the mean plane-of-the-sky field strength in a map region. To address this question, five DCF equations are applied to an idealized cloud map. Its polarization angles have a normal distribution with dispersion σ _θ , and its density and velocity dispersion have negligible variation. Each DCF equation specifies a global field strength B _DCF and a distribution of local DCF estimates. The “most-likely” DCF field strength B _ml is the distribution mode, for which a correction factor β _ml ≡ B _ml / B _DCF is calculated analytically. For each equation, β _ml < 1, indicating that B _DCF is a biased estimator of B _ml . The values of β _ml are β _ml ≈ 0.7 when ${B}_{\mathrm{DCF}}\propto {\sigma }_{\theta }^{-1}$ due to turbulent excitation of Alfvénic MHD waves, and β _ml ≈ 0.9 when ${B}_{\mathrm{DCF}}\propto {\sigma }_{\theta }^{-1/2}$ due to non-Alfvénic MHD waves. These statistical correction factors β _ml have partial agreement with correction factors ${\beta }_{\mathrm{sim}}$ obtained from MHD simulations. The relative importance of the statistical correction is estimated by assuming that each simulation correction has both a statistical and a physical component. Then the standard, structure function, and original DCF equations appear most accurate because they require the least physical correction. Their relative physical correction factors are 0.1, 0.3, and 0.4 on a scale from 0 to 1. In contrast, the large-angle and parallel- δ B equations have physical correction factors 0.6 and 0.7. These results may be useful in selecting DCF equations, within model limitations.