Source author record

Alex Lazarian

Alex Lazarian appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

55works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

55 published item(s)

preprint2026arXiv

Unveiling the 3D structure of the central molecular zone from stellar kinematics and photometry: The 50 and 20 km/s clouds

The central molecular zone (CMZ), surrounding the Galactic centre, is the largest reservoir of dense molecular gas in the Galaxy. Despite its relative proximity, the 3D structure of the CMZ remains poorly constrained, primarily due to projection effects. We aim to constrain the line-of-sight location of two molecular clouds in the CMZ -- the 50 and 20 km/s clouds -- and to investigate their possible physical connection using stellar kinematics and photometry. This study serves as a pilot for future applications across the full CMZ. We estimated the line-of-sight position of the clouds by analysing stellar kinematics, stellar densities, and stellar populations towards the cloud regions and a control field. We find an absence of westward moving stars in the cloud regions, which indicates that they lie on the near side of the CMZ. This interpretation is supported by the stellar density distributions. The similar behaviour observed in the two clouds, as well as in the region between them (the ridge), suggests that they are located at comparable distances and are physically linked. We also identified an intermediate-age stellar population (2-7 Gyr) in both regions, consistent with that observed on the near side of the CMZ. We estimated the line-of-sight distances at which the clouds and the ridge become kinematically detectable (i.e. where the proper motion component parallel to the Galactic plane differs from that of the control field at the 3 sigma level) by converting their measured proper motions parallel to the Galactic plane using a theoretical model of the stellar distribution. We find that the 50 and 20 km/s clouds are located at $43\pm8$ pc and $56\pm11$ pc from Sgr A*, respectively, and that the ridge lies at $56\pm11$ pc; this supports the idea that the clouds are physically connected through the ridge.

preprint2022arXiv

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

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.

preprint2022arXiv

Cosmic ray streaming in the turbulent interstellar medium

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Turbulent reconnection acceleration

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.

preprint2022arXiv

Turbulent universal galactic Kolmogorov velocity cascade over 6 decades

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.

preprint2021arXiv

Shock acceleration with oblique and turbulent magnetic fields

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.

preprint2021arXiv

Velocity Gradients: Magnetic Field Tomography towards the Supernova Remnant W44

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.

preprint2020arXiv

3D Turbulent Reconnection: Theory, Tests and Astrophysical Implications

Magnetic reconnection, topological change in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence magnetic lines constantly reconnect everywhere at all scales, making magnetic reconnection an intrinsic part of turbulent cascade. We note that this is inconsistent with the usual practice of regarding magnetic lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies, starting with Lazarian & Vishniac (1999), demonstrated that 3D turbulence makes magnetic reconnection fast and that these two processes are intrinsically connected. We discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence and demonstrate that in the presence of turbulence the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. This justifies an MHD-like treatment of magnetic reconnection at scales much larger than the relevant plasma scales. We discuss numerical and observational evidences supporting the turbulent reconnection model. In particular, we show that tearing reconnection is suppressed in 3D and, unlike the 2D case, the 3D reconnection induces turbulence that makes reconnection independent of resistivity. We show that turbulent reconnection dramatically affects the key astrophysical processes, e.g., star formation, turbulent dynamo, acceleration of cosmic rays. We provide criticism of the concept of "reconnection-mediated turbulence" and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity, and why the latter has fatal conceptual flaws.

preprint2020arXiv

Kelvin-Helmholtz versus Tearing Instability: What Drives Turbulence in Stochastic Reconnection?

Over the last few years it became clear that turbulent magnetic reconnection and magnetized turbulence are inseparable. It was not only shown that reconnection is responsible for violating the frozen-in condition in turbulence, but also that stochastic reconnection in 3D generates turbulence by itself. The actual mechanism responsible for this driving is still unknown. Processes such tearing mode or Kelvin-Helmholtz, among other plasma instabilities, could generate turbulence from irregular current sheets. We address the nature of driving mechanism for this process and consider a relative role of tearing and Kelvin-Helmholtz instabilities for the process of turbulence generation. In particular, we analyze the conditions for development of these two instabilities within three-dimensional reconnection regions. We show that both instabilities can excite turbulence fluctuations in reconnection regions. However, tearing mode has relatively slow growth rate, and at later times it becomes partially suppressed by transverse to the current sheet component of magnetic field, generated during the growth of turbulent fluctuations. On the contrary, the Kelvin-Helmholtz instability establishes quickly in the outflow region, and at later times it dominates the turbulence generation comparing to the contribution from tearing mode. Our results demonstrate that the tearing instability is subdominant compared to the the Kelvin-Helmholtz instability in terms of generation of turbulence in the 3D reconnection layers and therefore the self-driven reconnection is turbulent reconnection with tearing instability being important only at the initial stage of the reconnection.

preprint2020arXiv

Magnetic Field Morphology in Interstellar Clouds with the Velocity Gradient Technique

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.

preprint2020arXiv

Nonlinear turbulent dynamo during gravitational collapse

Via amplification by turbulent dynamo, magnetic fields can be potentially important for the formation of the first stars. To examine the dynamo behavior during the gravitational collapse of primordial gas, we extend the theory of nonlinear turbulent dynamo to include the effect of gravitational compression. The relative importance between dynamo and compression varies during contraction, with the transition from dynamo- to compression-dominated amplification of magnetic fields with the increase of density. In the nonlinear stage of magnetic field amplification with the scale-by-scale energy equipartition between turbulence and magnetic fields, reconnection diffusion of magnetic fields in ideal magnetohydrodynamic (MHD) turbulence becomes important. It causes the violation of flux-freezing condition and accounts for (a) the small growth rate of nonlinear dynamo, (b) the weak dependence of magnetic energy on density during contraction, (c) the saturated magnetic energy, and (d) the large correlation length of magnetic fields. The resulting magnetic field structure and the scaling of magnetic field strength with density are radically different from the expectations of flux-freezing.

preprint2020arXiv

Obtaining magnetic field strength using differential measure approach and velocity channel maps

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.

preprint2020arXiv

Studying the local magnetic field and anisotropy of magnetic turbulence by synchrotron polarization derivative

