Researcher profile

A. Lazarian

A. Lazarian contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
24works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

Decomposing Magnetic Fields in Three Dimensions over the Central Molecular Zone

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Multi-scale Magnetic Fields in the Central Molecular Zone: Inference from the Gradient Technique

The central molecular zone (CMZ) plays an essential role in regulating the nuclear ecosystem of our Galaxy. To get an insight into the magnetic fields of the CMZ, we employ the Gradient Technique (GT), which is rooted in the anisotropy of magnetohydrodynamic turbulence. Our analysis is based on the data of multiple wavelengths, including molecular emission lines, radio 1.4 GHz continuum image, and Herschel 70 $μ$m image, as well as ionized [Ne II] and Paschen-alpha emissions. The results are compared with the observations of Planck 353 GHz and High-resolution Airborne Wideband Camera Plus (HWAC+) 53 $μ$m polarized dust emissions. We map the orientation of the magnetic field at multiple wavelengths across the central molecular zone, including close-ups of the Radio Arc and Sagittarius A West regions, on multi scales from $\sim$ 0.1 pc to 10 pc. The magnetic fields towards the central molecular zone traced by GT are globally compatible with the polarization measurements, accounting for the contribution from the galactic foreground and background. This correspondence suggests that the magnetic field and turbulence are dynamically crucial in the galactic center. We find that the magnetic fields associated with the Arched filaments and the thermal components of the Radio Arc are in good agreement with the HAWC+ polarization. Our measurement towards the non-thermal Radio Arc reveals the poloidal magnetic field components in the galactic center. For Sagittarius A West region, we find a great agreement between the GT measurement using [Ne II] emission and HWAC+ 53 $μ$m observation. We use GT to predict the magnetic fields associated with ionized Paschen-alpha gas down to scales of 0.1 pc.

preprint2022arXiv

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

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$.

preprint2022arXiv

Superdiffusion of Cosmic Rays in Compressible Magnetized Turbulence

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$.

preprint2021arXiv

A Synergy of the Velocity Gradients Technique and the Probability Density Functions for Identifying Gravitational Collapse in Self-Absorbing Media

The Velocity Gradients Technique (VGT) and the Probability Density Functions (PDFs) of mass density are tools to study turbulence, magnetic fields, and self-gravity in molecular clouds. However, self-absorption can significantly make the observed intensity different from the column density structures. In this work, we study the effects of self-absorption on the VGT and the intensity PDFs utilizing three synthetic emission lines of CO isotopologs $^{12}$CO (1-0), $^{13}$CO (1-0), and C$^{18}$O (1-0). We confirm that the performance of VGT is insensitive to the radiative transfer effect. We numerically show the possibility of constructing 3D magnetic fields tomography through VGT. We find that the intensity PDFs change their shape from the pure log-normal to a distribution that exhibits a power-law tail depending on the optical depth for supersonic turbulence. We conclude the change of CO isotopologs' intensity PDFs can be independent of self-gravity, which makes the intensity PDFs less reliable in identifying gravitational collapsing regions. We compute the intensity PDFs for a star-forming region NGC 1333 and find the change of intensity PDFs in observation agrees with our numerical results. The synergy of VGT and the column density PDFs confirms that the self-gravitating gas occupies a large volume in NGC 1333.

preprint2021arXiv

Diffusion of large-scale magnetic fields by reconnection in MHD turbulence

The rate of magnetic field diffusion plays an essential role in several astrophysical plasma processes. It has been demonstrated that the omnipresent turbulence in astrophysical media induces fast magnetic reconnection, which consequently leads to large-scale magnetic flux diffusion at a rate independent of the plasma microphysics. This process is called ``reconnection diffusion&#39;&#39; (RD) and allows for the diffusion of fields which are dynamically important. The current theory describing RD is based on incompressible magnetohydrodynamic (MHD) turbulence. In this work, we have tested quantitatively the predictions of the RD theory when magnetic forces are dominant in the turbulence dynamics (Alfvénic Mach number $M_A < 1$). We employed the \textsc{Pencil Code} to perform numerical simulations of forced MHD turbulence, extracting the values of the diffusion coefficient $η_{RD}$ using the Test-Field method. Our results are consistent with the RD theory ($η_{RD} \sim M_A^{3}$ for $M_A < 1$) when turbulence approaches the incompressible limit (sonic Mach number $M_S \lesssim 0.02$), while for larger $M_S$ the diffusion is faster ($η_{RD} \sim M_A^{2}$). This work shows for the first time simulations of compressible MHD turbulence with the suppression of the cascade in the direction parallel to the mean magnetic field, which is consistent with incompressible weak turbulence theory. We also verified that in our simulations the energy cascading time does not follow the scaling with $M_A$ predicted for the weak regime, in contradiction with the RD theory assumption. Our results generally support and expand the RD theory predictions.

