Showing 1–50 of 448 results for author: Lazarian, A

Effects of Radiative Transfer on the Observed Anisotropy in MHD Turbulent Molecular Simulations

Authors: D. Hernández-Padilla, A. Esquivel, A. Lazarian, P. F. Velázquez, J. Cho

Abstract: We study the anisotropy of centroid and integrated intensity maps with synthetic observations. We perform post-process radiative transfer including the optically thick regime that was not covered in Hernández-Padilla et al. (2020). We consider the emission in various CO molecular lines, that range from optically thin to optically thick ($\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, $\mathrm{C^{18}O}$, a… ▽ More We study the anisotropy of centroid and integrated intensity maps with synthetic observations. We perform post-process radiative transfer including the optically thick regime that was not covered in Hernández-Padilla et al. (2020). We consider the emission in various CO molecular lines, that range from optically thin to optically thick ($\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, $\mathrm{C^{18}O}$, and $\mathrm{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. △ Less

Submitted 21 June, 2024; v1 submitted 19 June, 2024; originally announced June 2024.

Exploring magnetic fields in merging galaxy: combining polarization and velocity gradient in the Centaurus Galaxy

Authors: Quynh Lan Nguyen, Yue Hu, Alex Lazarian

Abstract: In this study, we apply the Velocity Gradient Technique (VGT) to the merging Centaurus galaxy. We compare gradient maps derived from the PHANGS-ALMA survey using CO emission lines with magnetic field tracings from dust polarization data obtained via the HAWC+ instrument. Our analysis reveals a strong correspondence between the directions indicated by these two tracers across most of the galactic i… ▽ More In this study, we apply the Velocity Gradient Technique (VGT) to the merging Centaurus galaxy. We compare gradient maps derived from the PHANGS-ALMA survey using CO emission lines with magnetic field tracings from dust polarization data obtained via the HAWC+ instrument. Our analysis reveals a strong correspondence between the directions indicated by these two tracers across most of the galactic image. Specifically, we identify jet regions as areas of anti-alignment, consistent with previous reports that gradients tend to rotate 90 degrees in outflow regions. Statistically, we find that the alignment of magnetic fields, as revealed by polarization, is most accurate in regions with the highest signal-to-noise ratios. Our findings underscore the utility of velocity gradients as a valuable complementary tool for probing magnetic fields and dynamical processes in merging galaxies. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 7 pages, 4 figures

Gradient Technique Theory: Tracing magnetic field and obtaining magnetic field strength

Authors: Alex Lazarian, Ka Ho Yuen, Dmitri Pogosyan

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 Alfven Mach number. We provide an analytical theory of gradient measurements and quantify the… ▽ More 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 Alfven 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 Alfven 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 Davis-Chandrasekhar-Fermi way of obtaining the magnetic field strength. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 31 pages, 11 figures, submitted to ApJ

Investigating Turbulence Effects on Magnetic Reconnection Rates Through High-Resolution Three-Dimensional Resistive Magnetohydrodynamical Simulations

Authors: Giovani H. Vicentin, Grzegorz Kowal, Elisabete M. de Gouveia Dal Pino, Alex Lazarian

Abstract: We investigate the impact of turbulence on magnetic reconnection through high-resolution 3D MHD simulations, covering Lundquist numbers from $S=10^3$ to $10^6$. Our simulations introduce small-scale perturbations injected into the system until $t=0.1\, t_A$. Even after the initial perturbation ceases, turbulence grows and persists, leading to sustained high reconnection rates of… ▽ More We investigate the impact of turbulence on magnetic reconnection through high-resolution 3D MHD simulations, covering Lundquist numbers from $S=10^3$ to $10^6$. Our simulations introduce small-scale perturbations injected into the system until $t=0.1\, t_A$. Even after the initial perturbation ceases, turbulence grows and persists, leading to sustained high reconnection rates of $V_\text{rec}/V_A \sim 0.03-0.08$, which surpass those generated solely by resistive tearing modes (plasmoids) observed in 2D and 3D PIC/MHD simulations by factors of 5 to 6. Our findings reveal that the achieved reconnection rates align with observations in solar phenomena and those reported in 3D MHD global simulations of accretion flows and relativistic jets. Notably, our simulations exhibit a steady-state reconnection rate coinciding with the full development of turbulence, showcasing the robustness and persistence of the reconnection in a turbulent environment. We establish the independence of the reconnection rate from the Lundquist number, consistent with the theory of turbulent reconnection proposed by Lazarian and Vishniac. Our results demonstrate a mild dependence of $V_\text{rec}$ on the plasma-$β$ parameter, the ratio between the thermal and magnetic pressure, illustrating a decrease from $0.036$ to $0.028$ as $β$ increases from $2$ to $64$ for simulations with $S=10^5$. Finally, we explore the influence of the magnetic Prandtl number ($\text{Pr}_m=ν/η$) on the reconnection rate and find this influence to be negligible during the turbulent regime across the range tested, from $\text{Pr}_m=1$ to $60$. These findings carry significant implications for research into reconnection and particle acceleration within realistic magnetized space and astrophysical environments. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: 22 pages, 21 figures

Cosmic Ray Diffusion in the Turbulent Interstellar Medium: Effects of Mirror Diffusion and Pitch Angle Scattering

Authors: Lucas Barreto-Mota, Elisabete M. de Gouveia Dal Pino, Siyao Xu, Alexandre Lazarian

Abstract: Cosmic rays (CRs) interact with turbulent magnetic fields in the intestellar medium, generating nonthermal emission. After many decades of studies, the theoretical understanding of their diffusion in the ISM continues to pose a challenge. This study numerically explores a recent prediction termed "mirror diffusion" and its synergy with traditional diffusion mechanism based on gyroresonant scatteri… ▽ More Cosmic rays (CRs) interact with turbulent magnetic fields in the intestellar medium, generating nonthermal emission. After many decades of studies, the theoretical understanding of their diffusion in the ISM continues to pose a challenge. This study numerically explores a recent prediction termed "mirror diffusion" and its synergy with traditional diffusion mechanism based on gyroresonant scattering. Our study combines 3D MHD simulations of star-forming regions with test particle simulations to analyze CR diffusion. We demonstrate the significance of mirror diffusion in CR diffusion parallel to the magnetic field, when the mirroring condition is satisfied. Our results support the theoretical expectation that the resulting particle propagation arising from mirror diffusion in combination with much faster diffusion induced by gyroresonant scattering resembles a Levy-flight-like propagation. Our study highlights the necessity to reevaluate the diffusion coefficients traditionally adopeted in the ISM based on gyroresonant scattering alone. For instance, our simulations imply a diffusion coefficient $\sim10^{27}cm^2/s$ for particles with a few hundred TeV within regions spanning a few parsecs around the source. This estimate is in agreement with gamma-ray observations, which shows the relevance of our results for understanding of diffuse gamma-ray emission in star-forming regions. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Probing Three-Dimensional Magnetic Fields: III -- Synchrotron Emission and Machine Learning

Authors: Yue Hu, Alex Lazarian

Abstract: Synchrotron observation serves as a fundamental tool for studying magnetic fields in various astrophysical settings, yet its ability to unveil three-dimensional (3D) magnetic fields-including plane-of-the-sky orientation, inclination angle relative to the line of sight, and magnetization-remains largely underexplored. Inspired by the latest insights into anisotropic magnetohydrodynamic (MHD) turbu… ▽ More Synchrotron observation serves as a fundamental tool for studying magnetic fields in various astrophysical settings, yet its ability to unveil three-dimensional (3D) magnetic fields-including plane-of-the-sky orientation, inclination angle relative to the line of sight, and magnetization-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, carrying detailed information about 3D magnetic fields. 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 model is physically interpretable and the CNN is capable of reconstructing 3D magnetic field topology and assessing magnetization. In addition, we test our methodology against noise and resolution effects. We show that this CNN-based approach maintains a high degree of robustness in tracing 3D magnetic fields, even when the low spatial frequencies of the synchrotron image are absent. This renders the method particularly suitable for application to interferometric data lacking single-dish measurements. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 12 pages, 8 figures, submitted to ApJ

Magnetic Field of Molecular Gas Measured with the Velocity Gradient Technique II: Curved Magnetic Field in kpc-Scale Bubble of NGC\,628

Authors: Mengke Zhao, Jianjun Zhou, Willem A. Baan, Yue Hu, A. Lazarian, Xindi Tang, Jarken Esimbek, Yuxin He, Dalei Li, Weiguang Ji, Zhengxue Chang, Kadirya Tursun

Abstract: 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 (VGT) on CO spectroscopic data from the ALMA-PHANGS, the magnetic field of NGC\,628 is measured at the scale of 191\,pc ($\sim$ 4\,$''$). The large-scale magnetic field is oriented parallel to the spiral arms and curves ar… ▽ More 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 (VGT) on CO spectroscopic data from the ALMA-PHANGS, the magnetic field of NGC\,628 is measured at the scale of 191\,pc ($\sim$ 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 (JWST). Twenty-one 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 OIII 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. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 15 pages, 7 figures, Accepted by ApJ

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

On the properties and implications of collapse-driven MHD turbulence

Authors: Enrique Vázquez-Semadeni, Yue Hu, Siyao Xu, Rubén Guerrero-Gamboa, Alex Lazarian

Abstract: We numerically investigate the driving of MHD turbulence by gravitational contraction using simulations of an initially spherical, magnetically supercritical cloud core with initially transonic and trans-Alfvénic turbulence. We perform a Helmholtz decomposition of the velocity field, and investigate the evolution of its solenoidal and compressible parts, as well as of the velocity component along… ▽ More We numerically investigate the driving of MHD turbulence by gravitational contraction using simulations of an initially spherical, magnetically supercritical cloud core with initially transonic and trans-Alfvénic turbulence. We perform a Helmholtz decomposition of the velocity field, and investigate the evolution of its solenoidal and compressible parts, as well as of the velocity component along the gravitational acceleration vector, a proxy for the infall component of the velocity field. We find that: 1) In spite of being supercritical, the core first contracts to a sheet perpendicular to the mean field, and the sheet itself collapses. 2) The solenoidal component of the turbulence remains at roughly its initial level throughout the simulation, while the compressible component increases continuously. This implies that turbulence does {\it not} dissipate towards the center of the core. 3) The distribution of simulation cells in the $B$-$ρ$ plane occupies a wide triangular region at low densities, bounded below by the expected trend for fast MHD waves ($B \propto ρ$, applicable for high local Alfvénic Mach number $\Ma$) and above by the trend expected for slow waves ($B \sim$ constant, applicable for low local $\Ma$). At high densities, the distribution follows a single trend $B \propto ρ^{\gamef}$, with $1/2 < \gamef < 2/3$, as expected for gravitational compression. 4) The measured mass-to-magnetic flux ratio $λ$ increases with radius $r$, due to the different scalings of the mass and magnetic flux with $r$. At a fixed radius, $λ$ increases with time due to the accretion of material along field lines. 5) The solenoidal energy fraction is much smaller than the total turbulent component, indicating that the collapse drives the turbulence mainly compressibly, even in directions orthogonal to that of the collapse. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Resubmitted to MNRAS after first set of reviewer's recommendations. Comments welcome