Due to the inevitable accumulation of the observed information in the direction of the line of sight, it is difficult to measure the local magnetic field of MHD turbulence. However, the correct understanding of the local magnetic field is a prerequisite for reconstructing the Galactic 3D magnetic field. We study how to reveal the local magnetic field direction and the eddy anisotropy on the basis of the statistics of synchrotron polarization derivative with respect to the squared wavelength $dP/dλ^2$. In the low frequency and strong Faraday rotation regime, we implement numerical simulations in the combination of multiple statistic techniques, such as structure function, quadrupole ratio modulus, spectral correlation function, correlation function anisotropy and spatial gradient techniques. We find that (1) statistic analysis of $dP/dλ^2$ indeed reveals the anisotropy of underlying MHD turbulence, the degree of which increases with the increase of the radiation frequency; (2) the synergy of both correlation function anisotropy and gradient calculation of $dP/dλ^2$ enables the measurement of the local magnetic field direction.

preprint2020arXiv

Trapping of cosmic rays in MHD turbulence

Astrophysical plasmas are turbulent and magnetized. The interaction between cosmic rays (CRs) and magnetohydrodynamic (MHD) turbulence is a fundamental astrophysical process. Based on the current understanding of MHD turbulence, we revisit the trapping of CRs by magnetic mirrors in the context of MHD turbulence. In compressible MHD turbulence, isotropic fast modes dominate both trapping and gyroresonant scattering of CRs. The presence of trapping significantly suppresses the pitch-angle scattering and the spatial diffusion of CRs along the magnetic field. The resulting parallel diffusion coefficient has a weaker dependence on CR energy at higher energies. In incompressible MHD turbulence, the trapping by pseudo-Alfvén modes dominates over the gyroresonant scattering by anisotropic Alfvén and pseudo-Alfvén modes at all pitch angles and prevents CRs from diffusion.

preprint2020arXiv

Turbulence in a self-gravitating molecular cloud core

Externally driven interstellar turbulence plays an important role in shaping the density structure in molecular clouds. Here we study the dynamical role of internally driven turbulence in a self-gravitating molecular cloud core. Depending on the initial conditions and evolutionary stages, we find that a self-gravitating core in the presence of gravity-driven turbulence can undergo constant, decelerated, and accelerated infall, and thus has various radial velocity profiles. In the gravity-dominated central region, a higher level of turbulence results in a lower infall velocity, a higher density, and a lower mass accretion rate. As an important implication of this study, efficient reconnection diffusion of magnetic fields against the gravitational drag naturally occurs due to the gravity-driven turbulence, without invoking externally driven turbulence.

preprint2020arXiv

Velocity Gradient in the Presence of Self-Gravity: Identifying Gravity-induced Inflow and Determining Collapsing Stage

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.

preprint2019arXiv

Intensity Gradients Technique: Synergy with Velocity Gradients and Polarization Studies

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.

preprint2019arXiv

Predictions of Cosmic Microwave Background Foregrounds Dust Polarization Using Velocity Gradients

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.

preprint2019arXiv

Tracing magnetic fields by the synergies of synchrotron emission gradients

This paper studies how to employ synchrotron emission gradient techniques to reveal the properties of the magnetic field within the interstellar media. Based on data cubes of three-dimensional numerical simulations of magnetohydrodynamic turbulence, we explore spatial gradients of synchrotron emission diagnostics to trace the direction of the magnetic field. According to our simulations, multifarious diagnostics for synchrotron emission can effectively determine the potential direction of projected magnetic fields. Applying the synergies of synchrotron diagnostic gradients to the archive data from the Canadian Galactic Plane Survey, we find that multifarious diagnostic techniques make consistent predictions for the Galactic magnetic field directions. With the high-resolution data presently available from Low Frequency Array for radio astronomy and those in the future from the Square Kilometer Array, the synergies of synchrotron emission gradients are supposed to perform better in tracing the actual direction of interstellar magnetic fields, especially in the low-frequency Faraday rotation regime where traditional synchrotron polarization measure fails.

preprint2016arXiv

Extending Velocity Channel Analysis for Studying Turbulence Anisotropies

We extend the velocity channel analysis (VCA), introduced by Lazarian & Pogosyan, of the intensity fluctuations in the velocity slices of position-position-velocity (PPV) spectroscopic data from Doppler broadened lines to study statistical anisotropy of the underlying velocity and density that arises in a turbulent medium from the presence of magnetic field. In particular, we study analytically how the anisotropy of the intensity correlation in the channel maps changes with the thickness of velocity channels. In agreement with the earlier VCA studies we find that the anisotropy in the thick channels reflects the anisotropy of the density field, while the relative contribution of density and velocity fluctuations to the thin velocity channels depends on the density spectral slope. We show that the anisotropies arising from Alfven, slow and fast magnetohydrodynamical modes are different, in particular, the anisotropy in PPV created by fast modes is opposite to that created by Alfven and slow modes, and this can be used to separate their contributions. We successfully compare our results with the recent numerical study of the PPV anisotropies measured with synthetic observations. We also extend our study to the medium with self-absorption as well as to the case of absorption lines. In addition, we demonstrate how the studies of anisotropy can be performed using interferometers.

preprint2016arXiv

Studying Magnetohydrodynamic Turbulence with Synchrotron Polarization Dispersion

We test a new technique of studying magnetohydrodynamic (MHD) turbulence suggested by Lazarian \& Pogosyan, using synthetic synchrotron polarization observations. This paper focuses on a one-point statistics, which is termed the polarization frequency analysis, that is characterized by the variance of polarized emission as a function of the square of wavelengths along a single line of sight. We adopt a ratio $η$ of the standard deviation of the line-of-sight turbulent magnetic field to the line-of-sight mean magnetic field to depict the level of turbulence. When this ratio is either large ($η\gg1$), which characterizes a turbulent field dominated region, or small ($η\lesssim0.2$), which characterizes a mean field dominated region, we obtain the polarization variance $\left<P^2\right>\proptoλ^{-2}$ and $\left<P^2\right>\proptoλ^{-2-2m}$, respectively. At small $η$, i.e., the mean field dominated region, we successfully recover the turbulent spectral index by the polarization variance. We find that our simulations agree well with the theoretical prediction of Lazarian \& Pogosyan. With existing and upcoming data cubes from the Low-Frequency Array for Radio astronomy (LOFAR) and Square Kilometer Array (SKA), this new technique can be applied to study the magnetic turbulence in the Milky Way and other galaxies.

preprint2015arXiv