preprint2021arXiv

Large-scale magnetic field in the Monoceros OB-1 East molecular cloud

We study the large-scale magnetic field structure and its interplay with the gas dynamics in the Monoceros OB1 East molecular cloud. We combine observations of dust polarised emission from the Planck telescope and CO molecular line emission observations from the Taeduk Radio Astronomy Observatory 14-metre telescope. We calculate the strength of the plane-of-the-sky magnetic field using a modified Chandrasekhar-Fermi method and estimate mass over flux ratios in different regions of the cloud. We use the comparison of the velocity and intensity gradients of the molecular line observations with the polarimetric observations to trace dynamically active regions. The molecular complex shows an ordered large-scale plane-of-the-sky magnetic field structure. In the Northern part, it is mostly orientated along the filamentary structures while the Southern part shows at least two regions with distinct magnetic field orientations. We find that in the Northern filaments the magnetic field is unlikely to provide support against fragmentation at large scales. Our analysis reveals a shock region in the Northern part of the complex right in-between two filamentary clouds which were previously suggested to be in collision. Moreover, the shock seems to extend farther towards the Western part of the complex. In the Southern part, we find that either the magnetic field guides the accretion of interstellar matter towards the cloud or it was dragged by the matter towards the densest regions. The large-scale magnetic field in Monoceros OB-1 East molecular clouds is tightly connected to the global structure of the complex and, in the Northern part, it seems to be dominated by gravity and turbulence, while in the Southern part it influences the structuring of matter.

preprint2020arXiv

Advancing the Velocity Gradient Technique: Using Gradient Amplitudes and handling thermal broadening

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&#39;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.

preprint2020arXiv

Alignment and rotational disruption of dust

We reveal a deep connection between alignment of dust grains by RAdiative torques (RATs) and MEchanical Torques (METs) and rotational disruption of grains introduced by Hoang et al. (2019). The disruption of grains happens if they have attractor points corresponding to high angular momentum (high-J). We introduce {\it fast disruption} for grains that are directly driven to the high-J attractor on a timescale of spin-up, and {\it slow disruption} for grains that are first moved to the low-J attractor and gradually transported to the high-J attractor by gas collisions. The enhancement of grain magnetic susceptibility via iron inclusions expands the parameter space for high-J attractors and increases percentage of grains experiencing the disruption. The increase in the magnitude of RATs or METs can increase the efficiency of fast disruption, but counter-intuitively, decreases the effect of slow disruption by forcing grains towards low-J attractors, whereas the increase in gas density accelerates disruption by faster transporting grains to the high-J attractor. We also show that disruption induced by RATs and METs depends on the angle between the magnetic field and the anisotropic flow. We find that pinwheel torques can increase the efficiency of {\it fast disruption} but may decrease the efficiency of {\it slow disruption} by delaying the transport of grains from the low-J to high-J attractors via gas collisions. The selective nature of the rotational disruption opens a possibility of observational testing of grain composition as well as physical processes of grain alignment.

preprint2020arXiv

Cosmic ray transport in starburst galaxies

