Magnetic fields in the intracluster medium (ICM) affect the structure and the evolution of galaxy clusters. However, their properties are largely unknown, and measuring magnetic fields in galaxy clusters is challenging, especially on large-scales outside of individual radio sources. In this work, we probe the plane-of-the-sky orientation of magnetic fields in clusters using the intensity gradients. The technique is a branch of the Gradient Technique (GT) that employs emission intensity maps from turbulent gas. We utilize the Chandra X-ray images of the Perseus, M 87, Coma, and A2597 galaxy clusters, and the VLA radio observations of the synchrotron emission from Perseus. We find that the fields predominantly follow the spiral arms in Perseus, which is in agreement with numerical simulations. The GT-predicted magnetic field shows signatures of magnetic draping around rising bubbles driven by supermassive black hole (SMBH) feedback in the centers of cool-core clusters, as well as draping around substructures merging with the Coma cluster. We calculate the mean-field orientation with respect to the radial direction in these clusters. The orientation is consistent with the prediction of heat-flux-driven buoyancy instability (HBI) in cluster cool cores, and possibly the magnetothermal instability (MTI) in cluster outskirts. There is a broad agreement between the magnetic field of Perseus predicted using the X-ray and radio data. Further numerical studies and better future observations with higher resolution and the larger effective area will help reduce the uncertainties of this method.

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.

We study the efficiency of grain alignment by radiative torques (RATs) for an ensemble of irregular grains. The grains are modeled as ensembles of oblate and prolate spheroids, deformed as Gaussian random ellipsoids, and their scattering interactions are solved using numerically exact methods. We define the fraction of the grains that both rotate fast and demonstrate perfect alignment with grain long axes perpendicular to the magnetic field. We demonstrate that for typical interstellar conditions the degree of alignment arising from the RAT mechanism is significantly larger than that arising from the Davis-Greenstein (DG) mechanism based on paramagnetic relaxation. We quantify a factor related to the efficacy of alignment and show that it is related to a $q_\mathrm{max}$-factor of analytical model (AMO) of the RAT theory. Our results indicate that the RAT alignment can potentially be sufficiently strong and to explain observations even if grains do not have magnetic inclusions.

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

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\'{e}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 $\eta_{RD}$ using the Test-Field method. Our results are consistent with the RD theory ($\eta_{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 ($\eta_{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.

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.

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

We provide a procedure for identification of dominating compressible and Alfvenic MHD modes or isotropic turbulence in synchrotron emission polarization maps of Galactic objects. The results for the region of North Galactic Pole, Orion molecular cloud complex and the star-forming complex Cygnus X are presented.

The first-year results from the Parker Solar Probe (PSP) reveal a gradual decrease of F-coronal dust from distances of $D=0.166-0.336$ AU (or the inner elongations of $\sim 9.22- 18.69~R_{\odot}$) to the Sun, which cannot be explained by the dust sublimation scenario that predicts a dust-free-zone $\lesssim 4-5R_{\odot}$. In this paper, we attempt to explain this puzzle using our newly introduced mechanism of dust destruction so-called Radiative Torque Disruption (RATD) mechanism. We demonstrate that RATD rapidly breaks large grains into nanoparticles so that they can be efficiently destroyed by nonthermal sputtering induced by bombardment of energetic protons from slow solar winds. This joint effect extends the dust-free-zone established by thermal sublimation to $R_{dfz}\sim 8R_{\odot}$, which is called a new dust-free-zone. Beyond this new dust-free-zone, we find that the dust mass decreases gradually from $R\sim 42R_{\odot}$ toward the Sun due to partial removal of nanodust by nonthermal sputtering. This feature can successfully reproduce the gradual decrease of the F-corona between $19-9R_{\odot}$ observed by the PSP. Finally, the RATD mechanism can efficiently produce nanoparticles usually observed in the inner solar system.

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

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

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

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

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

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

Starburst galaxies are efficient $\gamma$-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 $\gamma$-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 $\gamma$-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 $\gamma$-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.

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.

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

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

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

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

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

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

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.

Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) -high-energy charged particles accelerated in supernovae and young stars - can impact the efficiency of wind driving. The streaming instability limits the speed at which they can escape. However, in the presence of turbulence, the streaming instability is subject to suppression that depends on the magnetization of turbulence given by its Alfv\'en Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfv\'enic streaming of CRs enhances galactic winds, in the present paper we take into account a realistic model of streaming suppression. We perform three-dimensional magnetohydrodynamic simulations of a section of a galactic disk and find that turbulent damping dependent on local magnetization of turbulent interstellar medium (ISM) leads to more spatially extended gas and CR distributions compared to the earlier streaming calculations, and that scale-heights of these distributions increase for stronger turbulence. Our results indicate that the star formation rate increases with the level of turbulence in the ISM. We also find that the instantaneous wind mass loading is sensitive to local streaming physics with the mass loading dropping significantly as the strength of turbulence increases.

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

Alignment of dust grains in astrophysical environments results in the polarization of starlight as well as the polarization of radiation emitted by dust. We demonstrate the advances in grain alignment theory allow the use of linear and circular polarization to probe not only the magnetic field, but also dust composition, the dust environment, etc. We revisit the process of grain alignment by Radiative Torques (RATs) and focus on constraining magnetic susceptibility of grains via observations. We discuss the possibility of observational testing of the magnetic properties of grains as the alignment changes from being in respect to the magnetic field to being in respect to the radiation direction. This opens both a possibility of constraining the uncertain parameters of the RATs theory and provides a new way of measuring magnetic fields in interstellar medium and circumstellar regions. We provide a detailed discussion of the precession induced both by the magnetic field and the anisotropic radiation and revisit a number of key processes related to magnetic response of the grains. We consider various effects that increase the rate of magnetic relaxation both in silicate and carbonaceous grains. In particular, we find a new relaxation process related to the change of the amplitude of internal magnetization within a wobbling triaxial grain and identify a range of grain sizes in which this effect can dominate the internal alignment of angular momentum within grain axes. We show that these relaxation processes significantly change the dynamics of grains in the presence of RATs. We apply our analysis for observed grain alignment in special environments to put constraints on the enhanced magnetic properties of dust grains in the cloud near supernovae, in cometary coma, and protoplanetary disks.

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

We explore how the velocity gradient technique (VGT) can be applied to absorbing media in the case of$^{13}$CO 2-1 emission. The VGT is a new way to trace magnetic fields in the plane of the sky using only spectroscopic observations. We apply the VGT to magnetohydrodynamic turbulence simulations that have been post-processed to include $^{13}$CO 2-1 emission and we calculate the velocity centroid gradients. We find that the velocity centroid gradients trace the projected magnetic field in media with different $^{13}$CO abundances, densities and optical depths. Our study opens up the possibility of using velocity centroid gradients to trace magnetic fields in molecular clouds using 13CO emission.

We apply Velocity Channel Gradients (VChGs) and Reduced Velocity Centroids Gradients (RVCGs) to H I data (the LAB survey), obtaining the plane-of-sky component (POS) of the magnetic field as a function of the relative velocity. Assuming a galactic rotation curve, we transformed the relative velocities to distances, constructing the first map of the POS magnetic field at every point of the celestial sphere as a function of the distance. To test the accuracy of our 3D distribution, we used a set of stars with known distances from the stellar polarization catalog. We compared the polarization directions that we obtain with the VChGs and RVCGs against the starlight polarization directions. We find a good correspondence between the derived magnetic field and measured polarization directions, which testifies to the accuracy of this new way of probing the 3D galactic POS magnetic field structure.

We extend the use of the Tsallis statistic to measure the differences in gas dynamics relative to the mean magnetic field present from natural eddy-type motions existing in magnetohydrodynamical (MHD) turbulence. The variation in gas dynamics was estimated using the Tsallis parameters on the incremental probability distribution function of the observables (intensity and velocity centroid) obtained from compressible MHD simulations. We find that the Tsallis statistic is susceptible to the anisotropy produce by the magnetic field, even when anisotropy is percent the Tsallis statistic can be use to determine MHD parameters such as the Sonic Mach number. We quantize the goodness of the Tsallis parameters using the coefficient of determination to measure the differences in the gas dynamics. These parameters also determine the level of magnetization and compressibility of the medium. To further simulate realistic spectroscopic observational data we introduced smoothing, noise, and cloud boundaries to the MHD simulations.

This study is motivated by recent observations on ubiquitous interstellar density filaments and guided by modern theories of compressible magnetohydrodynamic (MHD) turbulence. The interstellar turbulence shapes the observed density structures. As the fundamental dynamics of compressible MHD turbulence, perpendicular turbulent mixing of density fluctuations entails elongated density structures aligned with the local magnetic field, accounting for low-density parallel filaments seen in diffuse atomic and molecular gas. The elongation of low-density parallel filaments depends on the turbulence anisotropy. When taking into account the partial ionization, we find that the minimum width of parallel filaments in the cold neutral medium and molecular clouds is determined by the neutral-ion decoupling scale perpendicular to magnetic field. In highly supersonic MHD turbulence in molecular clouds, both low-density parallel filaments due to anisotropic turbulent mixing and high-density filaments due to shock compression exist.

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

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

Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

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

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

The small-scale turbulent dynamo is an important process contributing to the cosmic magnetization. In partially ionized astrophysical plasmas, the dynamo growth of magnetic energy strongly depends on the coupling state between ions and neutrals and the ion-neutral collisional damping effect. A new damping stage of turbulent dynamo in a weakly ionized medium was theoretically predicted by Xu \& Lazarian (2016). By carrying out a 3D two-fluid dynamo simulation, here we for the first time numerically confirmed the physical conditions and the linear-in-time growth of magnetic field strength of the damping stage of dynamo. The dynamo-amplified magnetic field has a characteristic length as the damping scale, which increases with time and can reach the injection scale of turbulence after around eight largest eddy-turnover times given sufficiently low ionization fraction and weak initial magnetic field. Due to the weak coupling between ions and neutrals, most turbulent energy carried by neutrals cannot be converted to the magnetic energy, resulting in a relatively weak magnetic field at the end of dynamo. This result has important implications for the growth of magnetic fields in the partially ionized interstellar medium and shock acceleration of Galactic cosmic rays.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

The radiative torque (RAT) mechanism is the most promising way of explaining observed polarization arising from aligned grains. We explore the efficiency of the grain alignment by an anisotropic radiation flow for an extensive ensemble of grain shapes. We calculate the distribution of the ratios of the amplitudes of the two major components of the RATs, that is an essential parameter that enters the theory of RAT alignment in Lazarian & Hoang (2007, LH07). While this distribution is different for different classes of grain shapes that we considered, the most probable values of the parameter are centered in the range of $q^{max}\sim 0.5-1.5$. The functional form from RATs calculated is in good agreement with the analytical model (AMO). We find that the RAT efficiency scales as $(\lambda/a)^{-3}$ for $\lambda\gg a$ as previously found in LH07. This increases the power of predictions obtained with the RAT theory. We also confirm that superparamagnetic inclusions are necessary in achieving high degrees of alignment, and constrain the parameter space describing the requirements for achieving these alignment degrees.

Astrophysical fluids are turbulent, magnetized and frequently partially ionized. As an example of astrophysical turbulence, the interstellar turbulence extends over a remarkably large range of spatial scales and participates in key astrophysical processes happening in different ranges of scales. A significant progress has been achieved in the understanding of the magnetohydrodynamic (MHD) turbulence since the turn of the century, and this enables us to better describe turbulence in magnetized and partially ionized plasmas. In fact, the modern revolutionized picture of the MHD turbulence physics facilitates the development of various theoretical domains, including the damping process for dissipating MHD turbulence and the dynamo process for generating MHD turbulence with many important astrophysical implications. In this paper, we review some important findings from our recent theoretical works to demonstrate the interconnection between the properties of MHD turbulence and those of turbulent dynamo in a partially ionized gas. We also briefly exemplify some new tentative studies on how the revised basic processes influence the associated outstanding astrophysical problems in, such as, magnetic reconnection, cosmic ray scattering, magnetic field amplification in both the early and the present-day universe.

As a fundamental astrophysical process, the scattering of particles by turbulent magnetic fields has its physical foundation laid by the magnetohydrodynamic (MHD) turbulence theory. In the framework of the modern theory of MHD turbulence, we derive a generalized broadened resonance function by taking into account both the magnetic fluctuations and nonlinear decorrelation of turbulent magnetic fields arising in MHD turbulence, and we specify the energy range of particles for the dominance of different broadening mechanisms. The broadened resonance allows for scattering of particles beyond the energy threshold of the linear resonance. By analytically determining the pitch-angle diffusion coefficients for transit time damping (TTD) with slow and fast modes, we demonstrate that the turbulence anisotropy of slow modes suppresses their scattering efficiency. Furthermore, we quantify the dependence of the relative importance between slow and fast modes in TTD scattering on (i) particle energy, (ii) plasma $\beta$ (the ratio of gas pressure to magnetic pressure), and (iii) damping of MHD turbulence, and we also provide the parameter space for the dominance of slow modes. To exemplify its applications, we find that among typical partially ionized interstellar phases, in the warm neutral medium slow and fast modes have comparable efficiencies in TTD scattering of cosmic rays. For low-energy particles, e.g., sub-Alfv\'{e}nic charged grains, we show that slow modes always dominate TTD scattering.

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

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

The spectral slope of the magnetic energy in supernova remnants (SNRs) can be obtained by analysis of spatial two-point correlation functions of synchrotron intensities. This method has been originally developed for the analysis of magnetic field structure in diffuse interstellar medium and applied when the geometry of the emission region is simple and known. In this paper, applying this correlation analysis to Tycho's SNR, for which the synchrotron emission region is known to be a spherical shell, we find that the magnetic energy spectrum shows the Kolmogorov-like scaling. Our results can be explained by turbulence developed downstream of the shock front via Richtmyer-Meshkov instability or an amplification of upstream magnetic field induced by cosmic rays. They could be discriminated by future observations with a sub arcsecond resolution such as Square Kilometer Array.

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

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

We investigate how to recover the spectral properties of underlying magnetohydrodynamic (MHD) turbulence using fluctuation statistics of synchrotron polarization radiation, based on the synthetic observations. Taking spatially coincident, separated, and compounded synchrotron emission and Faraday rotation regions into account, we extract the power spectrum of synchrotron polarization intensities integrated along the line of sight. Our results demonstrate that in the short wavelength range, the power spectra reflect fluctuation statistics of the perpendicular component of turbulent magnetic fields, and the spectra at long wavelengths reveal the fluctuation of the Faraday rotation density, which is a product of the parallel component of magnetic field and thermal electron density. We find that our numerical results (in the case of spatially coincident regions) are in agreement with the analytical prediction in Lazarian \& Pogosyan (2016), and this theoretical prediction is applicable to more complicated settings, i.e., the spatially separated and compounded regions. We simulate telescopic observations that incorporate the effects of telescope angular resolution and noise, and find that statistics of underlying MHD turbulence can be recovered successfully. We expect that the technique can be applied to a variety of astrophysical environments, with existing synchrotron data cubes and a large number of forthcoming data sets from such as the LOw-Frequency Array for Radio astronomy (LOFAR), the Square Kilometer Array (SKA) and the Five-hundred-meter Aperture Spherical radio Telescope (FAST).

We study the effect of anisotropic radiation illumination on the alignment of polycyclic aromatic hydrocarbons (PAHs) and report that cross-sectional mechanism of alignment earlier considered in terms of gas-grain interactions can also be efficient for the photon-grain interaction. We demonstrate this by first calculating the angle-dependence rotational damping and excitation coefficients by photon absorption followed by infrared emission. We then calculate the degree of PAH alignment for the different environments and physical parameters, including the illumination direction, ionization fraction, and magnetic field strength. For the reflection nebula (RN) conditions with unidirectional radiation field, we find that the degree of alignment tends to increase with increasing the angle $\psi$ between the illumination direction and the magnetic field, as a result of the decrease of the cross-section of photon absorption with $\psi$. We calculate the polarization of spinning PAH emission using the obtained degree of alignment for the different physical parameters, assuming constant grain temperatures. We find that the polarization of spinning PAH emission from RN can be large, between $5-20~\%$ at frequencies $\nu > 20$ GHz, whereas the polarization is less than $3~\%$ for photodissociation regions (PDRs). In realistic conditions, the polarization is expected to be lower due to grain temperature fluctuations and magnetic field geometry. The polarization for the diffuse cold neutral medium (CNM) is rather low, below $1~\%$ at $\nu>20$ GHz, consistent with observations by WMAP and Planck. Our results demonstrate that the RNe are the favored environment to observe the polarization of spinning dust emission as well as polarized mid-IR emission from PAHs.

