Understanding how star formation is regulated requires studying the energy balance between turbulence, magnetic fields, stellar feedback, and gravity within molecular clouds. However, identifying the transition region where the gravity takes over remains elusive. Recent studies of the Velocity Gradient Technique (VGT), which is an advanced tool for magnetic field studies, reveal that the gradients of spectroscopic observables change their directions by 90 degrees with respect to the magnetic fields in the regions of gravitational collapse. In this study, we perform 3D MHD numerical simulations. We observe that star formation successfully proceeds in strongly magnetized and fully ionized media. We confirm that the self-gravity induces the change of gradients' orientation and high gradients' amplitude. We explore two ways of identifying collapsing self-gravitating regions through the double-peak feature in the histogram of gradients' orientation and the curvature of gradients. We show that velocity gradients' morphology and amplitude can be synthetically used to trace the convergent inflows. By comparing with the column density Probability Density Functions (N-PDFs) method, we show that VGT is a powerful new tool for studying the gas dynamics and tracing magnetic field in star-forming regions. By analogy with VGT, we extend the Intensity Gradient Technique (IGT) to locate the gravitational collapsing region and shocks. We demonstrate the synergy of VGT and IGT can determine the collapsing stages in a star-forming region.

Magnetic field lines in interstellar media have a rich morphology, which could be characterized by geometrical parameters such as curvature and torsion. In this paper, we explore the statistical properties of magnetic field line curvature $\kappa$ in compressible magnetized turbulence. We see that both the mean and standard deviation of magnetic field line curvature obey power-law relations to the magnetization. Moreover, the power-law tail of the curvature probability distribution function is also proportional to the Alfvenic Mach number. We also explore whether the curvature method could be used in the field-tracing Velocity Gradient Technique. In particular, we observe that there is a relation between the mean and standard deviation of the curvature probed by velocity gradients to $M_A$. Finally we discuss how curvature is contributed by different MHD modes in interstellar turbulence, and suggests that the eigenvectors of MHD modes could be possibly represented by the natural Fernet-Serrat frame of the magnetic field lines. We discuss possible theoretical and observational applications of the curvature technique, including the extended understanding on a special length scale that characterize the importance of magnetic field curvature in driving MHD turbulence, and how it could be potentially used to study self-gravitating system.

The recent development of the Velocity Gradient Technique allows observers to map magnetic field orientations and magnetization using the direction of velocity gradients. Aside from the directions, amplitudes of velocity gradients also contain valuable information about the underlying properties of magneto-hydrodynamic (MHD) turbulence. In this paper, we explore what physical information is contained in the amplitudes of velocity gradients and discuss how this information can be used to diagnose properties of turbulence in both diffuse and self-gravitating interstellar media. We identify the relations between amplitudes of both intensity and velocity centroid gradients and the sonic Mach number $M_s$ and they are consistent with the theory's predictions. We test the robustness of the method and discuss how to utilize the amplitudes of gradients into self-gravitating media. To extend the velocity gradient technique we also discuss the usage of amplitude method to Position-Position Velocity (PPV) space as a possible way to retrieve the velocity channel maps before the contamination of thermal broadening. We discuss that the Velocity Gradient Technique with these advancements could potentially give a significantly more accurate statistical insight into the properties magnetized turbulence.

The column density probability density functions (PDFs) of the mass density is a tool to study turbulence and self-gravity in molecular clouds. Self-absorption significantly makes the intensity PDFs different from the column density ones. In this work, we study the effects of self-absorption on the intensity PDFs utilizing three synthetic emission lines of CO isotopologs 12CO(1-0), 13CO(1-0), and C18O(1-0). We find that for supersonic turbulence, the intensity PDFs change their shape from the log-normal to a distribution that has a power-law part depending on the optical depth. We compute the intensity PDFs for massive star-forming region Vela C and find the change of intensity PDFs in observation agrees with our numerical results. We identify the gravitational collapsing regions on Ophiuchus molecular cloud through the PDFs of H2 column density. We see that the change of CO isotopologs' intensity PDFs is independent of self-gravity, which makes the intensity PDFs less reliable in identifying gravitational collapsing regions using emission lines.