Sub-Alfvenic Turbulence: Magnetic to Kinetic Energy Ratio, Modification of Weak Cascade and Implications for Magnetic Field Strength Measurement

Authors: Alex Lazarian, Ka Wai HO, Ka Ho Yuen, Ethan Vishniac

Abstract: We study the properties of sub-Alfvenic magnetohydrodynamic (MHD) turbulence, i.e., turbulence with Alfven Mach number $M_A=V_L/V_A<1$, where $V_L$ is the velocity at the injection scale and $V_A$ is the Alfven velocity. We demonstrate that weak turbulence can have different regimes depending on whether it is driven by velocity or magnetic fluctuations. If the turbulence is driven by isotropic bul… ▽ More We study the properties of sub-Alfvenic magnetohydrodynamic (MHD) turbulence, i.e., turbulence with Alfven Mach number $M_A=V_L/V_A<1$, where $V_L$ is the velocity at the injection scale and $V_A$ is the Alfven velocity. We demonstrate that weak turbulence can have different regimes depending on whether it is driven by velocity or magnetic fluctuations. If the turbulence is driven by isotropic bulk forces, i.e. velocity-driven, in an incompressible conducting fluid, we predict that the kinetic energy is $M_A^{-2}$ times larger than the energy of magnetic fluctuations. This effect arises from the long parallel wavelength tail of the forcing, which excites modes with $k_\|/k_\perp < M_A$. We also predict that as the turbulent cascade reaches the strong regime the energy of slow modes exceeds the energy of Alfven modes by a factor $M_A^{-1}$. These effects are absent if the turbulence is magnetically driven at the injection scale. We confirm these predictions with numerical simulations. As the assumption of magnetic and kinetic energy equipartition is at the core of the Davis-Chandrasekhar-Fermi (DCF) approach to measuring magnetic field strength in sub-Alfvenic turbulence, we conclude that the DCF technique is not universally applicable. In particular, we suggest that the dynamical excitation of long azimuthal wavelength modes in the galactic disk may compromise the use of the DCF technique. We discuss alternative expressions that can be used to obtain magnetic field strength from observations. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 22 pages, 10 figures, Submitted to ApJ

Velocity Gradient and Stellar Polarization: Magnetic Field Tomography towards the L1688 Cloud

Authors: Tyler Schmaltz, Yue Hu, A. Lazarian

Abstract: Magnetic fields are a defining yet enigmatic aspect of the interstellar medium (ISM), with their three-dimensional mapping posing a substantial challenge. In this study, we harness the innovative Velocity Gradient Technique (VGT), underpinned by magnetohydrodynamic (MHD) turbulence theories, to elucidate the magnetic field structure by applying it to the atomic neutral hydrogen (HI) emission line… ▽ More Magnetic fields are a defining yet enigmatic aspect of the interstellar medium (ISM), with their three-dimensional mapping posing a substantial challenge. In this study, we harness the innovative Velocity Gradient Technique (VGT), underpinned by magnetohydrodynamic (MHD) turbulence theories, to elucidate the magnetic field structure by applying it to the atomic neutral hydrogen (HI) emission line and the molecular tracer $^{12}$CO. We construct the tomography of the magnetic field in the low-mass star-forming region L1688, utilizing two approaches: (1) VGT-HI combined with the Galactic rotational curve, and (2) stellar polarization paired with precise star parallax measurements. Our analysis reveals that the magnetic field orientations deduced from stellar polarization undergo a distinct directional change in the vicinity of L1688, providing evidence that the misalignment between VGT-HI and stellar polarization stems from the influence of the molecular cloud's magnetic field on the polarization of starlight. When comparing VGT-$^{12}$CO to stellar polarization and Planck polarization data, we observe that VGT-$^{12}$CO effectively reconciles the misalignment noted with VGT-HI, showing statistical alignment with Planck polarization measurements. This indicates that VGT-$^{12}$CO could be integrated with VGT-HI, offering vital insights into the magnetic fields of molecular clouds, thereby enhancing the accuracy of our 3D magnetic field reconstructions. △ Less

Submitted 12 January, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 12 pages, 12 figures, accepted for publication in MNRAS

Detection of large-scale synchrotron radiation from the molecular envelope of the Sgr B cloud complex at the Galactic center

Authors: F. Yusef-Zadeh, M. Wardle, R. Arendt, J. W. Hewitt, Y. Hu, A. Lazarian, N. Kassim, S. Hyman, I. Heywood

Abstract: We present highly sensitive measurements taken with MeerKAT at 1280 MHz as well as archival GBT, MWA and VLA images at 333, 88 and 74 MHz. We report the detection of synchrotron radio emission from the infrared dark cloud (IRDC) associated with the halo of the Sgr B complex on a scale of ~60 pc. A strong spatial correlation between low-frequency radio continuum emission and dense molecular gas, co… ▽ More We present highly sensitive measurements taken with MeerKAT at 1280 MHz as well as archival GBT, MWA and VLA images at 333, 88 and 74 MHz. We report the detection of synchrotron radio emission from the infrared dark cloud (IRDC) associated with the halo of the Sgr B complex on a scale of ~60 pc. A strong spatial correlation between low-frequency radio continuum emission and dense molecular gas, combined with spectral index measurements, indicates enhanced synchrotron emission by cosmic-ray electrons. Correlation of the FeI 6.4 keV Kalpha line and synchrotron emission provides compelling evidence that the low energy cosmic-ray electrons are responsible for producing the Kalpha line emission. The observed synchrotron emission within the halo of the Sgr B cloud complex has mean spectral index alpha -1+/-1 gives the magnetic field strength ~100 muG for cloud densities nH = 10^4-10^5 cm-3 and estimate cosmic-ray ionization rates between 10^-13 and 10^-14 s^-1. Furthermore, the energy spectrum of primary cosmic-ray electrons is constrained to be E^-3 +/-1 for typical energies of few hundred MeV. The extrapolation of this spectrum to higher energies is consistent with X-ray and gamma-ray emission detected from this cloud. These measurements have important implications on the role that high cosmic-ray electron fluxes at the Galactic center play in production of radio synchrotron emission, the FeI Kalpha line emission at 6.4 keV and ~GeV gamma-ray emission throughout the central molecular zone (CMZ). △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 10 pages, 4 figures, MN (in press)

Probing Three-Dimensional Magnetic Fields: II -- An Interpretable Convolutional Neural Network

Authors: Yue Hu, A. Lazarian, Yan Wu, Chengcheng Fu

Abstract: Observing 3D magnetic fields, including orientation and strength, within the interstellar medium is vital but notoriously difficult. However, recent advances in our understanding of anisotropic magnetohydrodynamic (MHD) turbulence demonstrate that MHD turbulence and 3D magnetic fields leave their imprints on the intensity features of spectroscopic observations. Leveraging these theoretical framewo… ▽ More Observing 3D magnetic fields, including orientation and strength, within the interstellar medium is vital but notoriously difficult. However, recent advances in our understanding of anisotropic magnetohydrodynamic (MHD) turbulence demonstrate that MHD turbulence and 3D magnetic fields leave their imprints on the intensity features of spectroscopic observations. Leveraging these theoretical frameworks, we propose a novel Convolutional Neural Network (CNN) model to extract this embedded information, enabling the probe of 3D magnetic fields. This model examines not only the plane-of-the-sky magnetic field orientation ($φ$), but also the magnetic field's inclination angle ($γ$) relative to the line-of-sight, and the total magnetization level (M$_A^{-1}$) of the cloud. We train the model using synthetic emission lines of $^{13}$CO (J = 1 - 0) and C$^{18}$O (J = 1 - 0), generated from 3D MHD simulations that span conditions from sub-Alfvénic to super-Alfvénic molecular clouds. Our tests confirm that the CNN model effectively reconstructs the 3D magnetic field topology and magnetization. The median uncertainties are under $5^\circ$ for both $φ$ and $γ$, and less than 0.2 for M$_A$ in sub-Alfvénic conditions (M$_A\approx0.5$). In super-Alfvénic scenarios (M$_A\approx2.0$), they are under $15^\circ$ for $φ$ and $γ$, and 1.5 for M$_A$. We applied this trained CNN model to the L1478 molecular cloud. Results show a strong agreement between the CNN-predicted magnetic field orientation and that derived from Planck 353 GHz polarization data. The CNN approach enabled us to construct the 3D magnetic field map for L1478, revealing a global inclination angle of $\approx76^\circ$ and a global M$_A$ of $\approx1.07$. △ Less