Observational Diagnostics of Self-Gravitating MHD Turbulence in Giant Molecular Clouds

We study the observable signatures of self-gravitating MHD turbulence by applying the probability density functions (PDFs) and the spatial density power spectrum to synthetic column density maps. We find that there exists three characterizable stages of the evolution of the collapsing cloud which we term "early," "intermediate," and "advanced." At early times, i.e. $t<0.15t_{ff}$, the column density has a power spectral slope similar to nongravitating supersonic turbulence and a lognormal distribution. At an intermediate stage, i.e. $0.15t_{ff}< t \leq 0.35t_{ff}$, there exists signatures of the prestellar cores in the shallower PDF and power spectrum power law slopes. The column density PDF power law tails at these times have line of sight averaged slopes ranging from -2.5 to -1.5 with shallower values belonging to simulations with lower magnetic field strength. The density power spectrum slope becomes shallow and can be characterized by $P(k)=A_1k^{β_2}e^{-k/k_c}$, where $A_1$ describes the amplitude, $k^{β_2}$ describes the classical power law behavior and the scale $k_c$ characterizes the turn over from turbulence dominated to self-gravity dominated. At advanced stages of collapse, i.e. $\approx t>0.35t_{ff}$, the power spectral slope is positive valued, and a dramatic increase is observed in the PDF moments and the Tsallis incremental PDF parameters, which gives rise to deviations between PDF-sonic Mach number relations. Finally, we show that the imprint of gravity on the density power spectrum can be replicated in non-gravitating turbulence by introducing a delta-function with amplitude equivalent to the maximum valued point in a given self-gravitating map. We find that the turbulence power spectrum restored through spatial filtering of the high density material.

preprint2015arXiv

Principal Component Analysis studies of turbulence in optically thick gas

In this work we investigate the Principal Component Analysis (PCA) sensitivity to the velocity power spectrum in high opacity regimes of the interstellar medium (ISM). For our analysis we use synthetic Position-Position-Velocity (PPV) cubes of fractional Brownian motion (fBm) and magnetohydrodynamics (MHD) simulations, post processed to include radiative transfer effects from CO. We find that PCA analysis is very different from the tools based on the traditional power spectrum of PPV data cubes. Our major finding is that PCA is also sensitive to the phase information of PPV cubes and this allows PCA to detect the changes of the underlying velocity and density spectra at high opacities, where the spectral analysis of the maps provides the universal -3 spectrum in accordance with the predictions of Lazarian \& Pogosyan (2004) theory. This makes PCA potentially a valuable tool for studies of turbulence at high opacities provided that the proper gauging of the PCA index is made. The later, however, we found to be not easy, as the PCA results change in an irregular way for data with high sonic Mach numbers. This is in contrast to synthetic Brownian noise data used for velocity and density fields that show monotonic PCA behavior. We attribute this difference to the PCA's sensitivity to Fourier phase information.

preprint2015arXiv

Spectrum and Anisotropy of Turbulence from Multi-Frequency Measurement of Synchrotron Polarization

We consider turbulent synchrotron emitting media that also exhibits Faraday rotation and provide a statistical description of synchrotron polarization fluctuations. In particular, we consider these fluctuations as a function of the spatial separation of the direction of measurements and as a function of wavelength for the same line-of-sight. On the basis of our general analytical approach, we introduce several measures that can be used to obtain the spectral slopes and correlation scales of both the underlying magnetic turbulence responsible for emission and the spectrum of the Faraday rotation fluctuations. We show the synergetic nature of these measures and discuss how the study can be performed using sparsely sampled interferometric data. We also discuss how additional characteristics of turbulence can be obtained, including the turbulence anisotropy, the three dimensional direction of the mean magnetic field. We consider both cases when the synchrotron emission and Faraday rotation regions coincide and when they are spatially separated. Appealing to our earlier study in Lazarian and Pogosyan (2012) we explain that our new results are applicable to a wide range of spectral indexes of relativistic electrons responsible for synchrotron emission. We expect wide application of our techniques both with existing synchrotron data sets as well as with big forthcoming data sets from LOFAR and SKA.

preprint2015arXiv

The Phase Coherence of Interstellar Density Fluctuations

Studies of MHD turbulence often investigate the Fourier power spectrum to provide information on the nature of the turbulence cascade. However, the Fourier power spectrum only contains the Fourier amplitudes and rejects all information regarding the Fourier phases. Here we investigate the utility of two statistical diagnostics for recovering information on Fourier phases in ISM density data: the averaged amplitudes of the bispectrum and the phase coherence index (PCI), a new phase technique for the ISM. We create 3D density and 2D column density maps using a set of simulations of isothermal ideal MHD turbulence with a wide range of sonic and Alfvénic Mach numbers. We find that the bispectrum averaged along different angles with respect to either the $k_1$ or $k_2$ axis is primarily sensitive to the sonic Mach number while averaging the bispectral amplitudes over different annuli is sensitive to both the sonic and Alfvénic Mach numbers. The PCI of density suggests that the most correlated phases occur in supersonic sub-Alfvénic turbulence and also near the numerical dissipation regime. This suggests that non-linear interactions with correlated phases are strongest in shock dominated regions, in agreement with findings from the solar wind. Additionally, our results are particularly encouraging as they suggests the phase information contained in the bispectrum and PCI can be used to find parameters of turbulence in column density maps.

preprint2015arXiv

The Turbulence Velocity Power Spectrum of Neutral Hydrogen in the Small Magellanic Cloud

We present the results of the Velocity Coordinate Spectrum (VCS) technique to calculate the velocity power spectrum of turbulence in the Small Magellanic Cloud (SMC) in 21cm emission. We have obtained a velocity spectral index of -3.85 and an injection scale of 2.3 kpc. The spectral index is steeper than the Kolmogorov index which is expected for shock-dominated turbulence which is in agreement with past works on the SMC gas dynamics. The injection scale of 2.3 kpc suggests that tidal interactions with the Large Magellanic Cloud are the dominate driver of turbulence in this dwarf galaxy. This implies turbulence maybe driven by multiple mechanisms in galaxies in addition to supernova injection and that galaxy-galaxy interactions may play an important role.

preprint2015arXiv