Magnetic fields, while ubiquitous in many astrophysical environments, are challenging to measure observationally. Based on the properties of anisotropy of eddies in magnetized turbulence, the Velocity Gradient Technique is a method synergistic to dust polarimetry that is capable of tracing plane-of-the-sky magnetic field, measuring the magnetization of interstellar media and estimating the fraction of gravitational collapsing gas in molecular clouds using spectral line observations. In this paper, we apply this technique to five low-mass star-forming molecular clouds in the Gould Belt and compare the results to the magnetic-field orientation obtained from polarized dust emission. We find the estimates of magnetic field orientations and magnetization for both methods are statistically similar. We estimate the fraction of collapsing gas in the selected clouds. By means of the Velocity Gradient Technique, we also present the plane-of-the-sky magnetic field orientation and magnetization of the Smith cloud, for which dust polarimetry data are unavailable.

We introduce two new ways of obtaining the strength of plane-of-sky (POS) magnetic field by simultaneous use of spectroscopic Doppler-shifted lines and the information on magnetic field direction. The latter can be obtained either through polarization measurements or using the velocity gradient technique. We show the advantages that our techniques have compared to the traditional Davis-Chandrasekhar-Fermi (DCF) technique of estimating magnetic field strength from observations. The first technique that we describe in detail employs structure functions of velocity centroids and structure functions of Stokes parameters. We provide analytical expressions for obtaining magnetic field strength from observational data. We successfully test our results using synthetic observations obtained with results of MHD turbulence simulations. We measure velocity and magnetic field fluctuations at small scales using two, three and four point structure functions and compare the performance of these tools. We show that, unlike the DCF, our technique is capable of providing the detailed distribution of POS magnetic field and it can measure magnetic field strength in the presence of both velocity and magnetic field distortions arising from external shear and self-gravity. The second technique applies the velocity gradient technique to velocity channel maps in order to obtain the Alfven Mach number and uses the amplitudes of the gradients to obtain the sonic Mach number. The ratio of these two Mach numbers provides the intensity of magnetic field in the region contributing to the emission in the channel map. We test the technique and discuss obtaining the 3D distribution of POS galactic Magnetic field with it. We discuss the application of the second technique to synchrotron data.

The observations of fluctuations in the cosmic microwave background provide information about primordial inhomogeneities in the universe. However, the B-mode polarization of the inflationary gravitational wave is contaminated by the Galactic foreground polarized radiation arising from dust aligned by interstellar magnetic fields. To trace magnetic fields we use the Velocity Gradient Technique, which employs a modern understanding of the nature of magneto-hydrodynamics turbulent motions. In this paper, we combine the VGT with the Principal Component Analysis (PCA) to improve the accuracy of magnetic field tracing. We apply the VGT-PCA to the high-resolution neutral hydrogen data from the GALFA-HI survey to predict the polarization of dust. We report that the predicted directions of dust polarization provide good correspondence with those reported by Planck 353GHz, the alignment measure between the two AM=0.79. We show that our results statistically agree with the Planck polarization in terms of magnetic field tracing. We find that the variation of dust emission efficiency across the sky is small. Using our maps of predicted polarization we calculate the ratio of the E and B modes and show that BB/EE=0.53, which is similar to the result from Plank polarization.

Magnetic fields are ubiquitous in the interstellar medium but are notoriously difficult to study through observation. Making use of the advances in our understanding of MHD turbulence and turbulent reconnection, the Velocity Gradient Technique (VGT) was suggested and successfully applied to study magnetic fields utilizing spectroscopic data. Applying the tools developed for VGT to intensity statistics, we introduce the Intensity Gradients Technique (IGT) as a complementary tool that can be used synergistically with VGT. In this paper, we apply IGT to a diffuse HI region selected from the GALFA-HI survey and compare the intensity gradient maps with those obtained using velocity gradients as well as Planck polarization measurements. We demonstrate the possibility of using IGT and VGT for both studying the magnetic field and identifying shocks in the diffuse interstellar medium. We also explore the ability of IGT in locating self-gravitating regions and calculating Alfvenic Mach number, both alone and in combination with VGT and polarimetry. We compare IGT with the Histogram of Relative Orientation (HRO), which utilizes intensity gradients to characterize the relative orientation of column density structures and local magnetic fields.

Probing magnetic fields in Giant Molecular Clouds is often challenging. Fortunately, recently simulations show that analysis of velocity gradients (the Velocity Gradient Technique) can be used to map out the magnetic field morphology of different physical layers within molecular clouds when applied CO isotopologues with different optical depths. Here, we test the effectiveness of the Velocity Gradient Technique in reconstructing the magnetic field structure of the molecular cloud Vela C, employing seven chemical tracers that have different optical depths, i.e. 12CO, 13CO, C18O, CS, HNC, HCO+, and HCN. Our results show good correspondence between the magnetic field morphology inferred from velocity gradients using these different molecular tracers and the magnetic field morphology inferred from BLASTPol polarization observations. We also explore the possibility of using a combination of velocity gradients for multiple chemical tracers to explain the structure of the magnetic field in molecular clouds. We search for signatures of gravitational collapse in the alignment of the velocity gradients and magnetic field and conclude that collapsing regions constitute a small fraction of the cloud.

We employ synthetic observations obtained with MHD simulations to study how to trace the distribution of turbulent magnetic fields using the synchrotron polarization gradient techniques suggested by Lazarian \& Yuen (2018). Both synchrotron polarization gradients and its derivative gradients with regard to the squared wavelength $\lambda^2$ are used to explore the probing ability of the technique in magnetic fields from sub- to super-Alfv{\'e}nic turbulence. We focus on studies that involve multi-frequency measurements in the presence of strong Faraday rotation and show the ways of how to recover the projected mean magnetic fields in the plane of the sky and the local magnetic fields within a tomographic slice. We conclude that the new techniques can successfully reconstruct the 3D magnetic field within our Milky Way and other galaxies. This paper opens an avenue for applying our new techniques to a large number of data cubes such as those from the Low Frequency Array for Radio astronomy and the Square Kilometer Array.

A recent publication by Clark et.al (2019, CX19) uses both GALFA-HI observational data and numerical simulations to address the nature of intensity fluctuations in Position-Position-Velocity (PPV) space. The study questions the validity and applicability of the statistical theory of PPV space fluctuations formulated in Lazarian & Pogosyan (2000, LP00) to HI gas and concludes that $"${\it a significant reassessment of many observational and theoretical studies of turbulence in HI}$"$. This implies that dozens of papers that used LP00 theory to explore interstellar turbulence as well as the ongoing research based on LP00 theory are in error. This situation motivates the urgency of our public response. In our Comment we explain why we believe the criticism in CX19 is based on the incorrect understanding of the LP00 theory. In particular, we illustrate that the correlation between PPV slices and dust emissions in CX19 does not properly reveal the relative importance of velocity and density fluctuations in velocity channel maps. While CX19 provides an explanation of the change of the spectral index with respect to the thickness of PPV slice based on the two-phase nature of H1 gas, we failed to see any observational support for this idea. On the contrary, we show that the observations both in two-phase HI and one phase CO show similar results. Moreover, the observed change is in good agreement with LP00 predictions and spectral indexes of velocity and density spectra that are obtained following LP00 procedures are in good agreement with the numerically confirmed expectations of compressible MHD turbulence theory. In short, we could not find any justification of the criticism of LP00 theory that is provided in CX19. On the contrary, our analysis testifies that both available observational and numerical data agree well with the predictions of LP00 theory.

Observations of synchrotron polarization at multiple frequencies in the presence of Faraday rotation can provide a way to reconstruct the 3D magnetic field distribution. In this paper we compare the well known Faraday Tomography (FT, Burn 1966) technique to the new approach named Synchrotron Polarization Gradients (SPG, Lazarian & Yuen 2018b). We compare the strengths and limitations of the two techniques, and describe their synergy. In particular, we show that in the situations when FT technique fails, e.g. due to insufficient frequency coverage, the SPG can still trace the 3D structure of magnetic field.

Probing magnetic fields in self-gravitating molecular clouds are generally difficult even with the use of the polarimetry. Based on the properties of magneto-hydrodynamic (MHD) turbulence and turbulent reconnection, Velocity Gradient Technique (VGT) provides a new way in tracing magnetic field orientation and strength based on the spectroscopic data. Our study tests the applicability of VGT in various molecular tracers, e.g. 12CO, 13CO, and C18O. By inspecting synthetic molecular line maps of CO isotopologue generated through radiative transfer calculations, we show that the VGT method can be successfully applied in probing the magnetic field direction in the diffuse interstellar medium as well as in self-gravitating molecular clouds.

Magnetohydrodynamic(MHD) turbulence displays velocity anisotropies which reflect the direction of the magnetic field. This anisotropy has led to the development of a number of statistical techniques for studying magnetic fields in the interstellar medium. In this paper, we review and compare three techniques that use radio position-position-velocity data for determining magnetic field strength and morphology : the correlation function anisotropy (CFA), Principal Component Analysis of Anisotropies (PCAA), and the more recent Velocity Gradient Technique (VGT). We compare these three techniques and suggest improvements to the CFA and PCAA techniques to increase their accuracy and versatility. In particular, we suggest and successfully implement a much faster way of calculating non-periodic correlation functions for the CFA. We discuss possible improvements to the current implementation of the PCAA. We show the advantages of the VGT in terms of magnetic field tracing and stress the complementary nature with the other two techniques.

Recent developments of the Velocity Gradient Technique (VGT) show that the velocity gradients provide a reliable tracing of magnetic field direction in turbulent plasmas. In this paper, we explore the ability of velocity gradients to measure the magnetization of interstellar medium. We demonstrate that the distribution of velocity gradient orientations provides a reliable estimation of the magnetization of the media. In particular, we determine the relation between Alfvenic Mach number $M_A$ in the range of $M_A \in [0.2,1.7]$ and properties of the velocity gradient distribution, namely, with the dispersion of velocity gradient orientation as well as with the peak to base ratio of the amplitudes. We apply our technique for a selected GALFA-HI region and find the results consistent with the expected behavior of $M_A$. Using 3D MHD simulations we successfully compare the results with our new measure of magnetization that is based on the dispersion of starlight polarization. We demonstrate that, combined with the velocity dispersion along the line of sight direction, our technique is capable to delivering the magnetic field strength. The new technique opens a way to measure magnetization using other gradient measures such as synchrotron intensity gradients (SIGs) and synchrotron polarization gradients (SPGs).

We describe a new technique for probing galactic and extragalactic 2D and 3D magnetic field distribution using gradients of polarized synchrotron emission. The fluctuations of magnetic field are elongated along the ambient magnetic field. Therefore, the field variations are maximal perpendicular to the B-field. This allows tracing B-field with synchrotron polarization gradients (SPGs). We demonstrate that the Faraday depolarization allows to map 3D B-field structure by. The depolarization ensures that the polarization gradients sample the regions close to the observer with the sampling depth controlled by the frequency of radiation. We also analyze the B-field properties along the line-of-sight by applying the gradient technique to the wavelength derivative of synchrotron polarization. This Synchrotron Derivative Polarization Gradients (SDPGs) technique can recover the 3D vectors of the underlying B-fields. The new techniques are different from the Faraday tomography as they provide a way to map the 3D distribution of B-fields components perpendicular to the line of sight. In addition, we find that the alignment of gradients of polarization with the synchrotron polarization can be used to separate the contribution of the foreground from the polarization of cosmological origin. We notice that the same alignment is also present for the dust polarization.

Tracing of the magnetic field with Velocity Gradient Technique (VGT) allows observers to probe magnetic field directions with spectroscopic data. In this paper, we employ the method of Principal Component Analysis (PCA) to extract the spectroscopic information most valuable for VGT. By using synthetic observation data from numerical simulations, we show that PCA acts in a way similar to spatial filtering along the velocity axis. We study both subsonic and supersonic simulations and show that with the PCA filtering the tracing of magnetic fields by the VGT is significantly improved. Using 21 cm GALFA data, we demonstrate that the PCA filtering improves the alignment of the velocity gradients and the Planck dust polarization.

On the basis of the modern understanding of MHD turbulence, we propose a new way of using synchrotron radiation, namely using synchrotron intensity gradients for tracing astrophysical magnetic fields. We successfully test the new technique using synthetic data obtained with the 3D MHD simulations and provide the demonstration of the practical utility of the technique by comparing the directions of magnetic field that are obtained with PLANCK synchrotron intensity dats to the directions obtained with PLANCK synchrotron polarization data. We demonstrate that the synchrotron intensity gradients (SIGs) can reliably trace magnetic field in the presence of noise and can provide detailed maps of magnetic-field directions. We also show that the SIGs are relatively robust for tracing magnetic fields while the low spacial frequencies of the synchrotron image are removed. This makes the SIGs applicable to tracing of magnetic fields using interferometric data with single dish measurement absent. We discuss the synergy of using the SIGs together with synchrotron polarization in order to find the actual direction of the magnetic field and quantify the effects of Faraday rotation as well as with other ways of studying astrophysical magnetic fields. We test our method in the presence of noise and the resolution effects. We stress the complementary nature of the studies using the SIG technique and those employing the recently-introduced velocity gradient techniques that traces the magnetic fields using spectroscopic data.

We identify velocity channel map intensities as a new way to trace magnetic fields in turbulent media. This work makes use of both of the modern theory of MHD turbulence predicting that the magnetic eddies are aligned with the local direction of magnetic field, and also the theory of the spectral line Position-Position-Velocity (PPV) statistics describing how the velocity and density fluctuations are mapped into the PPV space. In particular, we use the fact that, the fluctuations of intensity of thin channel maps are mostly affected by the turbulent velocity, while the thick maps are dominated by density variations. We study how contributions of the fundamental MHD modes affect the Velocity Channel Gradients (VChGs) and demonstrate that the VChGs arising from Alfven and slow modes are aligned perpendicular to the local direction of magnetic field, while the VChGs produced by the fast mode are aligned parallel to the magnetic field. The dominance of the Alfven and slow modes in the interstellar media will therefore allow a reliable magnetic field tracing using the VChGs. We explore ways of identifying self-gravitating regions that do not require polarimetric information. In addition, we also introduce a new measure termed "Reduced Velocity Centroids" (RVCGs) and compare its abilities with the VChGs. We employed both measures to the GALFA 21cm data and successfully compared thus magnetic field directions with the {\it PLANCK} polarization. The applications of the suggested techniques include both tracing magnetic field in diffuse interstellar media and star forming regions, as well as removing the galactic foreground in the framework of cosmological polarization studies.

This study proceeds with the development of the technique employing velocity gradients that were identified in (\cite{GL17}, henceforth GL17) as a means of probing magnetic field in interstellar media. We demonstrate a number of practical ways on improving the accuracy of tracing magnetic fields in diffuse interstellar media using velocity centroid gradients (VCGs). Addressing the magnetic field tracing in super-Alfvenic turbulence we introduce the procedure of filtering low spatial frequencies, that enables magnetic field tracing in the situations when the kinetic energy of turbulent plasmas dominate its magnetic energy. We propose the synergic way of of using VCGs together with intensity gradients (IGs), synchrotron intensity gradients (SIGs) as well as dust polarimetry. We show that while the IGs trace magnetic field worse than the VCGs, the deviations of the angle between the IGs and VCGs trace the shocks in diffuse media. Similarly the perpendicular orientation of the VCGs and the SIGs or to the dust polarimetry data trace the regions of gravitational collapse. We demonstrate the utility of combining the VCGs, IGs and polarimetry using GALFA HI and Planck polarimetry data. We also provide an example of synergy of the VCGs and the SIGs using the HI4PI full-sky HI survey together with the Planck synchrotron intensity data.

The advancement of our understanding of MHD turbulence opens ways to develop new techniques to probe magnetic fields. In MHD turbulence, the velocity gradients are expected to be perpendicular to magnetic fields and this fact was used by Gonsalvez-Casanova & Lazarian to introduce a new technique to trace magnetic fields using velocity centroid gradients. The latter can be obtained from spectroscopic observations. We apply the technique to GALFA HI survey data and compare the directions of magnetic fields obtained with our technique with the direction of magnetic fields obtained using PLANCK polarization. We find excellent correspondence between the two ways of magnetic field tracing, which is obvious via visual comparison and through measuring of the statistics of magnetic field fluctuations obtained with the polarization data and our technique. This suggests that the velocity centroid gradients can provide a reliable way of measuring of the foreground magnetic field fluctuations and thus provide a new way of separating foreground and CMB polarization signals.

Most molecular clouds are filamentary or elongated. Among those forming low-mass stars, their long axes tend to be either parallel or perpendicular to the large-scale (10-100 pc) magnetic field (B-field) in the surrounding inter cloud medium. This arises because, along the dynamically dominant B-fields, the competition between self-gravity and turbulent pressure will shape the cloud to be elongated either perpendicular or parallel to the fields. Recent study also suggested that, on the scales of 0.1-0.01 pc, fields are dynamically important within cloud cores forming massive stars. But whether the core field morphologies are inherited from the inter cloud medium or governed by cloud turbulence is under vigorous debate, so is the role played by B-fields in cloud fragmentation at 10 - 0.1 pc scales. Here we report B-field maps covering 100-0.01 pc scales inferred from polarimetric observations of a massive-star forming region, NGC 6334. First, the main filament also lies perpendicular to the ambient field. NGC 6334 hosts young star-forming sites where fields are not severely affected by stellar feedback, and their directions do not change significantly over the entire scale range. This means that the fields are dynamically important. At various scales, we find that the hourglass-shaped field lines are pinched where the gas column density peaks and the field strength is proportional to the 0.4-power of the density. We conclude that B-fields play a crucial role in the fragmentation of NGC 6334.