Submitted 5 December, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

Magnetic Fields and Velocity Gradients in L1551: The Role of Stellar Feedback

Authors: Mingrui Liu, Yue Hu, Alex Lazarian

Abstract: Magnetic fields play a crucial role in star formation, yet tracing them becomes particularly challenging, especially in the presence of outflow feedback in protostellar systems. We targeted the star-forming region L1551, notable for its apparent outflows, to investigate the magnetic fields. These fields were probed using polarimetry observations from the Planck satellite at 353 GHZ/849 $μ$m, the S… ▽ More Magnetic fields play a crucial role in star formation, yet tracing them becomes particularly challenging, especially in the presence of outflow feedback in protostellar systems. We targeted the star-forming region L1551, notable for its apparent outflows, to investigate the magnetic fields. These fields were probed using polarimetry observations from the Planck satellite at 353 GHZ/849 $μ$m, the SOFIA/HAWC+ measurement at 214 $μ$m, and the JCMT/SCUPOL 850 $μ$m survey. Consistently, all three measurements show that the magnetic fields twist towards the protostar IRS 5. Additionally, we utilized the Velocity Gradients Technique (VGT) on the $^{12}$CO (J = 1-0) emission data to distinguish the magnetic fields directly associated with the protostellar outflows. These were then compared with the polarization results. Notably, in the outskirts of the region, these measurements generally align. However, as one approaches the center of IRS 5, the measurements tend to yield mostly perpendicular relative orientations. This suggests that the outflows might be dynamically significant from a scale of approximately $\sim0.2$~pc, causing the velocity gradient to change direction by 90 degrees. Furthermore, we discovered that the polarization fraction $p$ and the total intensity $I$ correlate as $p \propto I^{-α}$. Specifically, $α$ is approximately $1.044\pm0.06$ for SCUPOL and around $0.858\pm0.15$ for HAWC+. This indicates that the outflows could significantly impact the alignment of dust grains and magnetic fields in the L1551 region. △ Less

Submitted 22 March, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 11pages, 19 figures, accepted for publication in MNRAS

Particle acceleration in self-driven turbulent reconnection

Authors: Jian-Fu Zhang, Siyao Xu, Alex Lazarian, Grzegorz Kowal

Abstract: The theoretical prediction that magnetic reconnection spontaneously drives turbulence has been supported by magnetohydrodynamic (MHD) and kinetic simulations. While reconnection with externally driven turbulence is accepted as an effective mechanism for particle acceleration, the acceleration during the reconnection with self-driven turbulence is studied for the first time in this work. By using h… ▽ More The theoretical prediction that magnetic reconnection spontaneously drives turbulence has been supported by magnetohydrodynamic (MHD) and kinetic simulations. While reconnection with externally driven turbulence is accepted as an effective mechanism for particle acceleration, the acceleration during the reconnection with self-driven turbulence is studied for the first time in this work. By using high-resolution 3D MHD simulations of reconnection with self-generated turbulence, we inject test particles into the reconnection layer to study their acceleration process. We find that the energy gain of the particles takes place when they bounce back and forth between converging turbulent magnetic fields. The particles can be efficiently accelerated in self-driven turbulent reconnection with the energy increase by about 3 orders of magnitude in the range of the box size. The acceleration proceeds when the particle gyroradii become larger than the thickness of the reconnection layer. We find that the acceleration in the direction perpendicular to the local magnetic field dominates over that in the parallel direction. The energy spectrum of accelerated particles is time-dependent with a slope that evolves toward -2.5. Our findings can have important implications for particle acceleration in high-energy astrophysical environments. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: 13 pages, 5 Figures. Accepted for publication in Journal of High Energy Astrophysics

Mirror acceleration of cosmic rays in a high-$β$ medium

Authors: Alex Lazarian, Siyao Xu

Abstract: In a weakly compressible high-$β$ medium, pitch-angle scattering and the associated scattering acceleration of cosmic rays (CRs) by anisotropic Alfvén and slow modes of magnetohydrodynamic (MHD) turbulence is inefficient. To tap the energy from magnetic compressions for efficient particle acceleration, a diffusion mechanism that can effectively confine particles in space without causing their trap… ▽ More In a weakly compressible high-$β$ medium, pitch-angle scattering and the associated scattering acceleration of cosmic rays (CRs) by anisotropic Alfvén and slow modes of magnetohydrodynamic (MHD) turbulence is inefficient. To tap the energy from magnetic compressions for efficient particle acceleration, a diffusion mechanism that can effectively confine particles in space without causing their trapping or pitch-angle isotropization is needed. We find that the mirror diffusion in MHD turbulence recently identified in Lazarian and Xu (2021) satisfies all the above conditions and serves as a promising diffusion mechanism for efficient acceleration of CRs via their stochastic non-resonant interactions with magnetic compressions/expansions. The resulting mirror acceleration is dominated by the slow-mode eddies with their lifetime comparable to the mirror diffusion time of CRs. Consequently, we find that the acceleration time of mirror acceleration is independent of the spatial diffusion coefficient of CRs. The mirror acceleration brings new life for the particle acceleration in a weakly compressible/incompressible medium and has important implications for studying CR re-acceleration in the high-$β$ intracluster medium. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 10 pages, 2 figures, submitted to ApJ

Synchrotron Intensity Gradient Revealing Magnetic Fields in Galaxy Clusters

Authors: Yue Hu, C. Stuardi, A. Lazarian, G. Brunetti, A. Bonafede, Ka Wai Ho

Abstract: Magnetic fields and their dynamical interplay with matter in galaxy clusters contribute to the physical properties and evolution of the intracluster medium. However, the current understanding of the origin and properties of cluster magnetic fields is still limited by observational challenges. In this article, we map the magnetic fields at hundreds-kpc scales of five clusters RXC J1314.4 -2515, Abe… ▽ More Magnetic fields and their dynamical interplay with matter in galaxy clusters contribute to the physical properties and evolution of the intracluster medium. However, the current understanding of the origin and properties of cluster magnetic fields is still limited by observational challenges. In this article, we map the magnetic fields at hundreds-kpc scales of five clusters RXC J1314.4 -2515, Abell 2345, Abell 3376, MCXC J0352.4 -7401, and El Gordo using the innovative synchrotron intensity gradient technique in conjunction with high-resolution radio observations from JVLA and MeerKAT. We demonstrate that magnetic field orientation of radio relics derived from synchrotron intensity gradients is in very good agreement with that obtained with synchrotron polarization. Most important, synchrotron intensity gradients is not limited by Faraday depolarization in the cluster central regions and allows us to map magnetic fields in the radio halos of RXC J1314.4 -2515 and El Gordo. We find that magnetic fields in radio halos exihibit a preferential direction along the major merger axis and show turbulent structures at higher angular resolution. Results are consistent with expectations from numerical simulations which predict turbulent magnetic fields in cluster mergers that are stirred and amplified by matter motions. △ Less

Submitted 9 February, 2024; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 25 pages, 14 figures, accepted

Damping of MHD Turbulence in A Partially Ionized Medium

Authors: Yue Hu, Siyao Xu, Lev Arzamasskiy, James M. Stone, A. Lazarian

Abstract: The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions+neutrals) simulations generated with the AthenaK code. Specifically, we examine the velocity, densi… ▽ More The coupling state between ions and neutrals in the interstellar medium plays a key role in the dynamics of magnetohydrodynamic (MHD) turbulence, but is challenging to study numerically. In this work, we investigate the damping of MHD turbulence in a partially ionized medium using 3D two-fluid (ions+neutrals) simulations generated with the AthenaK code. Specifically, we examine the velocity, density, and magnetic field statistics of the two-fluid MHD turbulence in different regimes of neutral-ion coupling. Our results demonstrate that when ions and neutrals are strongly coupled, the velocity statistics resemble those of single-fluid MHD turbulence. Both the velocity structures and kinetic energy spectra of ions and neutrals are similar, while their density structures can be significantly different. With an excess of small-scale sharp density fluctuations in ions, the density spectrum in ions is shallower than that of neutrals. When ions and neutrals are weakly coupled, the turbulence in ions is more severely damped due to the ion-neutral collisional friction than that in neutrals, resulting in a steep kinetic energy spectrum and density spectrum in ions compared to the Kolmogorov spectrum. We also find that the magnetic energy spectrum basically follows the shape of the kinetic energy spectrum of ions, irrespective of the coupling regime. In addition, we find large density fluctuations in ions and neutrals and thus spatially inhomogeneous ionization fractions. As a result, the neutral-ion decoupling and damping of MHD turbulence take place over a range of length scales. △ Less

Submitted 10 November, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 17 pages, 14 figures, accepted for publication in MNRAS

Nature of Striation in 21 cm Channel Maps: Velocity Caustics

Authors: Yue Hu, A. Lazarian, D. Alina, D. Pogosyan, Ka Wai Ho