Turbulent Reconnection and Its Implications

Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes in magnetized plasmas. In most astrophysical environments the Reynolds numbers are large and therefore the transition to turbulence is inevitable. This turbulence must be taken into account for any theory of magnetic reconnection, since the initially laminar configurations can transit to the turbulence state, what is demonstrated by 3D high resolution numerical simulations. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (1999) reconnection model and present numerical evidence supporting the model and demonstrate that it is closely connected to the concept of Richardson diffusion and compatible with the Lagrangian dynamics of magnetized fluids. We point out that the Generalized Ohm's Law, that accounts for turbulent motion, predicts the subdominance of the microphysical plasma effects for a realistically turbulent media. We show that on of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a corner stone of most theories dealing with magnetized plasmas and therefore its change induces fundamental shifts in accepted paradigms like turbulent reconnection entailing the diffusion process that is essential for understanding star formation. We argue, that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares predicted by turbulent reconnection and relate them to solar flares and gamma ray bursts. We analyze solar observations, measurements in the solar wind or heliospheric current sheet, and show their correspondence with turbulent reconnection predictions. Finally, we discuss 1st Order Fermi acceleration as a natural consequence of the turbulent reconnection.

preprint2014arXiv

Cosmic Ray Acceleration in Magnetic Reconnection Sites

Cosmic Ray (CR) acceleration still challenges the researchers. Fast particles may be accelerated in astrophysical environments by a variety of processes. Acceleration in magnetic reconnection sites in particular, has lately attracted the attention of researchers not only for its potential importance in the solar system context, but also in other astrophysical environments, like compact stellar sources, AGNs and GRBs, and even in diffusive media like the ISM and the IGM, especially when the environment is magnetically dominated. In this talk we review this process and also present three-dimensional collisional MHD simulations with the injection of thousands of test particles showing from the evolution of their energy spectrum that they can be efficiently accelerated by reconnection through a first-order Fermi process within large scale magnetic current sheets (especially when local turbulence is present which makes reconnection fast and the acceleration layer thicker).

preprint2014arXiv

Filaments in the southern giant lobe of Centaurus A: constraints on nature and origin from modelling and GMRT observations

We present results from imaging of the radio filaments in the southern giant lobe of Centaurus A using data from GMRT observations at 325 and 235 MHz, and outcomes from filament modelling. The observations reveal a rich filamentary structure, largely matching the morphology at 1.4 GHz. We find no clear connection of the filaments to the jet. We seek to constrain the nature and origin of the vertex and vortex filaments associated with the lobe and their role in high-energy particle acceleration. We deduce that these filaments are at most mildly overpressured with respect to the global lobe plasma showing no evidence of large-scale efficient Fermi I-type particle acceleration, and persist for ~ 2-3 Myr. We demonstrate that the dwarf galaxy KK 196 (AM 1318-444) cannot account for the features, and that surface plasma instabilities, the internal sausage mode and radiative instabilities are highly unlikely. An internal tearing instability and the kink mode are allowed within the observational and growth time constraints and could develop in parallel on different physical scales. We interpret the origin of the vertex and vortex filaments in terms of weak shocks from transonic MHD turbulence or from a moderately recent jet activity of the parent AGN, or an interplay of both.

preprint2014arXiv

MHD Turbulence, Turbulent Dynamo and Applications

MHD Turbulence is common in many space physics and astrophysics environments. We first discuss the properties of incompressible MHD turbulence. A well-conductive fluid amplifies initial magnetic fields in a process called small-scale dynamo. Below equipartition scale for kinetic and magnetic energies the spectrum is steep (Kolmogorov -5/3) and is represented by critically balanced strong MHD turbulence. In this paper we report the basic reasoning behind universal nonlinear small-scale dynamo and the inertial range of MHD turbulence. We measured the efficiency of the small-scale dynamo $C_E=0.05$, Kolmogorov constant $C_K=4.2$ and anisotropy constant $C_A=0.63$ for MHD turbulence in high-resolution direct numerical simulations. We also discuss so-called imbalanced or cross-helical MHD turbulence which is relevant for in many objects, most prominently in the solar wind. We show that properties of incompressible MHD turbulence are similar to the properties of Alfvénic part of MHD cascade in compressible turbulence. The other parts of the cascade evolve according to their own dynamics. The slow modes are being cascaded by Alfvénic modes, while fast modes create an independent cascade. We show that different ways of decomposing compressible MHD turbulence into Alfvén, slow and fast modes provide consistent results and are useful in understanding not only turbulent cascade, but its interaction with fast particles.

preprint2014arXiv

Opacity Broadening of $^{13}$CO Linewidths and its Effect on the Variance-Sonic Mach Number Relation

We study how the estimation of the sonic Mach number ($M_s$) from $^{13}$CO linewidths relates to the actual 3D sonic Mach number. For this purpose we analyze MHD simulations which include post-processing to take radiative transfer effects into account. As expected, we find very good agreement between the linewidth estimated sonic Mach number and the actual sonic Mach number of the simulations for optically thin tracers. However, we find that opacity broadening causes $M_s$ to be overestimated by a factor of ~ 1.16 - 1.3 when calculated from optically thick $^{13}$CO lines. We also find that there is a dependency on the magnetic field: super-Alfvénic turbulence shows increased line broadening as compared with sub-Alfvénic turbulence for all values of optical depth for supersonic turbulence. Our results have implications for the observationally derived sonic Mach number--density standard deviation ($σ_{ρ/<ρ>}$) relationship, $σ^2_{ρ/<ρ>}=b^2M_s^2$, and the related column density standard deviation ($σ_{N/<N>}$) sonic Mach number relationship. In particular, we find that the parameter b, as an indicator of solenoidal vs. compressive driving, will be underestimated as a result of opacity broadening. We compare the $σ_{N/<N>}-M_s$ relation derived from synthetic dust extinction maps and $^{13}$CO linewidths with recent observational studies and find that solenoidally driven MHD turbulence simulations have values of $σ_{N/<N>}$ which are lower than real molecular clouds. This may be due to the influence of self-gravity which should be included in simulations of molecular cloud dynamics.

preprint2014arXiv

Studies of Anomalous Microwave Emission (AME) with the SKA