Starburst galaxies are efficient $γ$-ray producers, because their high supernova rates generate copious cosmic ray (CR) protons, and their high gas densities act as thick targets off which these protons can produce neutral pions and thence $γ$-rays. In this paper we present a first-principles calculation of the mechanisms by which CRs propagate through such environments, combining astrochemical models with analysis of turbulence in weakly ionised plasma. We show that CRs cannot scatter off the strong large-scale turbulence found in starbursts, because efficient ion-neutral damping prevents such turbulence from cascading down to the scales of CR gyroradii. Instead, CRs stream along field lines at a rate determined by the competition between streaming instability and ion-neutral damping, leading to transport via a process of field line random walk. This results in an effective diffusion coefficient that is nearly energy-independent up to CR energies of $\sim 1$ TeV. We apply our computed diffusion coefficient to a simple model of CR escape and loss, and show that the resulting $γ$-ray spectra are in good agreement with the observed spectra of the starbursts NGC 253, M82, and Arp 220. In particular, our model reproduces these galaxies&#39; relatively hard GeV $γ$-ray spectra and softer TeV spectra without the need for any fine-tuning of advective escape times or the shape of the CR injection spectrum.

preprint2020arXiv

Curvature of magnetic field lines in compressible magnetized turbulence: Statistics, magnetization predictions, gradient curvature, modes and self-gravitating media

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 $κ$ 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.

preprint2020arXiv

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

preprint2020arXiv

Modeling of Galactic Foreground Polarization with Velocity Gradients

The detection of primordial B-mode polarization is still challenging due to the relatively low amplitude compared to the galactic foregrounds. To remove the contribution from the foreground, a comprehensive picture of the galactic magnetic field is indispensable. The Velocity Gradient Technique (VGT) is promising in tracing magnetic fields based on the modern understanding of the magneto-hydrodynamic turbulence. In this work, we apply VGT to a HI region containing an intermediate velocity cloud and a local velocity cloud, which are distinguishable in position-position-velocity space. We show that VGT gives an excellent agreement with the Planck polarization and stellar polarization. We confirm the advantages of VGT in constructing the 3D galactic magnetic field.

preprint2020arXiv

Onset of turbulent fast magnetic reconnection observed in the solar atmosphere

Fast magnetic reconnection powers explosive events throughout the universe, from gamma-ray bursts to solar flares. Despite its importance, the onset of astrophysical fast reconnection is the subject of intense debate and remains an open question in plasma physics. Here we report high-cadence observations of two reconnection-driven solar microflares obtained by the \textit{Interface Region Imaging Spectrograph} that show persistent turbulent flows preceding flaring. The speeds of these flows are comparable to the local sound speed initially, suggesting the onset of fast reconnection in a highly turbulent plasma environment. Our results are in close quantitative agreement with the theory of turbulence-driven reconnection as well as with numerical simulations in which fast magnetic reconnection is induced by turbulence.

preprint2020arXiv

Study Turbulence and Probe Magnetic Field Using Gradients Technique: Application to HI-to-H2 Transition Regions

The atomic-to-molecular (HI-to-H$_2$) transition in photodissociation regions (PDRs) has been investigated over the last several decades through analytic and numerical modeling. However, classical PDR models typically assume uniform density gas, ignoring the turbulent nature of the interstellar medium. Recently, Bialy et al. (2017b, 2019) have presented a theoretical framework for studying the HI-to-H$_2$ in a realistic turbulent medium with a non-homogeneous density structure. Here we extend these turbulent-chemical models to explore the possibility of tracing the magnetic field direction in turbulent PDRs using the Gradients Technique. We utilize both subsonic and supersonic magnetohydrodynamic numerical simulations for chemical HI/H$_2$ balance calculations. We confirm that the density fluctuations induced by turbulence can disperse the distribution of H$_2$ and HI fraction. We find that the energy spectrum of moment maps gets shallower when the sonic Mach number MS increases. We explore the ability in magnetic field tracing of gradients of higher-order velocity centroids and compare their performance with that of traditional velocity centroid gradients (VCGs) and with intensity gradients (IGs). We find that the velocity gradients of the second-order centroids (VC$_2$Gs) are more accurate than VCGs and IGs in probing the magnetic field orientation.

preprint2020arXiv

Velocity Centroids Anisotropy and the Signature of different MHD Modes in the Turbulent ISM