Abstract: The alignment of striated intensity structures in thin neutral hydrogen (HI) spectroscopic channels with Galactic magnetic fields has been observed. However, the origin and nature of these striations are still debatable. Some studies suggest that the striations result solely from real cold-density filaments without considering the role of turbulent velocity fields, i.e., the velocity caustics effe… ▽ More The alignment of striated intensity structures in thin neutral hydrogen (HI) spectroscopic channels with Galactic magnetic fields has been observed. However, the origin and nature of these striations are still debatable. Some studies suggest that the striations result solely from real cold-density filaments without considering the role of turbulent velocity fields, i.e., the velocity caustics effect in shaping the channel's intensity distribution. To determine the relative contribution of density and velocity in forming the striations in channel maps, we analyze synthetic observations of channel maps obtained with simulations that represent realistic magnetized multi-phase HI. We vary the thickness of the channel maps and apply the Velocity Decomposition Algorithm to separate the velocity and density contributions. In parallel, we analyze GALFA HI observations and compare the results. Our analysis shows that the thin channels are dominated by velocity contribution, and velocity caustics mainly generate the HI striations. We show that velocity caustics can cause a correlation between unsharp-masked HI structures and far-infrared emission. We demonstrate that the linear HI fibers revealed by the Rolling Hough Transform (RHT) in thin velocity channels originate from velocity caustics. As the thickness of channel maps increases, the relative contribution of density to fluctuations in channel maps also increases. As a result, more RHT-detected fibers tend to be perpendicular to the magnetic field. Conversely, the alignment with the magnetic field is the most prominent in thin channels. We conclude that similar to the Velocity Channel Gradients (VChGs) approach, RHT traces magnetic fields through the analysis of velocity caustics in thin channel maps. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: 26 pages, 18 figures, accepted for publication in MNRAS

The origin of dust polarization in the Orion Bar

Authors: Valentin J. M. Le Gouellec, B-G Andersson, Archana Soam, Thiébaut Schirmer, Joseph M. Michail, Enrique Lopez-Rodriguez, Sophia Flores, David T. Chuss, John E. Vaillancourt, Thiem Hoang, Alex Lazarian

Abstract: 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 to this phenomenon, some aspects of this alignment mechanism still need to be quantitatively tested. One… ▽ More 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 to 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 multi-wavelength 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 occur 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 supra-thermal 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 supra-thermal rotation with the magnetic field, and potentially be subject to rotational disruption depending on their tensile strength. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Cosmic ray propagation in turbulent magnetic fields

Authors: Alex Lazarian, Siyao Xu, Yue Hu

Abstract: Propagation of cosmic rays (CRs) in turbulent and magnetized astrophysical media is a long-standing problem that requires both understanding of the properties of turbulent magnetic fields and their interaction with energetic particles. This review focuses on selected recent theoretical findings made based on the progress in understanding and simulating magnetohydrodynamic (MHD) turbulence. In part… ▽ More Propagation of cosmic rays (CRs) in turbulent and magnetized astrophysical media is a long-standing problem that requires both understanding of the properties of turbulent magnetic fields and their interaction with energetic particles. This review focuses on selected recent theoretical findings made based on the progress in understanding and simulating magnetohydrodynamic (MHD) turbulence. In particular, we address the problem of perpendicular and parallel propagation of CRs and identify the conditions when the perpendicular propagation is superdiffusive and diffusive. For the parallel diffusion, we discuss the problems of the traditionally used diffusion mechanism arising from pitch angle scattering and the possible solutions provided by the recently identified ``mirror diffusion" in the presence of turbulent magnetic mirrors. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: 21 pages, 4 figures, submitted to Frontiers in Astronomy and Space Sciences

Outlook on Magnetohydrodynamical Turbulence and its Astrophysical Implications

Authors: Elena Popova, Alexandre Lazarian

Abstract: Magnetohydrodynamical (MHD) turbulence is ubiquitous in magnetized astrophysical plasmas, and it radically changes a great variety of astrophysical processes. In this review, we give the concept of MHD turbulence and explain the origin of its scaling. We consider the implications of MHD turbulence to various problems: dynamo in different types of stars, flare activity, solar and stellar wind from… ▽ More Magnetohydrodynamical (MHD) turbulence is ubiquitous in magnetized astrophysical plasmas, and it radically changes a great variety of astrophysical processes. In this review, we give the concept of MHD turbulence and explain the origin of its scaling. We consider the implications of MHD turbulence to various problems: dynamo in different types of stars, flare activity, solar and stellar wind from different stars, propagation of cosmic rays, and star formation. We also discuss how the properties of MHD turbulence provide a new way of tracing magnetic fields in interstellar and intracluster media. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: arXiv admin note: text overlap with arXiv:astro-ph/0312046 by other authors

Mapping the Galactic Magnetic Field Orientation and Strength in Three Dimensions

Authors: Yue Hu, Alex Lazarian

Abstract: The mapping of the Galactic Magnetic Field (GMF) in three dimensions is essential to comprehend various astrophysical processes that occur within the Milky Way. This study endeavors to map the GMF by utilizing the latest MM2 technique, the Velocity Gradient Technique (VGT), the Column Density Variance Approach, and the GALFA-H I survey of Neutral Hydrogen (H I) emission. The MM2 and VGT methods re… ▽ More The mapping of the Galactic Magnetic Field (GMF) in three dimensions is essential to comprehend various astrophysical processes that occur within the Milky Way. This study endeavors to map the GMF by utilizing the latest MM2 technique, the Velocity Gradient Technique (VGT), the Column Density Variance Approach, and the GALFA-H I survey of Neutral Hydrogen (H I) emission. The MM2 and VGT methods rely on an advanced understanding of magnetohydrodynamics turbulence to determine the magnetic field strength and orientation respectively. The H I emission data, combined with the Galactic rotational curve, gives us the distribution of H I gas throughout the Milky Way. By combining these two techniques, we map the GMF orientation and strength, as well as the Alfvén Mach number $M_{\rm A}$ in 3D for a low galactic latitude ($b<30^{\rm o}$) region close to the Perseus Arm. The analysis of column density variance gives the sonic Mach number $M_{\rm s}$ distribution, The results of this study reveal the sub-Alfvénic and subsonic (or trans-sonic) nature of the H I gas. The variation of mean $M_{\rm A}$ along the line-of-sight approximately ranges from 0.6 to 0.9, while that of mean $M_{\rm s}$ is from 0.2 to 1.5. The mean magnetic field strength varies from ~0.5 $μ$G to ~2.5 $μ$G exhibiting a decreasing trend towards the Galaxy's outskirt. This work provides a new avenue for mapping the GMF, especially the magnetic field strength, in 3D. We discuss potential synergetic applications with other approaches. △ Less

Submitted 30 June, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

Tracing of Magnetic field with gradients: Sub-Sonic Turbulence

Authors: Ka Wai Ho, Alex Lazarian

Abstract: Recent development of the velocity gradient technique shows the {\toreferee capability} of the technique in the way of tracing magnetic fields morphology in diffuse interstellar gas and molecular clouds. In this paper, we perform the numerical systemic study of the performance of velocity and synchrotron gradient for a wide range of magnetization in the sub-sonic environment. Addressing the studie… ▽ More Recent development of the velocity gradient technique shows the {\toreferee capability} of the technique in the way of tracing magnetic fields morphology in diffuse interstellar gas and molecular clouds. In this paper, we perform the numerical systemic study of the performance of velocity and synchrotron gradient for a wide range of magnetization in the sub-sonic environment. Addressing the studies of magnetic field in atomic hydrogen, we also study the formation of velocity caustics in the spectroscopic channel maps in the presence of the thermal broadening. We show that the velocity caustics can be recovered when applied to the Cold Neutral Medium (CNM) and the Gradient Technique (GT) can reliably trace magnetic fields there. Finally, we discuss the changes of the anisotropy of observed structure functions when we apply to the analysis the procedures developed within the framework of GT studies. △ Less

Submitted 31 January, 2023; originally announced January 2023.

Comments: 11 pages, 7 figures

Journal ref: MNRAS 2023

Anomalous compressible mode generation by global frame projections of pure Alfven mode

Authors: Ka Ho Yuen, Huirong Yan, Alex Lazarian

Abstract: Alfven wave is the single most important physical phenomenon of magneto-hydrodynamic turbulence and has far-reaching impact to almost all studies related to astrophysical magnetic field. Yet the restoration of the Alfven wave fluctuations from a given magnetic field, aka the local Alfven wave problem, is never properly addressed in literature albeit its importance. Previous works model the Alfven… ▽ More Alfven wave is the single most important physical phenomenon of magneto-hydrodynamic turbulence and has far-reaching impact to almost all studies related to astrophysical magnetic field. Yet the restoration of the Alfven wave fluctuations from a given magnetic field, aka the local Alfven wave problem, is never properly addressed in literature albeit its importance. Previous works model the Alfven wave fluctuation as the perturbation along a straight-line, constant magnetic field. However, Lazarian & Pogosyan (2012) suggested that the decomposition of Alfven wave along a straight line, aka. the global frame decomposition, has a factor of discrepancy to the true local Alfven wave fluctuation. Here we provide a geometric interpretation on how the local Alfven wave is related to the global frame through the use of vector frame formulation. We prove both analytically and numerically that the local frame Alfven wave is an orthogonal transformation of that of the global frame and related by the local Alfvenic Mach number. In other words, when we observe Alfven wave in the global frame of reference, some of the Alfven wave will be mistaken as compressible waves. The importance of frame choices have a far-reaching impact to the analytical studies of MHD turbulence. Combining the frame formalism and the new techniques we can have accurate measurement to some of the fundamental turbulence properties like the inclination angle of mean magnetic field relative to the line of sight. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 18 pages, 12 figures, accepted by MNRAS for publication

Turbulent reconnection acceleration

Authors: Siyao Xu, Alex Lazarian