In this chapter, we will outline the scientific motivation for studying Anomalous Microwave Emission (AME) with the SKA. AME is thought to be due to electric dipole radiation from small spinning dust grains, although thermal fluctuations of magnetic dust grains may also contribute. Studies of this mysterious component would shed light on the emission mechanism, which then opens up a new window onto the interstellar medium (ISM). AME is emitted mostly in the frequency range $\sim 10$--100\,GHz, and thus the SKA has the potential of measuring the low frequency side of the AME spectrum, particularly in band 5. Science targets include dense molecular clouds in the Milky Way, as well as extragalactic sources. We also discuss the possibility of detecting rotational line emission from Poly-cyclic Aromatic Hydrocarbons (PAHs), which could be the main carriers of AME. Detecting PAH lines of a given spacing would allow for a definitive identification of specific PAH species.

preprint2013arXiv

Characterizing the Turbulent Properties of the Starless Molecular Cloud MBM16

We investigate turbulent properties of the non-star-forming, translucent molecular cloud, MBM16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (HI) observations obtained by the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a single power-law function, with a slope ranging from -3.3 to -3.7 and being consistent over the velocity range of MBM16 for a fixed velocity channel thickness. However, the slope varies significantly with the velocity slice thickness, suggesting that both velocity and density contribute to HI intensity fluctuations. By using this variation we estimate the slope of 3D density fluctuations in MBM16 to be -3.7\pm0.2. This is significantly steeper than what has been found for HI in the Milky Way plane, the Small Magellanic Cloud, or the Magellanic Bridge, suggesting that interstellar turbulence in MBM16 is driven on scales >17 pc and that the lack of stellar feedback could be responsible for the steep power spectrum.

preprint2013arXiv

Mass entrainment and turbulence-driven acceleration of ultra-high energy cosmic rays in Centaurus A

Observations of the FR I radio galaxy Centaurus A in radio, X-ray and gamma-ray bands provide evidence for lepton acceleration up to several TeV and clues about hadron acceleration to tens of EeV. Synthesising the available observational constraints on the physical conditions and particle content in the jets, inner lobes and giant lobes of Centaurus A, we aim to evaluate its feasibility as an ultra-high-energy cosmic-ray source. We apply several methods of determining jet power and affirm the consistency of various power estimates of ~ 1 x 10^43 erg s^-1. Employing scaling relations based on previous results for 3C 31, we estimate particle number densities in the jets, encompassing available radio through X-ray observations. Our model is compatible with the jets ingesting ~ 3 x 10^21 g s^-1 of matter via external entrainment from hot gas and ~ 7 x 10^22 g s^-1 via internal entrainment from jet-contained stars. This leads to an imbalance between the internal lobe pressure available from radiating particles and magnetic field, and our derived external pressure. Based on knowledge of the external environments of other FR I sources, we estimate the thermal pressure in the giant lobes as 1.5 x 10^-12 dyn cm^-2, from which we deduce a lower limit to the temperature of ~ 1.6 x 10^8 K. Using dynamical and buoyancy arguments, we infer ~ 440-645 Myr and ~ 560 Myr as the sound-crossing and buoyancy ages of the giant lobes respectively, inconsistent with their spectral ages. We re-investigate the feasibility of particle acceleration via stochastic processes in the lobes, placing new constraints on the energetics and on turbulent input to the lobes. The same 'very hot' temperatures that allow self-consistency between the entrainment calculations and the missing pressure also allow stochastic UHECR acceleration models to work.

preprint2013arXiv

Observational Diagnostics for Two-Fluid Turbulence in Molecular Clouds As Suggested by Simulations

We present high resolution simulations of two-fluid (ion-neutral) MHD turbulence with resolutions as large as 512^3. The simulations are supersonic and mildly sub-Alfvenic, in keeping with the conditions present in molecular clouds. Such turbulence is thought to influence star formation processes in molecular clouds because typical cores form on length scales that are comparable to the dissipation scales of this turbulence in the ions. The simulations are motivated by the fact that recent studies of isophotologue lines in molecular clouds have found significant differences in the linewidth-size relationship for neutral and ion species. The goals of this paper are to explain those observations using simulations and analytic theory, present a new set of density-based diagnostics by drawing on similar diagnostics that have been obtained by studying single-fluid turbulence, and show that our two-fluid simulations play a vital role in reconciling alternative models of star formation. The velocity-dependent diagnostics display a complementarity with the density-dependent diagnostics. We find that the linewidth-size relationships show a prominent difference between ions and neutrals when the line of sight is orthogonal to the mean field. This is because the MHD waves in the ions differ from hydrodynamic waves in the neutrals. We also find that the density probability distribution functions (PDFs) show prominent differences between the ions and neutrals when the line of sight is parallel to the mean field. This is because the velocity fluctuations in the ions tend to produce column density fluctuations along field lines. When the magnetic field makes an angle to the line of sight, both observable differences should be visible. These diagnostics should be easy for observers to test. This analysis assumes optically thin lines and a mean magnetic field that is uniform in direction within a cloud.

preprint2012arXiv

Particle Acceleration in Turbulence and Weakly Stochastic Reconnection

Fast particles are accelerated in astrophysical environments by a variety of processes. Acceleration in reconnection sites has attracted the attention of researchers recently. In this letter we analyze the energy distribution evolution of test particles injected in three dimensional (3D) magnetohydrodynamic (MHD) simulations of different magnetic reconnection configurations. When considering a single Sweet-Parker topology, the particles accelerate predominantly through a first-order Fermi process, as predicted in previous work (de Gouveia Dal Pino & Lazarian, 2005) and demonstrated numerically in Kowal, de Gouveia Dal Pino & Lazarian (2011). When turbulence is included within the current sheet, the acceleration rate, which depends on the reconnection rate, is highly enhanced. This is because reconnection in the presence of turbulence becomes fast and independent of resistivity (Lazarian & Vishniac, 1999; Kowal et al., 2009) and allows the formation of a thick volume filled with multiple simultaneously reconnecting magnetic fluxes. Charged particles trapped within this volume suffer several head-on scatterings with the contracting magnetic fluctuations, which significantly increase the acceleration rate and results in a first-order Fermi process. For comparison, we also tested acceleration in MHD turbulence, where particles suffer collisions with approaching and receding magnetic irregularities, resulting in a reduced acceleration rate. We argue that the dominant acceleration mechanism approaches a second order Fermi process in this case.

preprint2012arXiv

Reconnection Studies Under Different Types of Turbulence Driving

We study a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian and Vishniac (1999) using three-dimensional direct numerical simulations. The model has been already successfully tested in Kowal et al. (2009) confirming the dependencies of the reconnection speed $V_{rec}$ on the turbulence injection power $P_{inj}$ and the injection scale $l_{inj}$ expressed by a constraint $V_{rec} \sim P_{inj}^{1/2} l_{inj}^{3/4}$ and no observed dependency on Ohmic resistivity. In Kowal et al. (2009), in order to drive turbulence, we injected velocity fluctuations in Fourier space with frequencies concentrated around $k_{inj}=1/l_{inj}$, as described in Alvelius (1999). In this paper we extend our previous studies by comparing fast magnetic reconnection under different mechanisms of turbulence injection by introducing a new way of turbulence driving. The new method injects velocity or magnetic eddies with a specified amplitude and scale in random locations directly in real space. We provide exact relations between the eddy parameters and turbulent power and injection scale. We performed simulations with new forcing in order to study turbulent power and injection scale dependencies. The results show no discrepancy between models with two different methods of turbulence driving exposing the same scalings in both cases. This is in agreement with the Lazarian and Vishniac (1999) predictions. In addition, we performed a series of models with varying viscosity $ν$. Although Lazarian and Vishniac (1999) do not provide any prediction for this dependence, we report a weak relation between the reconnection speed with viscosity, $V_{rec}\simν^{-1/4}$.

preprint2011arXiv

Fast Magnetic Reconnection and Spontaneous Stochasticity

Magnetic field-lines in astrophysical plasmas are expected to be frozen-in at scales larger than the ion gyroradius. The rapid reconnection of magnetic flux structures with dimensions vastly larger than the gyroradius requires a breakdown in the standard Alfvén flux-freezing law. We attribute this breakdown to ubiquitous MHD plasma turbulence with power-law scaling ranges of velocity and magnetic energy spectra. Lagrangian particle trajectories in such environments become "spontaneously stochastic", so that infinitely-many magnetic field-lines are advected to each point and must be averaged to obtain the resultant magnetic field. The relative distance between initial magnetic field lines which arrive to the same final point depends upon the properties of two-particle turbulent dispersion. We develop predictions based on the phenomenological Goldreich & Sridhar theory of strong MHD turbulence and on weak MHD turbulence theory. We recover the predictions of the Lazarian & Vishniac theory for the reconnection rate of large-scale magnetic structures. Lazarian & Vishniac also invoked "spontaneous stochasticity", but of the field-lines rather than of the Lagrangian trajectories. More recent theories of fast magnetic reconnection appeal to microscopic plasma processes that lead to additional terms in the generalized Ohm's law, such as the collisionless Hall term. We estimate quantitatively the effect of such processes on the inertial-range turbulence dynamics and find them to be negligible in most astrophysical environments. For example, the predictions of the Lazarian-Vishniac theory are unchanged in Hall MHD turbulence with an extended inertial range, whenever the ion skin depth $δ_i$ is much smaller than the turbulent integral length or injection-scale $L_i.$

preprint2011arXiv

Low-Mach-number turbulence in interstellar gas revealed by radio polarization gradients

The interstellar medium of the Milky Way is multi-phase, magnetized and turbulent. Turbulence in the interstellar medium produces a global cascade of random gas motions, spanning scales ranging from 100 parsecs to 1000 kilometres. Fundamental parameters of interstellar turbulence such as the sonic Mach number (the speed of sound) have been difficult to determine because observations have lacked the sensitivity and resolution to directly image the small-scale structure associated with turbulent motion. Observations of linear polarization and Faraday rotation in radio emission from the Milky Way have identified unusual polarized structures that often have no counterparts in the total radiation intensity or at other wavelengths, and whose physical significance has been unclear. Here we report that the gradient of the Stokes vector (Q,U), where Q and U are parameters describing the polarization state of radiation, provides an image of magnetized turbulence in diffuse ionized gas, manifested as a complex filamentary web of discontinuities in gas density and magnetic field. Through comparison with simulations, we demonstrate that turbulence in the warm ionized medium has a relatively low sonic Mach number, M_s <~ 2. The development of statistical tools for the analysis of polarization gradients will allow accurate determinations of the Mach number, Reynolds number and magnetic field strength in interstellar turbulence over a wide range of conditions.

preprint2011arXiv

Magnetic Flux Transport by turbulent reconnection in astrophysical flows

The role of MHD turbulence in astrophysical environments is still highly debated. An important question that permeates this debate is the transport of magnetic flux. This is particularly important, for instance, in the context of star formation. When clouds collapse gravitationally to form stars, there must be some magnetic flux transport. otherwise the new born stars would have magnetic fields several orders of magnitude larger than the observed ones. Also, the magnetic flux that is dragged in the late stages of the formation of a star can remove all the rotational support from the accretion disk that grows around the protostar. The efficiency of the mechanism which is often invoked to allow the transport of magnetic fields in the different stages of star formation, namely, the ambipolar diffusion, has been lately put in check. We here discuss an alternative mechanism for magnetic flux transport which is based on turbulent fast magnetic reconnection. We review recent results obtained from 3D MHD numerical simulations that indicate that this mechanism is very efficient for decoupling and transport magnetic flux from the inner denser regions to the outskirts of collapsing clouds in the different stages of star formation. We also discuss this mechanism in the context of dynamo processes and speculate that it can play a role both in the solar dynamo and in accretion disk dynamo processes.

preprint2011arXiv

Particle Acceleration by Magnetic Reconnection in AGNs and in the IGM

There is no single mechanism by which fast particles are accelerated in astrophysical environments, and it is now recognized that the data require a rich variety of different mechanisms operating under different conditions. The mechanisms discussed in the literature include varying magnetic fields in compact sources, stochastic processes in turbulent environments, and acceleration behind shocks. An alternative, much less explored mechanism so far, involves particle acceleration within magnetic reconnection sites. In this work, we explore this mechanism in the AGN framework and show that particles are efficiently accelerated through a first-order Fermi process and have an exponential growth of energy. We also address briefly the propagation of cosmic rays (CRs) in the intergalactic medium (IGM). Since the latter is a collisionless environment, kinetic effects must be considered which will affect the turbulent magnetic field distribution and therefore, the CR propagation.

preprint2011arXiv

Statistical Description of Synchrotron Intensity Fluctuations: Studies of Astrophysical Magnetic Turbulence

We provide a theoretical description of synchrotron fluctuations arising from magnetic turbulence. We derive an expression that relates the correlation of synchrotron fluctuations for an arbitrary index of relativistic electrons with the correlations of squared fluctuations of magnetic field component perpendicular to the line of sight. The latter correlations we study assuming that the turbulence is axisymmetric. We obtain general relations valid for an arbitrary model of magnetic turbulence and analyse the relations for particular example of magnetic turbulence that is supported by numerical simulations. We predict that the synchrotron intensity fluctuations are anisotropic with larger correlation present along the direction of magnetic field. This anisotropy is dominated by the quadrupole component with the ratio between quadrupole and monopole parts being sensitive to the compressibility of underlying turbulence. Our work opens avenues for quantitative studies of magnetic turbulence in our galaxy and beyond using synchrotron emission. It also outlines the directions of how synchrotron foreground emission can be separated from cosmological signal, i.e. from CMB or highly redshifted HI emission. For the sake of completeness we also provide the expressions for the synchrotron polarization (Stocks parameters and their combinations) for the model of axisymmetric magnetic turbulence.

preprint2011arXiv

Statistical Tools of Interstellar Turbulence: Bridging the Gap Between Numerics and Observations

MHD Turbulence is a critical component of the current paradigms of star formation, particle transport, magnetic reconnection and evolution of the ISM. Progress on this difficult subject is made via numerical simulations and observational studies. However, due to limitations of resolution, scale discrepancies, and complexity of the observations, the best approach for connecting numerics to observations is not always obvious. Here we advocate for a approach that invokes statistical techniques to understand the underlying physics of turbulent astrophysical systems. The wealth of numerical and observational data calls for new statistical tools to be developed in order to study turbulence in the interstellar medium. We briefly review some of the recently developed statistics that focus on characterizing gas compressibility and magnetization and their uses to interstellar studies.

preprint2011arXiv

Statistical Tools of ISM turbulence

MHD Turbulence is a critical component of the current paradigms of star formation, particle transport, magnetic reconnection and evolution of the ISM, to name just a few. Progress on this difficult subject is made via numerical simulations and observational studies, however in order to connect these two, statistical methods are required. This calls for new statistical tools to be developed in order to study turbulence in the interstellar medium. Here we briefly review some of the recently developed statistics that focus on characterizing gas compressibility and magnetization and their uses to interstellar studies.

preprint2011arXiv

Velocity Anisotropy as a Diagnostic of the Magnetization of the Interstellar Medium and Molecular clouds

We use a set of magnetohydrodynamics (MHD) simulations of fully-developed (driven) turbulence to study the anisotropy in the velocity field that is induced by the presence of the magnetic field. In our models we study turbulence characterized by sonic Mach numbers M_s from 0.7 to 7.5, and Alfven Mach numbers M_A from 0.4 to 7.7. These are used to produce synthetic observations (centroid maps) that are analyzed. To study the effect of large scale density fluctuations and of white noise we have modified the density fields and obtained new centroid maps, which are analyzed. We show that restricting the range of scales at which the anisotropy is measured makes the method robust against such fluctuations. We show that the anisotropy in the structure function of the maps reveals the direction of the magnetic field for M_A \lesssim 1.5, regardless of the sonic Mach number. We found that the degree of anisotropy can be used to determine the degree of magnetization (i.e. M_A) for M_A \lesssim 1.5. To do this, one needs an additional measure of the sonic Mach number and an estimate of the LOS magnetic field, both feasible by other techniques, offering a new opportunity to study the magnetization state of the interstellar medium.

preprint2010arXiv

Velocity Field of Compressible MHD Turbulence: Wavelet Decomposition and Mode Scalings

We study compressible MHD turbulence, which holds key to many astrophysical processes, including star formation and cosmic ray propagation. To account for the variations of the magnetic field in the strongly turbulent fluid we use wavelet decomposition of the turbulent velocity field into Alfven, slow and fast modes, which presents an extension of the Cho & Lazarian (2003) decomposition approach based on Fourier transforms. The wavelets allow to follow the variations of the local direction of magnetic field and therefore improve the quality of the decomposition compared to the Fourier transforms which are done in the mean field reference frame. For each resulting component we calculate spectra and two-point statistics such as longitudinal and transverse structure functions, as well as, higher order intermittency statistics. In addition, we perform the Helmholtz-Hodge decomposition of the velocity field into the incompressible and compressible parts and analyze these components. We find that the turbulence intermittency is different for different components and we show that the intermittency statistics depend on whether the phenomenon was studied in the global reference frame related to the mean magnetic field or it was studied in the frame defined by the local magnetic field. The dependencies of the measures we obtained are different for different components of velocity, for instance, we show that while the Alfven mode intermittency changes marginally with the Mach number the intermittency of the fast mode is substantially affected by the change.

preprint2009arXiv

Characterizing Magnetohydrodynamic Turbulence in the Small Magellanic Cloud

We investigate the nature and spatial variations of turbulence in the Small Magellanic Cloud (SMC) by applying several statistical methods on the neutral hydrogen (HI) column density image of the SMC and a database of isothermal numerical simulations. By using the 3rd and 4th statistical moments we derive the spatial distribution of the sonic Mach number (M_s) across the SMC. We find that about 90% of the HI in the SMC is subsonic or transonic. However, edges of the SMC `bar' have M_s=4 and may be tracing shearing or turbulent flows. Using numerical simulations we also investigate how the slope of the spatial power spectrum depends on both sonic and Alfven Mach numbers. This allows us to gauge the Alfven Mach number of the SMC and conclude that its gas pressure dominates over the magnetic pressure. The super-Alfvenic nature of the HI gas in the SMC is also highlighted by the bispectrum, a three-point correlation function which characterizes the level of non-Gaussianity in wave modes. We find that the bispectrum of the SMC HI column density displays similar large-scale correlations as numerical simulations, however it has localized enhancements of correlations. In addition, we find a break in correlations at a scale of 160 pc. This may be caused by numerous expanding shells of a similar size.

preprint2009arXiv

On the role of magnetic reconnection in jet/accretion disk systems

The most accepted model for jet production is based on the magneto-centrifugal acceleration out off an accretion disk that surrounds the central source (Blandford & Payne, 1982). This scenario, however, does not explain, e.g., the quasi-periodic ejection phenomena often observed in different astrophysical jet classes. de Gouveia Dal Pino & Lazarian (2005) (hereafter GDPL) have proposed that the large scale superluminal ejections observed in microquasars during radio flare events could be produced by violent magnetic reconnection (MR) episodes. Here, we extend this model to other accretion disk systems, namely: active galactic nuclei (AGNs) and young stellar objects (YSOs), and also discuss its role on jet heating and particle acceleration.

preprint2008arXiv

Grain alignment induced by radiative torques: effects of internal relaxation of energy and complex radiation fields

Earlier studies of grain alignment dealt mostly with interstellar grains that have strong internal relaxation of energy which aligns grain axis of maximum moment of inertia with respect to grain's angular momentum. In this paper, we study the alignment by radiative torques for large irregular grains, e.g., grains in accretion disks, for which internal relaxation is subdominant. We use both numerical calculations and the analytical model of a helical grain introduced by us earlier. We demonstrate that grains in such a regime exhibit more complex dynamics. In particular, if initially the grain axis of maximum moment of inertia makes a small angle with angular momentum, then radiative torques can align the grain axis of maximum moment of inertia with angular momentum, and both axis of maximum moment of inertia and angular momentum are aligned with the magnetic field when attractors with high angular momentum (high-J attractors) are available. For the alignment without high-J attractors, beside the earlier studied attractors with low angular momentum (low-J attractors), there appears new low-J attractors. The former and later cases correspond to the alignment with long axes perpendicular and parallel to the angular momentum, respectively. In addition, we study the alignment of grains in the presence of strong internal relaxation, but induced not by a radiation beam as in earlier studies, instead, induced by a complex radiation field, that can be decomposed into dipole and quadrupole components. We find that in this situation, the parameter space $q^{max}$, for the existence of high-$J$ attractors is more extended, which entails higher degrees of polarization expected. Our obtained results are useful for modeling polarization arising from aligned grains in molecular clouds and accretion disks.

preprint2008arXiv

Wave Decay in MHD Turbulence

We present a model for nonlinear decay of the weak wave in three-dimensional incompressible magnetohydrodynamic (MHD) turbulence. We show that the decay rate is different for parallel and perpendicular waves. We provide a general formula for arbitrarily directed waves and discuss particular limiting cases known in the literature. We test our predictions with direct numerical simulations of wave decay in three-dimensional MHD turbulence, and discuss the influence of turbulent damping on the development of linear instabilities in the interstellar medium and on other important astrophysical processes.

preprint2001arXiv

Constraints on the Acceleration of Ultra-High-Energy Cosmic Rays in Accretion-Induced Collapse Pulsars

(Abridged) We have recently proposed that the ultra-high energy cosmic rays (UHECRs) observed above the GZK limit could be mostly protons accelerated in reconnection sites just above the magnetosphere of newborn millisecond pulsars originated by accretion induced collapse (AIC-pulsars). Although the expected rate of AIC sources in our own Galaxy is very small, our estimates have shown that the observed total flux of UHECRs could be obtained from the integrated contribution from AIC-pulsars of the whole distribution of galaxies located within a distance which is unaffected by the GZK cutoff ($\sim 50 $ Mpc). We presently examine the potential acceleration mechanisms in the magnetic reconnection site and find that first-order Fermi acceleration cannot provide sufficient efficiency (due to synchrotron losses). This leaves the one-shot acceleration via an induced electric field within the reconnection region as the only viable process for UHECR acceleration. We formulate the constraints on both the magnetic field topology and strength in order to accelerate the particles and allow them to freely escape from the system. Under fast reconnection, we find that AIC-pulsars with surface B-fields $10^{12} $ G $< B_{\star} \lesssim $ $10^{15}$ G and spin periods 1 ms $\lesssim P_{\star} < $ 60 ms, are able to accelerate particles to energies $\geq 10^{20} $ eV, but the magnetic field just above the Alfvén surface must be predominantly toroidal for the particles to be allowed to escape from the acceleration zone without being deflected. Synchrotron losses bring potentially important constraints on the B-field geometry of $any$ UHECR accelerators involving compact sources with strong magnetic fields.

preprint2001arXiv

Dust Polarization From Starlight Data

We present a statistical analysis of the interstellar medium (ISM) polarization from the largest compilation available of starlight data, which comprises ~ 5500 stars. The measured correlation between the mean polarization degree and extinction indicates that ISM dust grains are not fully aligned with the uniform component of the large-scale Galactic magnetic field. Moreover, we estimate the ratio of the uniform to the random plane-of-the-sky components of the magnetic field to be B_u/B_r = 0.8. From the analysis of starlight polarization degree and position angle we find that the magnetic field broadly follows Galactic structures on large-scales. On the other hand, the angular power spectrum C_l of the polarization degree for Galactic plane data is found to be consistent with a power-law, C_l ~ l^{-1.5} (where l = 180 deg/θis the multipole order), for angular scales θ> 10 arcmin. We argue that this data set can be used to estimate diffuse polarized emission at microwave frequencies.

preprint2001arXiv

Simulations of MHD Turbulence in a Strongly Magnetized Medium

We analyze 3D numerical simulations of driven incompressible magnetohydrodynamic (MHD) turbulence in a periodic box threaded by a moderately strong external magnetic field. We sum over nonlinear interactions within Fourier wavebands and find that the time scale for the energy cascade is consistent with the Goldreich-Sridhar model of strong MHD turbulence. Using higher order longitudinal structure functions we show that the turbulent motions in the plane perpendicular to the local mean magnetic field are similar to ordinary hydrodynamic turbulence while motions parallel to the field are consistent with a scaling correction which arises from the eddy anisotropy. We present the structure tensor describing velocity statistics of Alfvenic and pseudo-Alfvenic turbulence. Finally, we confirm that an imbalance of energy moving up and down magnetic field lines leads to a slow decay of turbulent motions and speculate that this imbalance is common in the interstellar medium where injection of energy is intermittent both in time and space.

preprint1995arXiv

Physics and Chemistry of Purcell's Alignment

Paramagnetic alignment of suprathermally rotating grains is discussed in view of recent progress in understanding subtle processes taking place over grain surface. It is shown that in typical ISM conditions, grains with surfaces of amorphous H$_{2}$O ice, defected silicate or polymeric carbonaceous material are likely to exhibit enhanced alignment, while those of aromatic carbonaceous material or graphite are not. The critical grain sizes and temperature for the Purcell's alignment are obtained and preferential alignment of large grains is established.