We present a model for atmospheric wind circulation in binary millisecond pulsar (MSP) companions, showing how the optical light curve (LC) and radial velocities (RV) are sensitive to the wind flow, causing LC orbital phase shifts and asymmetries, as observed for several `spider' MSP. Velocity widths of spectral lines offer additional opportunities for measuring surface wind speed. As examples, we fit optical data for a black widow pulsar J1959+2048 and a redback pulsar J2215+5135; the wind heating models (WH) are statistically strongly preferred over direct heating (DH) for both objects, although the latter is even better fit with a heated spot. In general, WH effects tend to increase the inferred orbital inclination $i$ and decrease the inferred companion center-of-mass radial velocity amplitude $K_\mathrm{c}$; both effects decrease the inferred neutron star mass. Even with such a decrease, we find large masses for the two neutron stars: $2.18\pm0.09M_\odot$ and $2.28^{+0.10}_{-0.09}M_\odot$, respectively (for the modest surface speeds fit from the bulk heat flow; supersonic photospheric winds can slightly change these values). These are among the highest masses known, and our improved modeling increases confidence that the results are important for understanding the dense matter equation of state.

We model millisecond pulsars winds colliding with radiatively-driven companion winds in black widow and redback systems. For the redbacks, the geometry of this intrabinary shock (IBS) is quite sensitive to the expected equatorial concentration in the pulsar outflow. We thus analytically extend IBS thin-shock models to $\sim \sin^{2n}\theta$ pulsar winds. We compute the synchrotron emission from such shocks, including the build-up and cooling of the particle population as it accelerates along the IBS. For reasonable parameters, this IBS flux dominates the binary emission in the X-ray band. The modeling shows subtle variation in spectrum across the IBS peak, accessible to sensitive X-ray studies. As example applications, we fit archival CXO/XMM data from the black widow pulsar J1959+2048 and the redback PSR J2339-0533, finding that the model reproduces well the orbital light curve profiles and energy spectra. The results show a very hard injected electron spectrum, indicating likely dominance by reconnection. The light curve fitting is sensitive to the geometric parameters, including the very important orbital inclination $i$. Coupled with optical fits of the companion star, such IBS X-ray light curve modeling can strongly constrain the binary geometry and the energetics of the MSP wind.

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.

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.

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.