Abstract: The ubiquitous turbulence in astrophysical plasmas is important for both magnetic reconnection and reconnection acceleration. We study the particle acceleration during fast 3D turbulent reconnection with reconnection-driven turbulence. Particles bounce back and forth between the reconnection-driven inflows due to the mirror reflection and convergence of strong magnetic fields. Via successive head-… ▽ More The ubiquitous turbulence in astrophysical plasmas is important for both magnetic reconnection and reconnection acceleration. We study the particle acceleration during fast 3D turbulent reconnection with reconnection-driven turbulence. Particles bounce back and forth between the reconnection-driven inflows due to the mirror reflection and convergence of strong magnetic fields. Via successive head-on collisions, the kinetic energy of the inflows is converted into the accelerated particles. Turbulence not only regulates the inflow speed but also introduces various inflow obliquities with respect to the local turbulent magnetic fields. As both the energy gain and escape probability of particles depend on the inflow speed, the spectral index of particle energy spectrum is not universal. We find it in the range from $\approx 2.5$ to $4$, with the steepest spectrum expected at a strong guide field, i.e. a small angle between the total incoming magnetic field and the guide field. Without scattering diffusion needed for confining particles, the reconnection acceleration can be very efficient at a large inflow speed and a weak guide field. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 16 pages, 8 figures, accepted for publication in ApJ

Studying the properties of compressible MHD turbulence by synchrotron fluctuation statistics

Authors: Ru-Yue Wang, Jian-Fu Zhang, Alex Lazarian, Hua-Ping Xiao, Fu-Yuan Xiang

Abstract: We study the observable properties of compressible MHD turbulence covering different turbulence regimes, based on synthetic synchrotron observations arising from 3D MHD numerical simulations. Using the synchrotron emissivity and intensity, we first explore how the cosmic ray spectral indices affect the measurements of turbulence properties by employing normalized correlation functions. We then stu… ▽ More We study the observable properties of compressible MHD turbulence covering different turbulence regimes, based on synthetic synchrotron observations arising from 3D MHD numerical simulations. Using the synchrotron emissivity and intensity, we first explore how the cosmic ray spectral indices affect the measurements of turbulence properties by employing normalized correlation functions. We then study how the anisotropy of synchrotron total and polarization intensities arising from three fundamental MHD modes vary with the viewing angle, i.e., the angle between the mean magnetic field and the line of sight. We employ the ratio of quadrupole moment to the monopole one (QM) for this purpose. Our numerical results demonstrate that: (1) the two-point correlation function of synchrotron statistics for the arbitrary cosmic ray spectral index is related to the special case of magnetic field index $γ=2$ in agreement with the analytical formulae provided by Lazarian \& Pogosyan (2012); (2) the anisotropy of synchrotron total and polarization intensities arising from Alfvén and slow modes increases with the increase of the viewing angle, while that of fast mode remains almost unchanged with the viewing angle; (3) the analytical formulae of synchrotron intensities for studying turbulence can be applied to describing statistics of polarization intensities, and the QM can be successfully used to recover turbulence anisotropy. This study validates Lazarian \& Pogosyan's analytical approach and opens a way to study turbulence from observations. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: Accepted by ApJ. 15 pages, 9 figures and 1 Table

Parallel and Perpendicular Alignments of Velocity Gradient and Magnetic Field observed in the Molecular Clouds L1478 and L1482

Authors: Tyler Schmaltz, Yue Hu, Alex Lazarian

Abstract: Star formation is a complex process that typically occurs in dense regions of molecular clouds mainly regulated by magnetic fields, magnetohydrodynamic (MHD) turbulence, and self-gravity. However, it remains a challenging endeavor to trace the magnetic field and determine regions of gravitational collapse where the star is forming. Based on the anisotropic properties of MHD turbulence, a new techn… ▽ More Star formation is a complex process that typically occurs in dense regions of molecular clouds mainly regulated by magnetic fields, magnetohydrodynamic (MHD) turbulence, and self-gravity. However, it remains a challenging endeavor to trace the magnetic field and determine regions of gravitational collapse where the star is forming. Based on the anisotropic properties of MHD turbulence, a new technique termed Velocity Gradient Technique (VGT) has been proposed to address these challenges. In this study, we apply the VGT to two regions of the giant California Molecular Cloud (CMC), namely, L1478 and L1482, and analyze the difference in their physical properties. We use the $^{12}$CO (J = 2 - 1), $^{13}$CO (J = 2 - 1), and C$^{18}$O (J = 2 - 1) emission lines observed with the Heinrich Hertz Submillimeter Telescope. We compare VGT results calculated in the resolutions of $3.3'$ and $10'$ to Planck polarization at 353 GHz and $10'$ to determine areas of MHD turbulence dominance and self-gravity dominance. We show that the resolution difference can introduce misalignment between the two measurements. We find the VGT-measured magnetic fields globally agree with that from Planck in L1478 suggesting self-gravity's effect is insignificant. The best agreement appears in VGT-$^{12}$CO. As for L1482, the VGT measurements are statistically perpendicular to the Planck polarization indicating the dominance of self-gravity. This perpendicular alignment is more significant in VGT-$^{13}$CO and VGT-C$^{18}$O. △ Less

Submitted 18 May, 2023; v1 submitted 22 October, 2022; originally announced October 2022.

Comments: 14 pages, 16 figures, accepted for publication in MNRAS

Characterizing three-dimensional magnetic field, turbulence, and self-gravity in the star-forming region L1688

Authors: Yue Hu, Alex Lazarian

Abstract: Interaction of three-dimensional magnetic fields, turbulence, and self-gravity in the molecular cloud is crucial in understanding star formation but has not been addressed so far. In this work, we target the low-mass star-forming region L1688 and use the spectral emissions of $^{12}$CO, $^{13}$CO, C$^{18}$O, and H I, as well as polarized dust emissions. To obtain the 3D direction of the magnetic f… ▽ More Interaction of three-dimensional magnetic fields, turbulence, and self-gravity in the molecular cloud is crucial in understanding star formation but has not been addressed so far. In this work, we target the low-mass star-forming region L1688 and use the spectral emissions of $^{12}$CO, $^{13}$CO, C$^{18}$O, and H I, as well as polarized dust emissions. To obtain the 3D direction of the magnetic field, we employ the novel polarization fraction analysis. In combining with the plane-of-the-sky (POS) magnetic field strength derived from the Davis-Chandrasekhar-Fermi (DCF) method and the new Differential Measure Analysis (DMA) technique, we present the first measurement of L1688's three-dimensional magnetic field, including its orientation and strength. We find that L1688's magnetic field has two statistically different inclination angles. The low-intensity tail has an inclination angle $\approx55^\circ$ on average, while that of the central dense clump is $\approx30^\circ$. We find the global mean value of total magnetic field strength is $B_{\rm tot}\approx135$ uG from DCF and $B_{\rm tot}\approx75$ uG from DMA. We use the velocity gradient technique (VGT) to separate the magnetic fields' POS orientation associated with L1688 and its foreground/background. The magnetic fields' orientations are statistically coherent. The probability density function of H$_2$ column density and VGT reveal that L1688 is potentially undergoing gravitational contraction at large scale $\approx1.0$ pc and gravitational collapse at small scale $\approx0.2$ pc. The gravitational contraction mainly along the magnetic field results in an approximate power-law relation $B_{\rm tot}\propto n_{\rm H}^{1/2}$ when volume density $n_{\rm H}$ is less than approximately $6.0\times10^3$ cm$^{-3}$. △ Less

Submitted 13 July, 2023; v1 submitted 20 October, 2022; originally announced October 2022.

Comments: 17 pages, 20 figures, accepted for publication in MNRAS

A Lunar Backup Record of Humanity

Authors: Carson Ezell, Alexandre Lazarian, Abraham Loeb

Abstract: The risk of a catastrophic or existential disaster for our civilization is increasing this century. A significant motivation for a near-term space settlement is the opportunity to safeguard civilization in the event of a planetary-scale disaster. A catastrophic event could destroy the significant cultural, scientific, and technological progress on Earth. However, early space settlements can preser… ▽ More The risk of a catastrophic or existential disaster for our civilization is increasing this century. A significant motivation for a near-term space settlement is the opportunity to safeguard civilization in the event of a planetary-scale disaster. A catastrophic event could destroy the significant cultural, scientific, and technological progress on Earth. However, early space settlements can preserve records of human activity by maintaining a backup data storage system. The backup can also store information about the events leading up to the disaster. The system would improve the ability of early space settlers to recover our civilization after collapse. We show that advances in laser communications and data storage enable the development of a data storage system on the lunar surface with a sufficient uplink data rate and storage capacity to preserve valuable information about the achievements of our civilization and the chronology of the disaster. △ Less

Submitted 7 November, 2022; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 12 pages, 2 figures; accepted for publication in the journal "Signals" (2022)

Multi-Phase Magnetic Fields in the Galaxy NGC 3627

Authors: Mingrui Liu, Yue Hu, A. Lazarian, Siyao Xu, Marian Soida

Abstract: Magnetic fields play an important role in the formation and evolution of a galaxy, but it is challenging to measure them by observation. Here we study the Seyfert galaxy NGC 3627's magnetic field orientations measured from the synchrotron polarization observed with the Very Large Array (VLA) and from the Velocity Gradient Technique (VGT) using spectroscopic data. The latter employs the magnetohydr… ▽ More Magnetic fields play an important role in the formation and evolution of a galaxy, but it is challenging to measure them by observation. Here we study the Seyfert galaxy NGC 3627's magnetic field orientations measured from the synchrotron polarization observed with the Very Large Array (VLA) and from the Velocity Gradient Technique (VGT) using spectroscopic data. The latter employs the magnetohydrodynamical (MHD) turbulence's anisotropy to probe the magnetic fields. Being applied to the CO (2-1) and H$α$ emission lines obtained from the PHANGS-ALMA and PHANGS-MUSE surveys, it reveals the magnetic field orientation globally consistent with the polarization. The agreement of the VGT-CO and polarization suggests that the magnetic fields associated with synchrotron emission also percolate through star-forming regions. The VGT-H$α$ measurement reveals the magnetic fields in the warm ionized medium that permeates the disk and disk's vicinity so that it exhibits less agreement with polarization. We find prominent radial fields measured by synchrotron polarization appear in the transition regions from the spiral arms to the galactic bar, while such morphology is less apparent in the VGT-CO and VGT-H$α$ measured magnetic fields. The radial fields suggest that the magnetic torque is important in removing orbiting gas' angular momentum. We notice that magnetic fields inferred from the dust polarization, VGT-CO, and synchrotron polarization are different in the east arm. We interpret this difference as arising from the fact that the three measurements are tracing the magnetic fields associated with pre-collision, the mixture of pre-collision and post-collision, and post-collision flows, respectively. △ Less

Submitted 28 November, 2022; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 12 pages, 8 figures, accepted for publication in MNRAS

Tracing Magnetic Fields with the Gradient Technique: Spatial Filtering Effect and Use of Interferometers

Authors: Yue Hu, Alex Lazarian

Abstract: 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-Alf… ▽ More 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-\ion{H}{1} survey. Our results show an increased alignment between the magnetic fields inferred from GT and the Planck 353 GHz polarization measurements. △ Less

Submitted 15 June, 2024; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 24 pages, 14 figures, accepted for publication in ApJ

Turbulent Magnetic Field Amplification by the Interaction of Shock Wave and Inhomogeneous Medium

Authors: Yue Hu, Siyao Xu, James M. Stone, Alex Lazarian

Abstract: Magnetic fields on the order of 100 $μ$G observed in young supernova remnants cannot be amplified by shock compression alone. To investigate the amplification caused by turbulent dynamo, we perform three-dimensional MHD simulations of the interaction between shock wave and inhomogeneous density distribution with a shallow spectrum in the preshock medium. The postshock turbulence is mainly driven b… ▽ More Magnetic fields on the order of 100 $μ$G observed in young supernova remnants cannot be amplified by shock compression alone. To investigate the amplification caused by turbulent dynamo, we perform three-dimensional MHD simulations of the interaction between shock wave and inhomogeneous density distribution with a shallow spectrum in the preshock medium. The postshock turbulence is mainly driven by the strongest preshock density contrast and follows the Kolmogorov scaling. The resulting turbulence amplifies the postshock magnetic field. The time evolution of the magnetic fields agrees with the prediction of the nonlinear turbulent dynamo theory in Xu & Lazarian (2016). When the initial weak magnetic field is perpendicular to the shock normal, the maximum amplification of the field's strength achieves a factor of $\approx200$, which is twice larger than that for a parallel shock. We find that the perpendicular shock exhibits a smaller turbulent Alfvén Mach number in the vicinity of the shock front than the parallel shock. However, the strongest magnetic field has a low volume filling factor and is limited by the turbulent energy due to the reconnection diffusion taking place in a turbulent and magnetized fluid. The magnetic field strength averaged along the $z$-axis is reduced by a factor $\gtrsim10$. We decompose the turbulent velocity and magnetic field into solenoidal and compressive modes. The solenoidal mode is dominant and evolves to follow the Kolmogorov scaling, even though the preshock density distribution has a shallow spectrum. When the preshock density distribution has a Kolmogorov spectrum, the turbulent velocity's compressive component increases. △ Less

Submitted 9 November, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 19 pages, 15 figures, accepted for publication in ApJ

Magnetic Field of Molecular Gas Measured with the Velocity Gradient Technique I. Orion A

Authors: Mengke Zhao, Jianjun Zhou, Yue Hu, A. Lazarian, Xindi Tang, Willem A. Baan, Jarken Esimbek, Yuxin He, Dalei Li, Weiguang Ji, Kadirya Tursun

Abstract: Magnetic fields play an important role in the evolution of molecular clouds and star formation. Using the Velocity Gradient Technique (VGT) model, we measured the magnetic field in Orion A using the 12CO, 13CO, and C18O (1-0) emission lines at a scale of 0.07 pc. The measured B-field shows an east-west orientation that is perpendicular to the integral shaped filament of Orion A at large scale. The… ▽ More Magnetic fields play an important role in the evolution of molecular clouds and star formation. Using the Velocity Gradient Technique (VGT) model, we measured the magnetic field in Orion A using the 12CO, 13CO, and C18O (1-0) emission lines at a scale of 0.07 pc. The measured B-field shows an east-west orientation that is perpendicular to the integral shaped filament of Orion A at large scale. The VGT magnetic fields obtained from 13CO and C18O are in agreement with the B-field that is measured from the Planck 353 GHz dust polarization at a scale of 0.55 pc. Removal of density effects by using a Velocity Decomposition Algorithm can significantly improve the accuracy of the VGT in tracing magnetic fields with the 12CO (1-0) line. The magnetic field strength of seven sub-clouds, OMC-1, OMC-2, OMC-3, OMC-4, OMC-5, L 1641-N, and NGC 1999 has also been estimated with the Davis-Chandrasekhar-Fermi (DCF) and MM2 technique, and these are found to be in agreement with previous results obtained from dust polarization at far-infrared and sub-millimeter wavelengths. At smaller scales, the VGT proves a good method to measure magnetic fields. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 15 pages, 6 figures, accepted to ApJ

Role of magnetic fields in fueling Seyfert nuclei

Authors: Yue Hu, Alex Lazarian, Rainer Beck, Siyao Xu

Abstract: Molecular gas is believed to be the fuel for star formation and nuclear activity in Seyfert galaxies. To explore the role of magnetic fields in funneling molecular gas into the nuclear region, measurements of the magnetic fields embedded in molecular gas are needed. By applying the new velocity gradient technique (VGT) to ALMA and PAWS's CO isotopolog data, we obtain the first detection of CO-asso… ▽ More Molecular gas is believed to be the fuel for star formation and nuclear activity in Seyfert galaxies. To explore the role of magnetic fields in funneling molecular gas into the nuclear region, measurements of the magnetic fields embedded in molecular gas are needed. By applying the new velocity gradient technique (VGT) to ALMA and PAWS's CO isotopolog data, we obtain the first detection of CO-associated magnetic fields in several nearby Seyfert galaxies and their unprecedented high-resolution magnetic field maps. The VGT-measured magnetic fields in molecular gas globally agree with those inferred from existing HAWC+ dust polarization and VLA synchrotron polarization. An overall good alignment between the magnetic fields traced by VGT-CO and by synchrotron polarization may support the correlation between star formation and cosmic ray generation. We find that the magnetic fields traced by VGT-CO have a significant radial component in the central regions of most Seyferts in our sample, where efficient molecular gas inflows or outflows may happen. In particular, we find local misalignment between the magnetic fields traced by CO and dust polarization within the nuclear ring of NGC 1097, and the former aligns with the central bar's orientation. This misalignment reveals different magnetic field configurations in different gas phases and may provide an observational diagnostic for the ongoing multi-phase fueling of Seyfert activity. △ Less

Submitted 27 October, 2022; v1 submitted 11 June, 2022; originally announced June 2022.

Comments: 24 pages, 14 figures

SOFIA observations of 30 Doradus: II -- Magnetic fields and large scale gas kinematics

Authors: Le Ngoc Tram, Lars Bonne, Yue Hu, Enrique Lopez-Rodriguez, Jordan A. Guerra, Pierre Lesaffre, Antoine Gusdorf, Thiem Hoang, Min-Young Lee, Alex Lazarian, B-G Andersson, Simon Coude, Archana Soam, William D. Vacca, Hyeseung Lee, Michael Gordon

Abstract: The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-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-fiel… ▽ More The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-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 R$\,$136. In addition, we use high spectral and angular resolution [\textsc{CII}] 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 sub-critical 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 R$\,$136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique (VGT) using [\textsc{CII}] and CO(2-1) is likely to confirm these results. △ Less

Submitted 12 December, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 26 pages, 17 figures, 3 tables, accepted to ApJ

Unveiling polarized emission from interstellar dust of the Large Magellanic Cloud with Planck

Authors: D. Alina, J. -Ph. Bernard, K. H. Yuen, A. Lazarian, A. Hughes, M. Iskakova, A. Akimkhan, A. Mukanova

Abstract: Polarization of interstellar dust emission is a powerful probe of dust properties and magnetic field structure. Yet studies of external galaxies are hampered by foreground dust contribution. The aim of this study is to separate the polarised signal from the Large Magellanic Cloud (LMC) from that of the Milky Way (MW) in order to construct a wide-field, spatially complete map of dust polarization u… ▽ More Polarization of interstellar dust emission is a powerful probe of dust properties and magnetic field structure. Yet studies of external galaxies are hampered by foreground dust contribution. The aim of this study is to separate the polarised signal from the Large Magellanic Cloud (LMC) from that of the Milky Way (MW) in order to construct a wide-field, spatially complete map of dust polarization using the Planck 353 GHz data. To estimate the foreground polarization direction, we used velocity gradients in HI spectral line data and assessed the performance of the output by comparing to starlight extinction polarization. We estimate the foreground intensity using dust properties derived previously from the Planck data and we assume the foreground polarization to be uniform and equal to the average of the MW around the galaxy. After foreground removal, the geometry of the plane-of-the-sky magnetic field in the LMC tends to follow the structure of the atomic gas in the LMC. This is notably the case along the molecular ridges extending south and south-east of the 30 Doradus star-forming complex, and along more diffuse southern arm extending towards the Small Magellanic Cloud. There is also an alignment between the magnetic field and the outer arm in the western part. The median polarization fraction in the LMC is slightly lower than that observed for the MW while the anti-correlation between the polarization angle dispersion function and the polarization fraction is slightly larger. Overall, polarization fraction distribution is similar to that observed in the MW. △ Less