The discovery of the first interstellar asteroid, 1I/2017 U1 (`Oumuamua), has opened a new era for research on interstellar objects. In this paper, we study the rotational dynamics of interstellar asteroids (ISAs) of irregular shapes moving through the interstellar gas. We find that regular mechanical torques resulting from the bombardment of gas flow on the irregular body could be important for the dynamics and destruction of ISAs. Mechanical torques can spin up the ISA, resulting in the breakup of the original ISA into small binary asteroids when the rotation rate exceeds the critical frequency. We find that the breakup timescale is short for ISAs of highly irregular shapes and low tensile strength. We apply our results to the first observed ISA, `Oumuamua, and suggest that its extreme elongated shape may originate from a reassembly of the binary fragments due to gravity along its journey in the interstellar medium. The tumbling of `Oumuamua could have been induced by rotational disruption due to mechanical torques. Finally, we discuss the survival possibility of high-velocity asteroids presumably formed by tidal disruption of planetary systems by the black hole at the Galactic center.

Recent Planck observations have revealed some of the important statistical properties of synchrotron and dust polarisation, namely, the $B$ to $E$ mode power and temperature-$E$ (TE) mode cross-correlation. In this paper, we extend our analysis in Kandel et al. (2017) that studied $B$ to $E$ mode power ratio for polarised dust emission to include TE cross-correlation and develop an analogous formalism for synchrotron signal, all using a realistic model of magnetohydrodynamical (MHD) turbulence. Our results suggest that the Planck results for both synchrotron and dust polarisation can be understood if the turbulence in the Galaxy is sufficiently sub-Alfv\'enic. We also show how $B$ to $E$ ratio as well as the TE cross-correlation can be used to study media magnetisation, compressibility, and level of density-magnetic field correlation.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016.

This study is motivated by recent observations on the ubiquitous interstellar density filaments and guided by the modern theory of magnetohydrodynamic (MHD) turbulence. The interstellar turbulence shapes the observed density structure. The perpendicular turbulent mixing, as the fundamental dynamics of MHD turbulence, naturally entails parallel filaments presented in both the diffuse medium and molecular clouds (MCs). The minimum width is determined by the perpendicular neutral-ion decoupling scale in partially ionized media. Differently, dense perpendicular filaments arise in highly supersonic turbulence in MCs as a result of shock compression. Their width specifically depends on the turbulence properties. We demonstrate that different alignments of filaments with respect to the magnetic field originate from the varying turbulence properties in the multi-phase interstellar medium.

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

We investigate roles of magnetic activity in the Galactic bulge region in driving large-scale outflows of size $\sim 10$ kpc. Magnetic buoyancy and breakups of channel flows formed by magnetorotational instability excite Poynting flux by the magnetic tension force. A three-dimensional global numerical simulation shows that the average luminosity of such \Alfvenic Poynting flux is $10^{40} - 10^{41}$ erg s$^{-1}$. We examine the energy and momentum transfer from the Poynting flux to the gas by solving time-dependent hydrodynamical simulations with explicitly taking into account low-frequency \Alfvenic waves of period of 0.5 Myr in a one-dimensional vertical magnetic flux tube. The \Alfvenic waves propagate upward into the Galactic halo, and they are damped through the propagation along meandering magnetic field lines. If the turbulence is nearly trans-Alfv\'{e}nic, the wave damping is significant, which leads to the formation of an upward propagating shock wave. At the shock front, the temperature $\gtrsim 5\times 10^6$ K, the density $\approx 6\times 10^{-4}$ cm$^{-3}$, and the outflow velocity $\approx 400-500$ km s$^{-1}$ at a height $\approx 10$ kpc, which reasonably explain the basic physical properties of the thermal component of the Fermi bubbles.

We consider a simple model for gamma-ray bursts induced by magnetic reconnection in turbulent media. The magnetic field in a jet is subject to kink instabilities, which distort the regular structure of the spiral magnetic field, drive turbulence, and trigger reconnection. The resulting reconnection takes place in a high Reynolds number medium, where turbulence is further enhanced and in turn accelerates the reconnection process. This boot-strap reconnection gives rise to bursts of reconnection events, through which the free energy of magnetic field is transformed into a gamma-ray burst. The efficiency of magnetic reconnection and magnetic energy dissipation is not constrained by microphysical properties of plasmas.

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

We study the alignment of irregular dust grains by mechanical torques due to the drift of grains through the ambient gas. We first calculate mechanical torques (MATs) resulting from specular reflection of gas atoms for seven irregular shapes: one shape of mirror symmetry, three highly irregular shapes (HIS), and three weakly irregular shapes (WIS). We find that the grain with mirror symmetry experiences negligible MATs due to its mirror-symmetry geometry. Three highly irregular shapes can produce strong MATs which exhibit some generic properties as radiative torques, while three weakly irregular shapes produce less efficient MATs. We then study grain alignment by MATs for the different angles between the drift velocity and the ambient magnetic field, for paramagnetic and superparamagnetic grains assuming efficient internal relaxation. We find that for HIS grains, MATs can align subsonically drifting grains in the same way as radiative torques, with low-J and high-J attractors. For supersonic drift, MATs can align grains with low-J and high-J attractors, analogous to radiative alignment by anisotropic radiation. We also show that the joint action of MATs and magnetic torques in grains with iron inclusions can lead to perfect MAT alignment. Our results point out the potential importance of MAT alignment for HIS grains predicted by the analytical model of Lazarian \& Hoang (2007b), although more theoretical and observational studies are required due to uncertainty in the shape of interstellar grains. We outline astrophysical environments where MAT alignment is potentially important.

Based on the new findings on the turbulent dynamo in \citet{XL16}, we examine the magnetic field amplification in the context of supernova remnants. Due to the strong ion-neutral collisional damping in the weakly-ionized interstellar medium, the dynamo in the preshock turbulence remains in the damping kinematic regime, which leads to {\it a linear-in-time growth of the magnetic field strength}. The resultant magnetic field structure enables effective diffusion upstream and shock acceleration of cosmic rays to energies above the "knee". Differently, the nonlinear dynamo in the postshock turbulence leads to {\it a linear-in-time growth of the magnetic energy} due to the turbulent magnetic diffusion. Given a weak initial field strength in the postshock region, the magnetic field saturates at a significant distance from the shock front as a result of the inefficiency of the nonlinear dynamo. This result is in a good agreement with existing numerical simulations and well explains the X-ray spots detected far behind the shock front.

In this paper, we report our detailed analysis of the new strong-coupling regime between Alfv\'en and fast modes in Poynting-dominated plasma turbulence, reported in our previous work Takamoto & Lazarian (2016), which is an important effect for many relativistic plasma phenomena, and calls for new theories of Poynting-dominated MHD turbulence. We performed numerical simulations of relativistic MHD turbulence in isothermal plasmas, and analyzed the ratio of fast to Alfv\'en mode energy. We found that the increase of the fast mode with the background $\sigma$-parameter can be observed even in isothermal plasma, showing that such a phenomena is universal in trans-Alfv\'enic turbulence in Poynting-dominated plasmas. To study the detailed energy conversion process, we also performed a series of simulations of decaying turbulence injecting pure Alfv\'en, fast, and slow modes, respectively, and investigated the development of the mode conversion from the each mode. We also found that the mode conversion between Alfv\'en and fast modes is nearly insensitive to the background temperature. Finally, we report a result of a simulation with initially fast mode dominated turbulence. It developed into a temporally strong-coupling regime, which is a strong evidence for the existence of our suggesting strong-coupling regime of fast and Alfv\'en modes. Our result suggests that the strong turbulence in Poynting-dominated plasma is very different from that in the non-relativistic plasma. It will also give an important guidance to studies of particle acceleration and non-thermal photon emission from Poynting-dominated plasma.

Recent Planck observations of dust polarization in the Galaxy have revealed that the power in $E$ mode is twice that in $B$ mode. Caldwell et. al. have formulated a theoretical model in the context of magnetohydrodynamic (MHD) turbulence and found it problematic to account for this result. In particular, they concluded that there is a very narrow range of theoretical parameters that could account for the observation. This poses a problem of whether the accepted description of MHD turbulence can apply to the interstellar medium. We revisit the problem and demonstrate that MHD turbulence corresponding to the high galactic latitudes range of Alfv\'en Mach numbers, i.e. $M_A\lesssim 0.5$, can successfully explain the available results for the $E$ to $B$ mode ratio.

We study the effect of dust absorption on the recovery velocity and density spectra as well as on the anisotropies of magnetohydrodynamic turbulence using the Velocity Channel Analysis (VCA), Velocity Coordinate Spectrum (VCS) and Velocity Centroids. The dust limits volume up to an optical depth of unity. We show that in the case of the emissivity proportional to the density of emitters, the effects of random density get suppressed for strong dust absorption intensity variations arise from the velocity fluctuations only. However, for the emissivity proportional to squared density, both density and velocity fluctuations affect the observed intensities. We predict a new asymptotic regime for the spectrum of fluctuations for large scales exceeding the physical depths to unit optical depth. The spectrum gets shallower by unity in this regime. In addition, the dust absorption removes the degeneracy resulted in the universal $K^{-3}$ spectrum of intensity fluctuations of self-absorbing medium reported by Lazarian \& Pogosyan. We show that the predicted result is consistent with the available HII region emission data. We find that for sub-Alfv\'enic and trans-Alfv\'enic turbulence one can get the information about both the magnetic field direction and the fundamental Alfv\'en, fast and slow modes that constitute MHD turbulence.

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

We apply the theory of radiative torque (RAT) alignment for studying protoplanetary disks around a T-Tauri star and perform 3D radiative transfer calculations to provide the expected maps of polarized radiation to be compared with observations, such as with ALMA. We revisit the issue of grain alignment for large grains expected in the protoplanetary disks and find that mm-sized grains at midplane do not align with the magnetic field as the Larmor precession timescale for such large grains becomes longer than the gaseous damping timescale. Hence, for these grains the RAT theory predicts that the alignment axis is determined by the grain precession with respect to the radiative flux. As a result, we expect that the polarization will be in the azimuthal direction for a face-on disk. It is also shown that if dust grains have superparamagnetic inclusions, magnetic field alignment is possible for (sub-)micron grains at the surface layer of disks, and this can be tested by mid-infrared polarimetric observations.

Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Within the standard Sweet-Parker model, this process would be too slow to explain observations (e.g. solar flares). In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen-in field lines and magnetic reconnection. In the presence of turbulence, the reconnection is independent of microscopic plasma properties, and may be much faster than previously thought, as proposed in Lazarian & Vishniac (1999) and tested in Kowal et al. (2009, 2012). However, the considered turbulence in the Lazarian-Vishniac model was imposed externally. In this work we consider reconnection-driven magnetized turbulence in realistic three-dimensional geometry initiated by stochastic noise. We demonstrate through numerical simulations that the stochastic reconnection is able to self-generate turbulence through interactions between its outflows. We analyze the statistical properties of velocity fluctuations using power spectra and anisotropy scaling, which demonstrates that the reconnection produces Kolmogorov-like turbulence, compatible with Goldreich-Sridhar (1995) model. Anisotropy statistics are, however, strongly affected by the dynamics of reconnection outflows. Once the broad turbulent region is formed, the typical anisotropy scaling $l_\parallel \propto l_\perp^{2/3}$ is formed, especially for high resolution models, were the broader range of scales is available. The decay of reconnection outflows to turbulent-like fluctuations, characterized by different anisotropy scalings, strongly depends on $\beta$ plasma. Moreover, the estimated reconnection rates are weakly dependent on the resolution, suggesting that no external processes are required to make reconnection fast.

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

The Breakthrough Starshot initiative aims to launch a gram-scale spacecraft to a speed of $v\sim 0.2$c, capable of reaching the nearest star system, $\alpha$ Centauri, in about 20 years. However, a critical challenge for the initiative is the damage to the spacecraft by interstellar gas and dust during the journey. In this paper, we quantify the interaction of a relativistic spacecraft with gas and dust in the interstellar medium. For gas bombardment, we find that damage by track formation due to heavy elements is an important effect. We find that gas bombardment can potentially damage the surface of the spacecraft to a depth of $\sim 0.1$ mm for quartz material after traversing a gas column of $N_{\rm H}\sim 2\times 10^{18}\rm cm^{-2}$ along the path to $\alpha$ Centauri, whereas the effect is much weaker for graphite material. The effect of dust bombardment erodes the spacecraft surface and produces numerous craters due to explosive evaporation of surface atoms. For a spacecraft speed $v=0.2c$, we find that dust bombardment can erode a surface layer of $\sim 0.5$ mm thickness after the spacecraft has swept a column density of $N_{\rm H}\sim 3\times 10^{17}\rm cm^{-2}$, assuming the standard gas-to-dust ratio of the interstellar medium. Dust bombardment also damages the spacecraft surface by modifying the material structure through melting. We calculate the equilibrium surface temperature due to collisional heating by gas atoms as well as the temperature profile as a function of depth into the spacecraft. Our quantitative results suggest methods for damage control, and we highlight possibilities for shielding strategies and protection of the spacecraft.

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.

This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath.

Strong Alfv\'enic turbulence develops eddy-like motions perpendicular to the local direction of magnetic fields. This local alignment induces velocity gradients perpendicular to the local direction of the magnetic field. We use this fact to propose a new technique of studying the direction of magnetic fields from observations, the Velocity Gradient Technique. We test our idea by employing the synthetic observations obtained via 3D MHD numerical simulations for different sonic and Alfv\'en Mach numbers. We calculate the velocity gradient, $\mathbf{\Omega}$, using the velocity centroids. We find that $\mathbf{\Omega}$ traces the projected magnetic field best for the synthetic maps obtained with sub-Alfv\'enic simulations providing good point-wise correspondence between the magnetic field direction and that of $\mathbf{\Omega}$. The reported alignment is much better than the alignment between the density gradients and the magnetic field and we demonstrated that it can be used to find the magnetic field strength using the Chandrasekhar-Fermi method. Our study opens a new way of studying magnetic fields using spectroscopic data.

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.

In this Letter, we report compressible mode effects on relativistic magnetohydrodynamic (RMHD) turbulence in Poynting-dominated plasmas using 3-dimensional numerical simulations. We decomposed fluctuations in the turbulence into 3 MHD modes (fast, slow, and Alfv\'en) following the procedure mode decomposition in (Cho & Lazarian 2002), and analyzed their energy spectra and structure functions separately. We also analyzed the ratio of compressible mode to Alfv\'en mode energy with respect to its Mach number. We found the ratio of compressible mode increases not only with the Alfv\'en Mach number but with the background magnetization, which indicates a strong coupling between the fast and Alfv\'en modes and appearance of a new regime of RMHD turbulence in Poynting-dominated plasmas where the fast and Alfv\'en modes strongly couples and cannot be separated, different from the non-relativistic MHD case. This finding will affect particle acceleration efficiency obtained by assuming Alfv\'enic critical balance turbulence, and will change the resulting photon spectra by non-thermal electrons.

We performed numerical calculations to test the suggestion by Desiati & Lazarian (2013) that the anisotropies of TeV cosmic rays may arise from their interactions with the heliosphere. For this purpose, we used a magnetic field model of the heliosphere and performed direct numerical calculations of particle trajectories. Unlike earlier papers testing the idea, we did not employ time-reversible techniques that are based on Liouville's theorem. We showed numerically that for scattering by the heliosphere the conditions of Liouville's theorem are not satisfied and the adiabatic approximation and time-reversibility of the particle trajectories are not valid. Our results indicate sensitivity to the magnetic structure of the heliospheric magnetic field, and we expect that this will be useful for probing this structure in future research.

We study statistical properties of synchrotron polarization emitted from media with magnetohydrodynamic (MHD) turbulence. We use both synthetic and MHD turbulence simulation data for our studies. We obtain the spatial spectrum and its derivative with respect to wavelength of synchrotron polarization arising from both synchrotron radiation and Faraday rotation fluctuations. In particular, we investigate how the spectrum changes with frequency. We find that our simulations agree with the theoretical predication in Lazarian \& Pogosyan (2016). We conclude that the spectrum of synchrotron polarization and it derivative can be very informative tools to get detailed information about the statistical properties of MHD turbulence from radio observations of diffuse synchrotron polarization. Especially, they are useful to recover the statistics of turbulent magnetic field as well as turbulent density of electrons. We also simulate interferometric observations that incorporate the effects of noise and finite telescope beam size, and demonstrate how we recover statistics of underlying MHD turbulence.

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 $\delta_{m}$ to produce high-J attractor points as functions of $q^{\max}$ and the anisotropic radiation angle relative to the magnetic field $\psi$. 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.

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

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

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

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.

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.

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.

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

We consider formation of accretion disks from a realistically turbulent molecular gas using 3D MHD simulations. In particular, we analyze the effect of the fast turbulent reconnection described by the Lazarian & Vishniac (1999) model for the removal of magnetic flux from a disk. With our numerical simulations we demonstrate how the fast reconnection enables protostellar disk formation resolving the so-called "magnetic braking catastrophe". In particular, we provide a detailed study of the dynamics of a 0.5 M$_\odot$ protostar and the formation of its disk for up to several thousands years. We measure the evolution of the mass, angular momentum, magnetic field, and turbulence around the star. We consider effects of two processes that strongly affect the magnetic transfer of angular momentum, both of which are based on turbulent reconnection: the first, "reconnection diffusion", removes the magnetic flux from the disk, the other involves the change of the magnetic field's topology, but does not change the absolute value of the magnetic flux through the disk. We demonstrate that for the first mechanism, turbulence causes a magnetic flux transport outward from the inner disk to the ambient medium, thus decreasing the coupling of the disk to the ambient material. A similar effect is achieved through the change of the magnetic field's topology from a split monopole configuration to a dipole configuration. We explore how both mechanisms prevent the catastrophic loss of disk angular momentum and compare both above turbulent reconnection mechanisms with alternative mechanisms from the literature.

We investigate how observations of synchrotron intensity fluctuations can be used to probe the sonic and Alfv\'enic 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\'enic 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.

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.

Galactic cosmic rays are believed to be mostly accelerated at supernova shocks. However, the interstellar magnetic field is too weak to efficiently accelerate galactic cosmic rays up to the highest energies, i.e. $10^{15}$ eV. A stronger magnetic field in the pre-shock region could provide the efficiency required. Bell's cosmic-ray nonresonant streaming instability has been claimed to be responsible for the amplification of precursor magnetic fields. However, an alternative mechanism has been proposed in which the cosmic-ray pressure gradient forms the shock precursor and drives turbulence, amplifying the magnetic field via the small-scale dynamo. A key ingredient for the mechanism to operate are the inhomogeneities present in the interstellar medium (ISM). These inhomogeneities are the consequence of turbulence. In this work we explore the magnetic field amplification in different ISM conditions through 3D MHD numerical simulations.

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.

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

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.

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

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

We report turbulence effects on magnetic reconnection in relativistic plasmas using 3-dimensional relativistic resistive magnetohydrodynamics simulations. We found reconnection rate became independent of the plasma resistivity due to turbulence effects similarly to non-relativistic cases. We also found compressible turbulence effects modified the turbulent reconnection rate predicted in non-relativistic incompressible plasmas; The reconnection rate saturates and even decays as the injected velocity approaches to the Alfv\'en velocity. Our results indicate the compressibility cannot be neglected when compressible component becomes about half of incompressible mode occurring when the Alfv\'en Mach number reaches about $0.3$. The obtained maximum reconnection rate is around $0.05$ to $0.1$, which will be able to reach around $0.1$ to $0.2$ if injection scales are comparable to the sheet length.

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.

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

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 $\gamma < 2$ by powerful radiation sources, which is lower than earlier estimates showing that $\gamma$ 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\mu$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.

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

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\'{e}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\'{e}nic turbulence, when the damping scale of Alfv\'{e}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\'{e}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.

In situ spacecraft data on the solar wind show events identified as magnetic reconnection with outflows and apparent "`$X$-lines" $10^{3-4}$ times ion scales. To understand the role of turbulence at these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic (MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the solar wind, however, with a hyperbolic flux-tube or "$X$-line" extending over integral length-scales.

Reflection nebulae--dense cores--illuminated by surrounding stars offer a unique opportunity to directly test our quantitative model of grain alignment based on radiative torques (RATs) and to explore new effects arising from additional torques. In this paper, we first perform detailed modeling of grain alignment by RATs for the IC 63 reflection nebula illuminated both by a nearby $\gamma$ Cas star and the diffuse interstellar radiation field. We calculate linear polarization $p_{\lambda}$ of background stars by radiatively aligned grains and explore the variation of fractional polarization ($p_{\lambda}/A_V$) with visual extinction $A_{V}$ across the cloud. Our results show that the variation of $p_{V}/A_{V}$ versus $A_{V}$ from the dayside of IC 63 to its center can be represented by a power-law ($p_{V}/A_{V}\propto A_{V}^{\eta}$) with different slopes depending on $A_{V}$. We find a shallow slope $\eta \sim- 0.1$ for $A_{V}< 3$ and a very steep slope $\eta\sim -2$ for $A_{V}> 4$. We then consider the effects of additional torques due to H$_{2}$ formation and model grain alignment by joint action of RATs and H$_2$ torques. We find that $p_{V}/A_{V}$ tends to increase with an increasing magnitude of H$_{2}$ torques. In particular, the theoretical predictions obtained for $p_{V}/A_{V}$ and peak wavelength $\lambda_{\max}$ in this case show an improved agreement with the observational data. Our results reinforce the predictive power of the RAT alignment mechanism in a broad range of environmental conditions and show the effect of pinwheel torques in environments with efficient H$_2$ formation. Physical parameters involved in H$_2$ formation may be constrained using detailed modeling of grain alignment combined with observational data. In addition, we discuss implications of our modeling for interpreting latest observational data by {\it Planck} and other ground-based instruments.

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

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

We investigate the nature of the Alfv\'enic 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\'en 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\'enic, 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).

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.

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.

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.

Magnetic reconnection is a process that changes magnetic field topology in highly conducting fluids. Traditionally, magnetic reconnection was associated mostly with solar flares. In reality, the process must be ubiquitous as astrophysical fluids are magnetized and motions of fluid elements necessarily entail crossing of magnetic frozen in field lines and magnetic reconnection. We consider magnetic reconnection in realistic 3D geometry in the presence of turbulence. This turbulence in most astrophysical settings is of pre-existing nature, but it also can be induced by magnetic reconnection itself. In this situation turbulent magnetic field wandering opens up reconnection outflow regions, making reconnection fast. We discuss Lazarian \& Vishniac (1999) model of turbulent reconnection, its numerical and observational testings, as well as its connection to the modern understanding of the Lagrangian properties of turbulent fluids. We show that the predicted dependences of the reconnection rates on the level of MHD turbulence make the generally accepted Goldreich \& Sridhar (1995) model of turbulence self-consistent. Similarly, we argue that the well-known Alfv\'en theorem on flux freezing is not valid for the turbulent fluids and therefore magnetic fields diffuse within turbulent volumes. This is an element of magnetic field dynamics that was not accounted by earlier theories. For instance, the theory of star formation that was developing assuming that it is only the drift of neutrals that can violate the otherwise perfect flux freezing, is affected and we discuss the consequences of the turbulent diffusion of magnetic fields mediated by reconnection.

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

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\mu$m) and very small ($a\sim 0.001\mu$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 $\lambda_{\max}<0.55\mu$m, which exhibits an excess UV polarization relative to the Serkowski law, compared to the typical case $\lambda_{\max}=0.55\mu$m. We interpret the correlation between the systematic increase of the UV polarization relative to maximum polarization (i.e. of $p(6\mu m^{-1})/p_{\max}$) with $\lambda_{\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 $\lambda_{\max}$ and $p(6\mu m^{-1})/p_{\max}$.[Abridged]

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

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$\alpha$ 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$\alpha$ data.

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

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

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\'enic 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 (\epsilon_+/\epsilon_-)^n \$ with n>2). These results are consistent with those obtained for imbalanced non-relativistic Alfv\'enic turbulence. This corresponds well to the earlier reported similarity of the relativistic and non-relativistic balanced magnetic turbulence.

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.

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

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\'enic 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\'enic 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\'enic clouds, ensuring that the detected magnetic flux removal is due to to the action of the RD rather than to numerical diffusivity.

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

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.

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

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.

Astrophysical plasmas can have parameters vastly different from the more studied laboratory and space plasmas. In particular, the magnetic fields can be the dominant component of the plasma, with energy-density exceeding the particle rest-mass energy density. Magnetic fields then determine the plasma dynamical evolution, energy dissipation and acceleration of non-thermal particles. Recent data coming from astrophysical high energy missions, like magnetar bursts and Crab nebula flares, point to the importance of magnetic reconnection in these objects. In this review we outline a broad spectrum of problems related to the astrophysical relevant processes in magnetically dominated relativistic plasmas. We discuss the problems of large scale dynamics of relativistic plasmas, relativistic reconnection and particle acceleration at reconnecting layers, turbulent cascade in force-fee plasmas. A number of astrophysical applications are also discussed.

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\'enic 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\'enic 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\'enic turbulence the spectral slopes are steeper with values near -3 which are similar to the very optically thick regime.

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.

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.

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 $\tau \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 $\sigma_{ln(\Sigma/\Sigma_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 $\tau \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.

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.

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.

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 ({\mu}) 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 {\mu} 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.

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.

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

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.

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 \rho^{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.

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 ($\sigma^2$) and the sonic Mach number (${\cal M}_s$). This is in part due to the fact that the $\sigma^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 $\sigma_{\Sigma/\Sigma_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 $\sigma_{\rho/\rho_0}^2$-${\cal M}_s$ trend but includes a scaling parameter $A$ such that: $\sigma_{ln(\Sigma/\Sigma_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 $\sigma^2$ for investigations of the relation between the sonic Mach number and the PDF variance in addition to existing PDF sonic Mach number relations.

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

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.

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.

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.

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

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

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.

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.

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

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.

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

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.

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.

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.

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.

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.

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.

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

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

Astrophysical fluids are generally turbulent and this preexisting turbulence must be taken into account for the models of magnetic reconnection which are attepmted to be applied to astrophysical, solar or heliospheric environments. In addition, reconnection itself induces turbulence which provides an important feedback on the reconnection process. In this paper we discuss both theoretical model and numerical evidence that magnetic reconnection gets fast in the approximation of resistive MHD. We consider the relation between the Lazarian & Vishniac turbulent reconnection theory and Lapenta's numerical experiments testifying of the spontaneous onset of turbulent reconnection in systems which are initially laminar.

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.

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

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.

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$\alpha$-correlated and 100 $\mu$m-correlated emission spectra.

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

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

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.

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

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.

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.

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.

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

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.

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.

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

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.

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.

We discuss a novel mechanism of dust acceleration which may dominate for particles smaller than $\sim0.1~\mu$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.

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.

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.

Present day cosmic microwave background (CMB) studies require more accurate removal of Galactic foreground emission. In this paper, we consider a way of filtering out the diffuse Galactic fluctuations on the basis of their statistical properties, namely, the power-law spectra of fluctuations. We focus on the statistical properties of two major Galactic foregrounds that arise from magnetized turbulence, namely, diffuse synchrotron emission and thermal emission from dust and describe how their power laws change with the Galactic latitude. We attribute this change to the change of the geometry of the emission region and claim that the universality of the turbulence spectrum provides a new way of removing Galactic foregrounds. We discuss and demonstrate how we can make use of our findings to remove Galactic foregrounds using a template of spatial fluctuations. In particular, we consider examples of spatial filtering of a foreground at small scales, when the separation into CMB signal and foregrounds is done at larger scales. We demonstrate that the new technique of spatial filtering of foregrounds may be promising for recovering the CMB signal in a situation when foregrounds are known at a scale different from the one under study. It can also improve filtering by combining measurements obtained at different scales.

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.

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.

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.

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

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.

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.

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.

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

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

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.

We use the data of Wisconsin H$\alpha$ 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$.

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.

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\'en Mach numbers of turbulence. This opens avenues for using Tsallis distributions to study the dynamical and perhaps magnetic states of interstellar gas.

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{\alpha}$-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.

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

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

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.

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.

We present arcsecond-resolution Submillimeter Array (SMA) polarimetric observations of the 880 um continuum emission from the protoplanetary disks around two nearby stars, HD 163296 and TW Hydrae. Although previous observations and theoretical work have suggested that a 2-3% polarization fraction should be common for the millimeter continuum emission from such disks, we detect no polarized continuum emission above a 3-sigma upper limit of 7 mJy in each arcsecond-scale beam, or <1% in integrated continuum emission. We compare the SMA upper limits with the predictions from the exploratory Cho & Lazarian (2007) model of polarized emission from T Tauri disks threaded by toroidal magnetic fields, and rule out their fiducial model at the ~10-sigma level. We explore some potential causes for this discrepancy, focusing on model parameters that describe the shape, magnetic field alignment, and size distribution of grains in the disk. We also investigate related effects like the magnetic field strength and geometry, scattering off of large grains, and the efficiency of grain alignment, including recent advances in grain alignment theory, which are not considered in the fiducial model. We discuss the impact each parameter would have on the data and determine that the suppression of polarized emission plausibly arises from rounding of large grains, reduced efficiency of grain alignment with the magnetic field, and/or some degree of magnetic field tangling (perhaps due to turbulence). A poloidal magnetic field geometry could also reduce the polarization signal, particularly for a face-on viewing geometry like the TW Hya disk. The data provided here offer the most stringent limits to date on the polarized millimeter-wavelength emission from disks around young stars.

We performed a series of high-resolution (up to 1024^3) direct numerical simulations of hydro and MHD strong turbulence. We found that for simulations with normal viscosity the slopes for spectra of MHD are similar, although slightly more shallower than for hydro simulations. However, for simulations with hyperviscosity the slopes were very different, for instance, the slopes for hydro simulations showed pronounced and well-defined bottleneck effect, while the MHD slopes were relatively much less affected. We believe that this is indicative of MHD strong turbulence being less local than Kolmogorov turbulence. This calls for revision of MHD strong turbulence models that assume local "as-in-hydro case" cascading. Nonlocality of MHD turbulence casts doubt on numerical determination of the slopes with currently available (512^3--1024^3) numerical resolutions, including simulations with normal viscosity. We also measure various so-called alignment effects and discuss their influence on the turbulent cascade.

In this paper we systematically study the spectrum and structure of incompressible MHD turbulence by means of high resolution direct numerical simulations. We considered both balanced and imbalanced (cross-helical) cases and simulated sub-Alfvenic as well as trans-Alfvenic turbulence. This paper extends numerics preliminarily reported in Beresnyak & Lazarian (2008). We confirm that driven imbalanced turbulence has a stationary state even for high degrees of imbalance. Our major finding is that the structure of the dominant and subdominant Alfvenic components are notably different. Using the most robust observed quantities, such as the energy ratio, we were able to reject several existing models of strong imbalanced turbulence.

We discuss a model of cosmic ray acceleration that accounts for the observations of anomalous cosmic rays by Voyager 1 and 2. The model appeals to fast magnetic reconnection rather than shocks as the driver of acceleration. The ultimate source of energy is associated with magnetic field reversals that occur in the heliosheath. It is expected that the magnetic field reversals will occur throughout the heliosheath, but especially near the heliopause where the flows slows down and diverge in respect to the interstellar wind and also in the boundary sector-in the heliospheric current sheet. While the First Order Fermi acceleration theory within reconnection layers is in its infancy, the available predictions do not contradict the available data on anomalous cosmic ray spectra measured by the spacecrafts. We argue that the Voyager data can be one of the first pieces of evidence favoring the acceleration within regions of fast magnetic reconnection, which we argue is a widely spread astrophysical process.

We present numerical simulations of electron magnetohydrodynamic (EMHD) and electron reduced MHD (ERMHD) turbulence. Comparing scaling relations, we find that both EMHD and ERMHD turbulence show similar spectra and anisotropy. We develop new techniques to study anisotropy of EMHD turbulence. Our detailed study of anisotropy of EMHD turbulence supports our earlier result of k_par ~ k_perp^(1/3) scaling. We find that the high-order statistics show a scaling that is similar to the She-Leveque scaling. We observe that the bispectra, which characterize the interaction of different scales within the turbulence cascade, are very different for EMHD and MHD turbulence. We show that both decaying and driven EMHD turbulence have the same statistical properties. We calculate the probability distribution functions (PDFs) of MHD and EMHD turbulence and compare them with those of interplanetary turbulence. We find that, as in the case of the solar wind, the PDFs of the increments of magnetic field strength in MHD and EMHD turbulence are well described by the Tsallis distribution. We discuss implications of our results for astrophysical situations, including the ADAFs and magnetic reconnection.

We study the effects of turbulence on magnetic reconnection using 3D numerical simulations. This is the first attempt to test a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian & Vishniac (1999). This model predicts that weak turbulence, generically present in most of astrophysical systems, enhances the rate of reconnection by reducing the transverse scale for reconnection events and by allowing many independent flux reconnection events to occur simultaneously. As a result the reconnection speed becomes independent of Ohmic resistivity and is determined by the magnetic field wandering induced by turbulence. To quantify the reconnection speed we use both an intuitive definition, i.e. the speed of the reconnected flux inflow, as well as a more sophisticated definition based on a formally derived analytical expression. Our results confirm the predictions of the Lazarian & Vishniac model. In particular, we find that Vrec Pinj^(1/2), as predicted by the model. The dependence on the injection scale for some of our models is a bit weaker than expected, i.e. l^(3/4), compared to the predicted linear dependence on the injection scale, which may require some refinement of the model or may be due to the effects like finite size of the excitation region. The reconnection speed was found to depend on the expected rate of magnetic field wandering and not on the magnitude of the guide field. In our models, we see no dependence on the guide field when its strength is comparable to the reconnected component. More importantly, while in the absence of turbulence we successfully reproduce the Sweet-Parker scaling of reconnection, in the presence of turbulence we do not observe any dependence on Ohmic resistivity, confirming that our reconnection is fast.

While the problem of grain alignment was posed more than 60 years ago the quantitative model of grain alignment that can account for the observed polarization arising from aligned grains has been formulated only recently. The quantitative predictions of the radiative torque mechanism, which is currently accepted as the dominant mechanism of grain alignment, open avenues to tracing magnetic fields in various astrophysical environments, including diffuse and dense interstellar gas, molecular clouds, circumstellar environments, accretion disks, comet tails, Zodiacal dust etc. At the same time, measurements of the absolute value of polarization and its variations can, in addition, provide unique information about the dust composition and dust environment. In the review I describe the analytical model describing well radiative torques acting on irregular grains and discuss how the alignment induced by radiative torques varies in the presence of superparamagnetic inclusions and pinwheel torques, e.g. arising from the molecular hydrogen formation over grain surface. I also describe observations that can establish whether grains are superparamagnetic and whether recoils from molecular hydrogen formations are powerful enough to give rise to substantial uncompensated torques. Answering to these questions should allow for reliable modeling of astrophysical polarization with numerous important applications, from accounting for dust contribution in Cosmic Microwave Background polarization studies to obtaining magnetic field strength using Chandrasekhar-Fermi technique.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of stars and stellar evolution. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of Galactic science. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of planetary systems and star formation. Community-based White Paper to Astro2010 in response to the call for such papers.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of extragalactic astrophysics. Community-based White Paper to Astro2010 in response to the call for such papers.

We discuss the outstanding issues of the interstellar medium which will depend on the application of knowledge from plasma physics. We particularly advocate attention to recent developments in experimental plasma physics, and urge that the astronomical community consider support of these experiments in the next decade.

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of cosmology. Community-based White Paper to Astro2010 in response to the call for such papers.

A combination of observation, theory, modeling, and laboratory plasma experiments provides a multifaceted approach to develop a much greater understanding of how magnetic fields arise in galactic settings and how these magnetic fields mediate important processes that affect the dynamics, distribution, and composition of galactic plasmas. An important emphasis below is the opportunity to connect laboratory experiments to astrophysics. This approach is especially compelling for the galactic neighborhood, where the distribution and character of magnetic fields can be observed with greater detail than what is possible elsewhere in the universe. The ability to produce laboratory plasmas with unparalleled accessibility permits an even greater level of detail to be assessed and exposed. Theory and modeling provide fundamental ways to understand important processes, and they act as the bridge to connect experimental validation to astronomical observations. In many cases the studies that utilize this approach can make use of existing laboratory facilities, resulting in a cost that is quite small compared to the cost of measurements in dedicated space missions.

The major questions relevant to star and planet formation are: What controls the rate, efficiency, spatial clustering, multiplicity, and initial mass function of star formation, now and in the past? What are the major feedback mechanisms through which star formation affects its environment? What controls the formation and orbital parameters of planets, especially terrestrial planets? These questions cannot be fully addressed without understanding key magnetohydrodynamics (MHD) and plasma physics processes. Although some of these basic problems have long been considered intractable, attacking them through a combination of laboratory experiment, theory, and numerical simulation is now feasible, and would be fruitful. Achieving a better understanding of these processes is critical to interpreting observations, and will form an important component of astrophysical models. These models in turn will serve as inputs to other areas of astrophysics, e.g. cosmology and galaxy formation.

We study the alignment of grains subject to both radiative torques and pinwheel torques while accounting for thermal flipping of grains. By pinwheel torques we refer to all systematic torques that are fixed in grain body axes, including the radiative torques arising from scattering and absorption of isotropic radiation. We discuss new types of pinwheel torques, which are systematic torques arising from infrared emission and torques arising from the interaction of grains with ions and electrons in hot plasma. We show that both types of torques are long-lived, i.e. may exist longer than gaseous damping time. We compare these torques with the torques introduced by E. Purcell, namely, torques due to H$_2$ formation, the variation of accommodation coefficient for gaseous collisions and photoelectric emission. Furthermore, we revise the Lazarian & Draine model for grain thermal flipping. We calculate mean flipping timescale induced by Barnett and nuclear relaxation for both paramagnetic and superparamagnetic grains, in the presence of stochastic torques associated with pinwheel torques, e.g. the stochastic torques arising from H$_2$ formation, and gas bombardment. We show that the combined effect of internal relaxation and stochastic torques can result in fast flipping for sufficiently small grains and, because of this, they get thermally trapped, i.e. rotate thermally in spite of the presence of pinwheel torques. For sufficiently large grains, we show that the pinwheel torques can increase the degree of grain alignment achievable with the radiative torques by increasing the magnitude of the angular momentum of low attractor points and/or by driving grains to new high attractor points.

We present an overview of the Velocity Coordinate Spectrum (VCS), a new technique for studying astrophysical turbulence that utilizes the line-of-sight statistics of Doppler-broadened spectral lines. We consider the retrieval of turbulence spectra from emission intensity observations of both high and low spatial resolution and find that the VCS allows one to study turbulence even when the emitting turbulent volume is not spatially resolved. This opens interesting prospects for using the technique for extragalactic research. VCS developed for spectral emission lines is applicable to absorption lines as well if the optical depth is used instead of intensity. VCS for absorption lines in point-source spectra benefit from effectively narrow beam and does not require dense sky coverage by sampling directions. Even strongly saturated absorption lines still carry the information about the small scale turbulence, albeit limited to the wings of a line. Combining different absorption lines one can develop tomography of the turbulence in the interstellar gas in all its complexity.

We review the use of velocity centroids statistics to recover information of interstellar turbulence from observations. Velocity centroids have been used for a long time now to retrieve information about the scaling properties of the turbulent velocity field in the interstellar medium. We show that, while they are useful to study subsonic turbulence, they do not trace the statistics of velocity in supersonic turbulence, because they are highly influenced by fluctuations of density. We show also that for sub-Alfv\'enic turbulence (both supersonic and subsonic) two-point statistics (e.g. correlation functions or power-spectra) are anisotropic. This anisotropy can be used to determine the direction of the mean magnetic field projected in the plane of the sky.

Alignment of dust by radiative torques (RATs) has proven to be the most promising mechanism to explain alignment in various astrophysical environments, from comet atmospheres to accretion disks, molecular clouds, and diffuse interstellar gas. We discuss some of the major advances, which include, first of all, formulating of the analytical model of RATs. This model was shown to reproduce well the torques acting on actual irregular dust grains and allowed studies of the parameter space for which the alignment happens with long axes perpendicular and parallel to the magnetic field. Such a study resulted in an important conclusion that, without any paramagnetic relaxation, the RAT alignment always happens for interstellar grains with long axes perpendicular to the magnetic field. We show that the gaseous bombardment in some cases increases the degree of alignment by knocking out grains from the positions of imperfect alignment when the grains rotate slowly to more stable positions of perfect alignment where grains rotate fast. In terms of pinwheel torques, important revisions have been made in the Lazarian and Draine model of grain flipping and thermal trapping. Those, however, do not change the major conclusion that very small grains (i.e. grain size smaller than ~0.03 micron) should be marginally aligned. Recent work made the RAT alignment a predictive theory which is ready for quantitative modeling of astrophysical polarization. We predict that the microwave emission from the Zodiacal dust presents an important contaminant, which should be included into foreground polarization templates.

Recent observation of 850 micron sub-mm polarization from T Tauri disks opens up the possibility of studying magnetic field structure within protostellar disks. The degree of polarization is around 3 % and the direction of polarization is perpendicular to the disk. Since thermal emission from dust grains dominates the spectral energy distribution at the sub-mm/FIR regime, dust grains are thought to be the cause of the polarization. We discuss grain alignment by radiation and we explore the efficiency of dust alignment in T Tauri disks. Calculations show that dust grains located far away from the Central proto-star are more efficiently aligned. In the presence of a regular magnetic field, the aligned grains produce polarized emission in sub-mm/FIR wavelengths. The direction of polarization is perpendicular to the local magnetic field direction. When we use a recent T Tauri disk model and take a Mathis-Rumpl-Nordsieck-type distribution with maximum grain size of 500-1000 $\mu$m, the degree of polarization is around 2-3 % level at wavelengths larger than $\sim100\mu$m. Our study indicates that multifrequency infrared polarimetric studies of protostellar disks can provide good insights into the details of their magnetic structure. We also provide predictions for polarized emission for disks viewed at different wavelengths and viewing angles.

Astrophysical fluids are generically turbulent, which means that frozen-in magnetic fields are, at least, weakly stochastic. Therefore realistic studies of astrophysical magnetic reconnection should include the effects of stochastic magnetic field. In the paper we discuss and test numerically the Lazarian & Vishniac (1999) model of magnetic field reconnection of weakly stochastic fields. The turbulence in the model is assumed to be subAlfvenic, with the magnetic field only slightly perturbed. The model predicts that the degree of magnetic field stochasticity controls the reconnection rate and that the reconnection can be fast independently on the presence or absence of anomalous plasma effects. For testing of the model we use 3D MHD simulations. To measure the reconnection rate we employ both the inflow of magnetic flux and a more sophisticated measure that we introduce in the paper. Both measures of reconnection provide consistent results. Our testing successfully reproduces the dependences predicted by the model, including the variations of the reconnection speed with the variations of the injection scale of turbulence driving as well as the intensity of driving. We conclude that, while anomalous and Hall-MHD effects in particular circumstances may be important for the initiation of reconnection, the generic astrophysical reconnection is fast due to turbulence, irrespectively of the microphysical plasma effects involved. This conclusion justifies numerical modeling of many astrophysical environments, e.g. interstellar medium, for which plasma-effect-based collisionless reconnection is not applicable.

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

Polarimeric maps have been used on the characterization of the magnetic field in molecular clouds. However, it is difficult to determine the 3-dimensional properties of these regions from the projected maps. In that case, numerical simulations can be used as benchmarks for polarimetric measurements, and evetually reveal more about the interplay of turbulence and the magnetic field lines. In this work we provide a number of MHD numerical simulations of turbulent molecular clouds and created their synthetic dust emission polarization maps, varying the direction of the observer. We determined the correlation of emission intensity and polarization degree for the simulated models. We were able to reproduce the decay on polarization degree at denser regions without any assumption regarding the properties of the dusty component. The anti-correlation arises from the simple cancelation of the polarization vectors along the line of sight. This effect is amplified within denser regions as the magnetic field configuration becomes more complex. We studied the probability distribution, the power spectrum and the structure function of the polarization angles. These statistical analysis revealed strong defferences depending on the turbulent regime (i.e. sub/supersonic and sub/super-Alfvenic). Therefore, these methods can be used on polarimetric observations to characterize the dynamics of molecular clouds. We also presented a modified Chandrashekhar-Fermi method to obtain the intensity of the local magnetic field. The proposed formulation showed no limitations regarding orientation or turbulent regime.

We discuss the model of magnetic field reconnection in the presence of turbulence introduced by us approximately ten years ago. The model does not require any plasma effects to be involved in order to make the reconnection fast. In fact, it shows that the degree of magnetic field stochasticity controls the reconnection. The turbulence in the model is assumed to be subAlfvenic, with the magnetic field only slightly perturbed. This ensures that the reconnection happens in generic astrophysical environments and the model does not appeal to any unphysical concepts, similar to the turbulent magnetic diffusivity concept, which is employed in the kinematic magnetic dynamo. The interest to that model has recently increased due to successful numerical testings of the model predictions. In view of this, we discuss implications of the model, including the first-order Fermi acceleration of cosmic rays, that the model naturally entails, bursts of reconnection, that can be associated with Solar flares, as well as, removal of magnetic flux during star-formation.

We study the effects of turbulence on magnetic reconnection using three-dimensional numerical simulations. This is the first attempt to test a model of fast magnetic reconnection proposed by Lazarian & Vishniac (1999), which assumes the presence of weak, small-scale magnetic field structure near the current sheet. 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. We performed a number of simulations to test the dependencies of the reconnection speed, defined as the ratio of the inflow velocity to the Alfven speed, on the turbulence power, the injection scale and resistivity. Our results show that turbulence significantly affects the topology of magnetic field near the diffusion region and increases the thickness of the outflow region. We confirm the predictions of the Lazarian & Vishniac model. In particular, we report the growth of the reconnection speed proportional to ~ V^2, where V is the amplitude of velocity at the injection scale. It depends on the injection scale l as ~ (l/L)^(2/3), where L is the size of the system, which is somewhat faster but still roughly consistent with the theoretical expectations. We also show that for 3D reconnection the Ohmic resistivity is important in the local reconnection events only, and the global reconnection rate in the presence of turbulence does not depend on it.

We present statistical analysis of diffuse Galactic synchrotron emission and polarized thermal emission from dust. Both Galactic synchrotron emission and polarized thermal emission from dust reflect statistics of magnetic field fluctuations and, therefore, Galactic turbulence. We mainly focus on the relation between observed angular spectra and underlying turbulence statistics. Our major findings are as follows. First, we find that magnetohydrodynamic (MHD) turbulence in the Galaxy can indeed explain diffuse synchrotron emission from high galactic latitude. Our model calculation suggests that either a one-component extended halo model or a two-component model, an extended halo component (scale height > 1kpc) plus a local component, can explain the observed angular spectrum of the synchrotron emission. However, discrete sources seem to dominate the spectrum for regions near the Galactic plane. Second, we study how star-light polarization is related with polarized emission from thermal dust. We also discuss the expected angular spectrum of polarized emission from thermal dust. Our model calculations suggest that C_l\propto l^{-11/3} for l > 1000 and a shallower spectrum for l < 1000.

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.

Grain alignment is a notoriously difficult problem, that is extremely rich in underlying physics. The long history of attempts theoretical handling of the problem resulted in rather sceptical approach to the theory on the part of some polarimetry practitioners. However, recently the theory has been very successful in accounting for the observational data. In view of that, I present a very concise discussion of the most promising mechanism of grain alignment, namely, radiative torque alignment. In particular, I discuss a new analytical model which excellently reproduces properties of radiative torques, as well as the application of the model to predicting the degree of alignment.

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)

Astrophysical fluids, including interstellar and interplanetary medium, are magnetized and turbulent. Their appearance, evolution, and overall properties are determined by the magnetic turbulence that stirs it. We argue that examining magnetic turbulence at a fundamental level is vital to understanding many processes. A point that frequently escapes the attention of researchers is that magnetic turbulence cannot be confidently understood only using "brute force" numerical approaches. In this review we illustrate this point on a number of examples, including intermittent heating of plasma by turbulence, interactions of turbulence with cosmic rays and effects of turbulence on the rate of magnetic reconnection. We show that the properties of magnetic turbulence may vary considerably in various environments, e.g. imbalanced turbulence in solar wind differs from balanced turbulence and both of these differ from turbulence in partially ionized gas. Appealing for the necessity of more observational data on magnetic fields, we discuss a possibility of studying interplanetary turbulence using alignment of Sodium atoms in the tail of comets.

Turbulent motions induce Doppler shifts of observable emission and absorption lines motivating studies of turbulence using precision spectroscopy. We provide the numerical testing of the two most promising techniques, Velocity Channel Analysis (VCA) and Velocity Coordinate Spectrum (VCS). We obtain an expression for the shot noise that the discretization of the numerical data entails and successfully test it. We show that the numerical resolution required for recovering the underlying turbulent spectrum from observations depend on the spectral index of velocity fluctuations, which makes low resolution testing misleading. We demonstrate numerically that, dealing with absorption lines, sampling of turbulence along just a dozen directions provides a high quality spectrum with the VCS technique.

Polarized microwave emission from dust is an important foreground that may contaminate polarized CMB studies unless carefully accounted for. Modeling of polarization from dust requires a quantitative understanding of grain alignment. I review the current status of grain alignment theory outlining recent advances in quantitative description of the alignment. In particular, I show that the grain-alignment theory is a predictive one, and its results nicely match observations. Those indicate that the most important process of alignment is related to radiative torques acting on irregular grains. The recently developed analytical model of radiative torque alignment has proven to be a very efficient tool for predicting the degree of grain alignment. We expect the alignment theory to further mature before CMBPol flight, which would ensure a better accounting for the dust-related polarization. At the same time, CMBPol should provide the additional testing of grain alignment, clarifying the reliability of polarimetry for tracing of magnetic field.

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.

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.

We consider stationary, forced, imbalanced, or cross-helical MHD Alfvenic turbulence where the waves traveling in one direction have higher amplitudes than the opposite waves. This paper is dedicated to so-called strong turbulence, which cannot be treated perturbatively. Our main result is that the anisotropy of the weak waves is stronger than the anisotropy of a strong waves. We propose that critical balance, which was originally conceived as a causality argument, has to be amended by what we call a propagation argument. This revised formulation of critical balance is able to handle the imbalanced case and reduces to old formulation in the balanced case. We also provide phenomenological model of energy cascading and discuss possibility of self-similar solutions in a realistic setup of driven turbulence.

We discuss a new technique of studying magnetic fields in diffuse astrophysical media, e.g. interstellar and intergalactic gas/plasma. This technique is based on the angular momentum alignment of atoms and ions in their ground or metastable states. As the life-time of atoms in such states is long, the alignment induced by anisotropic radiation is susceptible to weak magnetic fields ($1{\rm G}\gtrsim B\gtrsim0.1\mu$G). The alignment reveals itself in terms of the polarization of the absorbed and emitted light. A variety of atoms with fine or hyperfine splitting of the ground or metastable states exhibit the alignment and the resulting polarization degree in some cases exceeds 20%. We show that in the case of absorption the polarization direction is either parallel or perpendicular to magnetic field, while more complex dependencies emerge for the case of emission of aligned atoms. We show that the corresponding studies of magnetic fields can be performed with optical and UV polarimetry. A unique feature of these studies is that they can reveal the 3D orientation of magnetic field. In addition, we point out that the polarization of the radiation arising from the transitions between fine and hyperfine states of the ground level can provide yet another promising diagnostics of magnetic fields, including the magnetic fields in the Early Universe. We mention several cases of interplanetary, circumstellar and interstellar magnetic fields for which the studies of magnetic fields using ground state atomic alignment effect are promising.

We present an analysis of the distribution of H-alpha emission measures for the warm ionized medium (WIM) of the Galaxy using data from the Wisconsin H-Alpha Mapper (WHAM) Northern Sky Survey. Our sample is restricted to Galactic latitudes |b| > 10. We removed sightlines intersecting nineteen high-latititude classical H II regions, leaving only sightlines that sample the diffuse WIM. The distribution of EM sin |b| for the full sample is poorly characterized by a single normal distribution, but is extraordinarily well fit by a lognormal distribution, with <log EM sin |b|> = 0.146 +/- 0.001 and standard deviation 0.190 +/- 0.001. <log EM sin |b|> drops from 0.260 +/- 0.002 at Galactic latitude 10<|b|<30 to 0.038 +/- 0.002 at Galactic latitude 60<|b|<90. The distribution may widen slightly at low Galactic latitude. We compare the observed EM distribution function to the predictions of three-dimensional magnetohydrodynamic simulations of isothermal turbulence within a non-stratified interstellar medium. We find that the distribution of EM sin |b| is well described by models of mildy supersonic turbulence with a sonic Mach number of ~1.4-2.4. The distribution is weakly sensitive to the magnetic field strength. The model also successfully predicts the distribution of dispersion measures of pulsars and H-alpha line profiles. In the best fitting model, the turbulent WIM occupies a vertical path length of 400-500 pc within the 1.0-1.8 kpc scale height of the layer. The WIM gas has a lognormal distribution of densities with a most probable electron density n_{pk} = 0.03 cm^{-3}. We also discuss the implications of these results for interpreting the filling factor, the power requirement, and the magnetic field of the WIM.

In this paper, genus statistics have been applied to an HI column density map of the Small Magellanic Cloud in order to study its topology. To learn how topology changes with the scale of the system, we provide the study of topology for column density maps at varying resolution. To evaluate the statistical error of the genus we randomly reassign the phases of the Fourier modes while keeping the amplitudes. We find, that at the smallest scales studied ($40 {pc}\leq\lambda\leq 80 {pc}$) the genus shift is in all regions negative, implying a clump topology. At the larger scales ($110 {pc}\leq\lambda\leq 250 {pc}$) the topology shift is detected to be negative in 4 cases and positive (``swiss cheese'' topology) in 2 cases. In 4 regions there is no statistically significant topology shift at large scales.

We address the problem of particle acceleration in solar flares by fast modes which may be excited during the reconnection and undergo cascade and are subjected to damping. We extend the calculations beyond quasilinear approximation and compare the acceleration and scattering by transit time damping and gyroresonance interactions. We find that the acceleration is dominated by the so called transit time damping mechanism. We estimate the total energy transferred into particles, and show that our approach provides sufficiently accurate results We compare this rate with energy loss rate. Scattering by fast modes appears to be sufficient to prevent the protons from escaping the system during the acceleration. Confinement of electrons, on the other hand, requires the existence of plasma waves. Electrons can be accelerated to GeV energies through the process described here for solar flare conditions.

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

It is now well established that stellar winds of hot stars are fragmentary and that the X-ray emission from stellar winds has a strong contribution from shocks in winds. Chandra high spectral resolution observations of line profiles of O and B stars have shown numerous properties that had not been expected. Here we suggest explanations by considering the X-rays as arising from bow shocks that occur where the stellar wind impacts on spherical clumps in the winds. We use an accurate and stable numerical hydrodynamical code to obtain steady-state physical conditions for the temperature and density structure in a bow shock. We use these solutions plus analytic approximations to interpret some major X-ray features: the simple power-law distribution of the observed emission measure derived from many hot star X-ray spectra and the wide range of ionization stages that appear to be present in X-ray sources throughout the winds. Also associated with the adiabatic cooling of the gas around a clump is a significant transverse velocity for the hot plasma flow around the clumps, and this can help to understand anomalies associated with observed line widths, and the differences in widths seen in stars with high and low mass-loss rates. The differences between bow shocks and the planar shocks that are often used for hot stars are discussed. We introduce an ``on the shock'' (OTSh) approximation that is useful for interpreting the X-rays and the consequences of clumps in hot star winds and elsewhere in astronomy.

We continue our studies of atomic alignment in diffuse media, in particularly, in interstellar and circumstellar media, with the goal of developing new diagnostics of magnetic fields in these environments. We understand atomic alignment as alignment of atoms or ions in their ground state. Such atoms are sensitive to weak magnetic fields. In particular, we provide predictions of the polarization that arises from astrophysically important aligned atoms (ions) with fine structure of the ground level, namely, OI and SII and Ti II. Unlike our earlier papers which dealt with weak fields only, a substantial part of our current paper is devoted to the studies of atomic alignment when magnetic fields get strong enough to affect the emission from the excited level, i.e. with the regime when the magnetic splitting is comparable to the line-width. This is a regime of Hanle effect modified by the atomic alignment. Using an example of emission and absorption lines of SII ion we demonstrate how polarimetric studies can probe magnetic fields in circumstellar regions and accretion disks. In addition, we show that atomic alignment induced by anisotropic radiation can induce substantial variations of magnetic dipole transitions within the ground state, thus affecting abundance studies based on this emission. Moreover, the radio emission is polarized, provides a new way to study magnetic fields, e.g. at the epoch of Universe reionization.

We continue our work on developing techniques for studying turbulence with spectroscopic data. We show that Doppler-broadened absorption spectral lines, in particularly, saturated absorption lines, can be used within the framework of the earlier-introduced technique termed the Velocity Coordinate spectrum (VCS). The VCS relates the statistics of fluctuations along the velocity coordinate to the statistics of turbulence, thus it does not require spatial coverage by sampling directions in the plane of the sky. We consider lines with different degree of absorption and show that for lines of optical depth less than one, our earlier treatment of the VCS developed for spectral emission lines is applicable, if the optical depth is used instead of intensity. This amounts to correlating the logarithms of absorbed intensities. For larger optical depths and saturated absorption lines, we show, that the amount of information that one can use is, inevitably, limited by noise. In practical terms, this means that only wings of the line are available for the analysis. In terms of the VCS formalism, this results in introducing an additional window, which size decreases with the increase of the optical depth. As a result, strongly saturated absorption lines carry the information only about the small scale turbulence. Nevertheless, the contrast of the fluctuations corresponding to the small scale turbulence increases with the increase of the optical depth, which provides advantages for studying turbulence combining lines with different optical depths. Combining different absorption lines one can tomography turbulence in the interstellar gas in all its complexity.

We consider grains with superparamagnetic inclusions and report two new condensed matter effects that can enhance the internal relaxation of the energy of a wobbling grain, namely, superparamagnetic Barnett relaxation, as well as, an increase of frequencies for which nuclear relaxation becomes important. This findings extends the range of grain sizes for which grains are thermally trapped, i.e. rotate thermally, in spite of the presence of uncompensated pinwheel torques. In addition, we show that the alignment of dust grains by radiative torques gets modified for superparamagnetic grains, with grains obtaining perfect alignment with respect to magnetic fields as soon as the grain gaseous randomization time gets larger than that of paramagnetic relaxation. The same conclusion is valid for the mechanical alignment of helical grains. If observations confirm that the degrees of alignment are higher than radiative torques can produce alone, this may be a proof of the presence of superparamagentic inclusions.

Polarimetry is extensively used as a tool to trace the interstellar magnetic field projected on the plane of sky. Moreover, it is also possible to estimate the magnetic field intensity from polarimetric maps based on the Chandrasekhar-Fermi method. In this work, we present results for turbulent, isothermal, 3-D simulations of sub/supersonic and sub/super-Alfvenic cases. With the cubes, assuming perfect grain alignment, we created synthetic polarimetric maps for different orientations of the mean magnetic field with respect to the line of sight (LOS). We show that the dispersion of the polarization angle depends on the angle of the mean magnetic field regarding the LOS and on the Alfvenic Mach number. However, the second order structure function of the polarization angle follows the relation $SF \propto l^{\alpha}$, $\alpha$ being dependent exclusively on the Alfvenic Mach number. The results show an anti-correlation between the polarization degree and the column density, with exponent $\gamma \sim -0.5$, in agreement with observations, which is explained by the increase in the dispersion of the polarization angle along the LOS within denser regions. However, this effect was observed exclusively on supersonic, but sub-Alfvenic, simulations. For the super-Alfvenic, and the subsonic model, the polarization degree showed to be intependent on the column density. Our major quantitative result is a generalized equation for the CF method, which allowed us to determine the magnetic field strength from the polarization maps with errors $< 20%$. We also account for the role of observational resolution on the CF method.

Density is the turbulence statistics that is most readily available from observations. Different regimes of turbulence correspond to different density spectra. For instance, the viscosity-damped regime of MHD turbulence relevant, for instance, to partially ionized gas, can be characterized by shallow and very anisotropic spectrum of density. This spectrum can result in substantial variations of the column densities. Addressing MHD turbulence in the regime when viscosity is not important over the inertial range, we demonstrate with numerical simulations that it is possible to reproduce both the observed Kolmogorov spectrum of density fluctuations observed in ionized gas by measuring scintillations and more shallow spectra that are obtained through the emission measurements. We show that in supersonic turbulence the high density peaks dominate shallow isotropic spectrum, while the small-scale underlying turbulence that fills most of the volume has the Kolmogorov spectrum and demonstrates scale-dependent anisotropy. The limitations of the spectrum in studying turbulence induce searches of alternative statistics. We demonstrate that a measure called "bispectrum" may be a promising tool. Unlike spectrum, the bispectrum preserves the information about wave phases.

Numerical simulations shed light onto earlier not trackable problem of magnetohydrodynamic (MHD) turbulence. They allowed to test the predictions of different models and choose the correct ones. Inevitably, this progress calls for revisions in the picture of cosmic ray (CR) transport. It also shed light on the problems with the present day numerical modeling of CR. In this paper we focus on the analytical way of describing CR propagation and scattering, which should be used in synergy with the numerical studies. In particular, we use recently established scaling laws for MHD modes to obtain the transport properties for CRs. We include nonlinear effects arising from large scale trapping, to remove the 90 degree divergence. We determine how the efficiency of the scattering and CR mean free path depend on the characteristics of ionized media, e.g. plasma $\beta$, Coulomb collisional mean free path. Implications for particle transport in interstellar medium and solar corona are discussed. We also examine the perpendicular transport of CRs. Perpendicular transport depends on the comparison of parallel mean free path and the injection scale of the turbulence, as well as the Alfv\'enic Mach number. Normal turbulence does not allow subdiffusion unless there are slab waves. The critical scale below which subdiffusion applies is provided. These results can be used to compare with the numerical calculations, provided that these calculations use the structure of magnetic field which is consistent with the numerical studies of MHD turbulence.

Multi-wavelength imaging polarimetry at far-infrared wavelengths has proven to be an excellent tool for studying the physical properties of dust, molecular clouds, and magnetic fields in the interstellar medium. Although these wavelengths are only observable from airborne or space-based platforms, no first-generation instrument for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is presently designed with polarimetric capabilities. We study several options for upgrading the High-resolution Airborne Wideband Camera (HAWC) to a sensitive FIR polarimeter. HAWC is a 12 x 32 pixel bolometer camera designed to cover the 53 - 215 micron spectral range in 4 colors, all at diffraction-limited resolution (5 - 21 arcsec). Upgrade options include: (1) an external set of optics which modulates the polarization state of the incoming radiation before entering the cryostat window; (2) internal polarizing optics; and (3) a replacement of the current detector array with two state-of-the-art superconducting bolometer arrays, an upgrade of the HAWC camera as well as polarimeter. We discuss a range of science studies which will be possible with these upgrades including magnetic fields in star-forming regions and galaxies and the wavelength-dependence of polarization.

A detailed study of interstellar polarization efficiency toward molecular clouds is used to attempt discrimination between grain alignment mechanisms in dense regions of the ISM. Background field stars are used to probe polarization efficiency in quiescent regions of dark clouds, yielding a dependence on visual extinction well-represented by a power law. No significant change in this behavior is observed in the transition region between the diffuse outer layers and dense inner regions of clouds, where icy mantles are formed, and we conclude that mantle formation has little or no effect on the efficiency of grain alignment. Young stellar objects generally exhibit greater polarization efficiency compared with field stars at comparable extinctions, displaying enhancements by factors of up to 6. Of the proposed alignment mechanisms, that based on radiative torques appears best able to explain the data. The attenuated external radiation field accounts for the observed polarization in quiescent regions, and radiation from the embedded stars themselves may enhance alignment in the lines of sight to YSOs. Enhancements in polarization efficiency observed in the ice features toward several YSOs are of greatest significance, as they demonstrate efficient alignment in cold molecular clouds associated with star formation.

By applying recent results for the slab correlation time scale onto cosmic ray scattering theory, we compute cosmic ray parallel mean free paths within the quasilinear limit. By employing these results onto charged particle transport in the solar system, we demonstrate that much larger parallel mean free paths can be obtained in comparison to previous results. A comparison with solar wind observations is also presented to show that the new theoretical results are much closer to the observations than the previous results.

We consider the propagation of cosmic rays in turbulent magnetic fields. We use the models of magnetohydrodynamic turbulence that were tested in numerical simulations, in which the turbulence is injected on large scale and cascades to small scales. Our attention is focused on the models of the strong turbulence, but we also briefly discuss the effects that the weak turbulence and the slab Alfv\'enic perturbations can have. The latter are likely to emerge as a result of instabilities with in the cosmic ray fluid itself, e.g., beaming and gyroresonance instabilities of cosmic rays. To describe the interaction of cosmic rays with magnetic perturbations we develop a non-linear formalism that extends the ordinary Quasi-Linear Theory (QLT) that is routinely used for the purpose. This allows us to avoid the usual problem of 90 degree scattering and enable our computation of the mean free path of cosmic rays. We apply the formalism to the cosmic ray propagation in the galactic halo and in the Warm Ionized medium (WIM). In addition, we address the issue of the transport of cosmic rays perpendicular to the mean magnetic field and show that the issue of cosmic ray subdiffusion (i.e., propagation with retracing the trajectories backwards, which slows down the diffusion) is only important for restricted cases when the ambient turbulence is far from what numerical simulations suggest to us. As a result, this work provides formalism that can be applied for calculating cosmic ray propagation in a wide variety of circumstances.

We study the physical processes that affect the alignment of grains subject to radiative torques (RATs). To describe the action of RATs, we use the analytical model (AMO) of RATs introduced in Paper I. We focus our discussion on the alignment by anisotropic radiation flux with respect to magnetic field, which defines the axis of grain Larmor precession. Such an alignment does not invoke paramagnetic dissipation (i.e. Davis-Greenstein mechanism), but, nevertheless, grains tend to be aligned with long axes perpendicular to the magnetic field. When we account for thermal fluctuations within grain material, we show that for grains, which are characterized by a triaxial ellipsoid of inertia, the zero-$J$ attractor point obtained in our earlier study develops into a low-$J$ attractor point. We study effects of stochastic gaseous bombardment and show that gaseous bombardment can drive grains from low-$J$ to high-$J$ attractor points in cases when the high-$J$ attractor points are present. As the alignment of grain axes with respect to angular momentum is higher for higher values of $J$, counter-intuitively, gaseous bombardment can increase the degree of grain alignment in respect to the magnetic field. We also study the effects of torques induced by H$_2$ formation and show that they can change the value of angular momentum at high-$J$ attractor point, but marginally affect the value of angular momentum at low-$J$ attractor points. We compare the AMO results with those obtained using the direct numerical calculations of RATs acting upon irregular grains and validate the use of the AMO for realistic situations of RAT alignment.

We show that grains can be efficiently aligned by interacting with a subsonic gaseous flow. The alignment arises from grains having irregularities that scatter atoms with different efficiency in the right and left directions. The grains tend to align with long axes perpendicular to magnetic field, which corresponds to Davis-Greenstein predictions. Choosing conservative estimates, scattering efficiency of impinging atoms and conservative ``degree of helicity'', the alignment of helical grains is much more efficient than the Gold-type alignment processes.

We study scaling relations of compressible isothermal strongly magnetized turbulence using numerical simulations with resolution 512$^3$. We find a good correspondence of our results with the Fleck (1996) model of compressible hydrodynamic turbulence. In particular, we find that the density-weighted velocity, i.e. $\boldsymbol{u} \equiv \rho^{1/3} \boldsymbol{v}$, proposed in Kritsuk et al. (2007) obeys the Kolmogorov scaling, i.e. ${\cal E}_{u}(k)\sim k^{-5/3}$ for the high Mach number turbulence. Similarly, we find that the exponents of the third order structure functions for $\boldsymbol{u}$ stay equal to unity for all Mach numbers studied. The scaling of higher order correlations obeys the She-L\'{e}v\^{e}que (1994) scalings corresponding to the two-dimensional dissipative structures, and this result does not change with the Mach number either. In contrast to velocity $\boldsymbol{v}$ which exhibits different scaling parallel and perpendicular to the local magnetic field, the scaling of $\boldsymbol{u}$ is similar in both directions. In addition, we find that the peaks of density create a hierarchy in which both physical and column densities decrease with the scale in accordance to the Fleck (1996) predictions. This hierarchy can be related ubiquitous small ionized and neutral structures (SINS) in the interstellar gas. We believe that studies of statistics of the column density peaks can provide both consistency check for the turbulence velocity studies and insight into supersonic turbulence, when the velocity information is not available.

The structure of magnetic fields within protostellar disks may be studied via polarimetry provided that grains are aligned in respect to magnetic field within the disks. We explore alignment of dust grains by radiative torque in T Tauri disks and provide predictions for polarized emission for disks viewed at different wavelengths and viewing angles. We show that the alignment is especially efficient in outer part of the disks. In the presence of magnetic field, these aligned grains produce polarized emission in infrared wavelengths. We consider a simple model of an accretion disk and provide predictions for polarization that should be available to both instruments that do not resolve the disks and future instruments that will resolve the disks. As the surface magnetic field and the bulk magnetic field play different roles for the disk dynamics, we consider separately the contributions that arises from the surface areas of the disk and its interior. We find that the polarized emission drops for wavelengths shorter than $\sim 10 \mu m$. Between $\sim 10 \mu m$ and $\sim 100 \mu m$, the polarized emission is dominated by the emission from the surface layer of the disks and the degree of polarization can be as large as $\sim 10%$ for unresolved disks. The degree of polarization is around 2-3 % level at wavelengths larger than $\sim100\mu m$.

We attempt to get a physical insight into grain alignment processes by studying basic properties of radiative torques (RATs). For this purpose we consider a simple toy model of a helical grain that reproduces well the basic features of RATs. The model grain consists of a spheroidal body with a mirror attached at an angle to it. Being very simple, the model allows analytical description of RATs that act upon it. We show a good correspondence of RATs obtained for this model and those of irregular grains calculated by DDSCAT. Our analysis of the role of different torque components for grain alignment reveals that one of the three RAT components does not affect the alignment, but induces only for grain precession. The other two components provide a generic alignment with grain long axes perpendicular to the radiation direction, if the radiation dominates the grain precession, and perpendicular to magnetic field, otherwise. We study a self-similar scaling of RATs as a function of $\lambda/a_{eff}$. We show that the self-similarity is useful for studying grain alignment by a broad spectrum of radiation, i.e. interstellar radiation field.

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. One of the easiest ways to determine the magnetic field direction is via polarization of radiation resulting from extinction on or/and emission by aligned dust grains. Reliability of interpretation of the polarization maps in terms of magnetic fields depends on how well we understand the grain-alignment theory. Currently the grain-alignment theory is a predictive one, and its results nicely match observations. Among its predictions is a subtle phenomenon of radiative torques. after having stayed in oblivion for many years after its discovery, is currently viewed as the most powerful means of alignment. In this article, I shall review the basic physical processes involved in grain alignment, and the currently known mechanisms of alignment. I shall also discuss possible niches for different alignment mechanisms. I shall dwell on the importance of the concept of grain helicity for understanding of many properties of grain alignment, and shall demonstrate that rather arbitrarily shaped grains exhibit helicity when they interact with gaseous flows and radiative fluxes.

We study of the effect of turbulence on diffusion processes within magnetized medium. While we exemplify our treatment with heat transfer processes, our results are quite general and are applicable to different processes, e.g. diffusion of heavy elements. Our treatment is also applicable to describing the diffusion of cosmic rays arising from magnetic field wandering. In particular, we find that when the energy injection velocity is smaller than the Alfven speed the heat transfer is partially suppressed, while in the opposite regime the effects of turbulence depend on the intensity of driving. In fact, the scale $l_A$ at which the turbulent velocity is equal the Alfven velocity is a new important parameter. When the electron mean free path $\lambda$ is larger than $l_A$, the stronger the the turbulence, the lower thermal conductivity by electrons is. The turbulent motions, however, induces their own advective transport, that can provide effective diffusivity. For clusters of galaxies, we find that the turbulence is the most important agent for heat transfer. We also show that the domain of applicability of the subdiffusion concept is rather limited.

We study scaling relations of compressible strongly magnetized turbulence. We find a good correspondence of our results with the Fleck (1996) model of compressible hydrodynamic turbulence. In particular, we find that the density-weighted velocity, i.e. $u \equiv \rho^{1/3} v$, proposed in Kritsuk et al. (2007) obeys the Kolmogorov scaling, i.e. $E_{u}(k)\sim k^{-5/3}$ for the high Mach number turbulence. Similarly, we find that the exponents of the third order structure functions for $u$ stay equal to unity for the all the Mach numbers studied. The scaling of higher order correlations obeys the She-Leveque (1994) scalings corresponding to the two-dimensional dissipative structures, and this result does not change with the Mach number either. In contrast to $v$ which exhibits different scaling parallel and perpendicular to the local magnetic field, the scaling of $u$ is similar in both directions. In addition, we find that the peaks of density create a hierarchy in which both physical and column densities decrease with the scale in accordance to the Fleck (1996) predictions. This hierarchy can be related ubiquitous small ionized and neutral structures (SINS) in the interstellar gas. We believe that studies of statistics of the column density peaks can provide both consistency check for the turbulence velocity studies and insight into supersonic turbulence, when the velocity information is not available.

We attempt to explain the non-thermal emission arising from galaxy clusters as a result of the re-acceleration of electrons by compressible turbulence induced by cluster mergers. In our model intracluster medium (ICM) is represented by a high beta plasma in which turbulent motions are driven at large scales. The corresponding injection velocities are higher than the Alfven velocity. As a result, the turbulence is approximately isotropic up to the scale at which the turbulent velocity gets comparable with the Alfven velocity. Under the hypothesis that turbulence in the ICM is highly super- Alfvenic the magnetic field is passively advected and the field lines are bended on scales smaller than that of the classical, unmagnetized, ion-ion mean free path. This affects ion diffusion and the strength of the effective viscosity. Under these conditions the bulk of turbulence in hot (5-10 keV temperature) galaxy clusters is likely to be dissipated at collisionless scales via resonant coupling with thermal and fast particles. We use collisionless physics to derive the amplitude of the different components of the energy of the compressible modes, and review and extend the treatment of plasma damping in the ICM. We calculate the acceleration of both protons and electrons taking into account both TTD acceleration and non-resonant acceleration by large scale compressions. We find that relativistic electrons can be re-accelerated in the ICM up to energies of several GeV provided that the rms velocity of the compressible turbulent-eddies is (V_L/c_s)^2~0.15-0.3. We find that under typical conditions ~ 2-5 % of the energy flux of the cascading of compressible motions injected at large scales goes into the acceleration of fast particles and that this may explain the observed non-thermal emission from merging galaxy clusters.

Comet C/1999 S4 (LINEAR) was exceptional in many respects. Its nucleus underwent multiple fragmentations culminating in the complete disruption around July 20, 2000. We present circular polarization measurements along the cuts through the coma and nucleus of the comet during three separate observing runs, in June 28 - July 2, July 8 - 9, and July 21 - 22, 2000. The circular polarization was detected at a rather high level, up to 0.8%. The left-handed as well as right-handed polarization was observed over the coma with the left circularly polarized light systematically observed in the sunward part of the coma. During our observations the phase angle of the comet varied from 61 up to 122 deg., which allowed us to reveal variations of circular polarization with the phase angle. Correlation between the degree of circular polarization, visual magnitude, water production rate, and linear polarization of Comet C/1999 S4 (LINEAR) during its final fragmentation in July 2000 was found. The mechanisms that may produce circular polarization in comets and specifically in Comet C/1999 S4 (LINEAR) are discussed and some tentative interpretation is presented.

We perform three-dimensional (3D) compressible MHD simulations over many dynamical times for an extended range of sonic and Alfven Mach numbers and analyze the statistics of 3D density and 2D column density, which include probability distribution functions, spectra, skewness, kurtosis, She-Leveque exponents, and genus. In order to establish the relation between the statistics of the observables, i.e. column densities, and the underlying 3D statistics of density, we analyze the effects of cloud boundaries. We define the parameter space for 3D measures to be recovered from column densities. In addition, we show that for subsonic turbulence the spectra of density fluctuations are consistent with k^-7/3 in the case of a strong magnetic field and k^-5/3 in the case of a weak magnetic field. For supersonic turbulence we confirm the earlier findings of the shallow spectra of density and Kolmogorov spectra of the logarithm of density. We find that the intermittencies of the density and velocity are very different.

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 and this motivates studies of turbulence using precision spectroscopy. We provide high resolution numerical testing of the two promising techniques, namely, Velocity Channel Analysis and Velocity Coordinate Spectrum. We obtain an expression for the shot noise that the discretization of the numerical data entails and successfully test it. We show that numerical resolution required for recovering the underlying turbulent spectrum from observations depend on the spectral index of velocity fluctuations. Thus the low resolution testing may be misleading.

We analyze a technique of obtaining turbulence power spectrum using spectral line data along the velocity coordinate, which we refer to as Velocity Coordinate Spectrum (VCS). We formalize geometrical aspects of observation through a single factor, "geometric term". We find that all variety of particular observational configurations can be described using correspondent variants of this term, which we explicitly calculate. This allows us to obtain asymptotics for both parallel lines of sight and crossing lines of sight. The latter case is especially important for studies of turbulence within diffuse ISM in Milky Way. For verification of our results, we use direct calculation of VCS spectra, while the numerical simulations are presented in a companion paper.

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.

We show that atomic alignment presents a reliable way to study topology of astrophysical magnetic fields. The effect of atomic alignment arises from modulation of the relative population of the sublevels of atomic ground state pumped by anisotropic radiation flux. As such aligned atoms precess in the external magnetic field and this affects the properties of the polarized radiation arising from both scattering and absorption by the atoms. As the result the polarizations of emission and absorption lines depend on the 3D geometry of the magnetic field as well as the direction and anisotropy of incident radiation. We consider a subset of astrophysically important atoms with hyperfine structure. For emission lines we obtain the dependencies of the direction of linear polarization on the directions of magnetic field and the incident pumping radiation. For absorption lines we establish when the polarization is perpendicular and parallel to magnetic field. For both emission and absorption lines we find the dependence on the degree of polarization on the 3D geometry of magnetic field. We claim that atomic alignment provides a unique tool to study magnetic fields in circumstellar regions, AGN, interplanetary and interstellar medium. This tool allows studying of 3D topology of magnetic fields and establish other important astrophysical parameters. We consider polarization arising from both atoms in the steady state and also as they undergo individual scattering of photons. We exemplify the utility of atomic alignment for studies of astrophysical magnetic fields by considering a case of Na alignment in a comet wake.

Turbulence in an incompressible fluid with and without a magnetic field as well as moderately compressible MHD turbulence are compared. For the magnetohydrodynamic (MHD) models the probability distribution functions of the velocity components perpendicular to the external magnetic field are like the incompressible hydrodynamic (HD) model while those parallel to the field have a smaller range of velocities. The probability distribution functions of the transverse velocity increments for the MHD models decline slower than the incompressible HD model. The similarity of incompressible HD and both incompressible and compressible MHD turbulence persists over high order longitudinal structure function scaling exponents measured in the global reference frame as well as for motions perpendicular to the local mean field. In these two frames the longitudinal scaling exponents of both MHD models seem to follow theoretical incompressible HD dissipation structure predictions while the transverse scaling exponents of the incompressible MHD model seem to follow the predictions for incompressible MHD. In the local magnetic system the motions parallel to the local mean field for both MHD models are different from incompressible HD motions. The fields of the MHD simulations are decomposed into Alfvenic, fast, and slow modes and stuided separately. (abridged)

We report the first attempt to mode decomposition of the velocity field using wavelet transforms in the application to the magnetized turbulence analysis. We compare results of the separation with respect to the global and local mean magnetic field for models with a strong and weak external magnetic field. We show, that for models with a strong external magnetic field our results are relatively consistent for both methods. However, for superAlfv\'{e}nic turbulence with spatially variable direction of the local mean magnetic field the differences are substantial.

Grain alignment theory has reached the stage where quantitative predictions of the degree of alignment and its variations with optical depth are possible. With the goal of studying the effect of clumpiness on the sub-millimeter and far infrared polarization we computed the polarization due to alignment via radiative torques within clumpy models of cores and molecular clouds. Our models were based upon a highly inhomogeneous simulation of compressible magnetohydrodynamic turbulence. The P-I relations for our models reproduce those seen in observations. We show that the degree of polarization observed is extremely sensitive to the upper grain size cut-off, and is less sensitive to changes in the radiative anisotropy. Furthermore, despite a variety of dust temperatures along a single line of sight through our core and amongst dust grains of different sizes, the assumption of isothermality amongst the aligned grains does not introduce a significant error. Our calculations indicate that sub-mm polarization vectors can be reasonably good tracers for the underlying magnetic field structure, even for relatively dense clouds (A_V~10 to the cloud center). The current predictive power of the grain alignment theory should motivate future polarization observations using the next generation of multi-wavelength sub-mm polarimeters such as those proposed for SOFIA.

The relative population of the fine structure sublevels of an atom's ground state is affected by radiative transitions induced by an anisotropic radiation flux. This causes the alignment of atomic angular momentum. In terms of observational consequences for the interstellar and intergalactic medium, this results in the polarization of the absorption lines. In the paper we consider the conditions necessary for this effect and provide calculations of polarization from a few astrophysically important atoms and ions with multiple upper and lower levels for an arbitrary orientation of magnetic fields to the a) source of optical pumping, b) direction of observation, c) absorbed source. We also consider an astrophysically important ``degenerate'' case when the source of optical pumping coincides with the source of the absorbed radiation. We present analytical expressions that relate the degree of linear polarization and the intensity of absorption to the 3D orientation of the magnetic field with respect to the pumping source, the source of the absorbed radiation, and the direction of observations. We discuss how all these parameters can be determined via simultaneous observations of several absorption lines and suggest graphical means that are helpful in practical data interpretation. We prove that studies of absorption line polarization provide a unique tool to study 3D magnetic field topology in various astrophysical conditions.

The problems with explaining the Small Ionized and Neutral Structures (SINS) appealing to turbulence stem from inefficiency of the Kolmogorov cascade in creating large fluctuations at sufficiently small scales. However, other types of cascades are possible. When magnetic turbulence in a fluid with viscosity that is much larger than resistivity gets to a viscous damping scale, the turbulence does not vanish. Instead, it gets into a different new regime. Viscosity-damped turbulence produces fluctuations on the small scales. Magnetic fields sheared by turbulent motions by eddies not damped by turbulence create small scale filaments that are confined by the external plasma pressure. This creates small scale density fluctuations. In addition, extended current sheets create even stronger density gradients that accompany field reversals in the plane perpendicular to mean magnetic field. Those can be responsible for the SINS formation. This scenario is applicable to partially ionized gas. More studies of reconnection in the viscosity dominated regime are necessary to understand better the extend to which the magnetic reversals can compress the gas.

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.

It is most natural to assume that mysterious Small Ionized and Neutral Structures (SINS) indiffuse ISM arise from turbulence. There are two obvious problem with such an explanation, however. First of all, it is generally believed that at the small scales turbulence should be damped. Second, turbulence with Kolmogorov spectrum cannot be the responsible for the SINS. We consider, however, effects, that provide spectral index flatter than the Kolmogorov one and allow action at very small scales. These are the shocks that arise in high Mach number turbulence and transfer of energy to small scales by instabilities in cosmic rays. Our examples indicate that the origin of SINS may be discovered through systematic studies of astrophysical turbulence.

We study the scattering of low-energy Cosmic Rays (CRs) in a turbulent, compressive MHD fluid. We show that compressible MHD modes -- fast or slow waves with wave lengths smaller than CR mean free paths induce cyclotron instability in CRs. The instability feeds the new small-scale Alfvenic wave component with wave vectors mostly along magnetic field, which is not a part of the MHD turbulence cascade. This new component gives feedback on the instability through decreasing the CR mean free path. We show that the ambient turbulence fully suppresses the instability at large scales, while wave steepening constrains the amplitude of the waves at small scales. We provide the energy spectrum of the plane-parallel Alfvenic component and calculate mean free paths of CRs as a function of their energy. We find that for the typical parameters of turbulence in the interstellar medium and in the intercluster medium the new Alfvenic component provides the scattering of the low energy CRs that exceeds the direct resonance scattering by MHD modes. This solves the problem of insufficient scattering of low-energy CRs in the turbulent interstellar or intracluster medium that was reported in the literature.

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.

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.

Aligned non-spherical dust particles polarize starlight passing through the dust cloud. They also emit polarized far infrared and sub-mm radiation. Substantial progress in understanding of grain alignment theory makes the interpretation of the polarized radiation in terms of underlying magnetic fields much more reliable. I discuss a number of fundamental processes that affect grain alignment. In particular, I shall discuss how subtle effects related to nuclear spins of the atoms alter the dynamics of dust grains. I shall discuss how the theory explains the existing observational data and demonstrate when the polarization can and cannot be interpreted in terms of the underlying magnetic fields.

Intercluster medium is expected to be turbulent with turbulence being superAlfvenic at large scales. Magnetic fields substantially modify the turbulent cascade when the turbulence reaches the scales at which the fluctuation velocity gets less than the Alfven velocity. At those scales it is possible to consider three cascades, of fast, slow and Alfven modes with little energy exchange between them. As Alfvenic and slow modes are anisotropic they marginally scatter and accelerate cosmic rays, while fast modes dominate the processes. However, in the presence of cosmic rays the turbulence is modified as cosmic rays transfer the energy of compressible motions (i.e. slow and fast modes) from large scales to the scale of cosmic ray Larmor radius. This results in generation of a new small-scale Alfvenic component which is not a part of the ordinary MHD cascade. This component does scatter and accelerate cosmic rays. In addition, magnetic reconnection in turbulent medium accelerates cosmic rays. The complexity of the intracluster turbulence calls for observational studies. A new technique Velocity Coordinate Spectrum (VCS) is particularly promising for studies of velocity fluctuations with a new generation of X-ray observatories.

We discuss a new technique for studying astrophysical turbulence that utilizes the statistics of Doppler-broadened spectral lines. The technique relates the power Velocity Coordinate Spectrum (VCS), i.e. the spectrum of fluctuations measured along the velocity axis in Position-Position-Velocity (PPV) data cubes available from observations, to the underlying power spectra of the velocity/density fluctuations. Unlike the standard spatial spectra, that are function of angular wavenumber, the VCS is a function of the velocity wave number k_v ~ 1/v. We show that absorption affects the VCS to a higher degree for small k_v and obtain the criteria for disregarding the absorption effects for turbulence studies at large k_v. We consider the retrieval of turbulence spectra from observations for high and low spatial resolution observations and find that the VCS allows one to study turbulence even when the emitting turbulent volume is not spatially resolved. This opens interesting prospects for using the technique for extragalactic research. We show that, while thermal broadening interferes with the turbulence studies using the VCS, it is possible to separate thermal and non-thermal contributions. This allows a new way of determining the temperature of the interstellar gas using emission and absorption spectral lines.

We present models of turbulent mixing at the boundaries between hot (T~10^{6-7} K) and warm material (T~10^4 K) in the interstellar medium, using a three-dimensional magnetohydrodynamical code, with radiative cooling. The source of turbulence in our simulations is a Kelvin-Helmholtz instability, produced by shear between the two media. We found, that because the growth rate of the large scale modes in the instability is rather slow, it takes a significant amount of time (~1 Myr) for turbulence to produce effective mixing. We find that the total column densities of the highly ionized species (C IV, N V, and O VI) per interface (assuming ionization equilibrium) are similar to previous steady-state non-equilibrium ionization models, but grow slowly from log N ~10^{11} to a few 10^{12} cm^{-2} as the interface evolves. However, the column density ratios can differ significantly from previous estimates, with an order of magnitude variation in N(C IV)/N(O VI) as the mixing develops.

We propose collisionless damping of fast MHD waves as an important mechanism for the heating and acceleration of winds from rotating stars. Stellar rotation causes magnetic field lines anchored at the surface to form a spiral pattern and magneto-rotational winds can be driven. If the structure is a magnetically dominated, fast MHD waves generated at the surface can propagate almost radially outward and cross the field lines. The propagating waves undergo collisionless damping owing to interactions with particles surfing on magnetic mirrors that are formed by the waves themselves. The damping is especially effective where the angle between the wave propagation and the field lines becomes moderately large ($\sim 20$ to $80^{\circ}$). The angle tends naturally to increase into this range because the field in magneto-rotational winds develops an increasingly large azimuthal component. The dissipation of the wave energy produces heating and acceleration of the outflow. We show using specified wind structures that this damping process can be important in both solar-type stars and massive stars that have moderately large rotation rates. This mechanism can play a role in coronae of young solar-type stars which are rapidly rotating and show X-ray luminosities much larger than the sun. The mechanism could also be important for producing the extended X-ray emitting regions inferred to exist in massive stars of spectral type middle B and later.

The investigation of the statistical properties of maps of line centroids has been used for almost 50 years, but there is still no general agreement on their interpretation. We try to quantify which properties of underlying turbulent velocity fields can be derived from centroid velocity maps, and we test conditions under which the scaling behaviour of the centroid velocities matches the scaling of the three-dimensional velocity field. Using fractal cloud models we study systematically the relation between three-dimensional density and velocity fields and the statistical properties of the produced line centroid maps. We put special attention to cases with large density fluctuations resembling supersonic interstellar turbulence. Starting from the Delta-variance analysis we derive a new tool to compute the scaling behaviour of the three-dimensional velocity field from observed intensity and centroid velocity maps. We provide two criteria to decide whether the information from the centroid velocities directly reflects the properties of the underlying velocity field. Applying these criteria allows to understand the different results found so far in the literature on the interpretation of the statistics of velocity centroids. A new iteration scheme can be used to derive the three-dimensional velocity scaling from centroid velocity maps for arbitrary density and velocity fields, but it requires an accurate knowledge of the average density of the considered interstellar cloud.

Goldreich-Sridhar model of incompressible turbulence provides an elegant approach to describing strong MHD turbulence. It relies on the fact that interacting Alfvenic waves are independent and have random polarization. However, in case of strong interaction a spontaneous local assymetry can arise. We used direct numerical simulations to show that polarization alignment occurs and it grows larger at smaller scales. This effect would lead to a shallower spectrum and stronger anisotropy. Even small changes in these two properties will have important astrophysical consequencies, e.g. for the cosmic ray physics.

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-$\beta$) and gaseous pressure (high-$\beta$) 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.

(Abridged) We describe the cascade of plasma waves or turbulence injected, presumably by reconnection, at scales comparable to the size of a solar flare loop to scales comparable to particle gyroradii, and evaluate their damping by various mechanisms. We show that the classical viscous damping is unimportant for magnetically dominated or low beta plasmas and the primary damping mechanism is the collisionless damping by the background particles. We show that the damping rate is proportional to the total random momentum density of the particles. For solar flare conditions this means that in most flares, except the very large ones, the damping is dominated by thermal background electrons. For large flares one requires acceleration of essentially all background electrons into a nonthermal distribution so that the accelerated electrons can be important in the damping of the waves. In general, damping by thermal or nonthermal protons is negligible compared to that of electrons except for quasi-perpendicular propagating waves or for rare proton dominated flares with strong nuclear gamma-ray line emission. Using the rate for damping we determine the critical scale below which the damping becomes important and the spectrum of the turbulence steepens. This critical scale, however, has strong dependence on the angle of propagation with respect to the magnetic field direction. The waves can cascade down to very small scales, such as the gyroradii of the particles at small angles (quasi-parallel propagation) and possibly near 90 degree (quasi-perpendicular propagation) giving rise to a highly anisotropic spectral distribution.

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 $\beta \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\'en 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^{-\alpha_E}$, with $\alpha_E = 5/2, 2$, and corresponding synchrotron radio power-law spectra with spectral indices which are compatible with that observed during the flares ($ S_{\nu} \propto \nu^{-0.75, - 0.5}$).

We study alignment of grains by radiative torques. We found steep rise of radiative torque efficiency as grain size increases. This allows larger grains that are known to exist within molecular clouds to be aligned by the attenuated and reddened interstellar radiation field. In particular, we found that, even deep inside giant molecular clouds, e.g. at optical depths corresponding to less than Av of 10 large grains can still be aligned by radiative torque. This means that, contrary to earlier claims, far-infrared/submillimeter polarimetry provides a reliable tool to study magnetic fields of pre-stellar cores. Our results show that the grain size distribution is important for determining the relation between the degree of polarization and intensity.

In this presentation we review the link between the statistics of intensity fluctuations in spectral line data cubes with underlying statistical properties of turbulence in the interstellar medium. Both the formalism of Velocity Channel Analysis for optically thin lines and its extension to the lines with self-absorption is described. We demonstrate that by observing optically thin lines from cold gas in sufficiently narrow (thin) velocity channels one may recover the scaling of the stochastic velocities from turbulent cascade, in particular, Kolmogorov velocities give $K^{-2.7}$ contribution to the intensity power spectrum. Synthetically increasing the channel thickness separates out the underlying density inhomogeneities of the gas. Effects of self absorption, on the other hand, retain the velocity signature even for integrated lines. As a result, intensity fluctuations tend to show universal but featureless scaling of the power $\propto K^{-3}$ over the range of scales.

We use the results of magnetohydrodynamic (MHD) simulations to emulate spectroscopic observations and use maps of centroids to study their statistics. In order to assess under which circumstances the scaling properties of the velocity field can be retrieved from velocity centroids, we compare two point statistics (structure functions and power-spectra) of velocity centroids with those of the underlying velocity field and analytic predictions presented in a previous paper (Lazarian & Esquivel 2003). We tested a criterion for recovering velocity spectral index from velocity centroids derived in our previous work, and propose an approximation of the early criterion using only the variances of ``unnormalized'' velocity centroids and column density maps. It was found that both criteria are necessary, however not sufficient to determine if the centroids recover velocity statistics. Both criteria are well fulfilled for subsonic turbulence. We find that for supersonic turbulence with sonic Mach numbers > 2.5 centroids fail to trace the spectral index of velocity. Asymptotically, however, we claim that recovery of velocity statistics is always possible provided that the density spectrum is steep and the observed inertial range is sufficiently extended. In addition, we show that velocity centroids are useful for anisotropy studies and determining the direction of magnetic field, even if the turbulence is highly supersonic, but only if it is sub-Alfvenic. This provides a tool for mapping the magnetic field direction, and testing whether the turbulence is sub-Alfvenic or super-Alfvenic.

Most astrophysical fluids are turbulent and magnetized. Fluctuations of polarization provide a promising way to study astrophysical magnetic turbulence. We discuss polarization that arises from grains and atoms aligned in respect to magnetic field, describe the processes of alignment, explain when the alignment reflects the direction of magnetic field. We show that temporal fluctuations of interplanetary magnetic field can be studied using aligned sodium atoms.

Turbulent reconnection allows fast magnetic reconnection of astrophysical magnetic fields. This entails numerous astrophysical implications and opens new ways to approach long standing problems. I briefly discuss a model of turbulent reconnection within which the stochasticity of 3D magnetic field enables rapid reconnection through both allowing multiple reconnection events to take place simultaneously and by restricting the extension of current sheets. In fully ionized gas the model in Lazarian and Vishniac 99 predicts reconnection rates that depend only on the intensity of turbulence. In partially ionized gas a modification of the original model in Lazarian, Vishniac and Cho 04 predicts the reconnection rates that, apart from the turbulence intensity depend on the degree of ionization. In both cases the reconnection may be slow and fast depending on the level of turbulence in the system. As the result, the reconnection gets bursty, which provides a possible explanation to Solar flares and possibly to gamma ray busts. The implications of the turbulent reconnection model have not been yet studied in sufficient detail. I discuss first order Fermi acceleration of cosmic ray that takes place as the oppositely directed magnetic fluxes move together. This acceleration would work in conjunction with the second order Fermi acceleration that is caused by turbulence in the reconnection region. In partially ionized gas the stochastic reconnection enables fast removal of magnetic flux from star forming molecular clouds.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD modes to obtain the scattering frequency for cosmic rays. We account for the turbulence cutoff arising from both collisional and collisionless damping. We obtain the scattering rate and show that fast modes provide the dominant contribution to cosmic ray scattering for the typical interstellar conditions in spite of the fact that fast modes are subjected to damping. We determine how the efficiency of the scattering depends on the characteristics of ionized media, e.g. plasma $\beta$. We show that streaming instability is suppressed by the ambient MHD turbulence.

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 Velocity Centroids are dominated by density or velocity. We briefly discuss ways of going beyond power spectra by using of higher order correlations as well as genus analysis. We outline the relation between Spectral Correlation Functions and the statistics available through VCA and VCS.

We summarise basic properties of MHD turbulence. First, MHD turbulence is not so messy as it is believed. In fact, the notion of strong non-linear coupling of compressible and incompressible motions along MHD cascade is not tenable. Alfven, slow and fast modes of MHD turbulence follow their own cascades and exhibit degrees of anisotropy consistent with theoretical expectations. Second, the fast decay of turbulence is not related to the compressibility of fluid. Rates of decay of compressible and incompressible motions are very similar. Third, the properties of Alfven and slow modes are similar to their counterparts in the incompressible MHD. The properties of fast modes are similar to accoustic turbulence, which does require more studies. Fourth, the density at low Mach numbers and logarithm of density at higher Mach numbers exhibit Kolmogorov-type spectrum.

We study the statistics of density for supersonic turbulence in a medium with magnetic pressure larger than the gaseous pressure. This study is motivated by molecular cloud research. Our simulations exhibit clumpy density structures, which contrast increases with the Mach number. At 10 Machs densities of some clumps are three orders of magnitude higher than the mean density. These clumps give rise to flat and approximately isotropic density spectrum corresponding to the random distribution of clumps in space. We claim that the clumps originate from our random, isotropic turbulence driving. When the contribution from those clumps is suppressed by studying logarithm of density, the density statistics exhibit scale-dependent anisotropy consistent with the models where density structures arise from shearing by Alfv\'en waves. It is noteworthy that originally such models were advocated for the case of low-Mach, nearly incompressible turbulence.

We discuss the damping of strong magnetic turbulence downstream of the forward shock of young supernova remnants (SNR). We find that strong magnetic fields, that have been produced by the streaming instability in the upstream region of the shock, or by other kinetic instabilities at the shock, will be efficiently reduced, so the region of enhanced magnetic field strength would typically have a thickness of the order l_d=(10^{16}-10^{17}) cm. The non-thermal X-ray filaments observed in young SNR are thus likely limited by the magnetic field and not by the energy losses of the radiating electrons. Consequently the thickness of the filaments would not be a measure of the magnetic field strength and claims of efficient cosmic-ray acceleration on account of a run-away streaming instability appear premature.

Population of levels of the hyperfine and fine split ground state of an atom is affected by radiative transitions induced by anisotropic radiation flux. Such aligned atoms precess in the external magnetic field and this affects properties of polarized radiation arising from both scattering and absorption by atoms. As the result the degree of light polarization depends on the direction of the magnetic field. This provides a perspective tool for studies of astrophysical magnetic fields using optical and UV polarimetry. We provide calculations for several atoms and ions that can be used to study magnetic fields in interplanetary medium, interstellar medium, circumstellar regions and quasars.

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.

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.

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.

We discuss a few recent developments that are important for understanding of MHD turbulence. First, MHD turbulence is not so messy as it is usually believed. In fact, the notion of strong non-linear coupling of compressible and incompressible motions along MHD cascade is not tenable. Alfven, slow and fast modes of MHD turbulence follow their own cascades and exhibit degrees of anisotropy consistent with theoretical expectations. Second, the fast decay of turbulence is not related to the compressibility of fluid. Rates of decay of compressible and incompressible motions are very similar. Third, viscosity by neutrals does not suppress MHD turbulence in a partially ionized gas. Instead, MHD turbulence develops magnetic cascade at scales below the scale at which neutrals damp ordinary hydrodynamic motions. Forth, density statistics does not exhibit the universality that the velocity and magnetic field do. For instance, at small Mach numbers the density is anisotropic, but it gets isotropic at high Mach numbers. Fifth, the intermittency of magnetic field and velocity are different. Both depend on whether the measurements are done in local system of reference oriented along the local magnetic field or in the global system of reference related to the mean magnetic field.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for revisions in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD modes to obtain the scattering frequency for cosmic rays. Using quasilinear theory we calculate gyroresonance with MHD modes (Alfv\'{e}nic, slow and fast) and transit-time damping (TTD) by fast modes. We provide calculations of cosmic ray scattering for various phases of interstellar medium with realistic interstellar turbulence driving that is consistent with the velocity dispersions observed in diffuse gas. We account for the turbulence cutoff arising from both collisional and collisionless damping. We obtain analytical expressions for diffusion coefficients that enter Fokker-Planck equation describing cosmic ray evolution. We obtain the scattering rate and show that fast modes provide the dominant contribution to cosmic ray scattering for the typical interstellar conditions in spite of the fact that fast modes are subjected to damping. We determine how the efficiency of the scattering depends on the characteristics of ionized media, e.g. plasma $\beta$. We calculate the range of energies for which the streaming instability is suppressed by the ambient MHD turbulence.

We calculate the relative grain-grain motions arising from interstellar magnetohydrodynamic (MHD) turbulence. The MHD turbulence includes both fluid motions and magnetic fluctuations. While the fluid motions accelerate grains through hydro-drag, the electromagnetic fluctuations accelerate grains through resonant interactions. We consider both incompressive (Alfv\'{e}n) and compressive (fast and slow) MHD modes and use descriptions of MHD turbulence obtained in Cho & Lazarian (2002). Calculations of grain relative motion are made for realistic grain charging and interstellar turbulence that is consistent with the velocity dispersions observed in diffuse gas, including cutoff of the turbulence from various damping processes. We show that fast modes dominate grain acceleration, and can drive grains to supersonic velocities. Grains are also scattered by gyroresonance interactions, but the scattering is less important than acceleration for grains moving with sub-Alfv\'{e}nic velocities. Since the grains are preferentially accelerated with large pitch angles, the supersonic grains will be aligned with long axes perpendicular to the magnetic field. We compare grain velocities arising from MHD turbulence with those arising from photoelectric emission, radiation pressure and H$_{2}$ thrust. We show that for typical interstellar conditions turbulence should prevent these mechanisms from segregating small and large grains. Finally, gyroresonant acceleration is bound to preaccelerate grains that are further accelerated in shocks. Grain-grain collisions in the shock may then contribute to the overabundance of refractory elements in the composition of galactic cosmic rays.

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.

We present numerical studies of 3-dimensional electron magnetohydrodynamic (EMHD) turbulence. We investigate cascade timescale and anisotropy of freely decaying strong EMHD turbulence with zero electron skin depth. Cascade time scales with $k^{-4/3}$. Our numerical results clearly show scale-dependent anisotropy. We discuss that the observed anisotropy is consistent with $k_{\|} \propto k_{\perp}^{1/3}$, where $k_{\|}$ and $k_{\perp}$ are wave numbers parallel and perpendicular to (local) mean magnetic field, respectively.

Quantitative description of the statistics of intensity fluctuations within spectral line data cubes introduced in our earlier work is extended to the absorbing media. A possibility of extracting 3D velocity and density statistics from both integrated line intensity as well as from the individual channel maps is analyzed. We find that absorption enables the velocity effects to be seen even if the spectral line is integrated over frequencies. This regime that is frequently employed in observations is characterized by a non-trivial relation between the spectral index of velocities and the spectral index of intensity fluctuations. For instance when density is dominated by fluctuations at large scales, i.e. when correlations scale as r^{-\gamma}, \gamma<0, the intensity fluctuations exhibit a universal spectrum of fluctuations ~K^{-3} over a range of scales. When small scale fluctuations of density contain most of the energy, i.e. when correlations scale as r^{-\gamma}, \gamma>0, the resulting spectrum of the integrated lines depends on the scaling of the underlying density and scales as K^{-3+\gamma}. We show that if we take the spectral line slices that are sufficiently thin we recover our earlier results for thin slice data without absorption. As the result we extend the Velocity Channel Analysis (VCA) technique to optically thick lines enabling studies of turbulence in molecular clouds. In addition, the developed mathematical machinery enables a quantitative approach to solving other problems that involved statistical description of turbulence within emitting and absorbing gas.

Interstellar dust grains exhibit complex dynamics which is essential for understanding many key interstellar processes that involve dust, including grain alignment, grain growth, grain shattering etc. Grain rotational and translational motions are affected not only by gaseous collisions, but also by interactions with ions, photons, magnetic fields etc. Some of those interactions, e.g. interactions of ions with the dipole electric moment of dust grains, require the quantum nature of the process to be accounted for. Similarly, coupling of rotational and vibrational degrees of freedom in a grain happens due to relaxation processes, among which the process related to nuclear spins frequently is the dominant one. This coupling modifies substantially both the dynamics of rotational and translational motions by inducing grain flips. The flipping averages our certain systematic torques and forces that act on grains. As the rate of flipping is larger for smaller grains, these grains can be ``thermally trapped'', i.e. rotate at a thermal rate in spite of the presence of systematic torques. Moreover, a subset of small grains with high dipole moments may rotate at subthermal rates due to high damping arising from grain emission in microwave range of frequencies. Translational dynamics of grains is mostly dominated by grain interactions with magnetohydrodynamic turbulence. Efficient turbulent mixing of dust grains limits the degree to which grains of different sizes may be segregated in space.

Recent advances in understanding of the basic properties of compressible Magnetohydrodynamic (MHD) turbulence call for revisions of some of the generally accepted concepts. First, MHD turbulence is not so messy as it is usually believed. In fact, the notion of strong non-linear coupling of compressible and incompressible motions is not tenable. Alfven, slow and fast modes of MHD turbulence follow their own cascades and exhibit degrees of anisotropy consistent with theoretical expectations. Second, the fast decay of turbulence is not related to the compressibility of fluid. Rates of decay of compressible and incompressible motions are very similar. Third, viscosity by neutrals does not suppress MHD turbulence in a partially ionized gas. Instead, MHD turbulence develops magnetic cascade at scales below the scale at which neutrals damp ordinary hydrodynamic motions. The implications of those changes of MHD turbulence paradigm for molecular clouds require further studies. Those studies can benefit from testing of theoretical predictions using new statistical techniques that utilize spectroscopic data. We briefly discuss advances in development of tools using which the statistics of turbulent velocity can be recovered from observations.

One of the most informative techniques of studying magnetic fields in molecular clouds is based on the use of starlight polarization and polarized emission arising from aligned dust. How reliable the interpretation of the polarization maps in terms of magnetic fields is the issue that the grain alignment theory addresses. I briefly review basic physical processes involved in grain alignment.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for substantial revisions in the picture of cosmic ray transport. In this paper we use recently obtained scaling laws for MHD modes to calculate the scattering frequency for cosmic rays in the ISM. We consider gyroresonance with MHD modes (Alfvenic, slow and fast) and transit-time damping (TTD) by fast modes. We provide calculations of cosmic ray scattering for various phases of interstellar medium with realistic interstellar turbulence driving that is consistent with the velocity dispersions observed in diffuse gas. We account for the turbulence cutoff arising from both collisional and collisionless damping. We obtain analytical expressions for diffusion coefficients that enter Fokker-Planck equation describing cosmic ray evolution. We calculate the scattering rate and parallel spatial diffusion coefficients of cosmic rays for both Alfvenic and fast modes. We conclude that fast modes provides the dominant contribution to cosmic ray scattering for the typical interstellar conditions in spite of the fact that fast modes are subjected to damping. We show that the efficiency of the scattering depends on the plasma beta since it determines the damping of the fast modes. We also show that the streaming instability is modified in the presence of turbulence.

We consider stochastic reconnection in a magnetized, partially ionized medium. Stochastic reconnection is a generic effect, due to field line wandering, in which the speed of reconnection is determined by the ability of ejected plasma to diffuse away from the current sheet along magnetic field lines, rather than by the details of current sheet structure. We consider the limit of weak stochasticity, so that the mean magnetic field energy density is greater than either the turbulent kinetic energy density or the energy density associated with the fluctuating component of the field. We consider field line stochasticity generated through a turbulent cascade, which leads us to consider the effect of neutral drag on the turbulent cascade of energy. In a collisionless plasma, neutral particle viscosity and ion-neutral drag will damp mid-scale turbulent motions, but the power spectrum of the magnetic perturbations extends below the viscous cutoff scale. We give a simple physical picture of the magnetic field structure below this cutoff, consistent with numerical experiments. We provide arguments for the reemergence of the turbulent cascade well below the viscous cut-off scale and derive estimates for field line diffusion on all scales. We note that this explains the persistence of a single power law form for the turbulent power spectrum of the interstellar medium, from scales of tens of parsecs down to thousands of kilometers. We find that under typical conditions in the ISM stochastic reconnection speeds are reduced by the presence of neutrals, but by no more than an order of magnitude.

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.

We compare the properties of ordinary strong magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium with the recently discovered viscosity-damped regime. We focus on energy spectra, anisotropy, and intermittency. Our most surprising conclusion is that in ordinary strong MHD turbulence the velocity and magnetic fields show different high-order structure function scalings. Moreover this scaling depends on whether the intermittency is viewed in a global or local system of reference. This reconciles seemingly contradictory earlier results. On the other hand, the intermittency scaling for viscosity-damped turbulence is very different, and difficult to understand in terms of the usual phenomenological models for intermittency in turbulence. Our remaining results are in reasonable agreement with expectations. First, we find that our high resolution simulations for ordinary MHD turbulence show that the energy spectra are {\it compatible} with a Kolmogorov spectrum, while viscosity-damped turbulence shows a shallow $k^{-1}$ spectrum for the magnetic fluctuations. Second, a new numerical technique confirms that ordinary MHD turbulence exhibits Goldreich-Sridhar type anisotropy, while viscosity-damped MHD turbulence shows extremely anisotropic eddy structures. Finally, we show that many properties of incompressible turbulence for both the ordinary and viscosity-damped regimes carry over to the case of compressible turbulence.

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.

GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse.

We discuss a new mechanism of dust acceleration that acts in a turbulent magnetized medium. The magnetohydrodynamic (MHD) turbulence includes both fluid motions and magnetic fluctuations. We show that while the fluid motions bring about grain motions through the drag, the electromagnetic fluctuations, can accelerate grains through resonant interactions. In this paper we calculate grain acceleration by gyroresonance in Cold Neutral Medium. We consider both incompressible and compressible MHD modes. We show that fast modes dominate grain acceleration. For the parameters chosen, fast modes render supersonic velocities to grains that may shatter grains and enable efficient absorption of heavy elements. Since the grains are preferentially accelerated with large pitch angles, the supersonic grains get aligned with long axes perpendicular to the magnetic field.

We address the problem of studying interstellar turbulence using spectral line data. We find a criterion when the velocity centroids may provide trustworthy velocity statistics. To enhance the scope of centroids applications, we construct a measure that we term ``modified velocity centroids'' (MVCs) and derive an analytical solution that relates the 2D spectra of the modified centroids with the underlying 3D velocity spectrum. We test our results using synthetic maps constructed with data obtained through simulations of compressible magnetohydrodynamical (MHD) turbulence. We show that the modified velocity centroids (MVCs) are complementary to the the Velocity Channel Analysis (VCA) technique. Employed together, they make determining of the velocity spectral index more reliable and for wider variety of astrophysical situations.

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.

The radio - far infrared correlation is one of the tightest correlations found in astronomy. Many of the models explaining this correlation rely on the association of of global magnetic field strength with gas density. In this letter we put forward that the physical reason for this association lies within the processes of magnetohydrodynamic turbulence.

Using numerical methods, we systematically study in the framework of ideal MHD the effect of magnetic fields on heat transfer within a turbulent gas. We measure the rates of passive scalar diffusion within magnetized fluids and make the comparisons a) between MHD and hydro simulations, b) between different MHD runs with different values of the external magnetic field (up to the energy equipartition value), c) between thermal conductivities parallel and perpendicular to magnetic field. We do not find apparent suppression of diffusion rates by the presence of magnetic fields, which implies that magnetic fields do not suppress heat diffusion by turbulent motions.

In a previous work Lazarian and Pogosyan suggested a technique to extract velocity and density statistics, of interstellar turbulence, by means of analysing statistics of spectral line data cubes. In this paper we test that technique, by studying the effect of correlation between velocity and density fields, providing a systematic analysis of the uncertainties arising from the numerics, and exploring the effect of a linear shear. We make use of both compressible MHD simulations and synthetic data to emulate spectroscopic observations and test the technique. With the same synthetic spectroscopic data, we also studied anisotropies of the two point statistics and related those anisotropies with the magnetic field direction. This presents a new technique for magnetic field studies. The results show that the velocity and density spectral indices measured are consistent with the analytical predictions. We identified the dominant source of error with the limited number of data points along a given line of sight. We decrease this type of noise by increasing the number of points and by introducing Gaussian smoothing. We argue that in real observations the number of emitting elements is essentially infinite and that source of noise vanishes.

Magnetic reconnection, or the ability of the magnetic field lines that are frozen in plasma to change their topology, is a fundamental problem of magnetohydrodynamics (MHD). We briefly examine the problem starting with the well-known Sweet-Parker scheme, discuss effects of tearing modes, anomalous resistivity and the concept of hyperresistivity. We show that the field stochasticity by itself provides a way to enable fast reconnection even if, at the scale of individual turbulent wiggles, the reconnection happens at the slow Sweet-Parker rate. We show that fast reconnection allows efficient mixing of magnetic field in the direction perpendicular to the local direction of magnetic field. While the idea of stochastic reconnection still requires numerical confirmation, our numerical simulations testify that mixing motions perpendicular to the local magnetic field are up to high degree hydrodynamical. This suggests that the turbulent heat transport should be similar to that in non-magnetized turbulent fluid, namely, should have a diffusion coefficient \sim LV_L, where V_L is the amplitude of the turbulent velocity and L is the scale of the turbulent motions. We present numerical simulations which support this conclusion. The application of this idea to thermal conductivity in clusters of galaxies shows that this mechanism may dominate the diffusion of heat and may be efficient enough to prevent cooling flow formation.

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.

Understanding of the nature of gamma-ray bursts (GRBs) is one of the challenging problem facing the astrophysics community. The field of gamma ray bursts is a rapidly developing one and we expect that new missions, like SWIFT and GLAST will bring the field to a new quantitative level. The detail of the explosion, the formation of the fireball, its propagation, the generation of the shocks, the source of the magnetic fields and soft photons, the particle acceleration process, and the details of the radiation process are all still outstanding questions. The goal of this review is to attract the attention of the community to a number of physical processes (MHD turbulence and particle acceleration, magnetic reconnection, etc.) that seem to be very relevant for these sources.

Compressible turbulence, especially the magnetized version of it, traditionally has a bad reputation with researchers. However, recent progress in theoretical understanding of incompressible MHD as well as that in computational capabilities enabled researchers to obtain scaling relations for compressible MHD turbulence. We discuss scalings of Alfven, fast, and slow modes in both magnetically dominated (low $\beta$) and gas pressure dominated (high $\beta$) plasmas. We also show that the new regime of MHD turbulence below viscous cutoff reported earlier for incompressible flows persists for compressible turbulence. Our recent results show that this leads to density fluctuations. New understanding of MHD turbulence is likely to influence many key astrophysical problems.

Recent advances in understanding of magnetohydrodynamic (MHD) turbulence call for substantial revisions in our understanding of cosmic ray transport. In this paper we use gyroresonance recently obtained scaling laws for MHD modes to calculate the scattering frequency for cosmic rays in the ISM. We consider gyroresonance with MHD modes (Alfvenic, slow and fast) and transit-time damping (TTD) by fast modes. We conclude that the gyroresonance with fast modes is the dominant contribution to cosmic ray scattering for the typical interstellar conditions. In contrast to earlier studies, we find that Alfvenic and slow modes are inefficient because they are far from isotropy usually assumed.

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 Alfven, slow, and fast waves. We discuss the implications of this new insight into MHD turbulence for cosmic ray transport.

Aligned grains present a foreground for cosmic microwave emission studies. We review basic physical processes involved in grain alignment and discuss the niches for different alignment mechanisms. We show that mechanisms which were favored for decades do not look so promising right now, while the radiative torque mechanism ignored for more than 20 years looks quite attractive. We discuss alignment of small rotating grains that are thought to be responsible for 10-90 GHz anomalous foreground and polarization arising from magnetic grains.

Most astrophysical systems, e.g. stellar winds, the diffuse interstellar medium, molecular clouds, are magnetized with magnetic fields that influence almost all of their properties. One of the most informative techniques of magnetic field studies is based on the use of starlight polarization and polarized emission arising from aligned dust. How reliable the interpretation of the polarization maps in terms of magnetic fields is the issue that the grain alignment theory addresses. Although grain alignment is a problem of half a century standing, recent progress achieved in the field makes us believe that we are approaching the solution of this mystery. I review basic physical processes involved in grain alignment and discuss the niches for different alignment mechanisms. I show why mechanisms that were favored for decades do not look so promising right now, while the radiative torque mechanism ignored for more than 20 years looks so attractive. I define the observational tests and outline the circumstances when grain alignment theory predicts that new yet untapped information of magnetic field structure is available through polarimetry. In particular, I touch upon mapping magnetic fields in circumstellar regions, interplanetary space and in comet comae.

Neutron stars, asteroids, comets, cosmic-dust granules, spacecraft, as well as whatever other freely spinning body dissipate energy when they rotate about any axis different from principal. We discuss the internal-dissipation-caused relaxation of a freely precessing rotator towards its minimal-energy mode (mode that corresponds to the spin about the maximal-inertia axis). While the body nutates at some rate, the internal stresses and strains within the body oscillate at frequencies both higher and lower than this rate. The internal dissipation takes place mostly the second and higher harmonics. We discuss the application of our findings to asteroids. Regarding the comets, estimates show that the currently available angular resolution of spacecraft-based instruments makes it possible to observe wobble damping within year- or maybe even month-long spans of time. We also discuss cosmic-dust astrophysics; in particular, the role played by precession damping in the dust alignment. We show that this damping provides coupling of the grain's rotational and vibrational degrees of freedom; this entails occasional flipping of dust grains due to thermal fluctuations. During such a flip, grain preserves its angular momentum, but the direction of torques arising from H2 formation reverses. As a result, flipping grain will not rotate fast in spite of the action of uncompensated H2 formation torques. The grains get ``thermally trapped,'' and their alignment is marginal.

Astrophysical objects with negligible resistivity are often threaded by large scale magnetic fields. The generation of these fields is somewhat mysterious, since a magnetic field in a perfectly conducting fluid cannot change the flux threading a fluid element, or the field topology. Classical dynamo theory evades this limit by assuming that magnetic reconnection is fast, even for vanishing resistivity, and that the large scale field can be generated by the action of kinetic helicity. Both these claims have been severely criticized, and the latter appears to conflict with strong theoretical arguments based on magnetic helicity conservation and a series of numerical simulations. Here we discuss recent efforts to explain fast magnetic reconnection through the topological effects of a weak stochastic magnetic field component. We also show how mean-field dynamo theory can be recast in a form which respects magnetic helicity conservation, and how this changes our understanding of astrophysical dynamos. Finally, we comment briefly on why an asymmetry between small scale magnetic and velocity fields is necessary for dynamo action, and how it can arise naturally.

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

We discuss the relative grain motions due to MHD turbulence in interstellar medium. It has been known for decades that turbulent drag is an efficient way to induce grain relative motions. However, earlier treatments disregarded magnetic field and used Kolmogorov turbulence. Unlike hydro turbulence, MHD turbulence is anisotropic on small scales. Moreover, compressible modes are important for MHD and magnetic perturbations can directly interact with grains. We provide calculations of grain relative motion for realistic interstellar turbulence driving that is consistent with the velocity dispersions observed in diffuse gas and for realistic grain charging. We account for the turbulence cutoff arising from abmipolar drag. Our results on grain shattering are consistent with the customary accepted cutoff size. We obtain grain velocities for turbulence with parameters consistent with those in HI and dark clouds. These velocities are smaller than those in earlier papers, where MHD effects were disregarded. Finally, we consider grain velocities arising from photoelectric emission, radiation pressure and the thrust due to molecular hydrogen formation. These are lower than relative velocities induced by turbulence. We conclude that turbulence should prevent these mechanisms from segregating grains by size.

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.

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.

In astrophysical situations, e.g. in the interstellar medium (ISM), neutrals can provide viscous damping on scales much larger than the magnetic diffusion scale. Through numerical simulations, we have found that the magnetic field can have a rich structure below the dissipation cutoff scale. This implies that magnetic fields in the ISM can have structures on scales much smaller than parsec scales. Our results show that the magnetic energy contained in a wavenumber band is independent of the wavenumber and magnetic structures are intermittent and extremely anisotropic. We discuss the relation between our results and the formation of the tiny-scale atomic structure (TSAS).

Interstellar medium is turbulent and this induces relative motions of dust grains. We calculate relative velocities of charged grains in a partially ionized magnetized gas. We account for anisotropy of magnetohydrodynamic (MHD) turbulence, grain coupling with magnetic field, and the turbulence cutoff arising from the ambipolar drag. We obtain grain velocities for turbulence with parameters consistent with those in HI and dark clouds. Those velocities are smaller than those in earlier papers, where MHD effects were disregarded. Finally, we compare grain velocities arising from photoelectric emission, radiation pressure and thrusts from molecular hydrogen formation. We conclude that turbulence should prevent segregation of grains of different sizes.

Galactic HI is a gas that is coupled to magnetic field because of its fractional ionization. Many properties of HI are affected by turbulence. Recently, there has been a significant breakthrough on the theory of magnetohydrodynamic (MHD) turbulence. For the first time in the history of the subject we have a scaling model that is supported by numerical simulations. We review recent progress in studies of both incompressible and compressible turbulence. We also discuss the new regime of MHD turbulence that happens below the scale at which conventional turbulent motions get damped by viscosity. The viscosity in the case of HI is being produced by neutrals and truncates the turbulent cascade at parsec scales. We show that below this scale magnetic fluctuations with a shallow spectrum persist and point out to a possibility of the resumption of the MHD cascade after ions and neutrals decouple. We discuss the implications of the new insight into MHD turbulence for cosmic ray transport and grain dynamics.

HI data cubes are sources of unique information on interstellar turbulence. Doppler shifts due to supersonic motions contain information on turbulent velocity field which is otherwise difficult to obtain. However, the problem of separation of velocity and density fluctuations within HI data cubes is far from being trivial. Analytical description of the emissivity statistics of channel maps (velocity slices) in Lazarian & Pogosyan (2000) showed that the relative contribution of the density and velocity fluctuations depends on the thickness of the velocity slice. In particular, power-law assymptotics of the emissivity fluctuations change when the dispersion of the velocity at the scale under study becomes of the order of the velocity slice thickness (integrated width of the channel map). These results are the foundations of the Velocity-Channel Analysis (VCA) technique which allows to determine velocity and density statistics using 21-cm data cubes. The VCA has been successfully tested using data cubes obtained via compressible magnetohydrodynamic simulations and applied to Galactic and Magellanic Clouds data. As a tool it has become much more sophisticated recently when effects of absorption were accounted for. The systematic studies of vast 21-cm data sets to correlate the variations in the turbulence statistics with the astrophysical activity is on the agenda. This should allow to determine the interstellar energy injection mechanisms. Going beyond the VCA, we discuss other tools, namely, genus and anisotropy analysis. The first characterises the topology of HI, while the second provides magnetic field directions. We show a few applications of these new tools to HI data and MHD simulations.

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.

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

Polarized emission from dust is an important foreground that can hinder the progress in polarized CMB studies unless carefully accounted for. We discuss potential difficulties associated with the dust foreground, namely, the existence of different grain populations with very different emission/polarization properties and variations of the polarization yield with grain temperature. In this context we appeal for systematic studies of polarized dust emission as the means of dealing with this foreground.

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

We present a statistical analysis of Galactic interstellar polarization from the largest compilation available of starlight data. The data comprises ~ 9300 stars of which we have selected ~ 5500 for our analysis. We find a nearly linear growth of mean polarization degree with extinction. The amplitude of this correlation shows that interstellar grains are not fully aligned with the Galactic magnetic field, which can be interpreted as the effect of a large random component of the field. In agreement with earlier studies of more limited scope, we estimate the ratio of the uniform to the random plane-of-the-sky components of the magnetic field to be B_u/B_r = 0.8. Moreover, a clear correlation exists between polarization degree and polarization angle what provides evidence that the magnetic field geometry follows Galactic structures on large-scales. The angular power spectrum C_l of the starlight polarization degree for Galactic plane data (|b| < 10 deg) is consistent with a power-law, C_l ~ l^{-1.5} (where l ~ 180 deg/\theta is the multipole order), for all angular scales \theta > 10 arcmin. An investigation of sparse and inhomogeneous sampling of the data shows that the starlight data analyzed traces an underlying polarized continuum that has the same power spectrum slope, C_l ~ l^{-1.5}. Our findings suggest that starlight data can be safely used for the modeling of Galactic polarized continuum emission at other wavelengths.

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

In this paper we test the results of a recent analytical study by Lazarian and Pogosyan, on the statistics of emissivity in velocity channel maps, in the case of realistic density and velocity fields obtained from numerical simulations of MHD turbulence in the interstellar medium (ISM). To compensate for the lack of well-developed inertial ranges in the simulations due to the limited resolution, we apply a procedure for modifying the spectral slopes of the fields while preserving the spatial structures. We find that the density and velocity are moderately correlated in space and prove that the analytical results by Lazarian and Pogosyan hold in the case when these fields obey the fluid conservation equations. Our results imply that the spectra of velocity and density can be safely recovered from the position-position-velocity (PPV) data cubes available through observations, and confirm that the relative contributions of the velocity and density fluctuations to those of the emissivity depend on the velocity resolution used and on the steepness of the density spectral index. Furthermore, this paper supports previous reports that an interpretation of the features in the PPV data cubes as simple density enhancements (i.e., ``clouds'') can be often erroneous, as we observe that changes in the velocity statistics substantially modify the statistics of emissivity within the velocity data cubes.

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.

Observations of cosmic microwave background in the range 10-90 GHz have revealed an anomalous foreground component well correlated with 12 \mum, 60 \mum and 100 \mum emission from interstellar dust. As the recent cross-correlation analysis of WHAM H\alpha maps with the Tenerife 10 and 15 GHz maps supports an earlier conclusion that the emission does not arise from free-free radiation, the interstellar dust origin of it is left as the only suspect. Two competing models of this emission exist. The more favored at the moment is the spinning dust model, the other is the model that uses grains with strong magnetic response. In the spinning dust model the emission arises from rapid rotation of ultrasmall grains that have dipole moments, while in the other model magnetic grains emit due to thermal vibrations of magnetic dipoles. Both models predict the emission to be partially polarized and this emission can seriously interfere with the CMB polarization measurements. We discuss observational signatures that can be used to distinguish and eventually filter out the polarized component of the microwave dust radiation.

We show that the energy-level splitting arising from grain rotation ensures that paramagnetic dissipation acts at its maximum rate, i.e., the conditions for paramagnetic resonance are automatically fulfilled. We refer to this process as ``resonance relaxation''. The differences between the predictions of classical Davis-Greenstein relaxation and resonance relaxation are most pronounced for grains rotating faster than 1 GHz, i.e., in the domain where classical paramagnetic relaxation is suppressed. This mechanism can partially align even very small grains, resulting in linearly polarized microwave emission which could interfere with efforts to measure the polarization of the cosmic microwave background.

We here investigate the possibility that the ultra-high energy cosmic ray (UHECR) events observed above the GZK limit are mostly protons accelerated in reconnection sites just above the magnetosphere of newborn millisecond pulsars which are originated by accretion induced collapse (AIC). We show that AIC-pulsars with surface magnetic fields $10^{12} G < B_{\star} \lesssim 10^{15}$ G and spin periods $1 ms \lesssim P_{\star} < 60 ms$, are able to accelerate particles to energies $\geq 10^{20}$ eV. Because the expected rate of AIC sources in our Galaxy is very small (\sim 10^{-5} yr^{-1}), the corresponding contribution to the flux of UHECRs is neglegible, and the total flux is given by the integrated contribution from AIC sources produced by the distribution of galaxies located within the distance which is unaffected by the GZK cutoff ($\sim 50$ Mpc). We find that the reconnection efficiency factor needs to be $\xi \gtrsim 0.1$ in order to reproduce the observed flux of UHECRs.

We here investigate the possibility that the ultra-high energy cosmic ray (UHECR) events observed above the GZK limit are mostly protons accelerated in econnection sites just above the magnetosphere of newborn millisecond pulsars which are originated by accretion induced collapse (AIC).

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.

The distribution of atomic hydrogen in the Galactic plane is usually mapped using the Doppler shift of 21cm emission line. We calculate the emission spectrum in velocity slices of data (channel maps) and derive its dependence on the statistics of HI velocity and density fields. We find that (a) if the density spectrum is steep, i.e. n<-3, the short-wave asymptotics of the emissivity spectrum is dominated by velocity fluctuations; (b) the velocity fluctuations make the emission spectra shallower, provided that the data slices are sufficiently thin. In other words, turbulent velocity creates small scale structure that can erroneously be identified as clouds. The contribution of fluctuations in warm HI is suppressed relative to cold component when velocity channels used are narrower than warm HI thermal velocity and small angular scale fluctuations are measured. We calculate how emission spectra vary with the change of velocity slice thickness and show that the observational 21cm data is consistent with the explanation that intensity fluctuations within individual channel maps are generated by a turbulent velocity field. As the thickness of velocity slices increases density fluctuations get to dominate the emissivity. This allows to disentangle velocity and density statistics. Application of our technique to the Galactic and SMC data reveals spectra of density and velocity with the power law index close to -11/3. This is a Kolmogorov index, but the explanation of the spectrum appealing to the Kolmogorov-type cascade faces substantial difficulties. We generalize our treatment for the case of a statistical study of turbulence inside individual clouds. The mathematical machinery developed is applicable to other emission lines.

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.

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.

Observations of the cosmic microwave background have revealed a component of 10-60 GHz emission from the Galaxy which correlates with 100-140um emission from interstellar dust but has an intensity much greater than expected for the low-frequency tail of the "electric dipole vibrational" emission peaking at \~130um. This "anomalous emission" is more than can be accounted for by dust-correlated free-free emission. The anomalous emission could be due in part to magnetic dipole emission from thermal fluctuations of the magnetization within interstellar dust grains, but only if a substantial fraction of the Fe in interstellar dust resides in magnetic materials such as metallic iron or magnetite. The observed anomalous emission is probably due primarily to electric dipole radiation from spinning ultrasmall interstellar dust grains. This rotational emission is expected to be partially polarized. From the standpoint of minimizing confusion with non-CBR foregrounds, 60-120 GHz appears to be the optimal frequency window.

In this review, we intend to present the current knowledge of the polarized emission from thermal dust in our Galaxy. We show different methods to estimate the spatial distribution statistics of this emission in the lack of any data from the diffuse ISM, and compare it to the expected CMB polarized signal. We finally show how this contaminant could be efficiently removed from CMB maps using multi-frequency observations.

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.

Since the classical work by Purcell (1979) it has been generally accepted that most interstellar grains rotate suprathermally. Suprathermally rotating grains would be nearly perfectly aligned with the magnetic field by paramagnetic dissipation if not for ``crossovers'', intervals of low angular velocity resulting from reversals of the torques responsible for suprathermal rotation; during crossovers grains are susceptible to disalignment by random impulses. Lazarian and Draine (1997) identified thermal fluctuations within grain material as an important component of crossover dynamics. For grains of size less than 0.1 micron, these fluctuations ensure good correlation of angular momentum before and after crossover resulting in good alignment, in accord with observations of starlight polarization. In the present paper we discuss two new processes which are important for the dynamics of grains with a<0.1 micron. The first -- ``thermal flipping'' -- offers a way for small grains to bypass the period of greatly reduced angular momentum which would otherwise take place during a crossover, thereby enhancing the alignment of small grains. The second effect -- ``thermal trapping'' -- arises when thermal flipping becomes rapid enough to prevent the systematic torques from driving the grain to suprathermal rotation. This effect acts to reduce the alignment of small grains. The observed variation of grain alignment with grain size would then result from a combination of the thermal flipping process -- which suppresses suprathermal rotation of small grains -- and due to molecular hydrogen formation and starlight -- which drive large grains to suprathermal rotation rates.

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.

Motivated by a recent study by Lazarian and Draine, which showed that a high degree of grain alignment of the paramagnetic dust is achievable if the rates of internal relaxation are controlled by the Barnett relaxation process, we undertake a study of an alternative mechanism of internal dissipation, namely, the inelastic dissipation of energy in oblate dust grains. We find that deformations at double frequency that were disregarded in earlier studies dominate the inelastic relaxation. Our results indicate that for grains with 4:1 axis ratio, or/and grains formed via agglomeration inelastic relaxation dominates the Barnett relaxation within large (a>0.1 micron) grains even when they rotate at thermal rate. For grains with axis ratio less than 1:2 the inelastic relaxation is dominant when grains rotate suprathermally.

We propose that inward, subsonic flows arise from the local dissipation of turbulent motions in molecular clouds. Such "turbulent cooling flows" may account for recent observations of spatially extended inward motions towards dense cores. These pressure-driven flows may arise from various types of turbulence and dissipation mechanisms. For the example of MHD waves and turbulence damped by ion-neutral friction, sustained cooling flow requires that the outer gas be sufficiently turbulent, that the inner gas have marginal field-neutral coupling, and that this coupling decrease sufficiently rapidly with increasing density. These conditions are most likely met at the transition between outer regions ionized primarily by UV photons and inner regions ionized primarily by cosmic rays. If so, turbulent cooling flows can help form dense cores, with speeds faster than expected for ambipolar diffusion. Such motions could reduce the time needed for dense core formation and could precede and enhance the motions of star-forming gravitational infall.

Thermal fluctuations in the magnetization of interstellar grains will produce magnetic dipole emission at frequencies below ~100 GHz. We show how to calculate absorption and emission from small particles composed of magnetic materials. The Kramers-Kronig relations for a dusty medium are generalized to include the possibility of magnetic grains. The frequency-dependent magnetic permeability is discussed for candidate grain materials, including iron and magnetite. We calculate emission spectra for various interstellar grain candidates. While paramagnetic grains or magnetite grains cannot account for the observed "anomalous" emission from dust in the 14-90 GHz range, stronger magnetic dipole emission will result if a fraction of the grain material is ferromagnetic, as could be the case given the high Fe content of interstellar dust. The observed emission from dust near 90 GHz implies that not more than 5% of interstellar Fe is in the form of metallic iron grains or inclusions (e.g., in "GEMS"). However, we show that if most interstellar Fe is in a moderately ferromagnetic material, it could contribute a substantial fraction of the observed 14-90 GHz emission, perhaps comparable to the contribution from spinning ultrasmall dust grains. The two emission mechanisms can be distinguished by measuring the emission from dark clouds. The expected polarization of magnetic dipole emission is discussed

We discuss the role of tearing instabilities in magnetic reconnection. In three dimensions this instability leads to the formation of strong Alfvenic waves that remove plasma efficiently from the reconnection layer. As a result the instability proceeds at high rates while staying close to the linear regime. Our calculations show that for a resistive fluid the reconnection speed scales as the product of the Alfven speed V_A over the magnetic Reynolds number to the power -0.3. In the limit of vanishing resistivity, tearing modes proceed at a non-zero rate, driven by the electron inertia term, giving rise to a reconnection speed V_A (c/\omega_p L_x)^{3/5}, where \omega_p is the plasma frequency and L_x is the transverse scale of the reconnection layer. Formally this solves the problem of fast reconnection, but in practice this reconnection speed is small.

Understanding the properties of interstellar turbulence is a great intellectual challenge and the urge to solve this problem is partially motivated by a necessity to explain the star formation mystery. This review deals with a recently suggested inversion technique as applied to atomic hydrogen. This technique allows to determine 3D turbulence statistics through the variations of 21 cm intensity. We claim that a radio interferometer is an ideal tool for such a study as its visibility function is directly related to the statistics of galactic HI. Next, we show how galactic rotation curve can be used to study the turbulence slice by slice and relate the statistics given in galactic coordinates and in the velocity space. The application of the technique to HI data reveals a shallow spectrum of the underlying HI density that is not compatible with a naive Kolmogorov picture. We show that the random density corresponding to the found spectrum tends to form low contrast filaments that are elongated towards the observer.

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.

Spinning interstellar dust grains produce detectable rotational emission in the 10-100 GHz frequency range. We calculate the emission spectrum, and show that this emission can account for the ``anomalous'' Galactic background component which correlates with 100um thermal emission from dust. Implications for cosmic background studies are discussed.

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\'en 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.

A quantitative analysis of grain alignment in the filamentary dark cloud L1755 in Ophiuchus is presented. We show, that the observed decrease of the polarization to extinction ratio for the inner parts of this quiescent dark cloud can be explained as a result of the decrease of the efficiency of grain alignment.We make quantitative estimates of grain alignment efficiency for six mechanisms, involving grains with either thermal or suprathermal rotation, interacting with either magnetic field or gaseous flow. We also make semiquantitative estimates of grain alignment by radiative torques.We show that in conditions typical of dark cloud interiors, all known major mechanisms of grain alignment fail. All the studied mechanisms predict polarization at least an order of magnitude below the currently detectable levels of order of 1%. On the contrary, in the dark cloud environments where A_v<1 the grain alignment can be much more efficient. There the alignment of suprathermally rotating grains with superparamagentic inclusions, and possibly also radiative torques account for observed polarization. These results apply to L1755, which we model in detail and probably also to B216 and other similar dark clouds. Our study suggests an explanation for the difference in results obtained through polarimetry of background starlight and polarized thermal emission from the dust itself. We conjecture that the emission polarimetry selectively reveals aligned grains in the environment far from thermodynamic equilibrium as opposed to starlight polarization studie that probe the alignment of grains all the way along the line of sight, including the interiors of dark quiescent clouds, where no alignment is possible.

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.

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.

We present an exact formulation of the physics of Barnett relaxation. Our formulation is based on a realistic kinetic model of the relaxation mechanism which includes the alignment of the grain angular momentum in body coordinates by Barnett dissipation, disalignment by thermal fluctuations, and coupling of the angular momentum to the gas via gas damping. We solve the Fokker-Planck equation for the measure of internal alignment using numerical integration of the equivalent Langevin equation for Brownian rotation. The accuracy of our results is calibrated by comparing our numerical solutions with exact analytic results obtained for special cases.We describe an analytic approximation for the measure of alignment which fits our numerical results for cases of practical interest.

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.

We discuss the dynamics of dust grains subjected to uncompensated torques arising from H_2 formation. In particular, we discuss grain dynamics when a grain spins down and goes through a ``crossover''. As first pointed out by Spitzer & McGlynn (1979), the grain angular momentum before and after a crossover event are correlated, and the degree of this correlation affects the alignment of dust grains by paramagnetic dissipation. We calculate the correlation including the important effects of thermal fluctuations within the grain material. These fluctuations limit the degree to which the grain angular momentum J is coupled with the grain principal axis a_1 of maximal inertia. We show that this imperfect coupling of a_1 with J plays a critical role during crossovers and can substantially increase the efficiency of paramagnetic alignment for grains larger than 10^{-5} cm. As a result, we show that for reasonable choices of parameters, the observed alignment of $a\gtrsim 10^{-5}$ cm grains could be produced by paramagnetic dissipation in suprathermally rotating grains, if radiative torques due to starlight were not present. We also show that the efficiency of mechanical alignment in the limit of long alignment times is not altered by the thermal fluctuations in the grain material.

This paper presents a statistical explanation of filament formation in the galactic atomic hydrogen. Recently developed technique allows to determine the 3D spectrum of random HI density. We claim that even in the absence of dynamical factors the Gaussian field corresponding to the measured values of the spectrum of HI density should exhibit filamentary structure, the existence of which has long been claimed.

This paper provides a quantitative account of a recently introduced mechanism of mechanical alignment of suprathermally rotating grains. These rapidly rotating grains are essentially not susceptible to random torques arising from gas-grain collisions, as the timescales for such torques to have significant effect are orders of magnitude greater than the mean time between crossovers. Such grains can be aligned by gaseous torques during the short periods of crossovers and/or due to the difference in the rate at which atoms arrive at grain surface. The latter is a result of the difference in orientation of a grain in respect to the supersonic flow. This process, which we call cross-section alignment, is the subject of our present paper. We derive expressions for the measure of cross-section alignment for oblate grains and study how this measure depends upon the angle between the interstellar magnetic field and the gaseous flow and upon the grain shape.

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\'{e}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.

The frequency distributions of sizes of ``clouds" and ``clumps" within clouds are significantly flatter for extinction surveys than for CO spectral line surveys, even for comparable size ranges. A possible explanation is the blocking of extinction clouds by larger foreground clouds (occlusion), which should not affect spectral line surveys much because clouds are resolved in velocity space along a given line of sight. We present a simple derivation of the relation between the true and occluded size distributions, assuming clouds are uniformly distributed in space or the distance to a cloud comples is much greater than the size of the complex. Because the occlusion is dominated by the largest clouds, we find that occlusion does not affect the measured size distribution except for sizes comparable to the largest size, implying that occlusion is not responsible for the discrepancy if the range in sizes of the samples is large. However, we find that the range in sizes for many of the published observed samples is actually quite small, which suggests that occlusion does affect the extinction sample and/or that the discrepancy could arise from the different operational definitions and selection effects involved in the two samples. Size and mass spectra from an IRAS survey (Wood \etal\ 1994) suggest that selection effects play a major role in all the surveys. We conclude that a reliable determination of the ``true" size and mass spectra of clouds will require spectral line surveys with very high signal-to-noise and sufficient resolution and sampling to cover a larger range of linear sizes, as well as careful attention to selection effects.

Paramagnetic alignment of non-spherical dust grains rotating at thermal velocities is studied. The analytical solution is found for the alignment measure of oblate grains. Perturbative approach is used for solving the problem. It is shown that even the first approximation of the suggested iteration procedure provides the accuracy well within one percent of the expected measure of alignment. The results obtained are applicable both to paramagnetic and to superferromagnetic grains.

It is shown that mechanical alignment of grains can be efficient for grains rotating suprathermally, i.e. with kinetic energy substantially exceeding $k$ (the Boltzmann constant) over any temperature in the system. The paper studies suprathermal rotation caused by H$_{2}$ formation and the alignment that takes place due to crossover events. Gaseous bombardment in the course of a crossover as well as both gaseous friction and poisoning of active sites are shown to produce alignment. The first type of alignment happens due to the angular momentum deposited by a corpuscular flux with a grain, the second is caused by the change of the mean time back to crossover due to the interaction with a gaseous flux. We show that the two processes act as to decrease the grain cross section in respect to the flux and we find the Rayleigh reduction factor for the joint action of the two processes as well as the range of applicability of each of the processes. Our study indicates that mechanical alignment can be more widely spread than it is generally accepted.

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

Paramagnetic alignment of fractal suprathermally rotating grains is discussed. It is shown that if the concentration of H$_{2}$formation sites is low and resurfacing is active, fractal structure of grains enhances their alignment. Studying the influence of grain surface physics and chemistry on the alignment we found that there exist two critical values of grain sizes, and the alignment of smaller grains is expected to decrease as compared to the predictions of the Purcell's theory (1979). One of the critical sizes is due to intensified poisoning of active sites, while the other is a result of a finite value of the imaginary part of magnetic susceptibility in the limit of high angular velocities. On the contrary, if active sites completely cover grain surface, suprathermal rotation, and therefore efficient alignment, is possible only for a limited range of grain sizes.