Magnetic turbulence is anisotropic as the directions of motion are constrained by the magnetic field. Such anisotropy can be observed in velocity centroids obtained from spectroscopic observations. We use magnetohydrodynamics (MHD) simulations to produce synthetic spectroscopic observations (position-position-velocity data) and study the anisotropy in the structure function of velocity centroid maps. We decomposed the velocity in the simulations into Alfvén, slow and fast-modes and studied how each of them contribute to the observed anisotropy. We found that when the angle between the line of sight and the mean magnetic field is large the Alfvén-mode dominates the observed anisotropy, while for smaller angles the anisotropy is not large enough to be used to probe the magnetization of the media, and it is dominated by the slow-mode. Our results are in fair agreement with the theoretical predictions in Kandel et al.(2016,2017).

preprint2019arXiv

A comparison between Faraday Tomography and Synchrotron Polarization Gradients

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.

preprint2019arXiv

Producing Synthetic Maps of Dust Polarization using Velocity Channel Gradient Technique

In modern cosmology, many efforts have been put to detect primordial B-mode of cosmic microwave background (CMB) polarization from the gravitational waves generated during inflation. Considering the foreground dust contamination of microwave polarization maps, it is essential to obtain a precise prediction for polarization in dust emission. In this work, we show a new method to produce synthetic maps of dust polarization in magnetized turbulent ISM from more abundant high-resolution HI data. By using Velocity Channel Gradient (VChG) technique, we are able to predict both direction and degree of dust polarization by investigating spectroscopic HI information in position-position-velocity (PPV) space. We applied our approach to The Galactic Arecibo L-band feed Array HI (GALFA-HI) data, and find a good correspondence between synthesized maps and PLANCK&#39;s polarization measurements at 353 GHz.

preprint2010arXiv

Extending Big Power Law in the Sky with Turbulence Spectra from WHAM data

We use the data of Wisconsin H$α$ Mapper (WHAM) to test the hypothesis of whether the amplitudes and spectrum of density fluctuations measured by WHAM can be matched to the data obtained for interstellar scintillations and scattering. To do this, first of all, we adjusted the mean level of signal in the adjacent patches of the data. Then, assuming that the spectrum is Kolmogorov, we successfully matched the amplitudes of turbulence obtained from the WHAM data and the interstellar density fluctuations reported in the existing literature. As a result, we conclude that the existing data is consistent with the Kolmogorov cascade which spans from $10^6$ to $10^{17}$ $m$.

preprint2010arXiv

Reconnection Diffusion and Star Formation

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology or reconnect in highly conducting media. Recent progress in understanding fast magnetic reconnection in the presence of turbulence is reassuring that the magnetic field behavior in computer simulations and turbulent astrophysical environments is similar, as far as the magnetic reconnection is concerned. This makes it meaningful to perform MHD simulations of turbulent flows in order to understand the diffusion of magnetic field in astrophysical environments. These simulations support the concept of reconnection diffusion, which describes the ability of magnetic fields to get removed from magnetized clouds and cores in the process of star formation.

preprint2010arXiv

Turbulent magnetic reconnection in 2D and 3D

Magnetic field embedded in a perfectly conducting fluid preserves its topology for all time. Although ionized astrophysical objects, like stars and galactic disks, are almost perfectly conducting, they show indications of changes in topology, `magnetic reconnection&#39;, on dynamical time scales. Reconnection can be observed directly in the solar corona, but can also be inferred from the existence of large scale dynamo activity inside stellar interiors. Solar flares and gamma ray busts are usually associated with magnetic reconnection. Previous work has concentrated on showing how reconnection can be rapid in plasmas with very small collision rates. Here we present numerical evidence, based on three dimensional simulations, that reconnection in a turbulent fluid occurs at a speed comparable to the rms velocity of the turbulence, regardless of the value of the resistivity. In particular, this is true for turbulent pressures much weaker than the magnetic field pressure so that the magnetic field lines are only slightly bent by the turbulence. These results are consistent with the proposal by Lazarian and Vishniac (1999) that reconnection is controlled by the stochastic diffusion of magnetic field lines, which produces a broad outflow of plasma from the reconnection zone. This work implies that reconnection in a turbulent fluid typically takes place in approximately a single eddy turnover time, with broad implications for dynamo activity and particle acceleration throughout the universe. In contrast, the reconnection in 2D configurations in the presence of turbulence depends on resistivity, i.e. is slow.