Submitted 29 November, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Comments: 15 pages, 22 figures

Journal ref: MNRAS, 2022

Turbulent universal galactic Kolmogorov velocity cascade over 6 decades

Authors: Ka Ho Yuen, Ka Wai Ho, Chi Yan Law, Avi Chen, Alex Lazarian

Abstract: We report the evidence for the existence of the universal, continuous turbulent cascade of velocity fluctuations with Kolmogorov -5/3 slope spanning 6 orders of length scales, from $10^4$ pc down to $10^{-2}$ pc. This was achieved by applying our innovative technique of separating density and velocity fluctuations to a set of spectroscopic surveys featuring various galactic spectral lines. This un… ▽ More We report the evidence for the existence of the universal, continuous turbulent cascade of velocity fluctuations with Kolmogorov -5/3 slope spanning 6 orders of length scales, from $10^4$ pc down to $10^{-2}$ pc. This was achieved by applying our innovative technique of separating density and velocity fluctuations to a set of spectroscopic surveys featuring various galactic spectral lines. This unified velocity cascade involves different interstellar phases from diffuse galactic media to dense self-gravitating clouds and persists despite interstellar phase transitions. However, the turbulent density fluctuations do not show this universality as the value of the spectral slope changes in different interstellar phases. This agrees with the expectation of compressible turbulence theory and demonstrates that the density is only an indirect tracer of interstellar turbulence. We report that the density fluctuations for clouds and filaments that are preferentially parallel to magnetic fields exhibit the spectral slope of -2. The universal grand turbulent velocity cascade that is established in our paper has significant implications for key galactic physical processes, including star formation, cosmic ray transport etc. We anticipate our result to be a starting point for in vitro models of multiphase interstellar turbulence studies with a significant impact for modeling of spiral galaxies. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Magnetic field strength from turbulence theory (I): Using differential measure approach (DMA)

Authors: A. Lazarian, Ka Ho Yuen, Dmitri Pogosyan

Abstract: 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 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 fie… ▽ More 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 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, 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 the astrophysical conditions for which our technique can be applied, we consider the turbulence in both media with magnetic pressure larger than the gas pressure corresponding e.g. to molecular and the gas pressure larger than the magnetic pressure corresponding to the warm neutral medium. We provide general expressions for arbitrary admixture of Alfvén, slow and fast modes in these media and consider in detail the 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 (DMA), unlike the DCF, 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 large scale variations induces of non-turbulent nature. In parallel, our study uncover the deficiencies of the earlier DCF research. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: 53 pages, 32 figures, Accepted by ApJ for publications. Code: https://github.com/kyuen2/MHD_Mode

Probing Three-Dimensional Magnetic Fields: I -- Polarized Dust Emission

Authors: Yue Hu, A. Lazarian

Abstract: Polarized dust emission is widely used to trace the plane-of-the-sky (POS) component of interstellar magnetic fields in two dimensions. Its potential to access three-dimensional magnetic fields, including the inclination angle of the magnetic fields relative to the line-of-sight (LOS), is crucial for a variety of astrophysical problems. Based on the statistical features of observed polarization fr… ▽ More Polarized dust emission is widely used to trace the plane-of-the-sky (POS) component of interstellar magnetic fields in two dimensions. Its potential to access three-dimensional magnetic fields, including the inclination angle of the magnetic fields relative to the line-of-sight (LOS), is crucial for a variety of astrophysical problems. Based on the statistical features of observed polarization fraction and POS Alfvén Mach number $\overline{M_{\rm A}}_{,\bot}$ distribution, we present a new method for estimating the inclination angle. The magnetic field fluctuations raised by anisotropic magnetohydrodynamic (MHD) turbulence are taken into account in our method. By using synthetic dust emission generated from 3D compressible MHD turbulence simulations, we show that the fluctuations are preferentially perpendicular to the mean magnetic field. We find the inclination angle is the major agent for depolarization, while fluctuations of magnetic field strength and density have an insignificant contribution. We propose and demonstrate that the mean inclination angle over a region of interest can be calculated from the polarization fraction in a strongly magnetized reference position, where $\overline{M_{\rm A}}_{,\bot}^2\ll1$. We test and show that the new method can trace the 3D magnetic fields in sub-Alfvénic, trans-Alfvénic, and moderately super-Alfvénic conditions ($0.4\lesssim M_{\rm A}\lesssim1.2$). We numerically quantify that the difference between the estimated inclination angle and actual inclination angle ranges from 0 to $20^\circ$ with a median value of $\le10^\circ$. △ Less

Submitted 16 December, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 14 pages, 11 figures, accepted for publication in MNRAS

Correlation of velocity and density contributions to spectroscopic channel maps: Reality check on Kalberla et.al (2022)

Authors: Ka Ho Yuen, Ka Wai Ho, Alex Lazarian

Abstract: The existence of magnetized turbulence in the interstellar HI is well accepted. A number of techniques to obtain turbulence spectrum and magnetic field direction and strength have been developed and successfully applied to HI spectroscopic data. To better separate the imprints of density and velocity fluctuations to the channel maps, a new theory-based technique, the Velocity Decomposition Algorit… ▽ More The existence of magnetized turbulence in the interstellar HI is well accepted. A number of techniques to obtain turbulence spectrum and magnetic field direction and strength have been developed and successfully applied to HI spectroscopic data. To better separate the imprints of density and velocity fluctuations to the channel maps, a new theory-based technique, the Velocity Decomposition Algorithm (VDA,Yuen et.al 2021), has been created. The technique demonstrates that the intensity fluctuations are separated into a component pv that mostly arises from velocity fluctuations and pd that mostly arise from density fluctuations. The VDA helps to clarify the nature of the filamentary structure observed in channel maps. A recent publication (Kalberla et.al 2022,K22) claims that the application of VDA to HI4PI data provides negative correlation of pv and pd,which according to the authors invalidates the technique since it requires that pv and pd have zero correlation. However, the quantities pv and pd given by VDA are naturally orthogonal which can be trivially checked analytically or numerically. That means the correct application of the VDA to any data must provide zero correlation. This is the point that we clarify in this paper and search for the cause of the mistake in the application of the VDA in K22 that resulted in the erroneous conclusion. We prove analytically that by construction pv and pd are not correlated. We identify the likely mistake in the VDA expression that K22 used and reproduce their figures with the incorrect expression. We find that 14 out of 15 figures in K22 are invalid, and thus their criticism of the VDA is ill-founded and arises from their use of incorrect expressions. We conclude that the detrimental mistake that K22 made at their analysis completely invalidate their scientific claim that Y21 isnt compatible to observation. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 5 pages,3 figures, Submitted to A&A, Link of GitHub illustration, see https://github.com/kyuen2/kalberla_2022

Decomposing Magnetic Fields in Three Dimensions over the Central Molecular Zone

Authors: Yue Hu, A. Lazarian, Q. Daniel Wang

Abstract: Measuring magnetic fields in the interstellar medium and obtaining their distribution along line-of-sight is very challenging with the traditional techniques. The Velocity Gradient Technique (VGT), which utilizes anisotropy of magnetohydrodynamic (MHD) turbulence, provides an attractive solution. Targeting the central molecular zone (CMZ), we test this approach by applying the VGT to… ▽ More Measuring magnetic fields in the interstellar medium and obtaining their distribution along line-of-sight is very challenging with the traditional techniques. The Velocity Gradient Technique (VGT), which utilizes anisotropy of magnetohydrodynamic (MHD) turbulence, provides an attractive solution. Targeting the central molecular zone (CMZ), we test this approach by applying the VGT to $\rm ^{12}CO$ and $\rm ^{13}CO$ (J = 1-0) data cubes. We first used the SCOUSEPY algorithm to decompose the CO line emissions into separate velocity components, and then we constructed pseudo-Stokes parameters via the VGT to map the plane-of-the-sky magnetic fields in three-dimension. We present the decomposed magnetic field maps and investigate their significance. While the line-of-sight integrated magnetic field orientation is shown to be consistent with the polarized dust emission from the Planck survey at 353 GHz, individual velocity components may exhibit different magnetic fields. We present a scheme of magnetic field configuration in the CMZ based on the decomposed magnetic fields. In particular, we observe a nearly vertical magnetic field orientation in the dense clump near the Sgr B2 and a change in the outflow regions around the Sgr A*. Two high-velocity structures associated with an expanding ring in the CMZ show distinct swirling magnetic field structures. These results demonstrate the potential power of the VGT to decompose velocity or density-dependent magnetic structures. △ Less

Submitted 12 April, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 18 pages, 12 figures, accepted for publication in MNRAS

Damping of Alfvén waves in MHD turbulence and implications for cosmic ray streaming instability and galactic winds

Authors: Alex Lazarian, Siyao Xu

