Applying the Velocity Gradient Technique in NGC 1333: Comparison with Dust Polarization Observations

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the Velocity Gradient Technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC~1333 star-forming region ($\rm ^{12}CO$, $\rm ^{13}CO$, $\rm C^{18}O$, and $\rm N_{2}H^{+}$), and we compare the derived B-field properties with those inferred from 214 and 850~$\mu$m dust polarization observations of the region using SOFIA/HAWC+ and JCMT/POL-2, respectively. We estimate both the inclination angle and the 3D Alfv\'enic Mach Number $M_A$ from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelength, such as the plane-of-the-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM.