Source author record

A. Lazarian

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

138works
18topics
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

138 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'' (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'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' 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's polarization measurements at 353 GHz.

preprint2016arXiv

A unified model of grain alignment: radiative alignment of interstellar grains with magnetic inclusions

The radiative torque (RAT) alignment of interstellar grains with ordinary paramagnetic susceptibilities has been supported by earlier studies. The alignment of such grains depends on the so-called RAT parameter $q^{\max}$ that is determined by the grain shape. In this paper, we elaborate our model of RAT alignment for grains with enhanced magnetic susceptibility due to iron inclusions, such that RAT alignment is magnetically enhanced for which we term MRAT mechanism. Such grains can get aligned with high angular momentum at the so-called high-J attractor points, achieving a high degree of alignment. Using our analytical model of RATs we derive the critical value of the magnetic relaxation parameter $δ_{m}$ to produce high-J attractor points as functions of $q^{\max}$ and the anisotropic radiation angle relative to the magnetic field $ψ$. We find that if about $10\%$ of total iron abundance present in silicate grains are forming iron clusters, it is sufficient to produce high-J attractor points for all reasonable values of $q^{\max}$. To calculate the degree of grain alignment, we carry out numerical simulations of MRAT alignment by including stochastic excitations from gas collisions and magnetic fluctuations. We show that large grains can achieve perfect alignment when the high-J attractor point is present, regardless of the values of $q^{\max}$. Our obtained results pave the way for physical modeling of polarized thermal dust emission as well as magnetic dipole emission. We also find that millimeter-sized grains in accretion disks may be aligned with the magnetic field if they are incorporated with iron nanoparticles.

preprint2016arXiv

Cosmic Ray Small Scale Anisotropies and Local Turbulent Magnetic Fields

Cosmic ray anisotropy has been observed in a wide energy range and at different angular scales by a variety of experiments over the past decade. However, no comprehensive or satisfactory explanation has been put forth to date. The arrival distribution of cosmic rays at Earth is the convolution of the distribution of their sources and of the effects of geometry and properties of the magnetic field through which particles propagate. It is generally believed that the anisotropy topology at the largest angular scale is adiabatically shaped by diffusion in the structured interstellar magnetic field. On the contrary, the medium- and small-scale angular structure could be an effect of non-diffusive propagation of cosmic rays in perturbed magnetic fields. In particular, a possible explanation of the observed small-scale anisotropy observed at TeV energy scale, may come from the effect of particle scattering in turbulent magnetized plasmas. We perform numerical integration of test particle trajectories in low-$β$ compressible magnetohydrodynamic turbulence to study how the cosmic rays arrival direction distribution is perturbed when they stream along the local turbulent magnetic field. We utilize Liouville's theorem for obtaining the anisotropy at Earth and provide the theoretical framework for the application of the theorem in the specific case of cosmic ray arrival distribution. In this work, we discuss the effects on the anisotropy arising from propagation in this inhomogeneous and turbulent interstellar magnetic field.

preprint2016arXiv

Damping of Alfven waves by Turbulence and its Consequences: from Cosmic-Rays Streaming to Launching Winds

This paper considers turbulent damping of Alfven waves in magnetized plasmas. We identify two cases of damping, one related to damping of cosmic rays streaming instability, the other related to damping of Alfven waves emitted by a macroscopic wave source, e.g. stellar atmosphere. The physical difference between the two cases is that in the former case the generated waves are emitted in respect to the local direction of magnetic field, in the latter in respect to the mean field. The scaling of damping is different in the two cases. We the regimes of turbulence ranging from subAlfvenic to superAlfvenic we obtain analytical expressions for the damping rates and define the ranges of applicability of these expressions. Describing the damping of the streaming instability, we find that for subAlfvenic turbulence the range of cosmic ray energies influenced by weak turbulence is unproportionally large compared to the range of scales that the weak turbulence is present. On the contrary, the range of cosmic ray energies affected by strong Alfvenic turbulence is rather limited. A number of astrophysical applications of the process ranging from launching of stellar and galactic winds to propagation of cosmic rays in galaxies and clusters of galaxies is considered. In particular, we discuss how to reconcile the process of turbulent damping with the observed isotropy of the Milky Way cosmic rays.

preprint2016arXiv

Grain Alignment: Role of Radiative Torques and Paramagnetic Relaxation

Polarization arising from aligned dust grains presents a unique opportunity to study magnetic fields in the diffuse interstellar medium and molecular clouds. Polarization from circumstellar regions, accretion disks and comet atmospheres can also be related to aligned dust.To reliably trace magnetic fields quantitative theory of grain alignment is required. Formulating the theory that would correspond to observations was one of the longstanding problems in astrophysics. Lately this problem has been successfully addressed, and in this review we summarize some of the most important theoretical advances in the theory of grain alignment by radiative torques (RATs) that act on realistic irregular dust grains. We discuss an analytical model of RATs and the ways to make RAT alignment more efficient, e.g. through paramagnetic relaxation when grains have inclusions with strong magnetic response. For very small grains for which RAT alignment is inefficient, we also discuss paramagnetic relaxation and a process termed resonance relaxation. We provide an extensive analysis of the observational tests of grain alignment theory.

preprint2016arXiv

Limits on the ions temperature anisotropy in turbulent intracluster medium

Turbulence in the weakly collisional intracluster medium of galaxies (ICM) is able to generate strong thermal velocity anisotropies in the ions (with respect to the local magnetic field direction), if the magnetic moment of the particles is conserved in the absence of Coulomb collisions. In this scenario, the anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very different statistical behaviour from the standard MHD one and is unable to amplify seed magnetic fields, in disagreement with previous cosmological MHD simulations which are successful to explain the observed magnetic fields in the ICM. On the other hand, temperature anisotropies can also drive plasma instabilities which can relax the anisotropy. This work aims to compare the relaxation rate with the growth rate of the anisotropies driven by the turbulence. We employ quasilinear theory to estimate the ions scattering rate due to the parallel firehose, mirror, and ion-cyclotron instabilities, for a set of plasma parameters resulting from AMHD simulations of the turbulent ICM. We show that the ICM turbulence can sustain only anisotropy levels very close to the instabilities thresholds. We argue that the AMHD model which bounds the anisotropies at the marginal stability levels can describe the Alfvenic turbulence cascade in the ICM.

preprint2016arXiv

Radio Synchrotron Fluctuation Statistics as a Probe of Magnetized Interstellar Turbulence

We investigate how observations of synchrotron intensity fluctuations can be used to probe the sonic and Alfvénic Mach numbers of interstellar turbulence, based on mock observations performed on simulations of magnetohydrodynamic turbulence. We find that the structure function slope, and a diagnostic of anisotropy that we call the integrated quadrupole ratio modulus, both depend on the Alfvénic Mach number. However, these statistics also depend on the orientation of the mean magnetic field in the synchrotron emitting region relative to our line of sight, and this creates a degeneracy that cannot be broken by observations of synchrotron intensity alone. We conclude that the polarization of synchrotron emission could be analyzed to break this degeneracy, and suggest that this will be possible with the Square Kilometre Array.

preprint2016arXiv

Stochastic reacceleration of relativistic electrons by turbulent reconnection: a mechanism for cluster-scale radio emission ?

In this paper we investigate a situation where relativistic particles are reaccelerated diffusing across regions of reconnection and magnetic dynamo in super-Alfvenic, incompressible large-scale turbulence. We present an exploratory study of this mechanism in the intra-cluster-medium (ICM). In view of large-scale turbulence in the ICM we adopt a reconnection scheme that is based on turbulent reconnection and MHD turbulence. In this case particles are accelerated and decelerated in a systematic way in reconnecting and magnetic-dynamo regions, respectively, and on longer time-scales undergo a stochastic process diffusing across these sites (similar to second-order Fermi). Our study extends on larger scales numerical studies that focused on the acceleration in and around turbulent reconnecting regions. We suggest that this mechanism may play a role in the reacceleration of relativistic electrons in galaxy clusters providing a new physical scenario to explain the origin of cluster-scale diffuse radio emission. Indeed differently from current turbulent reacceleration models proposed for example for radio halos this mechanism is based on the effect of large-scale incompressible and super-Alfvenic turbulence. In this new model turbulence governs the interaction between relativistic particles and magnetic field lines that diffuse, reconnect and are stretched in the turbulent ICM.

preprint2016arXiv

Study of velocity centroids based on the theory of fluctuations in position-position-velocity space

We study possibility of obtaining velocity spectra by studying turbulence in an optically thick medium using velocity centroids (VCs).We find that the regime of universal, i.e. independent of underlying turbulence statistics, fluctuations discovered originally within the velocity channel analysis (VCA) carries over to the statistics of VCs. In other words, for large absorptions the VC lose their ability to reflect the spectra of turbulence. Combining our present study with the earlier studies of centroids in Esquivel & Lazarian, we conclude that centroids are applicable for studies subsonic/transsonic turbulence for the range of scales that is limited by the absorption effects. We also consider VCs based on absorption lines and define the range of their applicability. We address the problem of analytical description of spectra and anisotropies of fluctuations that are available through studies using VC. We obtain spectra and anisotropy of VC fluctuations arising from Alfvén, slow and fast modes that constitute the compressible MHD cascade to address the issue of anisotropy of VC statistics, and show how the VC anisotropy can be used to find the media magnetization as well as to identify and separate contributions from Alfvén, slow and fast modes. Our study demonstrates that VCs are complementary to the tools provided by the VCA. In order to study turbulent volume for which the resolution of single dish telescopes is insufficient, we demonstrate how the studies of anisotropy can be performed using interferometers.

preprint2016arXiv

Superdiffusion of Cosmic Rays: Implications for Cosmic Ray Acceleration

Diffusion of cosmic rays (CRs) is the key process of understanding their propagation and acceleration. We employ the description of spatial separation of magnetic field lines in MHD turbulence in Lazarian & Vishniac (1999) to quantify the divergence of magnetic field on scales less than the injection scale of turbulence and show this divergence induces superdiffusion of CR in the direction perpendicular to the mean magnetic field. The perpendicular displacement squared increases, not as distance $x$ along magnetic field, which is the case for a regular diffusion, but as the $x^{3}$ for freely streaming CRs. The dependence changes to $x^{3/2}$ for the CRs propagating diffusively along magnetic field. In the latter case we show that it is important to distinguish the perpendicular displacement in respect to the mean field and to the local magnetic field. We consider how superdiffusion changes the acceleration of CRs in shocks and show how it decreases efficiency of the CRs acceleration in perpendicular shocks. We also demonstrate that in the case when small-scale magnetic field is being generated in the pre-shock region, an efficient acceleration can take place for the CRs streaming without collisions along magnetic loops.

preprint2016arXiv

Turbulent dynamo in a conducting fluid and partially ionized gas

By following the Kazantsev theory and taking into account both microscopic and turbulent diffusion of magnetic fields, we develop a unified treatment of the kinematic and nonlinear stages of turbulent dynamo, and study the dynamo process for a full range of magnetic Prandtl number Pm and ionization fractions. We find a striking similarity between the dependence of dynamo behavior on Pm in a conducting fluid and R (a function of ionization fraction) in partially ionized gas. In a weakly ionized medium, the kinematic stage is largely extended, including not only exponential growth but a new regime of dynamo characterized by linear-in-time growth of magnetic field strength, and the resulting magnetic energy is much higher than the kinetic energy carried by viscous-scale eddies. Unlike the kinematic stage, the subsequent nonlinear stage is unaffected by microscopic diffusion processes and has a universal linear-in-time growth of magnetic energy with the growth rate as a constant fraction $3/38$ of the turbulent energy transfer rate, showing good agreement with earlier numerical results. Applying the analysis to the first stars and galaxies, we find that the kinematic stage is able to generate a field strength only an order of magnitude smaller than the final saturation value. But the generation of large-scale magnetic fields can only be accounted for by the relatively inefficient nonlinear stage and requires longer time than the free-fall time. It suggests that magnetic fields may not have played a dynamically important role during the formation of the first stars.

preprint2015arXiv

Damping of MHD turbulence in partially ionized plasma: implications for cosmic ray propagation

We study the damping from neutral-ion collisions of both incompressible and compressible magnetohydrodynamic (MHD) turbulence in partially ionized medium. We start from the linear analysis of MHD waves applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are then used in determining the damping scales of MHD turbulence. The physical connection between the damping scale of MHD turbulence and cutoff boundary of linear MHD waves is investigated. Our analytical results are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and solar chromosphere. As a significant astrophysical utility, we introduce damping effects to propagation of cosmic rays in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyroresonance is identified.

preprint2015arXiv

Linewidth Differences of Neutrals and Ions Induced by MHD Turbulence

We address the problem of the difference of line widths of neutrals and ions observed from molecular clouds and explore whether this difference can arise from the effects of magnetohydrodynamic (MHD) turbulence acting on partially ionized gas. Among the three fundamental modes of MHD turbulence, we find fast modes do not contribute to linewidth differences, whereas slow modes can have an effect on different line widths for certain parameters. We focus on Alfvénic component because they contain most of the turbulent energy, and consider the damping of this component taking into account both neutral-ion collisions and neutral viscosity. We consider different regimes of turbulence corresponding to different media magnetizations and turbulent drivings. In the case of super-Alfvénic turbulence, when the damping scale of Alfvénic turbulence is below $l_A$, where $l_A$ is the injection scale of anisotropic GS95-type turbulence, the linewidth difference does not depend on the magnetic field strength. While for other turbulent regimes, the dependence is present. For instance, the difference between the squares of the neutral and ion velocity dispersions in strong sub-Alfvénic turbulence allows evaluation of magnetic field. We discuss earlier findings on the neutral-ion linewidth differences in the literature and compare the expressions for magnetic field we obtain with those published earlier.

preprint2015arXiv

On origin and destruction of relativistic dust and its implication for ultrahigh energy cosmic rays

Dust grains may be accelerated to relativistic speeds by radiation pressure, diffusive shocks, and other acceleration mechanisms. Such relativistic grains have been suggested as primary particles of ultrahigh energy cosmic rays (UHECRs). In this paper, we first revisit the problem of acceleration by radiation pressure and calculate maximum grain velocities achieved. We find that grains can be accelerated to relativistic speeds with Lorentz factor $γ< 2$ by powerful radiation sources, which is lower than earlier estimates showing that $γ$ could reach $\sim 10$. We then investigate different destruction mechanisms for relativistic grains traversing a variety of environments. In solar radiation, we find that the destruction by thermal sublimation and Coulomb explosions is important. We also quantify grain destruction due to electronic sputtering by ions and grain-grain collisions. Electronic sputtering is found to be rather inefficient, whereas the evaporation following grain-grain collisions is shown to be an important mechanism for which the $a \sim 0.01- 1μ$m grains would be destroyed after sweeping a gas column $N_{Coul}\sim 5\times 10^{19}-5\times 10^{20}cm^{-2}$. Relativistic dust in the interstellar medium and intergalactic medium (IGM) would be disrupted by Coulomb explosions due to collisional charging after traversing a gas column $N_{Coul} \sim 10^{17}cm^{-2}$ unless grain material is very strong. We show that photoelectric emission by optical and ultraviolet background radiation is also significant for the destruction of relativistic dust in the IGM. The obtained results indicate that relativistic dust from galaxies would be destroyed before reaching the Earth's atmosphere and unlikely to account for UHECRs.

preprint2015arXiv

Polarization of Magnetic Dipole Emission and Spinning Dust Emission from Magnetic Nanoparticles

Magnetic dipole emission (MDE) from interstellar magnetic nanoparticles is an important Galactic foreground in the microwave frequencies, and its polarization level may pose great challenges for achieving reliable measurements of cosmic microwave background (CMB) B-mode signal. To obtain theoretical constraints on the polarization of MDE, we first compute the degree of alignment of big silicate grains incorporated with magnetic inclusions. We find that, in realistic conditions of the interstellar medium, thermally rotating big grains with magnetic inclusions are weakly aligned and achieve {\it alignment saturation} when the magnetic alignment rate becomes much faster than the rotational damping rate. We then compute the degree of alignment for free-flying magnetic nanoparticles, taking into account various interaction processes of grains with the ambient gas and radiation field, including neutral collisions, ion collisions, and infrared emission. We find that the rotational damping by infrared emission can significantly decrease the degree of alignment of small particles from the saturation level, whereas the excitation by ion collisions can enhance the alignment of ultrasmall particles. Using the computed degrees of alignment, we predict the polarization level of MDE from free-flying magnetic nanoparticles to be rather low. Such a polarization level is within the upper limits measured for anomalous microwave emission (AME), which indicates that MDE from free-flying iron particles may not be ruled out as a source of AME. We also quantify spinning dust emission from free-flying iron nanoparticles with permanent magnetic moments and find that its emissivity is one order of magnitude lower than that from spinning polycyclic aromatic hydrocarbons (PAHs). Finally, we compute the polarization spectra of spinning dust emission from PAHs for the different interstellar magnetic fields.

preprint2015arXiv

Theory and Applications of Non-Relativistic and Relativistic Turbulent Reconnection

Realistic astrophysical environments are turbulent due to the extremely high Reynolds numbers. Therefore, the theories of reconnection intended for describing astrophysical reconnection should not ignore the effects of turbulence on magnetic reconnection. Turbulence is known to change the nature of many physical processes dramatically and in this review we claim that magnetic reconnection is not an exception. We stress that not only astrophysical turbulence is ubiquitous, but also magnetic reconnection itself induces turbulence. Thus turbulence must be accounted for in any realistic astrophysical reconnection setup. We argue that due to the similarities of MHD turbulence in relativistic and non-relativistic cases the theory of magnetic reconnection developed for the non-relativistic case can be extended to the relativistic case and we provide numerical simulations that support this conjecture. We also provide quantitative comparisons of the theoretical predictions and results of numerical experiments, including the situations when turbulent reconnection is self-driven, i.e. the turbulence in the system is generated by the reconnection process itself. We show how turbulent reconnection entails the violation of magnetic flux freezing, the conclusion that has really far reaching consequences for many realistically turbulent astrophysical environments. In addition, we consider observational testing of turbulent reconnection as well as numerous implications of the theory. The former includes the Sun and solar wind reconnection, while the latter include the process of reconnection diffusion induced by turbulent reconnection, the acceleration of energetic particles, bursts of turbulent reconnection related to black hole sources as well as gamma ray bursts. Finally, we explain why turbulent reconnection cannot be explained by turbulent resistivity or derived through the mean field approach.

preprint2014arXiv

Alfvenic Turbulence Beyond the Ambipolar Diffusion Scale

We investigate the nature of the Alfvénic turbulence cascade in two fluid MHD simulations in order to determine if turbulence is damped once the ion and neutral species become decoupled at a critical scale called the ambipolar diffusion scale (L$_{AD}$). Using mode decomposition to separate the three classical MHD modes, we study the second order structure functions of the Alfvén mode velocity field of both neutrals and ions in the reference frame of the local magnetic field. On scales greater than L$_{AD}$ we confirm that two fluid turbulence strongly resembles single fluid MHD turbulence. Our simulations show that the behavior of two fluid turbulence becomes more complex on scales less than L$_{AD}$. We find that Alfvenic turbulence can exist past L$_{AD}$ when the turbulence is globally super-Alfvénic, with the ions and neutrals forming separate cascades once decoupling has taken place. When turbulence is globally sub-Alfvenic and hence strongly anisotropic with a large separation between the parallel and perpendicular decoupling scales, turbulence is damped at L$_{AD}$. We also find that the power spectrum of the kinetic energy in the damped regime is consistent with a $k^{-4}$ scaling (in agreement with the predictions of Lazarian, Vishniac & Cho 2004).

preprint2014arXiv

Galactic interstellar turbulence in the southern sky seen through spatial gradients of the polarization vector

Radio synchrotron polarization maps of the Galaxy can be used to infer the properties of interstellar turbulence in the diffuse warm ionized medium (WIM). In this paper, we investigate the spatial gradient of linearly polarized synchrotron emission ($|\nabla\textbf{P}|/|\textbf{P}|$) as a tracer of turbulence, the relationship of the gradient to the sonic Mach number of the WIM, and changes in morphology of the gradient as a function of Galactic position in the southern sky. We use data from the S-band Polarization All Sky Survey (S-PASS) to image the spatial gradient of the linearly polarized synchrotron emission ($|\nabla \textbf{P}|/|\textbf{P}|$) of the entire southern sky at $2.3$~GHz. The spatial gradient of linear polarization reveals rapid changes of the density and magnetic fluctuations in the WIM due to magnetic turbulence as a function of Galactic position; we make comparisons of these data to ideal MHD numerical simulations. In order to constrain the sonic Mach number ($M_{s}$), we apply a high order moments analysis to the observations and to the simulated diffuse, isothermal ISM with ideal magneto-hydrodynamic turbulence. We find that polarization gradient maps reveal elongated structures, which we associate with turbulence in the ISM. Our analysis corroborates the view of a turbulent WIM in a transonic regime $M_{s}\lesssim 2$. Filamentary structures with typical widths down to the angular resolution are seen and the observed morphologies match closely with numerical simulations and in some cases H$α$ contours. The $|\nabla \textbf{P}|/|\textbf{P}|$ intensity is found to be approximately log-normal distributed. No systematic variations of the sonic Mach number are observed as a function of Galactic coordinates, which is consistent with turbulence in the WIM inferred by the analysis of H$α$ data.

preprint2014arXiv

Grain Alignment by Radiative Torques in Special Conditions and Implications

Grain alignment by radiative torques (RATs) has been extensively studied for various environment conditions, including interstellar medium, dense molecular clouds, and accretion disks, thanks to significant progress in observational, theoretical and numerical studies. In this paper, we explore the alignment by RATs and provide quantitative predictions of dust polarization for a set of astrophysical environments that can be tested observationally. We first consider the alignment of grains in the local interstellar medium and compare predictions for linear polarization by aligned grains with recent observational data for nearby stars. We then revisit the problem of grain alignment in accretions disks by taking into account the dependence of RAT alignment efficiency on the anisotropic direction of radiation fields relative to magnetic fields. Moreover, we study the grain alignment in interplanetary medium, including diffuse Zodiacal cloud and cometary comae, and calculate the degree of circular polarization (CP) of scattered light arising from single scattering by aligned grains. We also discuss a new type of grain alignment, namely the alignment with respect to the ambient electric field instead of the alignment with the magnetic field. We show that this type of alignment can allow us to reproduce the systematic features of CP observed across a cometary coma. Our findings suggest that polarized Zodiacal dust emission may be an important polarized foreground component, which should be treated carefully in cosmic microwave background experiments.

preprint2014arXiv

Measuring the Alfvenic Nature of the Interstellar Medium: Velocity Anisotropy Revisited

The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel & Lazarian method to estimate the Alfven Mach number using the structure function anisotropy in velocity centroid data from position-position-velocity maps. We utilize 3D magnetohydrodynamic (MHD) simulations of fully developed turbulence, with a large range of sonic and Alfvenic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfven Mach number dependency found in Esquivel & Lazarian (2011) might change when taking the second moment of the position-position-velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e. Alfven Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line-of-sight is parallel to up to ~45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvenic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques.

preprint2014arXiv

Paramagnetic alignment of small grains: a novel method for measuring interstellar magnetic fields

We present a novel method to measure the strength of interstellar magnetic fields based on ultraviolet (UV) polarization of starlight, which is in part produced by weakly aligned, small interstellar grains. We begin with calculating degrees of alignment of small (size $a\sim 0.01μ$m) and very small ($a\sim 0.001μ$m) grains in the interstellar magnetic field due to the Davis-Greenstein paramagnetic relaxation and resonance paramagnetic relaxation. We compute the degrees of paramagnetic alignment with the ambient magnetic field $B$ using Langevin equations. In this paper, we take into account various processes essential for the dynamics of small grains, including infrared (IR) emission, electric dipole emission, plasma drag and collisions with neutral and ionized species. We find that the alignment of small grains is necessary to reproduce the observed polarization in the UV, although the polarization arising from these small grains is negligible at the optical and IR wavelengths. Based on fitting theoretical models to observed extinction and polarization curves, we find that the best-fit model requires a higher degree of alignment of small grains for the case with the peak wavelength of polarization $λ_{\max}<0.55μ$m, which exhibits an excess UV polarization relative to the Serkowski law, compared to the typical case $λ_{\max}=0.55μ$m. We interpret the correlation between the systematic increase of the UV polarization relative to maximum polarization (i.e. of $p(6μm^{-1})/p_{\max}$) with $λ_{\max}^{-1}$ by appealing to the higher degree of alignment of small grains. We identify paramagnetic relaxation as the cause of the alignment of small grains and utilize the dependence of the degree of alignment on the magnetic field strength $B$ to suggest a new way to measure $B$ using the observable parameters $λ_{\max}$ and $p(6μm^{-1})/p_{\max}$.[Abridged]

preprint2014arXiv

Reconnection in turbulent astrophysical fluids

Magnetic reconnection is a fundamental process of magnetic field topology change. We analyze the connection of this process with turbulence which is ubiquitous in astrophysical environments. We show how Lazarian & Vishniac (1999) model of turbulent reconnection is connected to the experimentally proven concept of Richardson diffusion and discuss how turbulence violates the generally accepted notion of magnetic flux freezing. We note that in environments that are laminar initially turbulence can develop as a result of magnetic reconnection and this can result in flares of magnetic reconnection in magnetically dominated media. In particular, magnetic reconnection can initially develop through tearing, but the transition to the turbulent state is expected for astrophysical systems. We show that turbulent reconnection predictions corresponds to the Solar and solar wind data.

preprint2013arXiv

Astrophysical hydromagnetic turbulence

Recent progress in astrophysical hydromagnetic turbulence is being reviewed. The physical ideas behind the now widely accepted Goldreich--Sridhar model and its extension to compressible magnetohydrodynamic turbulence are introduced. Implications for cosmic ray diffusion and acceleration is being discussed. Dynamo-generated magnetic fields with and without helicity are contrasted against each other. Certain turbulent transport processes are being modified and often suppressed by anisotropy and inhomogeneities of the turbulence, while others are being produced by such properties, which can lead to new large-scale instabilities of the turbulent medium. Applications of various such processes to astrophysical systems are being considered.

preprint2013arXiv

Collapse of Turbulent Cores and Reconnection Diffusion

For a molecular cloud clump to form stars some transport of magnetic flux is required from the denser, inner regions to the outer regions of the cloud, otherwise this can prevent the collapse. Fast magnetic reconnection which takes place in the presence of turbulence can induce a process of reconnection diffusion (RD). Extending earlier numerical studies of reconnection diffusion in cylindrical clouds, we consider more realistic clouds with spherical gravitational potentials and also account for the effects of the gas self-gravity. We demonstrate that within our setup RD is efficient. We have also identified the conditions under which RD becomes strong enough to make an initially subcritical cloud clump supercritical and induce its collapse. Our results indicate that the formation of a supercritical core is regulated by a complex interplay between gravity, self-gravity, the magnetic field strength and nearly transonic and trans-Alfvénic turbulence, confirming that RD is able to remove magnetic flux from collapsing clumps, but only a few of them become nearly critical or supercritical, sub-Alfvénic cores, which is consistent with the observations. Besides, we have found that the supercritical cores built up in our simulations develop a predominantly helical magnetic field geometry which is also consistent with observations. Finally, we have evaluated the effective values of the turbulent reconnection diffusion coefficient and found that they are much larger than the numerical diffusion, especially for initially trans-Alfvénic clouds, ensuring that the detected magnetic flux removal is due to to the action of the RD rather than to numerical diffusivity.

preprint2013arXiv

Constraint on the Polarization of Electric Dipole Emission from Spinning Dust

Planck results have revealed that the electric dipole emission from polycyclic aromatic hydrocarbons (PAHs) is the most reliable explanation for anomalous microwave emission that interferes with cosmic microwave background (CMB) radiation experiments. The emerging question is to what extent this emission component contaminates to the polarized CMB radiation. We present constraints on polarized dust emission for the model of grain size distribution and grain alignment that best fits to observed extinction and polarization curves. Two stars with a prominent polarization feature at wavelength 2175 Angstrom, HD 197770 and HD 147933-4, are chosen for our study. For HD 197770, we find that the model with aligned silicate grains plus weakly aligned PAHs can successfully reproduce the 2175 Angstrom polarization feature; whereas, for HD 147933-4, we find that the alignment of only silicate grains can account for that feature. The alignment function of PAHs for the best-fit model to the HD 197770 data is employed to constrain polarized spinning dust emission. We find that the degree of polarization of spinning dust emission is about 1.6 percent at frequency ~ 3 GHz and declines to below 0.9 percent for frequency above 20 GHz. We also predict the degree of polarization of thermal dust emission at 353 GHz to be ~ 11 percent and 14 percent for the lines of sight to the HD 197770 and HD 147933-4 stars, respectively.

preprint2013arXiv

Hierarchical Structure of Magnetohydrodynamic Turbulence In Position-Position-Velocity Space

Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic Position-Position-Velocity (PPV) emission cubes of optically thin isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the physics of the gas, i.e. self-gravity and the global sonic and Alfvenic Mach number. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observed a strong sonic and Alfvenic dependency when we apply the the statistical moments (i.e. mean, variance, skewness, kurtosis) to the dendrogram distribution. Larger magnetic field and sonic Mach number correspond to larger values of the moments. Application of the dendrogram to 3D density cubes, also known as Position-Position-Position cubes (PPP), reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions in the direction of obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP.

preprint2013arXiv

Imbalanced Relativistic Force-Free Magnetohydrodynamic Turbulence

When magnetic energy density is much larger than that of matter, as in pulsar/black hole magnetospheres, the medium becomes force-free and we need relativity to describe it. As in non-relativistic magnetohydrodynamics (MHD), Alfvénic MHD turbulence in the relativistic limit can be described by interactions of counter-traveling wave packets. In this paper we numerically study strong imbalanced MHD turbulence in such environments. Here, imbalanced turbulence means the waves traveling in one direction (dominant waves) have higher amplitudes than the opposite-traveling waves (sub-dominant waves). We find that (1) spectrum of the dominant waves is steeper than that of sub-dominant waves, (2) the anisotropy of the dominant waves is weaker than that of sub-dominant waves, and (3) the dependence of the ratio of magnetic energy densities of dominant and sub-dominant waves on the ratio of energy injection rates is steeper than quadratic (i.e., \$b_+^2/b_-^2 \propto (ε_+/ε_-)^n \$ with n>2). These results are consistent with those obtained for imbalanced non-relativistic Alfvénic turbulence. This corresponds well to the earlier reported similarity of the relativistic and non-relativistic balanced magnetic turbulence.

preprint2013arXiv

Magnetic field amplification and evolution in turbulent collisionless MHD: an application to the ICM

The amplification and maintenance of the observed magnetic fields in the ICM are usually attributed to the turbulent dynamo action. This is generally derived employing a collisional MHD model. However, in the ICM the ion mean free path between collisions is of the order of the dynamical scales, thus requiring a collisionless MHD description. Unlike collisional MHD simulations, our study uses an anisotropic plasma pressure with respect to the direction of the local magnetic field, which brings the plasma within a parameter space where collisionless instabilities should take place. Within the adopted model these instabilities are contained at bay through the relaxation term of the pressure anisotropy which simulates the feedback of the mirror and firehose instabilities. This relaxation acts to get the plasma distribution function consistent with the empirical studies of collisionless plasmas. Our 3D numerical simulations of forced transonic turbulence motivated by modeling of the turbulent ICM are performed for different initial values of the magnetic field intensity, and different relaxation rates of the pressure anisotropy. We found that in the high beta plasma regime corresponding to the ICM conditions, a fast anisotropy relaxation rate gives results which are similar to the collisional-MHD model as far as the statistical properties of the turbulence are concerned. Also, the amplification of seed magnetic fields due to the turbulent dynamo action is similar to the collisional-MHD model. Our simulations that do not employ the anisotropy relaxation deviate significantly from the collisional-MHD results, and show more power at the small-scale fluctuations of both density and velocity representing the results of the instabilities. For these simulations the large scale fluctuations in the magnetic field are mostly suppressed and the turbulent dynamo fails in amplifying seed magnetic fields.

preprint2013arXiv

Magnetic Field Measurement with Ground State Alignment

Observational studies of magnetic fields are crucial. We introduce a process "ground state alignment" as a new way to determine the magnetic field direction in diffuse medium. The alignment is due to anisotropic radiation impinging on the atom/ion. The consequence of the process is the polarization of spectral lines resulting from scattering and absorption from aligned atomic/ionic species with fine or hyperfine structure. The magnetic field induces precession and realign the atom/ion and therefore the polarization of the emitted or absorbed radiation reflects the direction of the magnetic field. The atoms get aligned at their low levels and, as the life-time of the atoms/ions we deal with is long, the alignment induced by anisotropic radiation is susceptible to extremely weak magnetic fields ($1{\rm G}\gtrsim B\gtrsim 10^{-15}$G). In fact, the effects of atomic/ionic alignment were studied in the laboratory decades ago, mostly in relation to the maser research. Recently, the atomic effect has been already detected in observations from circumstellar medium and this is a harbinger of future extensive magnetic field studies. A unique feature of the atomic realignment is that they can reveal the 3D orientation of magnetic field. In this article, we shall review the basic physical processes involved in atomic realignment. We shall also discuss its applications to interplanetary, circumstellar and interstellar magnetic fields. In addition, our research reveals that the polarization of the radiation arising from the transitions between fine and hyperfine states of the ground level can provide a unique diagnostics of magnetic fields in the Epoch of Reionization.

preprint2013arXiv

Magnetic Reconnection in Turbulent Plasmas and Gamma Ray Bursts

We discuss how the model of magnetic reconnection in the presence of turbulence proposed inLazarian & Vishniac 1999 makes the reconnection rate independent either of resistivity or microscopic plasma effects, but determined entirely by the magnetic field line wandering induced by turbulence. We explain that the model accounts for both fast and slow regimes of reconnection and that this property naturally induces flares of reconnection in low beta plasma environments. In addition, we show that the model involves volume reconnection which can convert a substantial part of the energy into energetic particles. It is important that the reconnection induces an efficient acceleration of the first order Fermi type. Finally, we relate the properties of the reconnection with the observed properties of gamma ray bursts and provide evidence supporting the explanation of gamma ray bursts based on energy release via reconnection.

preprint2013arXiv

Reconnection Diffusion, Star Formation and Numerical Simulations

We consider fast magnetic reconnection that takes place within turbulent magnetic flux and show that the process results in diffusion of magnetic fields and matter, which we term reconnection diffusion. The process of reconnection diffusion is based on the model of 3D reconnection of weakly turbulent magnetic fields and is applicable to both fully ionized and partially ionized gas. The rate of reconnection diffusion does not depend on the level of ionization and therefore the usually employed ambipolar diffusion idea gets irrelevant for magnetic field transport in turbulent fluids. We claim that the reconnection diffusion process is a manifestation of the violation of flux conservation in highly conducting turbulent fluids. We discuss the consequences of reconnection diffusion for star formation and stress. We show that reconnection diffusion on large scales is independent of small scale magnetic field dynamics of magnetic fields. We conclude that numerical simulations correctly represents the diffusion of actual astrophysical magnetic fields in flows with substantially larger Lundquist numbers if these simulated regions regions are turbulent.

preprint2013arXiv

The Effects of Radiative Transfer on the PDFs of Molecular MHD Turbulence

We study the effects of radiative transfer on the Probability Distribution Functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied $^{13}$CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have $τ\approx 1$ best matching the PDF of the column density. We fit a 2D variance-sonic Mach number relationship to our logarithmic PDFs of the form $σ_{ln(Σ/Σ_0)}^2=A\times ln(1+b^2{\cal M}_s^2)$ and find that, for parameter $b=1/3$, parameter $A$ depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the $τ\approx 1$ case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always sub-thermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed and consequently, the observable PDFs, are less determined by turbulence and more-often determined by radiative transfer effects.

preprint2013arXiv

The Turbulence Power Spectrum in Optically Thick Interstellar Clouds

The Fourier power spectrum is one of the most widely used statistical tools to analyze the nature of magnetohydrodynamic turbulence in the interstellar medium. Lazarian & Pogosyan (2004) predicted that the spectral slope should saturate to -3 for an optically thick medium and many observations exist in support of their prediction. However, there have not been any numerical studies to-date testing these results. We analyze the spatial power spectrum of MHD simulations with a wide range of sonic and Alfvénic Mach numbers, which include radiative transfer effects of the $^{13}$CO transition. We confirm numerically the predictions of Lazarian & Pogosyan (2004) that the spectral slope of line intensity maps of an optically thick medium saturates to -3. Furthermore, for very optically thin supersonic CO gas, where the density or CO abundance values are too low to excite emission in all but the densest shock compressed gas, we find that the spectral slope is shallower than expected from the column density. Finally, we find that mixed optically thin/thick CO gas, which has average optical depths on order of unity, shows mixed behavior: for super-Alfvénic turbulence, the integrated intensity power spectral slopes generally follow the same trend with sonic Mach number as the true column density power spectrum slopes. However, for sub-Alfvénic turbulence the spectral slopes are steeper with values near -3 which are similar to the very optically thick regime.

preprint2012arXiv

A Technique for Foreground Subtraction in Redshifted 21 cm Observations

One of the main challenges for future 21 cm observations is to remove foregrounds which are several orders of magnitude more intense than the HI signal. We propose a new technique for removing foregrounds of the redshifted 21 cm observations. We consider multi-frequency interferometer observations. We assume that the 21 cm signals in different frequency channels are uncorrelated and the foreground signals change slowly as a function of frequency. When we add the visibilities of all channels, the foreground signals increase roughly by a factor of ~N because they are highly correlated. However, the 21 cm signals increase by a factor of ~\sqrt{N} because the signals in different channels contribute randomly. This enables us to obtain an accurate shape of the foreground angular power spectrum. Then, we obtain the 21-cm power spectrum by subtracting the foreground power spectrum obtained this way. We describe how to obtain the average power spectrum of the 21 cm signal.

preprint2012arXiv

Acceleration of Energetic Particles through Reconnection of Weakly Stochastic Magnetic Field

Astrophysical media are turbulent and therefore reconnection should be treated in the presence of pre-existing turbulence. We consider the model of fast magnetic reconnection in Lazarian & Vishniac (1999) which predicts that the rate of reconnection is controlled by the intensity and the injection scale of turbulent motions. We provide new evidence of successful testing of the model and argue that the model presents a generic set up for astrophysical reconnection events. We study particle acceleration that takes place in volumes of driven turbulence as well turbulent volumes in the presence of large scale reconnection. We show that in the latter case the acceleration is of the first order Fermi type thus supporting the model of acceleration proposed in Gouveia dal Pino & Lazarian (2005).

preprint2012arXiv

Anisotropy of TeV Cosmic Rays and the Outer Heliospheric Boundaries

Cosmic rays in the energy range from about 10's GeV to several 100's TeV are observed on Earth with an energy-dependent anisotropy of order 0.01-0.1%, and a consistent topology that appears to significantly change at higher energy. The nearest and most recent galactic cosmic ray sources might stochastically dominate the observation and possibly explain a change in orientation of the anisotropy as a function of energy. However, the diffusion approximation is not able to explain its non-dipolar structure and, in particular, the significant contribution of small angular scale features. Particle propagation within the mean free path in the local interstellar medium might have a major role in determining the properties of galactic cosmic rays, such as their arrival distribution. In particular, scattering on perturbations induced in the local interstellar magnetic field by the heliosphere wake, may cause a re-distribution of anisotropic cosmic rays below about 100 TeV toward the direction of the elongated heliotail and of the local interstellar magnetic field in the outer heliosphere. Such scattering processes are considered responsible of the observed TeV cosmic ray global anisotropy and its fine angular structure.

preprint2012arXiv

Astrophysical Reconnection and Particle Acceleration

Astrophysical reconnection takes place in a turbulent medium. The turbulence in most cases is pre-existing, not caused by the reconnection itself. The model of magnetic reconnection in Lazarian & Vishniac (1999) predicts that in the presence of turbulence the reconnection becomes fast, i.e. it is independent of resistivity, but dependent on the level of turbulence. Magnetic reconnection injects energy into plasmas through a turbulent outflow from the reconnection region and this outflow can enhance the level of turbulence creating bursts of reconnection. Magnetic reconnection in the presence of turbulence can accelerate energetic particles through the first order Fermi mechanism, as was discussed in Gouveia dal Pino & Lazarian (2005). We discuss new numerical results on particle acceleration in turbulent reconnection, compare the acceleration arising from turbulent reconnection to the acceleration of energetic particles in turbulent medium.

preprint2012arXiv

Cosmic Rays and Stochastic Magnetic Reconnection in the Heliotail

Galactic cosmic rays are believed to be generated by diffusive shock acceleration processes in Supernova Remnants, and the arrival direction is likely determined by the distribution of their sources throughout the Galaxy, in particular by the nearest and youngest ones. Transport to Earth through the interstellar medium is expected to affect the cosmic ray properties as well. However, the observed anisotropy of TeV cosmic rays and its energy dependence cannot be explained with diffusion models of particle propagation in the Galaxy. Within a distance of a few parsec, diffusion regime is not valid and particles with energy below about 100 TeV must be influenced by the heliosphere and its elongated tail. The observation of a highly significant localized excess region of cosmic rays from the apparent direction of the downstream interstellar flow at 1-10 TeV energies might provide the first experimental evidence that the heliotail can affect the transport of energetic particles. In particular, TeV cosmic rays propagating through the heliotail interact with the 100-300 AU wide magnetic field polarity domains generated by the 11 year cycles. Since the strength of non-linear convective processes is expected to be larger than viscous damping, the plasma in the heliotail is turbulent. Where magnetic field domains converge on each other due to solar wind gradient, stochastic magnetic reconnection likely occurs. Such processes may be efficient enough to re-accelerate a fraction of TeV particles as long as scattering processes are not strong. Therefore the fractional excess of TeV cosmic rays from the narrow region toward the heliotail direction traces sightlines with the lowest smearing scattering effects, that can also explain the observation of a harder than average energy spectrum.

preprint2012arXiv

Disc formation in turbulent cloud cores: is magnetic flux loss necessary to stop the magnetic braking catastrophe or not?

Recent numerical analysis of Keplerian disk formation in turbulent, magnetized cloud cores by Santos-Lima, de Gouveia Dal Pino, & Lazarian (2012) demonstrated that reconnection diffusion is an efficient process to remove the magnetic flux excess during the build up of a rotationally supported disk. This process is induced by fast reconnection of the magnetic fields in a turbulent flow. In a similar numerical study, Seifried et al. (2012) concluded that reconnection diffusion or any other non-ideal MHD effects would not be necessary and turbulence shear alone would provide a natural way to build up a rotating disk without requiring magnetic flux loss. Their conclusion was based on the fact that the mean mass-to-flux ratio (μ) evaluated over a spherical region with a radius much larger than the disk is nearly constant in their models. In this letter we compare the two sets of simulations and show that this averaging over large scales can mask significant real increases of μ in the inner regions where the disk is built up. We demonstrate that turbulence-induced reconnection diffusion of the magnetic field happens in the initial stages of the disk formation in the turbulent envelope material that is accreting. Our analysis is suggestive that reconnection diffusion is present in both sets of simulations and provides a simple solution for the "magnetic braking catastrophe" which is discussed in the literature in relation to the formation of protostellar accretion disks.

preprint2012arXiv

Fast Reconnection and Reconnection Diffusion: Implications for Star Formation

Fast reconnection of magnetic field in turbulent fluids allows magnetic field to change its topology and connections. As a result, the traditional concept of magnetic fields being frozen into the plasma is no longer applicable. The diffusion of plasmas and magnetic field is enabled by reconnection and therefore is termed "reconnection diffusion". We explore the consequences of reconnection diffusion for star formation. In the paper we explain the physics of reconnection diffusion both from macroscopic and microscopic points of view. We quantify the reconnection diffusion rate both for weak and strong MHD turbulence and address the problem of reconnection diffusion acting together with ambipolar diffusion. In addition, we provide a criterion for correctly representing the magnetic diffusivity in simulations of star formation. We show that the role of the plasma effects is limited to "breaking up lines" on small scales and does not affect the rate of reconnection diffusion. We address the existing observational results and demonstrate how reconnection diffusion can explain the puzzles presented by observations, in particular, the observed higher magnetization of cloud cores in comparison with the magnetization of envelopes. We also outline a possible set of observational tests of the reconnection diffusion concept and discuss how the application of the new concept changes our understanding of star formation and its numerical modeling. Finally, we outline the differences of the process of reconnection diffusion and the process of accumulation of matter along magnetic field lines that is frequently invoked to explain the results of numerical simulations

preprint2012arXiv

Heat Transfer and Reconnection Diffusion in Turbulent Magnetized Plasmas

It is well known that magnetic fields constrain motions of charged particles, impeding the diffusion of charged particles perpendicular to magnetic field direction. This modification of transport processes is of vital importance for a wide variety of astrophysical processes including cosmic ray transport, transfer of heavy elements in the interstellar medium, star formation etc. Dealing with these processes one should keep in mind that in realistic astrophysical conditions magnetized fluids are turbulent. In this review we single out a single transport process, namely, heat transfer and consider how it occurs in the presence of the magnetized turbulence. We show that the ability of magnetic field lines to constantly change topology and connectivity is at the heart of the correct description of the 3D magnetic field stochasticity in turbulent fluids. This ability is ensured by fast magnetic reconnection in turbulent fluids and puts forward the concept of reconnection diffusion at the core of the physical picture of heat transfer in astrophysical plasmas. Appealing to reconnection diffusion we describe the ability of plasma to diffuse between different magnetized eddies explaining the advection of the heat by turbulence. Adopting the structure of magnetic field that follows from the modern understanding of MHD turbulence, we also discuss thermal conductivity that arises as electrons stream along stochastic magnetic field lines. We compare the effective heat transport that arise from the two processes and conclude that in many astrophysically-motivated cased eddy advection of heat dominates. Finally, we discuss the concepts of sub and superdiffusion and show that the subdiffusion requires rather restrictive settings. At the same time, accelerated diffusion or superdiffusion of heat is possible on the scales less than the injection scale of the turbulence.

preprint2012arXiv

Magnetization of cloud cores and envelopes and other observational consequences of reconnection diffusion

Recent observational results for magnetic fields in molecular clouds reviewed by Crutcher (2012) seem to be inconsistent with the predictions of the ambipolar diffusion theory of star formation. These include the measured decrease in mass to flux ratio between envelopes and cores, the failure to detect any self-gravitating magnetically subcritical clouds, the determination of the flat PDF of the total magnetic field strengths implying that there are many clouds with very weak magnetic fields, and the observed scaling $B \propto ρ^{2/3}$ that implies gravitational contraction with weak magnetic fields. We consider the problem of magnetic field evolution in turbulent molecular clouds and discuss the process of magnetic field diffusion mediated by magnetic reconnection. For this process that we termed "reconnection diffusion" we provide a simple physical model and explain that this process is inevitable in view of the present day understanding of MHD turbulence. We address the issue of the expected magnetization of cores and envelopes in the process of star formation and show that reconnection diffusion provides an efficient removal of magnetic flux that depends only on the properties of MHD turbulence in the core and the envelope. As a result, the magnetic flux trapped during the collapse in the envelope is being released faster than the flux trapped in the core, resulting in much weaker fields in envelopes than in cores, as observed. We provide simple semi-analytical model calculations which support this conclusion and qualitatively agree with the observational results. We argue that magnetic reconnection provides a solution to the magnetic flux problem of star formation that agrees better with observations than the long-standing ambipolar diffusion paradigm.

preprint2012arXiv

Spinning Dust Emission from Wobbling Grains: Important Physical Effects and Implications

We review major progress on the modeling of electric dipole emission from rapidly spinning tiny dust grains, including polycyclic aromatic hydrocarbons (PAHs). We begin by summarizing the original model of spinning dust proposed by Draine and Lazarian and recent theoretical results improving the Draine and Lazarian model. The review is focused on important physical effects that were disregarded in earlier studies for the sake of simplicity and recently accounted for by us, including grain wobbling due to internal relaxation, impulsive excitation by single-ion collisions, the triaxiality of grain shape, charge fluctuations, and the turbulent nature of astrophysical environments. Implications of the spinning dust emission for constraining physical properties of tiny dust grains and environment conditions are discussed. We discuss the alignment of tiny dust grains and possibility of polarized spinning dust emission. Suggestions for constraining the alignment of tiny grains and polarization of spinning dust emission are also discussed.

preprint2012arXiv

The Column Density Variance-Sonic Mach Number Relationship

Although there are a wealth of column density tracers for both the molecular and diffuse interstellar medium, there are very few observational studies investigating the relationship between the density variance ($σ^2$) and the sonic Mach number (${\cal M}_s$). This is in part due to the fact that the $σ^2$-${\cal M}_s$ relationship is derived, via MHD simulations, for the 3D density variance only, which is not a direct observable. We investigate the utility of a 2D column density $σ_{Σ/Σ_0}^2$-${\cal M}_s$ relationship using solenoidally driven isothermal MHD simulations and find that the best fit follows closely the form of the 3D density $σ_{ρ/ρ_0}^2$-${\cal M}_s$ trend but includes a scaling parameter $A$ such that: $σ_{ln(Σ/Σ_0)}^2=A\times ln(1+b^2{\cal M}_s^2)$, where $A=0.11$ and $b=1/3$. This relation is consistent with the observational data reported for the Taurus and IC 5146 molecular clouds with $b=0.5$ and $A=0.16$ and $b=0.5$ and $A=0.12$, respectively. These results open up the possibility of using the 2D column density values of $σ^2$ for investigations of the relation between the sonic Mach number and the PDF variance in addition to existing PDF sonic Mach number relations.

preprint2012arXiv

Tracing Magnetic Fields with Ground State Alignment

Observational studies of magnetic fields are vital as magnetic fields play a crucial role in various astrophysical processes, including star formation, accretion of matter, transport processes (e.g., transport of heat), and cosmic rays. We identified a process "ground state alignment" as a new way to determine the magnetic field direction in diffuse medium. The alignment is due to anisotropic radiation impinging on the atom/ion, while the magnetic field induces precession and realign the atom/ion and therefore the polarization of the emitted or absorbed radiation reflects the direction of the magnetic field. The atoms get aligned at their low levels and, as the life-time of the atoms/ions we deal with is long, the alignment induced by anisotropic radiation is susceptible to weak magnetic fields ($1{\rm G}\gtrsim B\gtrsim 10^{-15}$G). Compared to the upper level Hanle effect, atomic realignment is most suitable for the studies of magnetic field in the diffuse medium, where magnetic field is relatively weak. In fact, the effects of atomic/ionic alignment, including the realignment in magnetic field, were studied in the laboratory decades ago, mostly in relation to the maser research. Recently, the atomic effect has been already detected in observations from circumstellar medium and this is a harbinger of future extensive magnetic field studies. A unique feature of the atomic realignment is that they can reveal the 3D orientation of magnetic field. In this article, we shall review the basic physical processes involved in atomic realignment and its applications to interplanetary, circumstellar and interstellar magnetic fields. In addition, our research reveals that the polarization of the radiation arising from the transitions between fine and hyperfine states of the ground level can provide a unique diagnostics of magnetic fields, including those in the Early Universe.

preprint2012arXiv

Turbulence, Magnetic Reconnection in Turbulent Fluids and Energetic Particle Acceleration

Turbulence is ubiquitous in astrophysics. It radically changes many astrophysical phenomena, in particular, the propagation and acceleration of cosmic rays. We present the modern understanding of compressible magnetohydrodynamic (MHD) turbulence, in particular its decomposition into Alfvén, slow and fast modes, discuss the density structure of turbulent subsonic and supersonic media, as well as other relevant regimes of astrophysical turbulence. All this information is essential for understanding the energetic particle acceleration that we discuss further in the review. For instance, we show how fast and slow modes accelerate energetic particles through the second order Fermi acceleration, while density fluctuations generate magnetic fields in pre-shock regions enabling the first order Fermi acceleration of high energy cosmic rays. Very importantly, however, the first order Fermi cosmic ray acceleration is also possible in sites of magnetic reconnection. In the presence of turbulence this reconnection gets fast and we present numerical evidence supporting the predictions of the Lazarian & Vishniac (1999) model of fast reconnection. The efficiency of this process suggests that magnetic reconnection can release substantial amounts of energy in short periods of time. As the particle tracing numerical simulations show that the particles can be efficiently accelerated during the reconnection, we argue that the process of magnetic reconnection may be much more important for particle acceleration than it is currently accepted. In particular, we discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays measured by Voyagers as well as the origin cosmic ray excess in the direction of Heliotail.

preprint2011arXiv

Acceleration of Small Dust Grains due to Charge Fluctuations

We consider the acceleration of very small dust grains including Polycyclic Aromatic Hydrocarbons (PAHs) arising from the electrostatic interactions of dust grains that have charge fluctuations in time due to charging events. We simulate the charge fluctuations of very small grains due to their sticking collisions with electrons and ions in plasma, and the emission of photoelectrons by UV photons using Monte Carlo method. We identify the acceleration induced by the charge fluctuations as the dominant acceleration mechanism of very small grains in the diffuse interstellar medium (ISM). We show that this acceleration mechanism is efficient for environments with low degree of ionization (i.e., large Debye length), where the charge fluctuations are slow but have a large amplitude. We also discuss the implications of the present mechanism for grain coagulation and shattering in the diffuse ISM, molecular clouds, and protoplanetary disks.

preprint2011arXiv

Cosmic Ray streaming from SNRs and gamma ray emission from nearby molecular clouds

High-energy gamma ray emission has been detected recently from supernovae remnants (SNRs) and their surroundings. The existence of molecular clouds near some of the SNRs suggests that the gamma rays originate predominantly from p-p interactions with cosmic rays accelerated at a closeby SNR shock wave. Here we investigate the acceleration of cosmic rays and the gamma ray production in the cloud self-consistently by taking into account the interactions of the streaming instability and the background turbulence both at the shock front and in the ensuing propagation to the clouds. We focus on the later evolution of SNRs, when the conventional treatment of the streaming instability is valid but the magnetic field is enhanced due to either Bell's current instability and/or due to the dynamo generation of magnetic field in the precursor region. We calculate the time dependence of the maximum energy of the accelerated particles. This result is then used to determine the diffusive flux of the runaway particles escaping the shock region, from which we obtain the gamma spectrum consistent with observations. Finally, we check the self-consistency of our results by comparing the required level of diffusion with the level of the streaming instability attainable in the presence of turbulence damping. The energy range of cosmic rays subject to the streaming instability is able to produce the observed energy spectrum of gamma rays.

preprint2011arXiv

Cosmic Ray transport in MHD turbulence: large and small scale interactions

Cosmic ray (CR) transport and acceleration is essential for many astrophysical problems, e.g., CMB foreground, ionization of molecular clouds and all high energy phenomena. Recent advances in MHD turbulence call for revisions in the paradigm of cosmic ray transport. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulation, in which turbulence is injected at large scale and cascades to to small scales. We shall address the issue of the transport of CRs, both parallel and perpendicular to the magnetic field and show that the issue of cosmic ray subdiffusion is only important for restricted cases when the ambient turbulence is far from that suggested by numerical simulations. Moreover, on scales less than injection scale of turbulence, CRs's transport becomes super-diffusive. We also shall discuss the nonlinear growth of kinetic gyroresonance instability of cosmic rays induced by large scale compressible turbulence. This gyroresonance of cosmic rays on turbulence is demonstrated an important scattering mechanism in addition to direct interaction with the compressible turbulence. The feedback of the instability on large scale turbulence cannot be neglected, and should be included in future simulations.

preprint2011arXiv

Cosmic Ray transport through gyroresonance instability in compressible turbulence

We study the nonlinear growth of kinetic gyroresonance instability of cosmic rays (CRs) induced by large scale compressible turbulence. This feedback of cosmic rays on turbulence was shown to induce an important scattering mechanism in addition to direct interaction with the compressible turbulence. The linear growth is bound to saturate due to the wave-particle interactions. By balancing increase of CR anisotropy via the large scale compression and its decrease via the wave-particle scattering, we find the steady state solutions. The nonlinear suppression due to the wave-particle scattering limit the energy range of CRs that can excite the instabilities and be scattered by the induced slab waves. The direct interaction with large scale compressible modes still appears to be the dominant mechanism for isotropization of high energy cosmic rays (> 100 GeV).

preprint2011arXiv

Dynamo in the Intra-Cluster Medium: Simulation of CGL-MHD Turbulent Dynamo

The standard magnetohydrodynamic (MHD) description of the plasma in the hot, magnetized gas of the intra-cluster (ICM) medium is not adequate because it is weakly collisional. In such collisionless magnetized gas, the microscopic velocity distribution of the particles is not isotropic, giving rise to kinetic effects on the dynamical scales. These kinetic effects could be important in understanding the turbulence, as so as the amplification and maintenance of the magnetic fields in the ICM. It is possible to formulate fluid models for collisonless or weakly collisional gas by introducing modifications in the MHD equations. These models are often referred as kinetic MHD (KMHD). Using a KMHD model based on the CGL-closure, which allows the adiabatic evolution of the two components of the pressure tensor (the parallel and perpendicular components with respect to the local magnetic field), we performed 3D numerical simulations of forced turbulence in order to study the amplification of an initially weak seed magnetic field. We found that the growth rate of the magnetic energy is comparable to that of the ordinary MHD turbulent dynamo, but the magnetic energy saturates in a level smaller than of the MHD case. We also found that a necessary condition for the dynamo works is to impose limits to the anisotropy of the pressure.

preprint2011arXiv

Evolution and lifetime of transient clumps in the turbulent interstellar medium

We study the evolution of dense clumps and provide argument that the existence of the clumps is not limited by the crossing time of the clump. We claim that the lifetimes of the clumps are determined by the turbulent motions on larger scale and predict the correlation of the clump lifetime and its column density. We use numerical simulations and successfully test this relation. In addition, we study the morphological asymmetry and the magnetization of the clumps as a function of their masses.

preprint2011arXiv

Interstellar Sonic and Alfvénic Mach Numbers and the Tsallis Distribution

In an effort to characterize the Mach numbers of ISM magnetohydrodynamic (MHD) turbulence, we study the probability distribution functions (PDFs) of patial increments of density, velocity, and magnetic field for fourteen ideal isothermal MHD simulations at resolution 512^3. In particular, we fit the PDFs using the Tsallis function and study the dependency of fit parameters on the compressibility and magnetization of the gas. We find that the Tsallis function fits PDFs of MHD turbulence well, with fit parameters showing sensitivities to the sonic and Alfven Mach numbers. For 3D density, column density, and position-position-velocity (PPV) data we find that the amplitude and width of the PDFs shows a dependency on the sonic Mach number. We also find the width of the PDF is sensitive to global Alfvenic Mach number especially in cases where the sonic number is high. These dependencies are also found for mock observational cases, where cloud-like boundary conditions, smoothing, and noise are introduced. The ability of Tsallis statistics to characterize sonic and Alfvenic Mach numbers of simulated ISM turbulence point to it being a useful tool in the analysis of the observed ISM, especially when used simultaneously with other statistical techniques.

preprint2011arXiv

Magnetohydrodynamic Simulations of Reconnection and Particle Acceleration: Three-Dimensional Effects

The magnetic fields can change their topology through a process known as magnetic reconnection. This process in not only important for understanding the origin and evolution of the large-scale magnetic field, but is seen as a possibly efficient particle accelerator producing cosmic rays mainly through the first order Fermi process. In this work we study the properties of particle acceleration in reconnection zones and show that the velocity component parallel to the magnetic field of test particles inserted in nearly non-resistive magnetohydrodynamic (MHD) domains of reconnection without including kinetic effects, such as pressure anisotropy, the Hall term, or anomalous effects, increases exponentially. Also, the acceleration of the perpendicular component is always possible in such models. We have found that within contracting magnetic islands or current sheets the particles accelerate predominantly through the first order Fermi process, as previously described, while outside the current sheets and islands the particles experience mostly drift acceleration due to magnetic fields gradients. Considering two dimensional MHD models without a guide field, we find that the parallel acceleration stops at some level. This saturation effect is however removed in the presence of an out-of-plane guide field or in three dimensional models. Therefore, we stress the importance of the guide field and fully three dimensional studies for a complete understanding of the process of particle acceleration in astrophysical reconnection environments.

preprint2011arXiv

MHD turbulence-Star Formation Connection: from pc to kpc scales

The transport of magnetic flux to outside of collapsing molecular clouds is a required step to allow the formation of stars. Although ambipolar diffusion is often regarded as a key mechanism for that, it has been recently argued that it may not be efficient enough. In this review, we discuss the role that MHD turbulence plays in the transport of magnetic flux in star forming flows. In particular, based on recent advances in the theory of fast magnetic reconnection in turbulent flows, we will show results of three-dimensional numerical simulations that indicate that the diffusion of magnetic field induced by turbulent reconnection can be a very efficient mechanism, especially in the early stages of cloud collapse and star formation. To conclude, we will also briefly discuss the turbulence-star formation connection and feedback in different astrophysical environments: from galactic to cluster of galaxy scales.

preprint2011arXiv

Power of Turbulent Reconnection: Star Formation, Acceleration of Cosmic Rays, Heat Transfer, Flares and Gamma Ray Bursts

Turbulence is ubiquitous in astrophysical fluids. Therefore it is necessary to study magnetic reconnection in turbulent environments. The model of fast turbulent reconnection proposed in Lazarian & Vishniac 1999 has been successfully tested numerically and it suggests numerous astrophysical implications. Those include a radically new possibility of removing magnetic field from collapsing clouds which we termed "reconnection diffusion", acceleration of cosmic rays within shrinking filaments of reconnected magnetic fields, flares of reconnection, from solar flares to much stronger ones which can account for gamma ray bursts. In addition, the model reveals a very intimate relation between magnetic reconnection and properties of strong turbulence, explaining how turbulent eddies can transport heat in magnetized plasmas. This is a small fraction the astrophysical implications of the quantitative insight into the fundamental process of magnetic reconnection in turbulent media.

preprint2011arXiv

Properties of Interstellar Turbulence from Gradients of Linear Radio Polarization Maps

Faraday rotation of linearly polarized radio signals provides a very sensitive probe of fluctuations in the interstellar magnetic field and ionized gas density resulting from magnetohydrodynamic (MHD) turbulence. We used a set of statistical tools to analyze images of the spatial gradient of linearly polarized radio emission ($|\nabla \textbf{P}|$) from the ISM for both observational data from a test image of the Southern Galactic Plane Survey (SGPS) and isothermal simulations of MHD turbulence. We compared the observational data with results of synthetic observations obtained with the simulations of 3D turbulence. Visually, in both data sets, a complex network of filamentary structures is seen. Our analysis shows that the filaments in the gradient can be produced by shocks as well as random fluctuations characterizing the non-differentiable field of MHD turbulence. The latter dominates for subsonic turbulence, while the former dominates for supersonic turbulence. In order to quantitatively characterize these differences we use the topology tool known as a genus curve as well as the moments of the image distribution. We find that higher values for the moments correspond to cases of $|\nabla \textbf{P}|$ with larger Mach numbers, but the strength of the dependency is connected to the telescope angular resolution. In regards to the topology, the supersonic filaments observed in $|\nabla \textbf{P}|$ have a positive genus shift, which indicates a "swisscheese" like topology, while the subsonic cases show a negative genus, indicating a "clump" like topology. In the case of the genus, the dependency on the telescope resolution is not as strong. The SGPS test region data has a distribution and morphology that matches subsonic to transsonic type turbulence, which independently confirms what is now expected for the WIM.

preprint2011arXiv

Relation of Astrophysical Turbulence and Magnetic Reconnection

Astrophysical fluids are generically turbulent and this must be taken into account for most transport processes. We discuss how the preexisting turbulence modifies magnetic reconnection and how magnetic reconnection affects the MHD turbulent cascade. We show the intrinsic interdependence and interrelation of magnetic turbulence and magnetic reconnection, in particular, that strong magnetic turbulence in 3D requires reconnection and 3D magnetic turbulence entails fast reconnection. We follow the approach in Eyink, Lazarian & Vishniac 2011 to show that the expressions of fast magnetic reconnection in Lazarian & Vishniac 1999 can be recovered if Richardson diffusion of turbulent flows is used instead of ordinary Ohmic diffusion. This does not revive, however, the concept of magnetic turbulent diffusion which assumes that magnetic fields can be mixed up in a passive way down to a very small dissipation scales. On the contrary, we are dealing the reconnection of dynamically important magnetic field bundles which strongly resist bending and have well defined mean direction weakly perturbed by turbulence. We argue that in the presence of turbulence the very concept of flux-freezing requires modification. The diffusion that arises from magnetic turbulence can be called reconnection diffusion as it based on reconnection of magnetic field lines. The reconnection diffusion has important implications for the continuous transport processes in magnetized plasmas and for star formation. In addition, fast magnetic reconnection in turbulent media induces the First order Fermi acceleration of energetic particles, can explain solar flares and gamma ray bursts. However, the most dramatic consequence of these developments is the fact that the standard flux freezing concept must be radically modified in the presence of turbulence.

preprint2011arXiv

Revisiting Acceleration of Charged Grains in Magnetohydrodynamic Turbulence

We study the acceleration of charged grains by magnetohydrodynamics (MHD) turbulence in the interstellar medium (ISM). We begin with revisiting gyroresonance acceleration by taking into account the fluctuations of grain guiding center along a uniform magnetic field (i.e. nonlinear theory--NLT). We calculate grain velocities due to gyroresonance by fast MHD modes using the NLT for different phases of the ISM, and compare with results obtained using quasi-linear theory (QLT). We find for the parameters applicable to the typical ISM phases that the fluctuations of grain guiding center reduce grain velocities by less than 15 percent, but they can be important for more special circumstances. We confirm that large grains can be accelerated to super-Alfvenic velocities through gyroresonance. For such super-Alfvenic grains, we investigate the effect of further acceleration via transit time damping (TTD) by fast modes. We find that due to the broadening of resonance condition in the NLT, the TTD acceleration is not only important for the cosines of grain pitch angle relative to the magnetic field mu>V_{A}/v, but also for mu<V_{A}/v where v is the grain velocity and V_{A} is the Alfven speed. We show that the TTD acceleration is dominant over the gyroresonance for large grains, and can increase substantially grain velocities induced by gyroresonance acceleration. We quantify another stochastic acceleration mechanism arising from low frequency Alfven waves. We discuss the range of applicability of the mechanisms and their implications.

preprint2011arXiv

Spinning Dust Emission: Effects of irregular grain shape, transient heating and comparison with WMAP results

Planck is expected to answer crucial questions on the early Universe, but it also provides further understanding on anomalous microwave emission. Electric dipole emission from spinning dust grains continues to be the favored interpretation of anomalous microwave emission. In this paper, we present a method to calculate the rotational emission from small grains of irregular shape with moments of inertia $I_{1}> I_{2}> I_{3}$. We show that a torque-free rotating irregular grain with a given angular momentum radiates at multiple frequency modes. The resulting spinning dust spectrum has peak frequency and emissivity increasing with the degree of grain shape irregularity, which is defined by $I_{1}:I_{2}:I_{3}$. We discuss how the orientation of dipole moment $\bmu$ in body coordinates affects the spinning dust spectrum for different regimes of internal thermal fluctuations. We show that the spinning dust emissivity for the case of strong thermal fluctuations is less sensitive to the orientation of $\bmu$ than in the case of weak thermal fluctuations. We calculate spinning dust spectra for a range of gas density and dipole moment. The effect of compressible turbulence on spinning dust emission intensity is investigated. We show that the emission intensity in a turbulent medium increases by a factor from 1.2-1.4 relative to that in a uniform medium, as sonic Mach number $M_{s}$ increases from 2-7. Finally, spinning dust parameters are constrained by fitting our improved model to five-year {\it Wilkinson Microwave Anisotropy Probe} cross-correlation foreground spectra, for both the H$α$-correlated and 100 $μ$m-correlated emission spectra.

preprint2011arXiv

The energy source of the filaments around the giant galaxy NGC1275

The brightest galaxy in the nearby Perseus cluster, NGC1275, is surrounded by a network of filaments. These were first observed through their Halpha emission but are now known to have a large molecular component with a total mass approaching 10^11Msun of gas. The filaments are embedded in hot intracluster gas and stretch over 80 kpc. They have an unusual low excitation spectrum which is well modelled by collisional heating and ionization by secondary electrons. Here we note that the surface radiative flux from the outer filaments is close to the energy flux impacting on them from particles in the hot gas. We propose that the secondary electrons within the cold filaments, which excite the observed submillimetre through UV emission, are due to the hot surrounding gas efficiently penetrating the cold gas through reconnection diffusion. Some of the soft X-ray emission seen from the filaments is then due to charge exchange, although this is insufficient to account for all the observed X-ray flux. The filaments are complex with multiphase gas. Interpenetration of hot and cold gas leads to the filaments growing in mass, at a rate of up to 100Msunpyr. The lack of soft X-ray cooling emission in cool core clusters is then due to the non-radiative cooling of hot gas on mixing with cold gas around and within the central galaxy.

preprint2011arXiv

The role of turbulent magnetic reconnection on the formation of rotationally supported protostellar disks

The formation of protostellar disks out of molecular cloud cores is still not fully understood. Under ideal MHD conditions, the removal of angular momentum from the disk progenitor by the typically embedded magnetic field may prevent the formation of a rotationally supported disk during the main protostellar accretion phase of low mass stars. This has been known as the magnetic braking problem and the most investigated mechanism to alleviate this problem and help removing the excess of magnetic flux during the star formation process, the so called ambipolar diffusion (AD), has been shown to be not sufficient to weaken the magnetic braking at least at this stage of the disk formation. In this work, motivated by recent progress in the understanding of magnetic reconnection in turbulent environments, we appeal to the diffusion of magnetic field mediated by magnetic reconnection as an alternative mechanism for removing magnetic flux. We investigate numerically this mechanism during the later phases of the protostellar disk formation and show its high efficiency. By means 3D MHD simulations, we show that this mechanism is able to transport magnetic flux to the outskirts of the disk progenitor at time scales compatible with the collapse, allowing the formation of a rotationally supported disk around the protostar of dimensions ~100 AU. Since MHD turbulence is expected to be present in protostellar disks, this is a natural mechanism for removing magnetic flux excess and allowing the formation of these disks. This mechanism dismiss the necessity of postulating a hypothetical increase of the Ohmic resistivity as discussed in the literature. Together with our earlier work which showed that magnetic flux removal from molecular cloud cores is very efficient, this work calls for reconsidering the relative role of AD for the processes of star and planet formation.

preprint2011arXiv

Turbulence, Reconnection and Cosmic Rays in Galaxy Clusters

Recent years have been marked by substantial changes in our understanding of magnetic turbulence and magnetic reconnection, which, in its turn induced better understanding of cosmic ray diffusion and acceleration. Current models of magnetized turbulence are no more ad hoc constructions, but numerically tested theories. In this very short review we summarize topics presented in two talks given at the conference and provide a brief sketch of the vast and rapidly developing field. We discuss how turbulence decreases the efficient mean free path of the particles in the collisionless plasmas in galaxy clusters and claim that this makes MHD turbulence description applicable to a wider range of scales. We discuss the properties of MHD turbulence and its relation to magnetic reconnection. Finally, we overview how turbulence induces particle acceleration via second order Fermi process and affects first order Fermi acceleration in shocks and reconnection regions.

preprint2010arXiv

Acceleration of primary and secondary particles in galaxy clusters by compressible MHD turbulence: from radio halos to gamma rays

Radio observations discovered large scale non thermal sources in the central Mpc regions of dynamically disturbed galaxy clusters (radio halos). The morphological and spectral properties of these sources suggest that the emitting electrons are accelerated by spatially distributed and gentle mechanisms, providing some indirect evidence for turbulent acceleration in the inter-galactic-medium (IGM). Only deep upper limits to the energy associated with relativistic protons in the IGM have been recently obtained through gamma and radio observations. Yet these protons should be (theoretically) the main non-thermal particle component in the IGM implying the unavoidable production, at some level, of secondary particles that may have a deep impact on the gamma ray and radio properties of galaxy clysters. Following Brunetti and Lazarian (2007), in this paper we consider the advances in the theory of MHD turbulence to develop a comprehensive picture of turbulence in the IGM and extend our previous calculations of particle acceleration by compressible MHD turbulence by considering self-consistently the reacceleration of both primary and secondary particles. Under these conditions we expect that radio to gamma ray emission is generated from galaxy clusters with a complex spectrum that depends on the dynamics of the thermal gas and Dark Matter. The non-thermal emission results in very good agreement with radio observations and with present constraints from hard X-ray and gamma ray observations. In our model giant radio halos are generated in merging (turbulent) clusters only. However, in case secondaries dominate the electron component in the IGM, we expect that the level of the Mpc-scale synchrotron emission in more relaxed clusters is already close to that of the radio upper limits derived by present observations of clusters without radio halos.

preprint2010arXiv

Acceleration of small astrophysical grains due to charge fluctuations

We discuss a novel mechanism of dust acceleration which may dominate for particles smaller than $\sim0.1~μ$m. The acceleration is caused by their direct electrostatic interactions arising from fluctuations of grain charges. The energy source for the acceleration are the irreversible plasma processes occurring on the grain surfaces. We show that this mechanism of charge-fluctuation-induced acceleration likely affects the rate of grain coagulation and shattering of the population of small grains.

preprint2010arXiv

Damping of MHD turbulence in partially ionized gas and the observed difference of velocities of neutrals and ions

Theoretical and observational studies on the turbulence of the interstellar medium developed fast in the past decades. The theory of supersonic magnetized turbulence, as well as the understanding of projection effects of observed quantities, are still in progress. In this work we explore the characterization of the turbulent cascade and its damping from observational spectral line profiles. We address the difference of ion and neutral velocities by clarifying the nature of the turbulence damping in the partially ionized. We provide theoretical arguments in favor of the explanation of the larger Doppler broadening of lines arising from neutral species compared to ions as arising from the turbulence damping of ions at larger scales. Also, we compute a number of MHD numerical simulations for different turbulent regimes and explicit turbulent damping, and compare both the 3-dimensional distributions of velocity and the synthetic line profile distributions. From the numerical simulations, we place constraints on the precision with which one can measure the 3D dispersion depending on the turbulence sonic Mach number. We show that no universal correspondence between the 3D velocity dispersions measured in the turbulent volume and minima of the 2D velocity dispersions available through observations exist. For instance, for subsonic turbulence the correspondence is poor at scales much smaller than the turbulence injection scale, while for supersonic turbulence the correspondence is poor for the scales comparable with the injection scale. We provide a physical explanation of the existence of such a 2D-3D correspondence and discuss the uncertainties in evaluating the damping scale of ions that can be obtained from observations. However, we show that the statistics of velocity dispersion from observed line profiles can provide the spectral index and the energy transfer rate of turbulence. Also, comparing two similar simulations with different viscous coefficients it was possible to constrain the turbulent cut-off scale. This may especially prove useful since it is believed that ambipolar diffusion may be one of the dominant dissipative mechanism in star-forming regions. In this case, the determination of the ambipolar diffusion scale may be used as a complementary method for the determination of magnetic field intensity in collapsing cores. We discuss the implications of our findings in terms of a new approach to magnetic field measurement proposed by Li & Houde (2008).

preprint2010arXiv

Diffusion of magnetic field and removal of magnetic flux from clouds via turbulent reconnection

The diffusion of astrophysical magnetic fields in conducting fluids in the presence of turbulence depends on whether magnetic fields can change their topology via reconnection 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 magnetic reconnection is concerned. Our studies of magnetic field diffusion in turbulent medium reveal interesting new phenomena. In the presence of gravity and turbulence, our 3D simulations show the decrease of the magnetic flux-to-mass ratio as the gaseous density at the center of the gravitational potential increases. We observe this effect both in the situations when we start with equilibrium distributions of gas and magnetic field and when we follow the evolution of collapsing dynamically unstable configurations. Thus the process of turbulent magnetic field removal should be applicable both to quasi-static subcritical molecular clouds and cores and violently collapsing supercritical entities. The increase of the gravitational potential as well as the magnetization of the gas increases the segregation of the mass and magnetic flux in the saturated final state of the simulations, supporting the notion that the reconnection-enabled diffusivity relaxes the magnetic field + gas system in the gravitational field to its minimal energy state. This effect is expected to play an important role in star formation, from its initial stages of concentrating interstellar gas to the final stages of the accretion to the forming protostar.

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

Fast Magnetic Reconnection and Energetic Particle Acceleration

Our numerical simulations show that the reconnection of magnetic field becomes fast in the presence of weak turbulence in the way consistent with the Lazarian and Vishniac (1999) model of fast reconnection. We trace particles within our numerical simulations and show that the particles can be efficiently accelerated via the first order Fermi acceleration. We discuss the acceleration arising from reconnection as a possible origin of the anomalous cosmic rays measured by Voyagers.

preprint2010arXiv

Improving the model of emission from spinning dust: effects of grain wobbling and transient spin-up

Observations continue to support the interpretation of the anomalous microwave foreground as electric dipole radiation from spinning dust grains as proposed by Draine and Lazarian (1998ab). In this paper we present a refinement of the original model by improving the treatment of a number of physical effects. First, we consider a disk-like grain rotating with angular velocity at an arbitrary angle with respect to the grain symmetry axis and derive the rotational damping and excitation coefficients arising from infrared emission, plasma-grain interactions and electric dipole emission. The angular velocity distribution and the electric dipole emission spectrum for grains is calculated using the Langevin equation, for cases both with and without fast internal relaxation. Our results show that, the peak emissivity of spinning dust, compared to earlier studies, increases by a factor of ~2 for the Warm Neutral Medium (WNM), the Warm Ionized Medium (WIM), the Cold Neutral Medium (CNM) and the Photodissociation Region (PDR), and by a factor ~4 for Reflection Nebulae (RN). The frequency at the emission peak also increases by factors ~1.4 to ~2 for these media. The increased emission and peak frequency result from the non-sphericity of grain shape and from the anisotropy in damping and excitation along directions parallel and perpendicular to the grain symmetry axis. Second, we provide a detailed numerical study including transient spin-up of grains by single-ion collisions. The impulses broaden the emission spectrum and increase the peak emissivity for the CNM, WNM and WIM. In addition, we present an improved treatment of rotational excitation and damping by infrared emission.

preprint2010arXiv

Magnetic Field Effects near the launching region of Astrophysical Jets

One of the fundamental properties of astrophysical magnetic fields is their ability to change topology through reconnection and in doing so, to release magnetic energy, sometimes violently. In this work, we review recent results on the role of magnetic reconnection and associated heating and particle acceleration in jet/accretion disk systems, namely young stellar objects (YSOs), microquasars, and active galactic nuclei (AGNs).

preprint2010arXiv

Magnetic reconnection as the cause of cosmic ray excess from the heliospheric tail

The observation of a broad excess of sub-TeV cosmic rays compatible with the direction of the heliospheric tail (Nagashima et al. 1998) and the discovery of two significant localized excess regions of multi-TeV cosmic rays by the MILAGRO collaboration (Abdo et al. 2008), also from the same region of the sky, have raised questions on their origin. In particular, the coincidence of the most significant localized region with the direction of the heliospheric tail and the small angular scale of the observed anisotropy (~ 10deg) is suggestive a local origin and of a possible connection to the low energy broad excess. Cosmic ray acceleration from magnetic reconnection in the magnetotail is proposed as a possible source of the energetic particles.

preprint2010arXiv

Numerical study of Cosmic Ray Diffusion in MHD turbulence

We study diffusion of Cosmic Rays (CRs) in turbulent magnetic fields using test particle simulations. Electromagnetic fields are produced in direct numerical MHD simulations of turbulence and used as an input for particle tracing, particle feedback on turbulence being ignored. Statistical transport coefficients from the test particle runs are compared with earlier analytical predictions. We find qualitative correspondence between them in various aspects of CR diffusion. In the incompressible case, that we consider in this paper, the dominant scattering mechanism occurs to be the non-resonant mirror interactions with the slow-mode perturbations. Perpendicular transport roughly agrees with being produced by magnetic field wandering.

preprint2010arXiv

Particle reacceleration by compressible turbulence in galaxy clusters: effects of reduced mean free path

Direct evidence for in situ particle acceleration mechanisms in the inter-galactic-medium (IGM) is provided by the diffuse Mpc--scale synchrotron emissions observed from galaxy clusters. It has been proposed that MHD turbulence, generated during cluster-cluster mergers, may be a source of particle reacceleration in the IGM. Calculations of turbulent acceleration must account self-consistently for the complex non--linear coupling between turbulent waves and particles. This has been calculated in some detail under the assumption that turbulence interacts in a collisionless way with the IGM. In this paper we explore a different picture of acceleration by compressible turbulence in galaxy clusters, where the interaction between turbulence and the IGM is mediated by plasma instabilities and maintained collisional at scales much smaller than the Coulomb mean free path. In this regime most of the energy of fast modes is channeled into the reacceleration of relativistic particles and the acceleration process approaches a universal behaviour being self-regulated by the back-reaction of the accelerated particles on turbulence itself. Assuming that relativistic protons contribute to several percent (or less) of the cluster energy, consistent with the FERMI observations of nearby clusters, we find that compressible turbulence at the level of a few percent of the thermal energy can reaccelerate relativistic electrons at GeV energies, that are necessary to explain the observed diffuse radio emission in the form of giant radio halos.

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

Scaling laws and diffuse locality of balanced and imbalanced MHD turbulence

The search of ways to generalize the theory of strong MHD turbulence for the case of non-zero cross-helicity (or energy imbalance) has attracted considerable interest recently. In our earlier publications we performed three-dimensional numerical simulations and showed that some of existing models are inconsistent with numerics. In this paper we focused our attention on low-imbalance limit and performed new high-resolution simulations. The results strongly suggest that in the limit of small imbalances we smoothly transition to a standard Goldreich-Sridhar (1995) balanced model. We also claim that Perez-Boldyrev (2009) model that predicts the same nonlinear timescale for both components due to so-called "dynamic alignment" strongly contradicts numerical evidence.

preprint2010arXiv

Turbulence in collisionless plasmas: statistical analysis from numerical simulations with pressure anisotropy

In the past years we have experienced an increasing interest in understanding of the physical properties of collisionless plasmas, mostly because of the large number of astrophysical environments, e.g. the intracluster medium (ICM), containing magnetic fields which are strong enough to be coupled with the ionized gas and characterized by densities sufficiently low to prevent the pressure isotropization with respect to the magnetic line direction. Under these conditions a new class of kinetic instabilities arises, such as firehose and mirror ones, which were extensively studied in the literature. Their role in the turbulence evolution and cascade process in the presence of pressure anisotropy, however, is still unclear. In this work we present the first statistical analysis of turbulence in collisionless plasmas using three dimensional double isothermal magnetohydrodynamical with the Chew-Goldberger-Low closure (CGL-MHD) numerical simulations. We study models with different initial conditions to account for the firehose and mirror instabilities and to obtain different turbulent regimes. We study the probability distribution functions, spectra, structure functions and anisotropy of density and velocity fluctuations. The results indicate that in some cases the instabilities significantly modifies the statistical properties of turbulence and even though preliminary and restricted to very specific conditions, show that the physical properties of turbulence in collisionless plasmas, as those found in the ICM, may be very different from what has been largely believed. Implications can range from interchange of energies to cosmic rays acceleration.

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', 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.

preprint2010arXiv

Velocity Spectrum for HI at High Latitudes

In this paper we present the results of the statistical analysis of high-latitude HI turbulence in the Milky Way. We have observed HI in the 21 cm line, obtained with the Arecibo L-Band Feed Array (ALFA) receiver at the Arecibo radio telescope. For recovering of velocity statistics we have used the Velocity Coordinate Spectrum (VCS) technique. In our analysis we have used direct fitting of the VCS model, as its asymptotic regimes are questionable for Arecibo's resolution and given the restrictions from thermal smoothing of the turbulent line. We have obtained a velocity spectral index $3.87 \pm 0.11$, an injection scale of $140 \pm 80$ pc, and an HI cold phase temperature of $52 \pm 11$ K. The spectral index is steeper than the Kolmogorov index and can be interpreted as being due to shock-dominated turbulence.

preprint2009arXiv

Reconnection in weakly stochastic B-fields in 2D

We study two dimensional turbulent magnetic reconnection in a compressible fluid in the gas pressure dominated limit. We use open boundary conditions and start from a Harris current sheet configuration with a uniform total pressure. A small perturbation to the vector potential initiates laminar reconnection at the Sweet-Parker rate, which is allowed to evolve for several dynamical times. Subsequently sub-Alfvenic turbulence is produced through random forcing at small wave numbers. The magnetic field topology near the current sheet is strongly affected by the turbulence. However, we find that the resulting reconnection speed depends on the resistivity. In contrast to previous results in three dimensions, we find no evidence for fast reconnection. The reconnection speed exhibits large variations but the time averages increase smoothly with the strength of the turbulence.

preprint2009arXiv

Tsallis statistics as a tool for studying interstellar turbulence

We used magnetohydrodynamic (MHD) simulations of interstellar turbulence to study the probability distribution functions (PDFs) of increments of density, velocity, and magnetic field. We found that the PDFs are well described by a Tsallis distribution, following the same general trends found in solar wind and Electron MHD studies. We found that the PDFs of density are very different in subsonic and supersonic turbulence. In order to extend this work to ISM observations we studied maps of column density obtained from 3D MHD simulations. From the column density maps we found the parameters that fit to Tsallis distributions and demonstrated that these parameters vary with the sonic and Alfvén Mach numbers of turbulence. This opens avenues for using Tsallis distributions to study the dynamical and perhaps magnetic states of interstellar gas.

preprint2009arXiv

Turbulence and the formation of filaments, loops and shock fronts in NGC 1275 in the Perseus Galaxy Cluster

NGC1275, the central galaxy in the Perseus cluster, is the host of gigantic hot bipolar bubbles inflated by AGN jets observed in the radio as Perseus A. It presents a spectacular $Hα$-emitting nebulosity surrounding NGC1275, with loops and filaments of gas extending to over 50 kpc. The origin of the filaments is still unknown, but probably correlates with the mechanism responsible for the giant buoyant bubbles. We present 2.5 and 3-dimensional MHD simulations of the central region of the cluster in which turbulent energy, possibly triggered by star formation and supernovae (SNe) explosions is introduced. The simulations reveal that the turbulence injected by massive stars could be responsible for the nearly isotropic distribution of filaments and loops that drag magnetic fields upward as indicated by recent observations. Weak shell-like shock fronts propagating into the ICM with velocities of 100-500 km/s are found, also resembling the observations. The isotropic outflow momentum of the turbulence slows the infall of the intracluster medium, thus limiting further starburst activity in NGC1275. As the turbulence is subsonic over most of the simulated volume, the turbulent kinetic energy is not efficiently converted into heat and additional heating is required to suppress the cooling flow at the core of the cluster. Simulations combining the MHD turbulence with the AGN outflow can reproduce the temperature radial profile observed around NGC1275. While the AGN mechanism is the main heating source, the supernovae are crucial to isotropize the energy distribution.

preprint2009arXiv

Turbulence-induced magnetic fields and the structure of Cosmic Ray modified shocks

We propose a model for Diffusive Shock Acceleration (DSA) in which stochastic magnetic fields in the shock precursor are generated through purely fluid mechanisms of a so-called small-scale dynamo. This contrasts with previous DSA models that considered magnetic fields amplified through cosmic ray streaming instabilities; i.e., either by way of individual particles resonant scattering in the magnetic fields, or by macroscopic electric currents associated with large-scale cosmic ray streaming. Instead, in our picture, the solenoidal velocity perturbations that are required for the dynamo to work are produced through the interactions of the pressure gradient of the cosmic ray precursor and density perturbations in the inflowing fluid. Our estimates show that this mechanism provides fast growth of magnetic field and is very generic. We argue that for supernovae shocks the mechanism is capable of generating upstream magnetic fields that are sufficiently strong for accelerating cosmic rays up to around 10^16 eV. No action of any other mechanism is necessary.

preprint2008arXiv

CMBPol Mission Concept Study: Foreground Science Knowledge and Prospects

We report on our knowledge of Galactic foregrounds, as well as on how a CMB satellite mission aiming at detecting a primordial B-mode signal (CMBPol) will contribute to improving it. We review the observational and analysis techniques used to constrain the structure of the Galactic magnetic field, whose presence is responsible for the polarization of Galactic emissions. Although our current understanding of the magnetized interstellar medium is somewhat limited, dramatic improvements in our knowledge of its properties are expected by the time CMBPol flies. Thanks to high resolution and high sensitivity instruments observing the whole sky at frequencies between 30 GHz and 850 GHz, CMBPol will not only improve this picture by observing the synchrotron emission from our galaxy, but also help constrain dust models. Polarized emission from interstellar dust indeed dominates over any other signal in CMBPol's highest frequency channels. Observations at these wavelengths, combined with ground-based studies of starlight polarization, will therefore enable us to improve our understanding of dust properties and of the mechanism(s) responsible for the alignment of dust grains with the Galactic magnetic field. CMBPol will also shed new light on observations that are presently not well understood. Morphological studies of anomalous dust and synchrotron emissions will indeed constrain their natures and properties, while searching for fluctuations in the emission from heliospheric dust will test our understanding of the circumheliospheric interstellar medium. Finally, acquiring more information on the properties of extra-Galactic sources will be necessary in order to maximize the cosmological constraints extracted from CMBPol's observations of CMB lensing. (abridged)

preprint2008arXiv

Density Studies of MHD Interstellar Turbulence: Statistical Moments, Correlations and Bispectrum

We present a number of statistical tools for obtaining studying turbulence in molecular clouds and diffuse interstellar medium. For our tests we used of three-dimensional 512 cube compressible MHD isothermal simulations performed for different sonic and Alfvenic Mach numbers. We introduce the bispectrum, a new tool for statistical studies of the interstellar medium which, unlike an ordinary power spectrum of turbulence, preserves the phase information of the stochastic field. We show that the bispectra of the 3D stochastic density field and of column densities, available from observations, are similar. We use the bispectrum technique to define the role of non-linear wave-wave interactions in the turbulent energy cascade. We also obtained the bispectrum function for density and column densities with varying magnetic field strength. Larger values of sonic Mach number result in increased correlations for modes with different wavenumbers. This effect becomes more evident with increasing magnetic field intensity. In addition to the bispectrum, we calculated the 3rd and 4th statistical moments of density and column density, namely, skewness and kurtosis, respectively. We found a strong dependence of skewness and kurtosis with the sonic Mach number. In particular, as this number increases, so does the asymmetry of the density distribution. We also studied the correlations of 2D column density with dispersion of velocities and magnetic field, as well as the correlations of 3D density with magentic and kinetic energy and Alfven Mach number for comparison. Our results show that column density is linearly correlated with magnetic field for high sonic Mach number. This trend is independent of the turbulent kinetic energy and can be used to characterize inhomogeneities of physical properties in low density clumps in the ISM.

preprint2008arXiv

Growth of Magnetic Fields Induced by Turbulent Motions

We present numerical simulations of driven magnetohydrodynamic (MHD) turbulence with weak/moderate imposed magnetic fields. The main goal is to clarify dynamics of magnetic field growth. We also investigate the effects of the imposed magnetic fields on the MHD turbulence, including, as a limit, the case of zero external field. Our findings are as follows. First, when we start off simulations with weak mean magnetic field only (or with small scale random field with zero imposed field), we observe that there is a stage at which magnetic energy density grows linearly with time. Runs with different numerical resolutions and/or different simulation parameters show consistent results for the growth rate at the linear stage. Second, we find that, when the strength of the external field increases, the equilibrium kinetic energy density drops by roughly the product of the rms velocity and the strength of the external field. The equilibrium magnetic energy density rises by roughly the same amount. Third, when the external magnetic field is not very strong (say, less than ~0.2 times the rms velocity when measured in the units of Alfven speed), the turbulence at large scales remains statistically isotropic, i.e. there is no apparent global anisotropy of order B_0/v. We discuss implications of our results on astrophysical fluids.

preprint2008arXiv

Obtaining Spectra of Turbulent Velocity from Observations

We discuss a long-standing problem of how turbulence can be studied using observations of Doppler broadened emission and absorption lines. The focus of the present review is on two new techniques, the Velocity-Channel Analysis (VCA), which makes use of the channel maps, and the Velocity Coordinate Spectrum (VCS), which utilizes the fluctuations measured along the velocity axis of the Position-Position Velocity (PPV) data cubes. Both techniques have solid theoretical foundations based on analytical calculations as well as on numerical testings. Among the two the VCS, which has been developed quite recently, has two unique features. First of all, it is applicable to turbulent volumes that are not spatially resolved. Second, it can be used with absorption lines that do not provide good spatial sampling of different lags over the image of the turbulent object. In fact, recent studies show that measurements of absorption line along less than 10 absorption directions are sufficient for a reliable recovering of the underlying spectrum of the turbulence. Moreover, both weak absorption lines and absorption lines in saturated regime can be used, which extends the applicability of the technique. Our comparison of the VCA and the VCS with a more traditional technique of Velocity Centroids shows that the former two techniques recover reliably the spectra of supersonic turbulence, while the Velocity Centroids may be used only for studying subsonic turbulence. We discuss spectra of astrophysical turbulence obtained with the VCA and the VCS techniques.

preprint2007arXiv

Cascading of Fast-Mode Balanced and Imbalanced Turbulence

We study the cascading of fast MHD modes in magnetically dominated plasma by performing one-dimensional (1D) dynamical simulations. We find that the cascading becomes more efficient as an angle between wave vector and underlying magnetic field increases and fast mode becomes more compressive. We also study imbalanced turbulence, in which wave packets propagating in one-direction have more power than those in the opposite direction. Unlike imbalanced \Alfvenic turbulence, the imbalanced fast mode turbulence shows faster cascading as the degree of imbalance increases. We find that the spectral index of the velocity and magnetic field, which are carried by the fast-mode turbulence, quickly reaches stationary value of -2. Thus we conclude that the dissipation of fast mode, at least in 1D case, happens not due to weak or strong turbulent cascading, but mostly due to nonlinear steepening. The density fluctuations, which are carried by slow-mode perturbation in the larger scale and entropy-mode perturbation in the smaller scale, depend on initial driving spectrum and a ratio of specific heat.

preprint2006arXiv

Enhancement and suppression of heat transfer by MHD turbulence

We study of the effect of turbulence on heat transfer within magnetized plasmas for energy injection velocities both larger and smaller that the Alfven speed. We find that in the latter regime the heat transfer is partially suppressed, while in the former regime the effects of turbulence depend on the intensity of driving. In fact, the scale l at which the turbulent velocity is equal the Alfven velocity is a new important parameter. When the electron mean free path is larger than l, the stronger the the turbulence, the lower thermal conductivity by electrons is. The turbulent motions, however, induces their own advective heat transport, which, for the parameters of intracluster medium (ICM) provides effective heat diffusivity that exceeds the classical Spitzer value.

preprint2006arXiv

Intermittency of Magnetohydrodynamic Turbulence: Astrophysical Perspective

Intermittency is an essential property of astrophysical fluids, which demonstrate an extended inertial range. As intermittency violates self-similarity of motions, it gets impossible to naively extrapolate the properties of fluid obtained computationally with relatively low resolution to the actual astrophysical situations. In terms of Astrophysics, intermittency affects turbulent heating, momentum transfer, interaction with cosmic rays, radio waves and many more essential processes. Because of its significance, studies of intermittency call for coordinated efforts from both theorists and observers. In terms of theoretical understanding we are still just scratching a surface of a very rich subject. We have some theoretically well justified models that are poorly supported by experiments, we also have She-Leveque model, which could be vulnerable on theoretical grounds, but, nevertheless, is well supported by experimental and laboratory data. I briefly discuss a rather mysterious property of turbulence called ``extended self-similarity'' and the possibilities that it opens for the intermittency research. Then I analyze simulations of MHD intermittency performed by different groups and show that their results do not contradict to each other. Finally, I discuss the intermittency of density, intermittency of turbulence in the viscosity-dominated regime as well as the intermittency of polarization of Alfvenic modes. The latter provides an attractive solution to account for a slower cascading rate that is observed in some of the numerical experiments. I conclude by claiming that a substantial progress in the field may be achieved by studies of the turbulence intermittency via observations.

preprint2006arXiv

Turbulence Spectra from Doppler-shifted Spectral Lines

Turbulence is a key element of the dynamics of astrophysical fluids, including those of interstellar medium, clusters of galaxies and circumstellar regions. Turbulent motions induce Doppler shifts of observable emission and absorption lines. In the review we discuss new techniques that relate the spectra of underlying velocity turbulence and spectra of Doppler-shifted lines. In particular, the Velocity-Channel Analysis (VCA) makes use of the channel maps, while the Velocity Coordinate Spectrum (VCS) utilizes the fluctuations measured along the velocity axis of the Position-Position Velocity (PPV) data cubes. Both techniques have solid foundations based on analytical calculations as well as on numerical testings. Among the two the VCS, which has been developed quite recently, has two advantages. First of all, it is applicable to turbulent volumes that are not spatially resolved. Second, it can be used with absorption lines that do not provide good spatial sampling of different lags over the image of turbulent object. In fact, numerical testing shows that measurements of Doppler shifted absorption lines over a few directions is sufficient for a reliable recovering of the underlying spectrum of the turbulence. Our comparison of the VCA and the VCS with a more traditional technique of Velocity Centroids, shows that the former two techniques recover reliably the spectra of supersonic turbulence, while the Velocity Centroids may have advantages for studying subsonic turbulence. In parallel with theoretical and numerical work on the VCA and the VCS, the techniques have been applied to spectroscopic observations. We discuss results on astrophysical turbulence obtained with the VCA and the VCS.

preprint2005arXiv

Particle Acceleration by MHD Turbulence

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions in the picture of particle acceleration. We make use of the recently established scaling of slow and fast MHD modes in strong and weak MHD turbulence to provide a systematic study of particle acceleration in magnetic pressure (low-$β$) and gaseous pressure (high-$β$) dominated plasmas. We consider the acceleration by large scale compressions in both slow and fast particle diffusion limits. We compare the results with the acceleration rate that arises from resonance scattering and Transit-Time Damping (TTD). We establish that fast modes accelerate particles more efficiently than slow modes. We find that particle acceleration by pitch-angle scattering and TTD dominates acceleration by slow or fast modes when the spatial diffusion rate is small. When the rate of spatial diffusion of particles is high, we establish an enhancement of the efficiency of particle acceleration by slow and fast modes in weak turbulence. We show that highly supersonic turbulence is an efficient agent for particle acceleration. We find that even incompressible turbulence can accelerate particles on the scales comparable with the particle mean free path.

preprint2005arXiv

The role of Violent Magnetic Reconnection on the Production of the Large Scale Superluminal Ejections of the Microquasar GRS 1915+105

We here propose that the large scale superluminal ejections observed in the galactic microquasar GRS 1915+105 during radio flare events are produced by violent magnetic reconnection episodes in the corona just above the inner edge of the magnetized accretion disk that surrounds the central $\sim 10 M_{\odot}$ black hole. The process occurs when a large scale magnetic field is established by turbulent dynamo in the inner disk region with a ratio between the gas$+$radiation and the magnetic pressures $β\simeq 1$, implying a magnetic field intensity of $\sim 7 \times 10^8$ G. During this process, substantial angular momentum is removed from the disk by the wind generated by the vertical magnetic flux therefore increasing the disk mass accretion to a value near (but below) the critical one ($\dot M \sim 10^{19} $ g s$^{-1}$). Part of the magnetic energy released by reconnection heats the coronal gas ($T_c \lesssim 5 \times 10^8 $ K) that produces a steep, soft X-ray spectrum with luminosity $L_X \simeq 10^{39} $ erg s$^{- 1}$, in consistency with observations. The remaining magnetic energy released goes to accelerate the particles to relativistic velocities ($v \sim v_A \sim c$, where $v_A$ is the Alfvén speed) in the reconnection site through first-order Fermi processes. In this context, two possible mechanisms have been examined which produce power-law electron distributions $N(E) \propto E^{-α_E}$, with $α_E = 5/2, 2$, and corresponding synchrotron radio power-law spectra with spectral indices which are compatible with that observed during the flares ($ S_ν \propto ν^{-0.75, - 0.5}$).

preprint2004arXiv

Alignment of Dust in Molecular Clouds

Polarimetry is one of the most informative techniques of studying magnetic fields in molecular clouds. How reliable the interpretation of the polarization maps in terms of magnetic fields is the issue that the grain alignment theory addresses. We show that grain alignment involves several processes acting simultaneously, but on different time-scales. We explain that rotating dust grains get substantial magnetic moment that allows them precess fast about magnetic field lines. As the result, grains preserve their orientation to magnetic field when the magnetic field direction fluctuates. We point out to the importance of internal alignment, i.e. the process forces grain axes to be aligned in respect to the grain angular momentum. We show that subtle quantum effects, in particular relaxation related to nuclear magnetic moments of atoms composing the grain, brings to live complex grain motions, e.g. flips. These flips substantially alter the dynamics of grain and limit the applicability of earlier theories that did not account for them. We also briefly review basic physical processes involved in the alignment of grain angular momentum in respect to interstellar magnetic field. We claim that the bulk of existing observational data is consistent with the radiative torque alignment mechanism. In particular, we show that large grains that are known to exist in the cores of molecular clouds may be aligned by the attenuated external interstellar radiation field.

preprint2004arXiv

Magnetic turbulence in compressible fluids

We present high resolution numerical simulations of compressible magnetohydrodynamic (MHD) turbulence. We concentrate on studies of spectra and anisotropy of velocity and density. We describe a technique of separating different magnetohydrodynamic (MHD) modes (slow, fast and Alfven) and apply it to our simulations. For the mildly supersonic case, velocity and density show anisotropy. However, for the highly supersonic case, we observe steepening of the slow mode velocity spectrum and isotropization of density. Our present studies show that the spectrum of density gets substantially {\it flatter} than the Kolmogorov one.

preprint2004arXiv

Obtaining Statistics of Turbulent Velocity from Astrophysical Spectral Line Data

Turbulence is a crucial component of dynamics of astrophysical fluids dynamics, including those of ISM, clusters of galaxies and circumstellar regions. Doppler shifted spectral lines provide a unique source of information on turbulent velocities. We discuss Velocity-Channel Analysis (VCA) and its offspring Velocity Coordinate Spectrum (VCS) that are based on the analytical description of the spectral line statistics. Those techniques are well suited for studies of supersonic turbulence. We stress that a great advantage of VCS is that it does not necessary require good spatial resolution. Addressing the studies of mildly supersonic and subsonic turbulence we discuss the criterion that allows to determine whether a traditional tool for such a research, namely, Velocity Centroids are dominated by density or velocity. We briefly discuss the use of higher order correlations as the means to study intermittency of turbulence. We discuss observational data available and prospects of the field.

preprint2004arXiv

Spurious contribution to CR scattering calculations

The quasilinear theory for cosmic ray propagation is a well known and widely accepted theory. In this paper, we discuss the different contributions to the pitch-angle Fokker-Planck coefficient from large and small scales for slab geometry using the damping model of dynamical turbulence. These examinations will give us a hint on the limitation range where quasilinear approximation is a good approximation.

preprint2003arXiv

Angular Spectra of Polarized Galactic Foregrounds

It is believed that magnetic field lines are twisted and bend by turbulent motions in the Galaxy. Therefore, both Galactic synchrotron emission and thermal emission from dust reflects statistics of Galactic turbulence. Our simple model of Galactic turbulence, motivated by results of our simulations, predicts that Galactic disk and halo exhibit different angular power spectra. We show that observed angular spectra of synchrotron emission are compatible with our model. We also show that our model is compatible with the angular spectra of star-light polarization for the Galactic disk. Finally, we discuss how one can estimate polarized microwave emission from dust in the Galactic halo using star-light polarimetry.

preprint2003arXiv

Compressible Magnetohydrodynamic Turbulence: mode coupling, scaling relations, anisotropy, new regime and astrophysical implications

We present numerical simulations and explore scalings and anisotropy of compressible magnetohydrodynamic (MHD) turbulence. Our study covers both gas pressure dominated (high beta) and magnetically dominated (low beta) plasmas at different Mach numbers. In addition, we present results for superAlfvenic turbulence and discuss in what way it is similar to the subAlfvenic turbulence. We describe a technique of separating different magnetohydrodynamic (MHD) modes (slow, fast and Alfven) and apply it to our simulations. We show that, for both high and low beta cases, Alfven and slow modes reveal the Kolmogorov spectrum (index=-5/3) and scale-dependent Goldreich-Sridhar anisotropy, while fast modes exhibit spectrum with index=-3/2 and isotropy. We discuss the statistics of density fluctuations arising from MHD turbulence at different regimes. Our findings entail numerous astrophysical implications ranging from cosmic ray propagation to gamma ray bursts and star formation. In particular, we show that the rapid decay of turbulence reported by earlier researchers is not related to compressibility and mode coupling in MHD turbulence. In addition, we show that magnetic field enhancements and density enhancements are marginally correlated. Addressing the density structure of partially ionized interstellar gas on AU scales, we show that the new regime of MHD turbulence that we earlier reported for incompressible flows persists for compressible turbulence and therefore may provide an explanation for those mysterious structures.

preprint2003arXiv

Generation of Compressible Modes in MHD Turbulence

Astrophysical turbulence is magnetohydrodynamic (MHD) in its nature. We discuss fundamental properties of MHD turbulence. In particular, we discuss the generation of compressible MHD waves by Alfvenic turbulence and show that this process is inefficient. This allows us to study the evolution of different types of MHD perturbations separately. We describe how to separate MHD fluctuations into 3 distinct families - Alfven, slow, and fast modes. We find that the degree of suppression of slow and fast modes production by Alfvenic turbulence depends on the strength of the mean field. We show that Alfven modes in compressible regime exhibit scalings and anisotropy similar to those in incompressible regime. Slow modes passively mimic Alfven modes. However, fast modes exhibit isotropy and a scaling similar to that of acoustic turbulence both in high- and low-beta plasmas. We show that our findings entail important consequences for theories of star formation, cosmic ray propagation, dynamics of dust, and gamma ray bursts. We anticipate many more applications of the new insight to MHD turbulence and expect more revisions of the existing paradigms of astrophysical processes as the field matures.

preprint2003arXiv

Microwave Emission from Aligned Dust

Polarized microwave emission from dust is an important foreground that may contaminate polarized CMB studies unless carefully accounted for. We discuss potential difficulties associated with this foreground, namely, the existence of different grain populations with very different emission/polarization properties and variations of the polarization yield with grain temperature. In particular, we discuss observational evidence in favor of rotational emission from tiny PAH particles with dipole moments, i.e. ``spinning dust'', and also consider magneto-dipole emission from strongly magnetized grains. We argue that in terms of polarization, the magneto-dipole emission may dominate even if its contribution to total emissivity is subdominant. Addressing polarized emission at frequencies larger than approsimately 100 GHz, we discuss the complications arising from the existence of dust components with different temperatures and possibly different alignment properties.

preprint2002arXiv

Compressible Sub-Alfvenic MHD turbulence in Low-beta Plasmas

We present a model for compressible sub-Alfvenic isothermal magnetohydrodynamic (MHD) turbulence in low-beta plasmas and numerically test it. We separate MHD fluctuations into 3 distinct families - Alfven, slow, and fast modes. We find that, production of slow and fast modes by Alfvenic turbulence is suppressed. As a result, Alfven modes in compressible regime exhibit scalings and anisotropy similar to those in incompressible regime. Slow modes passively mimic Alfven modes. However, fast modes show isotropy and a scaling similar to acoustic turbulence.

preprint2002arXiv

MHD Turbulence as a Foreground for CMB Studies

Measurements of intensity and polarization of diffuse Galactic synchrotron emission as well as starlight polarization reveal power law spectra of fluctuations. We show that these fluctuations can arise from magnetohydrodynamic (MHD) turbulence in the Galactic disk and halo. To do so we take into account the converging geometry of lines of sight for the observations when the observer is within the turbulent volume. Assuming that the intensity of turbulence changes along the line of sight, we get a reasonable fit to the observed synchrotron data. As for the spectra of polarized starlight we get a good fit to the observations taking into account the fact that the observational sample is biased toward nearby stars.

preprint2002arXiv

MHD Turbulence: Scaling Laws and Astrophysical Implications

Turbulence is the most common state of astrophysical flows. In typical astrophysical fluids, turbulence is accompanied by strong magnetic fields, which has a large impact on the dynamics of the turbulent cascade. Recently, there has been a significant breakthrough on the theory of magnetohydrodynamic (MHD) turbulence. For the first time we have a scaling model that is supported by both observations and numerical simulations. We review recent progress in studies of both incompressible and compressible turbulence. We compare Iroshnikov-Kraichnan and Goldreich-Sridhar models, and discuss scalings of Alfvén, slow, and fast waves. We also discuss the completely new regime of MHD turbulence that happens below the scale at which hydrodynamic turbulent motions are damped by viscosity. In the case of the partially ionized diffuse interstellar gas the viscosity is due to neutrals and truncates the turbulent cascade at $\sim$parsec scales. We show that below this scale magnetic fluctuations with a shallow spectrum persist and discuss the possibility of a resumption of the MHD cascade after ions and neutrals decouple. We discuss the implications of this new insight into MHD turbulence for cosmic ray transport, grain dynamics, etc., and how to test theoretical predictions against observations.

preprint2001arXiv

Inelastic Dissipation in Wobbling Asteroids and Comets

Asteroids and comets dissipate energy when they rotate about the axis different from the axis of the maximal moment of inertia. We show that the most efficient internal relaxation happens at the double frequency of body's tumbling. Therefore the earlier estimates that ignore double frequency input underestimate the internal relaxation in asteroids and comets. We show that the Earth seismological data may poorly represent acoustic properties of asteroids and comet as internal relaxation increases in the presence of moisture. At the same time owing to non-linearlity of inelastic relaxation small angle nutations can persist for very long time spans, but our ability to detect such precessions is limited by the resolution of the radar-generated images. Wobbling may provide valuable information on the composition and structure of asteroids and on their recent history of external impacts.

preprint2001arXiv

Velocity and density spectra of the Small Magellanic Cloud

This paper reports results on the statistical analysis of HI turbulence in the Small Magellanic Cloud (SMC). We use 21 cm channel maps, obtained with the Australia Telescope Compact Array and the Parkes telescope, and analyze the spectrum of observed intensity fluctuations as a function of the velocity slice thickness. We confirm predictions by Lazarian & Pogosyan (2000) on the change of the power law index and establish the spectra of 3-D density and velocity. The obtained spectral indices, -3.3 and -3.4, are slightly more shallow than the predictions for the Kolmogorov spectrum. This contrasts to the predictions for the shock-type spectra that are steeper than the Kolmogorov one. The nature of the energy injection in the SMC is unclear as no distinct energy injection scales are observed up to the entire scale of the SMC.

preprint2000arXiv

Fast Turbulent Reconnection

Reconnection is the process by which magnetic fields in a conducting fluid change their topology. This process is essential for understanding a wide variety of astrophysical processes, including stellar and galactic dynamos and astrophysical turbulence. To account for solar flares, solar cycles and the structure of the galactic magnetic field reconnection must be fast, propagating with a speed close to the Alfven speed. We show that the presence of a random magnetic field component substantially enhances the reconnection rate and enables fast reconnection, i.e. reconnection that does not depend on fluid resistivity. The enhancement of the reconnection rate is achieved via a combination of two effects. First of all, only small segments of magnetic field lines are subject to direct Ohmic annihilation. Thus the fraction of magnetic energy that goes directly into fluid heating goes to zero as fluid resistivity vanishes. However, the most important enhancement comes from the fact that unlike the laminar fluid case where reconnection is constrained to proceed line by line, the presence of turbulence enables many magnetic field lines to enter the reconnection zone simultaneously. A significant fraction of magnetic energy goes into MHD turbulence and this enhances reconnection rates through an increase in the field stochasticity. In this way magnetic reconnection becomes fast when field stochasticity is accounted for. As a consequence solar and galactic dynamos are also fast, i.e. do not depend on fluid resistivity.

preprint2000arXiv

Statistics of Turbulence from Spectral-Line Data Cubes

Emission in spectral lines can provide unique information on interstellar turbulence. Doppler shifts due to supersonic motions contain information on turbulent velocity field which is otherwise difficult to measure. However, the problem of separation of velocity and density fluctuations is far from being trivial. Using atomic hydrogen (HI) as a test case, I review techniques applicable to emission line studies with the emphasis on those that can provide information on the underlying power spectra of velocity and density. I show that recently developed mathematical machinery is promising for the purpose. Its application to HI shows that in cold neutral hydrogen the velocity fluctuations dominate the small scale structures observed in spectral-line data cubes and this result is very important for the interpretation of observational data, including the identification of clouds. Velocity fluctuations are shown to dominate the formation of small scale structures that can be erroneously identified as diffuse clouds. One may argue that the HI data is consistent with the Goldreich-Shridhar picture of magnetohydrodynamic turbulence, but the cascade from the scales of several kpc that this interpretation involves does not fit well in the current paradigm of energy injection. The issue whether magnetic field does make the turbulence anisotropic is still open, but if this is the case, I show that studies of emission lines can provide a reliable way of determining magnetic field direction. I discuss various techniques for studying interstellar turbulence using emission lines, e.g. spectral correlation functions, genus statistics and principal component analysis.

preprint1999arXiv

Nuclear Spin Relaxation within Interstellar Grains

We identify a new mechanism of internal dissipation of rotational kinetic energy in spinning dust grains, arising from the reorientation of nuclear angular momentum, e.g., spins of protons. Grain rotation induces magnetization of the nuclear spin system, with net alignment of nuclear spins parallel to the grain angular velocity. When the grain does not rotate around a principal axis, the nuclear magnetization vector precesses in grain body coordinates, resulting in dissipation of energy. The analogous process involving electron spins was discovered by Purcell and termed ``Barnett relaxation''. We revisit the physics of the ``Barnett relaxation'' process and correct the estimate for the Barnett relaxation rate. We show that nuclear relaxation can be orders of magnitude more important than Barnett relaxation. This finding implies that the processes of ``thermal flipping'' and ``thermal trapping'' are important for a broad range of grain sizes.

preprint1999arXiv

Reconnection in a Weakly Stochastic Field

We examine the effect of weak, small scale magnetic field structure on the rate of reconnection in a strongly magnetized plasma. This affects the rate of reconnection by reducing the transverse scale for reconnection flows, and by allowing many independent flux reconnection events to occur simultaneously. Allowing only for the first effect and using Goldreich and Sridhar's model of strong turbulence in a magnetized plasma with negligible intermittency, we find that the lower limit for the reconnection speed is the Alfven speed times the Lundquist number to the power (-3/16). The upper limit on the reconnection speed is typically a large fraction of Alfven speed. We argue that generic reconnection in turbulent plasmas will normally occur at close to this upper limit. The fraction of magnetic energy that goes directly into electron heating scales as Lundquist number to the power (-2/5) and the thickness of the current sheet scales as the Lundquist number to the power (-3/5). A significant fraction of the magnetic energy goes into high frequency Alfven waves. We claim that the qualitative sense of these conclusions, that reconnection is fast even though current sheets are narrow, is almost independent of the local physics of reconnection and the nature of the turbulent cascade. As the consequence of this the Galactic and Solar dynamos are generically fast, i.e. do not depend on the plasma resistivity.

preprint1998arXiv

Davis-Greenstein alignment of oblate spheroidal grains

We present extensive calculations on the efficiency of grain alignment by the Davis-Greenstein mechanism. We model the grains as oblate spheroids with arbitrary axis ratios. Our description of the grain dynamics includes (i) magnetic dissipation and the inverse process driven by thermal fluctuations in the grain magnetization; (ii) gas-grain collisions and thermal evaporation of molecules from the grain surface; (iii) the transformation of rotational energy into heat by the Barnett effect and the inverse process driven by thermal fluctuations; and (iv) rapid Larmor precession of the grain angular momentum about the interstellar magnetic field. For ordinary paramagnetic grains, we calculate the Rayleigh reduction factor, R, for more than 1000 combinations of the 3 dimensionless parameters which characterise the alignment. For superparamagnetic grains, we calculate R from an exact analytic solution for the relevant distribution function. We calibrate the accuracy of a recently-proposed perturbative approximation and show that it yields R values with a mean error of about 17 percent.

preprint1998arXiv

Electric Dipole Radiation from Spinning Dust Grains

We discuss the rotational excitation of small interstellar grains and the resulting electric dipole radiation from spinning dust. Attention is given to excitation and damping of rotation by: collisions with neutrals; collisions with ions; plasma drag; emission of infrared radiation; emission of microwave radiation; photoelectric emission; and formation of H_2 on the grain surface. We introduce dimensionless functions F and G which allow direct comparison of the contributions of different mechanisms to rotational drag and excitation. Emissivities are estimated for dust in different phases of the interstellar medium, including diffuse HI, warm HI, low-density photoionized gas, and cold molecular gas. Spinning dust grains can explain much, and perhaps all, of the 14-50 GHz background component recently observed in CBR studies. It should be possible to detect rotational emission from small grains by ground-based observations of molecular clouds.

preprint1997arXiv

Cosmic rays and grain alignment

The recent detection of interstellar polarization in the solid CO feature near 4.67 micron shows that CO-mantled grains can be aligned in cold molecular clouds. These observations conflict with a theory of grain alignment which attributes the polarization in molecular clouds to the effects of cosmic rays: according to this theory, oblate spheroidal grains with H_2O and CO_2-dominated ice mantles are spun up to suprathermal energies by molecular evaporation from cosmic ray impact sites but spin up does not occur for CO-mantled grains. Motivated by this conflict, we reexamine the effects of cosmic rays on the alignment of icy grains. We show that the systematic torques produced by cosmic rays are insufficient to cause suprathermal spin. In principle, the random torques due to cosmic rays can enhance the efficiency of Davis-Greenstein alignment by raising the grain rotational temperature. However, a significant enhancement would require cosmic ray fluxes 6--7 orders of magnitude larger than the flux in a typical cold cloud.

preprint1997arXiv

Gold alignment & internal dissipation

The measures of mechanical alignment were obtained for both prolate and oblate grains when their temperature is comparable with grain kinetic energy devided by k, the Boltzmann constant. For such grains, the alignment of angular momentum, J, with the axis of maximal inertia, a, is only partial. This substantially alters the alignment as compared with the results in Lazarian (1995) and Roberge, Hanany & Messinger (1996) obtained on the assumption of perfect alignment. We also describe the Gold alignment when the Barnett dissipation is suppressed and derive an analytical expression which relates the measure of alignment with parameters of grain nonsphericity and the direction of the gas - grain drift. This solution provides the lower limit for the alignment measure, while the upper limit is given by the analytics derived in Lazarian (1994). Using results of a recent study of incomplete internal relaxation in Lazarian & Roberge (1996), we find measures of alignment for the whole range of ratios of grain rotational energy to k over T_s, where T_s is the grain temperature. To describe alignment for mildly supersonic drifts, we suggest an analytical approach which provides good correspondence with the results of direct numerical simulations in Roberge, Hanany & Messinger (1995). We also extend our approach to account for the simultaneous action of the Gold and Davis-Greenstein mechanisms.

preprint1997arXiv

Paramagnetic alignment of thermally rotating dust

Paramagnetic alignment of thermally rotating oblate dust grains is studied analytically for finite ratios of grain to gas temperatures. For such ratios, the alignment of angular momentum J in respect to the grain axis of maximal inertia is only partial. We treat the alignment of J using perturbative methods and disentangle the problem of J alignment in grain body axes from that of J alignment in respect to magnetic field. This enables us to find the alignment of grain axes to magnetic field and thus relate our theory to polarimetric observations.Our present results are applicable to the alignment of both paramagnetic and superparamagnetic grains.

preprint1997arXiv

Reconnection in the ISM

We discuss the role of ambipolar diffusion for simple reconnection in a partially ionized gas, following the reconnection geometry of Parker and Sweet. When the recombination time is short the mobility and reconnection of the magnetic field is substantially enhanced as matter escapes from the reconnection region via ambipolar diffusion. Our analysis shows that in the interstellar medium it is the recombination rate that usually limits the rate of reconnection. We show that heating effects can reduce this rate by increasing the recombination time and raising the local ion pressure. In the colder parts of the ISM, when temperatures are or the order of 100K or less, we obtain a significant enhancement over the usual Sweet-Parker rate, but only in dense molecular clouds will the reconnection velocity exceed 0.001 times the Alfvén speed. The ratio of the ion orbital radius to the reconnection layer thickness is typically a few percent, except in dense molecular clouds where it can approach unity. We briefly discuss prospects for obtaining much faster reconnection speeds in astrophysical plasmas.

preprint1995arXiv

Alignment of suprathermally rotating grains

It is shown, that mechanical alignment can be efficient for suprathermally rotating grains, provided that they drift with supersonic velocities. Such a drift should be widely spread due to both Alfvénic waves and ambipolar diffusion. Moreover, if suprathermal rotation is caused by grain interaction with a radiative flux, it is shown, that mechanical alignment may be present even in the absence of supersonic drift. This means that the range of applicability of mechanical alignment is wider that it is generally accepted and it can rival the paramagnetic one. We also study the latter mechanism and reexamine the interplay between poisoning of active sites and desorption of molecules blocking the access to the active sites of H$_{2}$ formation in order to explain the observed poor alignment of small grains and good alignment of large grains. To have a more comprehensive picture of alignment we briefly discuss the alignment by radiation fluxes and caused by grain magnetic moments.