Abstract: Alfvénic component of MHD turbulence damps Alfvénic waves. The consequences of this effect are important for many processes, from cosmic ray (CR) propagation to launching outflows and winds in galaxies and other magnetized systems. We discuss the differences in the damping of the streaming instability by turbulence and the damping of a plane parallel wave. The former takes place in the system of r… ▽ More Alfvénic component of MHD turbulence damps Alfvénic waves. The consequences of this effect are important for many processes, from cosmic ray (CR) propagation to launching outflows and winds in galaxies and other magnetized systems. We discuss the differences in the damping of the streaming instability by turbulence and the damping of a plane parallel wave. The former takes place in the system of reference aligned with the local direction of magnetic field along which CRs stream. The latter is in the reference frame of the mean magnetic field and traditionally considered in plasma studies. We also compare the turbulent damping of streaming instability with ion-neutral collisional damping, which becomes the dominant damping effect at a sufficiently low ionization fraction. Numerical testing and astrophysical implications are also discussed. △ Less

Submitted 13 January, 2022; originally announced January 2022.

Comments: 19 pages, 3 figures, accepted for publication in Frontiers in Astronomy and Space Sciences

Cosmic ray streaming in the turbulent interstellar medium

Authors: Siyao Xu, Alex Lazarian

Abstract: We study the streaming instability of GeV$-100~$GeV cosmic rays (CRs) and its damping in the turbulent interstellar medium (ISM). We find that the damping of streaming instability is dominated by ion-neutral collisional damping in weakly ionized molecular clouds, turbulent damping in the highly ionized warm medium, and nonlinear Landau damping in the Galactic halo. Only in the Galactic halo, is th… ▽ More We study the streaming instability of GeV$-100~$GeV cosmic rays (CRs) and its damping in the turbulent interstellar medium (ISM). We find that the damping of streaming instability is dominated by ion-neutral collisional damping in weakly ionized molecular clouds, turbulent damping in the highly ionized warm medium, and nonlinear Landau damping in the Galactic halo. Only in the Galactic halo, is the streaming speed of CRs close to the Alfvén speed. Alfvénic turbulence plays an important role in both suppressing the streaming instability and regulating the diffusion of streaming CRs via magnetic field line tangling, with the effective mean free path of streaming CRs in the observer frame determined by the Alfvénic scale in super-Alfvénic turbulence. The resulting diffusion coefficient is sensitive to Alfvén Mach number, which has a large range of values in the multi-phase ISM. Super-Alfvénic turbulence contributes to additional confinement of streaming CRs, irrespective of the dominant damping mechanism. △ Less

Submitted 20 January, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 11 pages, 4 figures, accepted for publication in ApJ

Is the recently discovered large scale filamentary feature Cattail in the cold or unstable phase?

Authors: Ka Ho Yuen, Avi Chen, Ka Wai Ho, Alex Lazarian

Abstract: A recent publication (Li et al. 2021) discovered one of the largest filamentary neutral hydrogen features dubbed Cattail from high resolution FAST observations that might be a new galactic arm of our own Milky Way. However in the analysis, it was suggested that this neutral hydrogen feature is cold despite having 12km/s total linewidth. We evaluate the probability whether the Cattail is actually c… ▽ More A recent publication (Li et al. 2021) discovered one of the largest filamentary neutral hydrogen features dubbed Cattail from high resolution FAST observations that might be a new galactic arm of our own Milky Way. However in the analysis, it was suggested that this neutral hydrogen feature is cold despite having 12km/s total linewidth. We evaluate the probability whether the Cattail is actually cold neutral media via the newly developed Velocity Decomposition Algorithm (Yuen et al. 2021a) and Force Balancing Model (Ho et al. 2021a). We discovered that even with the inclusion of the galactic shear term, the feature is still at the unstable neutral media regime. Moreover, we also discover that the Cattail is two disjoint features in caustics space, suggesting that the Cattail might have two different turbulent systems. We check the spectra of the individual system separated via VDA to confirm this argument. We do not exclude the existence of smaller scale cold media being embedded within this structure. △ Less

Submitted 4 December, 2021; originally announced December 2021.

Superdiffusion of Cosmic Rays in Compressible Magnetized Turbulence

Authors: Yue Hu, A. Lazarian, Siyao Xu

Abstract: Owing to the complexity of turbulent magnetic fields, modeling the diffusion of cosmic rays is challenging. Based on the current understanding of anisotropic magnetohydrodynamic (MHD) turbulence, we use test particles to examine the cosmic rays' superdiffusion in the direction perpendicular to the mean magnetic field. By changing Alfven Mach number $M_A$ and sonic Mach number $M_S$ of compressible… ▽ More Owing to the complexity of turbulent magnetic fields, modeling the diffusion of cosmic rays is challenging. Based on the current understanding of anisotropic magnetohydrodynamic (MHD) turbulence, we use test particles to examine the cosmic rays' superdiffusion in the direction perpendicular to the mean magnetic field. By changing Alfven Mach number $M_A$ and sonic Mach number $M_S$ of compressible MHD simulations, our study covers a wide range of astrophysical conditions including subsonic warm gas phase and supersonic cold molecular gas. We show that freely streaming cosmic rays' perpendicular displacement increases as 3/2 to the power of the time traveled along local magnetic field lines. This power-law index changes to 3/4 if the parallel propagation is diffusive. We find that the cosmic rays' parallel mean free path decreases in a power-law relation of $M_A^{-2}$ in supersonic turbulence. We investigate the energy fraction of slow, fast, and Alfvenic modes and confirm the dominance of Alfvenic modes in the perpendicular superdiffusion. In particular, the energy fraction of fast mode, which is the main agent for pitch-angle scattering, increases with $M_A$ but is insensitive to $M_S \ge 2$. Accordingly, our results suggest that the suppressed diffusion in supersonic molecular clouds arises primarily due to the variations of $M_A$ instead of $M_S$. △ Less

Submitted 1 February, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

How the existence of unstable neutral media restricts the aspect ratio of cold neutral media?

Authors: Ka Wai Ho, Ka Ho Yuen, Alex Lazarian

Abstract: The ubiquity of very thin and lengthy cold neutral media (CNM) has been reported by multiple authors in the HI community. Yet, the reason of how the CNM can be so long and lengthy is still in debate. In this paper, we recognize a new type of instability due to the attractive nature of the pressure force in the unstable phase. We provide a new estimation of the average CNM filament aspect ratio wit… ▽ More The ubiquity of very thin and lengthy cold neutral media (CNM) has been reported by multiple authors in the HI community. Yet, the reason of how the CNM can be so long and lengthy is still in debate. In this paper, we recognize a new type of instability due to the attractive nature of the pressure force in the unstable phase. We provide a new estimation of the average CNM filament aspect ratio with the consideration of force balances at the phase boundary, which is roughly 5-20 in common CNM environment. We show that most of the cold filaments are less filamentary than what usually predicted via MHD turbulence theory or inferred from observations: The average length of CNM filament is roughly 1/2 of that in isothermal MHD turbulence with similar turbulence conditions. This suggests that the "cold filaments" that is identified in observations might not be in pressure equilibrium or generated via other mechanisms. △ Less

Submitted 8 February, 2023; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 10 pages, 7 figures, accepted by MNRAS

Shock acceleration with oblique and turbulent magnetic fields

Authors: Siyao Xu, Alex Lazarian

Abstract: We investigate shock acceleration in a realistic astrophysical environment with density inhomogeneities. The turbulence induced by the interaction of the shock precursor with upstream density fluctuations amplifies both upstream and downstream magnetic fields via the turbulent dynamo. The dynamo-amplified turbulent magnetic fields (a) introduce variations of shock obliquities along the shock face,… ▽ More We investigate shock acceleration in a realistic astrophysical environment with density inhomogeneities. The turbulence induced by the interaction of the shock precursor with upstream density fluctuations amplifies both upstream and downstream magnetic fields via the turbulent dynamo. The dynamo-amplified turbulent magnetic fields (a) introduce variations of shock obliquities along the shock face, (b) enable energy gain through a combination of shock drift and diffusive processes, (c) give rise to various spectral indices of accelerated particles, (d) regulate the diffusion of particles both parallel and perpendicular to the magnetic field, and (e) increase the shock acceleration efficiency. Our results demonstrate that upstream density inhomogeneities and dynamo amplification of magnetic fields play an important role in shock acceleration, and thus shock acceleration depends on the condition of the ambient interstellar environment. The implications on understanding radio spectra of supernova remnants are also discussed. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 17 pages, 12 figures, accepted for publication in ApJ

Velocity Gradients: Magnetic Field Tomography towards the Supernova Remnant W44

Authors: Mingrui Liu, Yue Hu, Alex Lazarian

Abstract: As a novel approach for tracing interstellar magnetic fields, the Velocity Gradient Technique (VGT) has been proven to be effective for probing magnetic fields in the diffuse interstellar medium (ISM). In this work, we verify the VGT in a broader context by applying the technique to a molecular cloud interacting with the supernovae remnant (SNR) W44. We probe the magnetic fields with the VGT using… ▽ More As a novel approach for tracing interstellar magnetic fields, the Velocity Gradient Technique (VGT) has been proven to be effective for probing magnetic fields in the diffuse interstellar medium (ISM). In this work, we verify the VGT in a broader context by applying the technique to a molecular cloud interacting with the supernovae remnant (SNR) W44. We probe the magnetic fields with the VGT using CO, $\rm HCO^+$, and H I emission lines and make a comparison with the Planck 353 GHZ dust polarization. We show that the VGT gives an accurate measurement that coheres with the Planck polarization especially in intense molecular gas emission regions. We further study the foreground's contribution to the polarization that results in a misalignment between the VGT and the Planck measurements in low-intensity molecular gas areas. We advance the VGT to achieve magnetic field tomography by decomposing the W44 into various velocity components. We show that W44's velocity component at $v\sim45$ km s$^{-1}$ exhibits the largest coverage and gives the best agreement with Planck polarization in terms of magnetic field orientation. △ Less

Submitted 24 December, 2021; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: 10 pages, 7 figures, accepted for publication in MNRAS