We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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.

We perform numerical simulations of hydrodynamic (HD) and magnetohydrodynamic (MHD) turbulence driven by compressive driving to study generation of solenoidal velocity component and small-scale magnetic field. We mainly focus on the effects of mean magnetic field ($B_0$) and the sonic Mach number ($M_s$). We also consider two different driving schemes in terms of correlation timescale of forcing vectors: a finite-correlated driving and a delta-correlated driving. The former has a longer correlation timescale of forcing vectors, which is comparable to large-eddy turnover time, than the latter. Our findings are as follows. First, when we fix the value of $B_0$, the level of solenoidal velocity component after saturation increases as $M_s$ increases. A similar trend is observed for generation of magnetic field when $B_0$ is small. Second, when we fix the value of $M_s$, HD and MHD simulations result in similar level of the solenoidal component when $B_0$ $\lesssim$ 0.2 (or Alfven Mach number of $\sim$ 5). However, the level increases when $B_0$ $\gtrsim$ 0.2. Roughly speaking, the magnetic energy density after saturation is a linearly increasing function of $B_0$ irrespective of $M_s$. Third, generation of solenoidal velocity component is not sensitive to numerical resolution, but that of magnetic energy density is mildly sensitive. Lastly, when initial conditions are same, the finite-correlated driving always produces more solenoidal velocity and small-scale magnetic field components than the delta-correlated driving. We additionally analyze the vorticity equation to understand why higher $M_s$ and $B_0$ yield larger quantity of the solenoidal velocity component.

We present numerical simulations of decaying hydrodynamic turbulence initially driven by solenoidal (divergence-free) and compressive (curl-free) driving. Most previous numerical studies for decaying turbulence assume an isothermal equation of state (EOS). Here we use a polytropic EOS, $P$ $\propto$ $\rho^{\gamma}$, with polytropic $\gamma$ ranging from 0.7 to 5/3. We mainly aim at determining the effects of polytropic $\gamma$ and driving schemes on the decay law of turbulence energy, E $\propto$ $t^{-\alpha}$. We additionally study probability density function (PDF) of gas density and skewness of the distribution in polytropic turbulence driven by compressive driving. Our findings are as follows. First of all, we find that even if polytropic $\gamma$ does not strongly change scaling relation of the decay law, the driving schemes weakly change the relation; in our all simulations, turbulence decays with $\alpha$ $\approx$ 1, but compressive driving yields smaller $\alpha$ than solenoidal driving at the same sonic Mach number. Second, we calculate compressive and solenoidal velocity components separately and compare their decay rates in turbulence initially driven by compressive driving. We find that the former decays much faster so that it ends up having a smaller fraction than the latter. Third, the density PDF of compressively driven turbulence with polytropic $\gamma$ $>$ 1 deviates from log-normal distribution: it has a power-law tail at low density as in the case of solenoidally driven turbulence. However, as it decays, the density PDF becomes approximately log-normal. We discuss why decay rates of compressive and solenoidal velocity components are different in compressively driven turbulence and astrophysical implication of our findings.

Dust polarization depends on the physical and mechanical properties of dust, as well as the properties of local environments. To understand how dust polarization varies with grain mechanical properties and the local environment, in this paper, we model the wavelength-dependence polarization of starlight and polarized dust emission by aligned grains by simultaneously taking into account grain alignment and rotational disruption by radiative torques (RATs). We explore a wide range of the local radiation field and grain mechanical properties characterized by tensile strength. We find that the maximum polarization and the peak wavelength shift to shorter wavelengths as the radiation strength $U$ increases due to the enhanced alignment of small grains. Grain rotational disruption by RATs tends to decrease the optical-near infrared polarization but increases the ultraviolet polarization of starlight due to the conversion of large grains into smaller ones. In particular, we find that the submillimeter (submm) polarization degree at $850~\mu \rm m$ ($P_{850}$) does not increase monotonically with the radiation strength or grain temperature ($T_{d}$), but it depends on the tensile strength of grain materials. Our physical model of dust polarization can be tested with observations toward star-forming regions or molecular clouds irradiated by a nearby star, which have higher radiation intensity than the average interstellar radiation field. Finally, we compare our predictions of the $P_{850}-T_{d}$ relationship with {\it Planck} data and find that the observed decrease of $P_{850}$ with $T_{d}$ can be explained when grain disruption by RATs is accounted for, suggesting that interstellar grains unlikely to have a compact structure but perhaps a composite one. The variation of the submm polarization with U (or $T_{d}$) can provide a valuable constraint on the internal structures of cosmic dust.

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.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

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.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We performed a molecular line survey of 82 high-mass protostellar objects in search for inflow signatures associated with high-mass star formation. Using the H$^{13}$CO$^+$ (1$-$0) line as an optically thin tracer, we detected a statistically significant excess of blue asymmetric line profiles in the HCO$^+$ (1$-$0) transition, but nonsignificant excesses in the HCO$^+$ (3$-$2) and H$_2$CO (2$_{12}$$-$1$_{11}$) transitions. The negative blue excess for the HCN (3$-$2) transition suggests that the line profiles are affected by dynamics other than inflow motion. The HCO$^+$ (1$-$0) transition thus seems to be the suitable tracer of inflow motions in high-mass star-forming regions, as previously suggested. We found 27 inflow candidates that have at least one blue asymmetric profile and no red asymmetric profile, and derived the inflow velocities to be 0.23$-$2.00 km s$^{-1}$ for 20 of them using a simple two-layer radiative transfer model. Our sample is divided into two groups in different evolutionary stages. The blue excess of the group in relatively earlier evolutionary stages was estimated to be slightly higher than that of the other in the HCO$^+$ (1$-$0) transition.

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.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

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

In many astrophysical systems, smoothly-varying large-scale variations coexist with small-scale fluctuations. For example, a large-scale velocity or density gradient can exist in molecular clouds that exhibit small-scale turbulence. In redshifted 21cm observations, we also have two types of signals - the Galactic foreground emissions that change smoothly and the redshifted 21cm signals that change fast in frequency space. Sometimes the large-scale variations make it difficult to extract information on small-scale fluctuations. We propose a simple technique to remove smoothly varying large-scale variations. Our technique relies on multi-point structure functions and can obtain the magnitudes of small-scale fluctuations. It can also help us to filter out large-scale variations and retrieve small-scale maps. We discuss applications of our technique to astrophysical observations.

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J = 1 - 0 lines, and dense cores in L1641 were observed in the 2.0-0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of 35 km/s and a precession period of 1600 years. The CO outflow length gives a dynamical timescale of 6300 years, and the protostar must be extremely young. The inferred binary separation of 6-70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of beta = 0.3-2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of beta = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices.

During the protostellar phase of stellar evolution, accretion onto the star is expected to be variable, but this suspected variability has been difficult to detect because protostars are deeply embedded. In this paper, we describe a sub-mm luminosity burst of the Class I protostar EC 53 in Serpens Main, the first variable found during our dedicated JCMT/SCUBA-2 monitoring program of eight nearby star-forming regions. EC 53 remained quiescent for the first 6 months of our survey, from February to August 2016. The sub-mm emission began to brighten in September 2016, reached a peak brightness of $1.5$ times the faint state, and has been decaying slowly since February 2017. The change in sub-mm brightness is interpreted as dust heating in the envelope, generated by a luminosity increase of the protostar of a factor of $\ge 4$. The 850~$\mu$m lightcurve resembles the historical $K$-band lightcurve, which varies by a factor of $\sim 6$ with a 543 period and is interpreted as accretion variability excited by interactions between the accretion disk and a close binary system. The predictable detections of accretion variability observed at both near-infrared and sub-mm wavelengths make the system a unique test-bed, enabling us to capture the moment of the accretion burst and to study the consequences of the outburst on the protostellar disk and envelope.

We determine the magnetic field strength in the OMC 1 region of the Orion A filament via a new implementation of the Chandrasekhar-Fermi method using observations performed as part of the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-Forming Region Observations (BISTRO) survey with the POL-2 instrument. We combine BISTRO data with archival SCUBA-2 and HARP observations to find a plane-of-sky magnetic field strength in OMC 1 of $B_{\rm pos}=6.6\pm4.7$ mG, where $\delta B_{\rm pos}=4.7$ mG represents a predominantly systematic uncertainty. We develop a new method for measuring angular dispersion, analogous to unsharp masking. We find a magnetic energy density of $\sim1.7\times 10^{-7}$ Jm$^{-3}$ in OMC 1, comparable both to the gravitational potential energy density of OMC 1 ($\sim 10^{-7}$ Jm$^{-3}$), and to the energy density in the Orion BN/KL outflow ($\sim 10^{-7}$ Jm$^{-3}$). We find that neither the Alfv\'{e}n velocity in OMC 1 nor the velocity of the super-Alfv\'{e}nic outflow ejecta is sufficiently large for the BN/KL outflow to have caused large-scale distortion of the local magnetic field in the $\sim$500-year lifetime of the outflow. Hence, we propose that the hour-glass field morphology in OMC 1 is caused by the distortion of a primordial cylindrically-symmetric magnetic field by the gravitational fragmentation of the filament and/or the gravitational interaction of the BN/KL and S clumps. We find that OMC 1 is currently in or near magnetically-supported equilibrium, and that the current large-scale morphology of the BN/KL outflow is regulated by the geometry of the magnetic field in OMC 1, and not vice versa.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We perform integrated spectroscopy of 24 Galactic globular clusters. Spectra are observed from one core radius for each cluster with a high wavelength resolution of ~2.0 A FWHM. In combination with two existing data sets from Puzia et al. (2002) and Schiavon et al. (2005), we construct a large database of Lick spectral indices for a total of 53 Galactic globular clusters with a wide range of metallicities, -2.4 < [Fe/H] < 0.1, and various horizontal-branch morphologies. The empirical index-to-metallicity conversion relationships are provided for the 20 Lick indices for the use of deriving metallicities for remote, unresolved stellar systems.

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.

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.

The Chandrasekhar-Fermi method is a powerful technique for estimating the strength of the mean magnetic field projected on the plane of the sky. In this paper, we present a technique for improving the Chandrasekhar-Fermi method, in which we take into account the averaging effect arising from independent eddies along the line of sight . In the conventional Chandrasekhar-Fermi method, the strength of fluctuating magnetic field divided by $\sqrt{4 \pi \bar{\rho}}$, where $\bar{\rho}$ is average density, is assumed to be comparable to the line-of-sight velocity dispersion. This however is not true when the driving scale of turbulence $L_f$, i.e. the outer scale of turbulence, is smaller than the size of the system along the line of sight $L_{los}$. In fact, the conventional Chandrasekhar-Fermi method over-estimates the strength of the mean plane-of-the-sky magnetic field by a factor of $\sim \sqrt{ L_{los}/L_f}$. We show that the standard deviation of centroid velocities divided by the average line-of-sight velocity dispersion is a good measure of $\sqrt{ L_{los}/L_f}$, which enables us to propose a modified Chandrasekhar-Fermi method.

Equatorial superrotation under zonally-symmetric thermal forcing is investigated in a setup close to that of the classic Held & Suarez (1994) setup. In contrast to the behaviour in the classic setup, a transition to equatorial superrotation occurs when the equator-to-pole surface equilibrium entropy gradient is weakened. Two factors contribute to this transition: 1) the reduction of breaking Rossby waves from the mid-latitude that decelerate the equatorial flow and 2) the presence of barotropic instability in the equatorial region, providing stirring to accelerate the equatorial flow. In the latter, Kelvin waves excited by instability near the equator generate and maintain the superrotation. However, the superrotation is unphysically enhanced if simulations are under-resolved and/or over-dissipated.

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 revisit the issue of sensitivity to initial flow and intrinsic variability in hot-Jupiter atmospheric flow simulations, originally investigated by Cho et al. (2008) and Thrastarson & Cho (2010). The flow in the lower region (~1 to 20 MPa) `dragged' to immobility and uniform temperature on a very short timescale, as in Liu & Showman (2013), leads to effectively a complete cessation of variability as well as sensitivity in three-dimensional (3D) simulations with traditional primitive equations. Such momentum (Rayleigh) and thermal (Newtonian) drags are, however, ad hoc for 3D giant planet simulations. For 3D hot-Jupiter simulations, which typically already employ strong Newtonian drag in the upper region, sensitivity is not quenched if only the Newtonian drag is applied in the lower region, without the strong Rayleigh drag: in general, both sensitivity and variability persist if the two drags are not applied concurrently in the lower region. However, even when the drags are applied concurrently, vertically-propagating planetary waves give rise to significant variability in the ~0.05 to 0.5 MPa region, if the vertical resolution of the lower region is increased (e.g. here with 1000 layers for the entire domain). New observations on the effects of the physical setup and model convergence in `deep' atmosphere simulations are also presented.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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 origin of magnetic fields in clusters of galaxies is still an unsolved problem, which is largely due to our poor understanding of initial seed magnetic fields. If the seed magnetic fields have primordial origins, it is likely that large-scale pervasive magnetic fields were present before the formation of the large-scale structure. On the other hand, if they were ejected from astrophysical bodies, they were highly localized in space at the time of injection. In this paper, using turbulence dynamo models for high magnetic Prandtl number fluids, we find constraints on the seed magnetic fields. The hydrodynamic Reynolds number based on the Spitzer viscosity in the intracluster medium (ICM) is believed to be less than $O(10^2)$, while the magnetic Reynolds number can be much larger than that. In this case, if the seed magnetic fields have primordial origins, they should be stronger than $O(10^{-11})$G, which is very close to the upper limit of $O(10^{-9})$G set by the cosmic microwave background (CMB) observations. On the other hand, if the seed magnetic fields were ejected from astrophysical bodies, any seed magnetic fields stronger than $O(10^{-9})$G can safely magnetize the intracluster medium. Therefore, it is less likely that primordial magnetic fields are the direct origin of present-day magnetic fields in the ICM.

Strong ionization on close-in extrasolar giant planets suggests that their atmospheres may be affected by ion drag and resistive heating arising from wind-driven electrodynamics. Recent models of ion drag on these planets, however, are based on thermal ionization only and do not include the upper atmosphere above the 1 mbar level. These models are also based on simplified equations of resistive MHD that are not always valid in extrasolar planet atmospheres. We show that photoionization dominates over thermal ionization over much of the dayside atmosphere above the 100 mbar level, creating an upper ionosphere dominated by ionization of H and He and a lower ionosphere dominated by ionization of metals such as Na, K, and Mg. The resulting dayside electron densities on close-in exoplanets are higher than those encountered in any planetary ionosphere of the solar system, and the conductivities are comparable to the chromosphere of the Sun. Based on these results and assumed magnetic fields, we constrain the conductivity regimes on close-in EGPs and use a generalized Ohm's law to study the basic effects of electrodynamics in their atmospheres. We find that ion drag is important above the 10 mbar level where it can also significantly alter the energy balance through resistive heating. Due to frequent collisions of the electrons and ions with the neutral atmosphere, however, ion drag is largely negligible in the lower atmosphere below the 10 mbar level for a reasonable range of planetary magnetic moments [abridged].

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.

Turbulence is ubiquitous in astrophysical fluids such as the interstellar medium (ISM) and the intracluster medium (ICM). In turbulence studies, it is customary to assume that fluid is driven on a single scale. However, in astrophysical fluids, there can be many different driving mechanisms that act on different scales. If there are multiple energy-injection scales, the process of energy cascade and turbulence dynamo will be different compared with the case of single energy-injection scale. In this work, we perform three-dimensional incompressible/compressible magnetohydrodynamic (MHD) turbulence simulations. We drive turbulence in Fourier space in two wavenumber ranges, 2\$\leq k \leq \sqrt{12}\$ (large-scale) and 15 \$\lesssim k \lesssim \$ 26 (small-scale). We inject different amount of energy in each range by changing the amplitude of forcing in the range. We present the time evolution of the kinetic and magnetic energy densities and discuss the turbulence dynamo in the presence of energy injections at two scales. We show how kinetic, magnetic and density spectra are affected by the two-scale energy injections and we discuss the observational implications. In the case \$\epsilon_L < \epsilon_S\$, where \$\epsilon_L\$ and \$\epsilon_S\$ are energy-injection rates at the large and small scales, respectively, our results show that even a tiny amount of large-scale energy injection can significantly change the properties of turbulence. On the other hand, when \$\epsilon_L \gtrsim \epsilon_S\$, the small-scale driving does not influence the turbulence statistics much unless \$\epsilon_L \sim \epsilon_S\$.

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ different numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and `cubed-sphere' grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should -- except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably different. When subject to a typical `hot-Jupiter'-like forcing, all five models show quantitatively different behavior -- although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

Color distributions of globular clusters (GCs) in most massive galaxies are bimodal. Assuming linear color-to-metallicity conversions, bimodality is viewed as the presence of merely two GC subsystems with distinct metallicities, which serves as a critical backbone of various galaxy formation theories. Recent studies, however, revealed that the color-metallicity relations (CMRs) often used to derive GC metallicities (e.g., CMRs of g-z, V-I and C-T1) are in fact inflected. Such inflection can create bimodal color distributions if the underlying GC metallicity spread is simply broad as expected from the hierarchical merging paradigm of galaxy formation. In order to test the nonlinear-CMR scenario for GC color bimodality, the u-band photometry is proposed because the u-related CMRs (e.g., CMRs of u-g and u-z) are theoretically predicted to be least inflected and most distinctive among commonly used optical CMRs. Here, we present Hubble Space Telescope (HST)/WFC3 F336W (u-band) photometry of the GC system in M84, a giant elliptical in the Virgo galaxy cluster. Combining the u data with the existing HST ACS/WFC g and z data, we find that the u-z and u-g color distributions are different from the g-z distribution in a very systematic manner and remarkably consistent with our model predictions based on the nonlinear-CMR hypothesis. The results lend further confidence to validity of the nonlinear-CMR scenario as an explanation for GC color bimodality. There are some GC systems showing bimodal spectroscopic metallicity, and in such systems the inflected CMRs often create stronger bimodality in the color domain.

Data from the Galileo Probe, collected during its descent into Jupiter's atmosphere, is used to obtain a vertical profile of the zonal wind from $\mathbf{\sim 0.5}$ bar (upper troposphere) to $\mathbf{\sim 0.1\, \mu{bar}}$ (lower thermosphere) at the probe entry site. This is accomplished by constructing a map of gravity wave Lomb-Scargle periodograms as a function of altitude. The profile obtained from the map indicates that the wind speed above the visible cloud deck increases with height to $\mathbf{\sim 150}$ m\,s$\mathbf{^{-1}}$ and then levels off at this value over a broad altitude range. The location of the turbopause, as a region of wide wave spectrum, is also identified from the map. In addition, a cross-equatorial oscillation of a jet, which has previously been linked to the quasi-quadrennial oscillation in the stratosphere, is suggested by the profile.

Turbulence is believed to play important roles in the origin of cosmic magnetism. While it is well known that turbulence can efficiently amplify a uniform or spatially homogeneous seed magnetic field, it is not clear whether or not we can draw a similar conclusion for a localized seed magnetic field. The main uncertainty is the rate of magnetic field diffusion on scales larger than the outer scale of turbulence. To measure the diffusion rate of magnetic field on those large scales, we perform a numerical simulation in which the outer scale of turbulence is much smaller than the size of the system. We numerically compare diffusion of a localized seed magnetic field and a localized passive scalar. We find that diffusion of the magnetic field can be much faster than that of the passive scalar and that turbulence can efficiently amplify the localized seed magnetic field. Based on the simulation result, we construct a model for fast diffusion of magnetic field. Our model suggests that a localized seed magnetic field can fill the whole system in (L_sys/L) times the large-eddy turnover time and that growth of the magnetic field stops in max(15, L_sys/L) times the large-eddy turnover time, where L_sys is the size of the system and $L$ is the driving scale. Our finding implies that, regardless of the shape of the seed field, fast magnetization is possible in turbulent systems, such as large-scale structure of the universe or galaxies.

Turbulence dynamo deals with amplification of a seed magnetic field in a turbulent medium and has been studied mostly for uniform or spatially homogeneous seed magnetic fields. However, some astrophysical processes (e.g. jets from active galaxies, galactic winds, or ram-pressure stripping in galaxy clusters) can provide localized seed magnetic fields. In this paper, we numerically study amplification of localized seed magnetic fields in a turbulent medium. Throughout the paper, we assume that driving scale of turbulence is comparable to the size of the system. Our findings are as follows. First, turbulence can amplify a localized seed magnetic field very efficiently. The growth rate of magnetic energy density is as high as that for a uniform seed magnetic field. This result implies that a magnetic field ejected from an astrophysical object can be a viable source of magnetic field in a cluster. Second, the localized seed magnetic field disperses and fills the whole system very fast. If turbulence in a system (e.g. a galaxy cluster or a filament) is driven at large scales, we expect that it takes a few large-eddy turnover times for magnetic field to fill the whole system. Third, growth and turbulence diffusion of a localized seed magnetic field are also fast in high magnetic Prandtl number turbulence. Fourth, even in decaying turbulence, a localized seed magnetic field can ultimately fill the whole system. Although the dispersal rate of magnetic field is not fast in purely decaying turbulence, it can be enhanced by an additional forcing.

We investigate baroclinic instability in flow conditions relevant to hot extrasolar planets. The instability is important for transporting and mixing heat, as well as for influencing large-scale variability on the planets. Both linear normal mode analysis and non-linear initial value calculations are carried out -- focusing on the freely-evolving, adiabatic situation. Using a high-resolution general circulation model (GCM) which solves the traditional primitive equations, we show that large-scale jets similar to those observed in current GCM simulations of hot extrasolar giant planets are likely to be baroclinically unstable on a timescale of few to few tens of planetary rotations, generating cyclones and anticyclones that drive weather systems. The growth rate and scale of the most unstable mode obtained in the linear analysis are in qualitative, good agreement with the full non-linear calculations. In general, unstable jets evolve differently depending on their signs (eastward or westward), due to the change in sign of the jet curvature. For jets located at or near the equator, instability is strong at the flanks -- but not at the core. Crucially, the instability is either poorly or not at all captured in simulations with low resolution and/or high artificial viscosity. Hence, the instability has not been observed or emphasized in past circulation studies of hot extrasolar planets.

Deep images of 10 early-type galaxies in low-density environments have been obtained with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The global properties of the globular cluster (GC) systems of the galaxies have been derived in order to investigate the role of the environment in galaxy formation and evolution. Using the ACS Virgo Cluster Survey (ACSVCS) as a high-density counterpart, the similarities and differences between the GC properties in high- and low-density environments are presented. We find a strong correlation of the GC mean colours and the degree of colour bimodality with the host galaxy luminosity in low-density environments, in good agreement with high-density environments. In contrast, the GC mean colours at a given host luminosity are somewhat bluer (\Delta(g-z) ~ 0.05) than those for cluster galaxies, indicating more metal-poor (\Delta[Fe/H] ~ 0.10-0.15) and/or younger (\Delta age > 2 Gyr) GC systems than those in dense environments. Furthermore, with decreasing host luminosity, the colour bimodality disappears faster, when compared to galaxies in cluster environments. Our results suggest that: (1) in both high- and low-density environments, the mass of the host galaxy has the dominant effect on GC system properties, (2) the local environment has only a secondary effect on the history of GC system formation, (3) GC formation must be governed by common physical processes across a range of environments.

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.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

One of the conundrums in extragalactic astronomy is the discrepancy in observed metallicity distribution functions (MDFs) between the two prime stellar components of early-type galaxies-globular clusters (GCs) and halo field stars. This is generally taken as evidence of highly decoupled evolutionary histories between GC systems and their parent galaxies. Here we show, however, that new developments in linking the observed GC colors to their intrinsic metallicities suggest nonlinear color-to-metallicity conversions, which translate observed color distributions into strongly-peaked, unimodal MDFs with broad metal-poor tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved field stars in nearby elliptical galaxies and those produced by chemical evolution models of galaxies. The GC MDF shape, characterized by a sharp peak with a metal-poor tail, indicates a virtually continuous chemical enrichment with a relatively short timescale. The characteristic shape emerges across three orders of magnitude in the host galaxy mass, suggesting a universal process of chemical enrichment among various GC systems. Given that GCs are bluer than field stars within the same galaxy, it is plausible that the chemical enrichment processes of GCs ceased somewhat earlier than that of field stellar population, and if so, GCs preferentially trace the major, vigorous mode of star formation events in galactic formation. We further suggest a possible systematic age difference among GC systems, in that the GC systems in more luminous galaxies are older. This is consistent with the downsizing paradigm of galaxies and supports additionally the similar nature shared by GCs and field stars. Our findings suggest that GC systems and their parent galaxies have shared a more common origin than previously thought, and hence greatly simplify theories of galaxy formation.

The optical color distributions of globular clusters (GCs) in most large elliptical galaxies are bimodal. Based on the assumed linear relationship between GC colors and their metallicities, the bimodality has been taken as evidence of two GC subsystems with different metallicities in each galaxy and led to a number of theories in the context of galaxy formation. More recent observations and modeling of GCs, however, suggest that the color-metallicity relations (CMRs) are inflected, and thus colors likely trace metallicities in a nonlinear manner. The nonlinearity could produce bimodal color distributions from a broad underlying metallicity spread, even if it is unimodal. Despite the far-reaching implications, whether CMRs are nonlinear and whether the nonlinearity indeed causes the color bimodality are still open questions. Given that the spectroscopic refinement of CMRs is still very challenging, we here propose a new photometric technique to probe the possible nonlinear nature of CMRs. In essence, a color distribution of GCs is a "projected" distribution of their metallicities. Since the form of CMRs hinges on which color is used, the shape of color distributions varies depending significantly on the colors. Among other optical colors, the u-band related colors (e.g., u - g and u - z) are theoretically predicted to exhibit significantly less inflected CMRs than other preferred CMRs (e.g., for g - z). As a case study, we performed the HST/WFPC2 archival u-band photometry for the M87 (NGC 4486) GC system with confirmed color bimodality. We show that the u-band color distributions are significantly different from that of g - z, and consistent with our model predictions. With more u-band measurements, this method will support or rule out the nonlinear-CMR scenario for the origin of GC color bimodality with high confidence. [Abridged]

The gaseous media of galaxy clusters and cosmic filaments, which constitute most of the baryonic matter in the universe, is highly dynamic. It is also probably turbulent, although the turbulence properties are poorly known. The gas is highly rarefied, essentially fully ionized plasma. Observational evidence suggests intracluster media (ICMs) are magnetized at some level. There are several possible origins for ICM seed fields; the observed fields are likely the result of turbulence in the ICM. We are engaged in a simulation study designed to understand in this context how very weak initial magnetic fields evolve in driven turbulence. We find that the magnetic fields eventually evolve towards equipartition levels with the vortical, solenoidal kinetic energy in the turbulence. As they do so the topology of the field structures transition from filamentary forms into ribbon-like structures in which the field orientations are laminated with vorticity structures.

Electron magnetohydrodynamics (EMHD) provides a fluid-like description of small-scale magnetized plasmas. An EMHD wave (also known as whistler wave) propagates along magnetic field lines. The direction of propagation can be either parallel or anti-parallel to the magnetic field lines. We numerically study propagation of 3-dimensional (3D) EMHD wave packets moving in one direction. We obtain two major results: 1. Unlike its magnetohydrodynamic (MHD) counterpart, an EMHD wave packet is dispersive. Because of this, EMHD wave packets traveling in one direction create opposite traveling wave packets via self-interaction and cascade energy to smaller scales. 2. EMHD wave packets traveling in one direction clearly exhibit inverse energy cascade. We find that the latter is due to conservation of magnetic helicity. We compare inverse energy cascade in 3D EMHD turbulence and 2-dimensional (2D) hydrodynamic turbulence.

We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and $8~\mu$m obtained with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, HST and ground-based $V, I, H$ and $K_s$ published observations, the range $0.5-10~\mu$m can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the dataset. Representative climate models were calculated by using a three-dimensional, pseudo-spectral general circulation model with idealised thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio calculated, linelist for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water and other molecules. No clear evidence of carbon monoxide and dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesised to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence of abundant methane, but further spectroscopic data are desirable to confirm this scenario.

Cluster media are dynamical, not static; observational evidence suggests they are turbulent. High-resolution simulations of the intracluster media (ICMs) and of idealized, similar media help us understand the complex physics and astrophysics involved. We present a brief overview of the physics behind ICM turbulence and outline the processes that control its development. High-resolution, compressible, isothermal MHD simulations are used to illustrate important dynamical properties of turbulence that develops in media with initially very weak magnetic fields. The simulations follow the growth of magnetic fields and reproduce the characteristics of turbulence. These results are also compared with full cluster simulations that have examined the properties of ICM turbulence.

We elucidate the interplay between Newtonian thermal relaxation and numerical dissipation, of several different origins, in flow simulations of hot extrasolar planet atmospheres. Currently, a large range of Newtonian relaxation, or "cooling", times (~10 days to ~1 hour) is used among different models and within a single model over the model domain. In this study we demonstrate that a short relaxation time (much less than the planetary rotation time) leads to a large amount of unphysical, grid-scale oscillations that contaminate the flow field. These oscillations force the use of an excessive amount of artificial viscosity to quench them and prevent the simulation from "blowing up". Even if the blow-up is prevented, such simulations can be highly inaccurate because they are either severely over-dissipated or under-dissipated, and are best discarded in these cases. Other numerical stability and timestep size enhancers (e.g., Robert-Asselin filter or semi-implicit time-marching schemes) also produce similar, but less excessive, damping. We present diagnostics procedures to choose the "optimal" simulation and discuss implications of our findings for modeling hot extrasolar planet atmospheres.

We report the detection of 24 micron variations from the planet-hosting upsilon Andromedae system consistent with the orbital periodicity of the system's innermost planet, upsilon And b. We find a peak-to-valley phase curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of constant surface brightness and assuming a planetary radius of 1.3 Jupiter radii gives a planetary temperature contrast of >900 K and an orbital inclination of >28 degrees. We further report the largest phase offset yet observed for an extrasolar planet: the flux maximum occurs ~80 degrees before phase 0.5. Such a large phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on earlier observations of this system in several important ways: (1) observations of a flux calibrator star demonstrate the MIPS detector is stable to 10^-4 on long timescales, (2) we note that the background light varies systematically due to spacecraft operations, precluding use of this background as a flux calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier observations of this system is consistent with our new result.

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.

We compare non-locality of interactions between different scales in hydrodynamic (HD) turbulence and magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium. We use 3-dimensional incompressible direct numerical simulations to evaluate non-locality of interactions. Our results show that non-locality in MHD turbulence is much more pronounced than that in HD turbulence. Roughly speaking, non-local interactions count for more than 10\% of total interactions in our MHD simulation on a grid of $512^3$ points. However, there is no evidence that non-local interactions are important in our HD simulation with the same numerical resolution. We briefly discuss how non-locality affects energy spectrum.

We use a general circulation model to study the three-dimensional (3-D) flow and temperature distributions of atmospheres on tidally synchronized extrasolar planets. In this work, we focus on the sensitivity of the evolution to the initial flow state, which has not received much attention in 3-D modeling studies. We find that different initial states lead to markedly different distributions-even under the application of strong forcing (large day-night temperature difference with a short "thermal drag time") that may be representative of close-in planets. This is in contrast with the results or assumptions of many published studies. In general, coherent jets and vortices (and their associated temperature distributions) characterize the flow, and they evolve differently in time, depending on the initial condition. If the coherent structures reach a quasi- stationary state, their spatial locations still vary. The result underlines the fact that circulation models are currently unsuitable for making quantitative predictions (e.g., location and size of a "hot spot") without better constrained, and well posed, initial conditions.

We study the effects of gravity waves, or g-modes, on hot extrasolar planets. These planets are expected to possess stably-stratified atmospheres, which support gravity waves. In this paper, we review the derivation of the equation that governs the linear dynamics of gravity waves and describe its application to a hot extrasolar planet, using HD209458 b as a generic example. We find that gravity waves can exhibit a wide range of behaviors, even for a single atmospheric profile. The waves can significantly accelerate or decelerate the background mean flow, depending on the difference between the wave phase and mean flow speeds. In addition, the waves can provide significant heating (~100 to ~1000 K per planetary rotation), especially to the region of the atmosphere above about 10 scale heights from the excitation region. Furthermore, by propagating horizontally, gravity waves provide a mechanism for transporting momentum and heat from the dayside of a tidally locked planet to its nightside. We discuss work that needs to be undertaken to incorporate these effects in current atmosphere models of extrasolar planets.

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 survey the basic principles of atmospheric dynamics relevant to explaining existing and future observations of exoplanets, both gas giant and terrestrial. Given the paucity of data on exoplanet atmospheres, our approach is to emphasize fundamental principles and insights gained from Solar-System studies that are likely to be generalizable to exoplanets. We begin by presenting the hierarchy of basic equations used in atmospheric dynamics, including the Navier-Stokes, primitive, shallow-water, and two-dimensional nondivergent models. We then survey key concepts in atmospheric dynamics, including the importance of planetary rotation, the concept of balance, and scaling arguments to show how turbulent interactions generally produce large-scale east-west banding on rotating planets. We next turn to issues specific to giant planets, including their expected interior and atmospheric thermal structures, the implications for their wind patterns, and mechanisms to pump their east-west jets. Hot Jupiter atmospheric dynamics are given particular attention, as these close-in planets have been the subject of most of the concrete developments in the study of exoplanetary atmospheres. We then turn to the basic elements of circulation on terrestrial planets as inferred from Solar-System studies, including Hadley cells, jet streams, processes that govern the large-scale horizontal temperature contrasts, and climate, and we discuss how these insights may apply to terrestrial exoplanets. Although exoplanets surely possess a greater diversity of circulation regimes than seen on the planets in our Solar System, our guiding philosophy is that the multi-decade study of Solar-System planets reviewed here provides a foundation upon which our understanding of more exotic exoplanetary meteorology must build.

In the framework of turbulence dynamo, flow motions amplify a weak seed magnetic field through the stretching of field lines. Although the amplification process has been a topic of active research, less attention has been paid to the length scales of magnetic field. In this paper, we described a numerical study on characteristic lengths of magnetic field in magnetohydrodynamic turbulence. We considered the case of very weak or zero mean magnetic field, which is applicable to the turbulence in the intergalactic space. Our findings are as follows. (1) At saturation, the peak of magnetic field spectrum occurs at $\sim L_0/2$, where $L_0$ is the energy injection scale, while the most energy containing scale is $\sim L_0/5$. The peak scale of spectrum of projected, two-dimensional field is $\sim L_0$. (2) During the stage of magnetic field amplification, the energy equipartition scale shows a power-law increase of $\sim t^{1.5}$, while the integral and curvature scales show a linear increase. The equipartition, integral, and curvature scales saturate at $\sim L_0$, $\sim 0.3L_0$, and $\sim 0.15L_0$, respectively. (3) The coherence length of magnetic field defined in the Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) is related to the integral scale. We presented a formula that expresses the standard deviation of RM, $\sigma_{RM}$, in terms of the integral scale and rms strength of the IGMF, and estimated that $\sigma_{RM}$ would be $\sim 100$ and $\sim$ a few rad m$^{-2}$ for clusters and filaments, respectively.

We studied Faraday rotation measure (RM) in turbulent media with the rms Mach number of unity, using isothermal, magnetohydrodynamic turbulence simulations. Four cases with different values of initial plasma beta were considered. Our main findings are as follows. (1) There is no strong correlation between the fluctuations of magnetic field strength and gas density. So the magnetic field strength estimated with RM/DM (DM is the dispersion measure) correctly represents the true mean strength of the magnetic field along the line of sight. (2) The frequency distribution of RMs is well fitted to the Gaussian. In addition, there is a good relation between the width of the distribution of RM/$\bar{\rm RM}$ ($\bar{\rm RM}$ is the average value of RMs) and the strength of the regular field along the line of sight; the width is narrower, if the field strength is stronger. We discussed the implications of our findings in the warm ionized medium where the Mach number of turbulent motions is around unity.

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 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 present the first 350 micron polarization measurement for the disk of the T Tauri star (TTS) DG Tau. The data were obtained using the SHARP polarimeter at the Caltech Submillimeter Observatory. We did not detect polarization but we provide constraints on the normalized Stokes parameters q and u. We derive information about the polarization spectrum by comparing our 350 micron result with an 850 micron polarization detection previously published for this source. This empirical information concerning the polarization spectrum disagrees with the predictions of a recent model for TTS disk polarization. We conclude, however, that adding more mass to the model disk would probably result in model polarization spectra that agree much better with the 350 and 850 micron polarimetric observations. We suggest that multiwavelength polarimetry of TTS disk emission may provide a promising method for probing the opacity of TTS disks.

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.

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 study internal extinction of late-type galaxies in the Sloan Digital Sky Survey. We find that the degree of internal extinction depends on both the concentration index c and K_s-band absolute magnitude M_K. We give simple fitting functions for internal extinction. In particular, we present analytic formulae giving the extinction-corrected magnitudes from the observed optical parameters. For example, the extinction-corrected r-band absolute magnitude can be obtained by M_{r,0}=-20.77 +(-1+\sqrt{1+4\Delta (M_{r,obs}+20.77+4.93\Delta)})/2\Delta, where \Delta = 0.236\{1.35(c-2.48)^2-1.14\}\log(a/b), c=R_{90}/R_{50} is the the concentration index, and a/b is the isophotal axis ratio of the 25 mag/arcsec^2 isophote in the i-band. The 1\sigma error in M_{r,0} is 0.21{\rm log}(a/b). The late-type galaxies with very different inclinations are found to trace almost the same sequence in the (u-r)-$M_r$ diagram when our prescriptions for extinction correction are applied. We also find that (u-r) color can be a third independent parameter that determines the degree of internal extinction.

The physical characterization of exoplanets will require to take spectra at several orbital positions. For that purpose, a direct imaging capability is necessary. Direct imaging requires an efficient stellar suppression mechanism, associated with an ultrasmooth telescope. We show that before future large space missions (interferometer, 4-8 m class coronograph, external occulter or Fresnel imager), direct imaging of giant planets and close-by super-Earth are at the cross-road of a high scientific interest and a reasonable feasibility. The scientific interest lies in the fact that super-Earths share common geophysical attributes with Earths. They already begin to be detected by radial velocity (RV) and, together with giant planets, they have a larger area than Earths, making them detectable with a 1.5-2 m class telescope in reflected light. We propose such a (space) telescope be a first step before large direct imaging missions.

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.

The Super-Earth Explorer is an Off-Axis Space Telescope (SEE-COAST) designed for high contrast imaging. Its scientific objective is to make the physico-chemical characterization of exoplanets possibly down to 2 Earth radii >. For that purpose it will analyze the spectral and polarimetric properties of the parent starlight reflected by the planets, in the wavelength range 400-1250 nm

We present the 3.5m SPICA space telescope, a proposed Japanese-led JAXA-ESA mission scheduled for launch around 2017. The actively cooled (<5 K), single aperture telescope and monolithic mirror will operate from ~3.5 to ~210 um and will provide superb sensitivity in the mid- and far-IR spectral domain (better than JWST at lambda > 18 um). SPICA is one of the few space missions selected to go to the next stage of ESA's Cosmic Vision 2015-2025 selection process. In this White Paper we present the main specifications of the three instruments currently baselined for SPICA: a mid-infrared (MIR) coronagraph (~3.5 to ~27 um) with photometric and spectral capabilities (R~200), a MIR wide-field camera and high resolution spectrometer (R~30,000), and a far-infrared (FIR ~30 to ~210 um) imaging spectrometer - SAFARI - led by a European consortium. We discuss their capabilities in the context of MIR direct observations of exo-planets (EPs) and multiband photometry/high resolution spectroscopy observations of transiting exo-planets. We conclude that SPICA will be able to characterize the atmospheres of transiting exo-planets down to the super-Earth size previously detected by ground- or space-based observatories. It will also directly detect and characterize Jupiter/Neptune-size planets orbiting at larger separation from their parent star (>5-10 AU), by performing quantitative atmospheric spectroscopy and studying proto-planetary and debris disks. In addition, SPICA will be a scientific and technological precursor for future, more ambitious, IR space missions for exo-planet direct detection as it will, for example, quantify the prevalence exo-zodiacal clouds in planetary systems and test coronographic techniques, cryogenic systems and lightweight, high quality telescopes. (abridged)

The nature and origin of turbulence and magnetic fields in the intergalactic space are important problems that are yet to be understood. We propose a scenario in which turbulent flow motions are induced via the cascade of the vorticity generated at cosmological shocks during the formation of the large scale structure. The turbulence in turn amplifies weak seed magnetic fields of any origin. Supercomputer simulations show that the turbulence is subsonic inside clusters/groups of galaxies, whereas it is transonic or mildly supersonic in filaments. Based on a turbulence dynamo model, we then estimate that the average magnetic field strength would be a few microgauss inside clusters/groups, approximately 0.1 microgauss around clusters/groups, and approximately 10 nanogauss in filaments. Our model presents a physical mechanism that transfers the gravitation energy to the turbulence and magnetic field energies in the large scale structure of the universe.

If ultra-high-energy cosmic rays (UHECRs) originate from extragalactic sources, understanding the propagation of charged particles through the magnetized large scale structure (LSS) of the universe is crucial in the search for the astrophysical accelerators. Based on a novel model of the turbulence dynamo, we estimate the intergalactic magnetic fields (IGMFs) in cosmological simulations of the formation of the LSS. Under the premise that the sources of UHECRs are strongly associated with the LSS, we consider a model in which protons with E >10^{19} eV are injected by sources that represent active galactic nuclei located inside clusters of galaxies. With the model IGMFs, we then follow the trajectories of the protons, while taking into account the energy losses due to interactions with the cosmic background radiation. For observers located inside groups of galaxies like ours, about 70% and 35% of UHECR events above 60 EeV arrive within ~15 degree and ~5 degree, respectively, of the source position with time delays of less than ~10^7 yr. This implies that the arrival direction of super-GZK protons might exhibit a correlation with the distribution of cosmological sources on the sky. In this model, nearby sources (within 10 - 20 Mpc) should contribute significantly to the particle flux above ~10^{20} eV.

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.

Turbulence is ubiquitous in Solar System planetary atmospheres. In hot Jupiter atmospheres, the combination of moderately slow rotation and thick pressure scale height may result in dynamical weather structures with unusually large, planetary-size scales. Using equivalent-barotropic, turbulent circulation models, we illustrate how such structures can generate a variety of features in the thermal phase curves of hot Jupiters, including phase shifts and deviations from periodicity. Such features may have been spotted in the recent infrared phase curve of HD 189733b. Despite inherent difficulties with the interpretation of disk-integrated quantities, phase curves promise to offer unique constraints on the nature of the circulation regime present on hot Jupiters.

We present results from a set of over 300 pseudospectral simulations of atmospheric circulation on extrasolar giant planets with circular orbits. The simulations are of high enough resolution (up to 341 total and sectoral modes) to resolve small-scale eddies and waves, required for reasonable physical accuracy. In this work, we focus on the global circulation pattern that emerges in a shallow, ``equivalent-barotropic'', turbulent atmosphere on both tidally synchronized and unsynchronized planets. A full exploration of the large physical and numerical parameter-space is performed to identify robust features of the circulation. For some validation, the model is first applied to Solar System giant planets. For extrasolar giant planets with physical parameters similar to HD209458b--a presumably synchronized extrasolar giant planet, representative in many dynamical respects--the circulation is characterized by the following features: 1) a coherent polar vortex that revolves around the pole in each hemisphere; 2) a low number--typically two or three--of slowly-varying, broad zonal (east-west) jets that form when the maximum jet speed is comparable to, or somewhat stronger than, those observed on the planets in the Solar System; and, 3) motion-associated temperature field, whose detectability and variability depend on the strength of the net heating rate and the global root mean square wind speed in the atmosphere. In many ways, the global circulation is Earth-like, rather than Jupiter-like. However, if extrasolar giant planets rotate faster and are not close-in (therefore not synchronized), their circulations become more Jupiter-like, for Jupiter-like rotation rates.

Hot Jupiters are new laboratories for the physics of giant planet atmospheres. Subject to unusual forcing conditions, the circulation regime on these planets may be unlike anything known in the Solar System. Characterizing the atmospheric circulation of hot Jupiters is necessary for reliable interpretation of the multifaceted data currently being collected on these planets. We discuss several fundamental concepts of atmospheric dynamics that are likely central to obtaining a solid understanding of these fascinating atmospheres. A particular effort is made to compare the various modeling approaches employed so far to address this challenging problem.

The propagation of UHECRs is affected by the intergalactic magnetic fields that were produced during the course of the large scale structure formation of the universe. We adopt a novel model where the large scale extragalactic magnetic fields (EGMF) are estimated from local dynamic properties of the gas flows in hydrodynamic simulations of a concordance LambdaCDM universe. With the model magnetic fields, we calculate the deflection angle, time delay and energy spectrum of protons with E > 10^{19} eV that are injected at cosmological sources and then travel through the large scale structure of the universe, losing the energy due to interactions with the cosmic background radiation. Implications of this study on the origin of UHECRs are discussed.

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

Physical conditions in the atmospheres of tidally-locked, slowly-rotating hot Jupiters correspond to dynamical circulation regimes with Rhines scales and Rossby deformation radii comparable to the planetary radii. Consequently, the large spatial scales of moving atmospheric structures could generate significant photospheric variability. Here, we estimate the level of thermal infrared variability expected in successive secondary eclipse depths, according to hot Jupiter turbulent ``shallow-layer'' models. The variability, at the few percent level or more in models with strong enough winds, is within the reach of Spitzer measurements. Eclipse depth variability is thus a valuable tool to constrain the circulation regime and global wind speeds in hot Jupiter atmospheres.

Recent Spitzer infrared measurements of hot Jupiter eclipses suggest that eclipse mapping techniques could be used to spatially resolve the day-side photospheric emission of these planets using partial occultations. As a first step in this direction, we simulate ingress/egress lightcurves for the three brightest known eclipsing hot Jupiters and evaluate the degree to which parameterized photospheric emission models can be distinguished from each other with repeated, noisy eclipse measurements. We find that the photometric accuracy of Spitzer is insufficient to use this tool effectively. On the other hand, the level of photospheric details that could be probed with a few JWST eclipse measurements could greatly inform hot Jupiter atmospheric modeling efforts. A JWST program focused on non-parametric eclipse map inversions for hot Jupiters should be actively considered.

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)

The star upsilon Andromeda is orbited by three known planets, the innermost of which has an orbital period of 4.617 days and a mass at least 0.69 that of Jupiter. This planet is close enough to its host star that the radiation it absorbs overwhelms its internal heat losses. Here we present the 24 micron light curve of this system, obtained with the Spitzer Space Telescope. It shows a clear variation in phase with the orbital motion of the innermost planet. This is the first demonstration that such planets possess distinct hot substellar (day) and cold antistellar (night) faces.

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.

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.

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.

We discuss atmosphere models of HD209458b in light of the recent day-side flux measurement of HD209458b's secondary eclipse by Spitzer-MIPS at 24 microns. In addition, we present a revised secondary eclipse IRTF upper limit at 2.2 microns which places a stringent constraint on the adjacent H2O absorption band depths. These two measurements are complementary because they are both shaped by H2O absorption and because the former is on the Wien tail of the planet's thermal emission spectrum and the latter is near the thermal emission peak. A wide range of models fit the observational data, confirming our basic understanding of hot Jupiter atmospheric physics. Although a range of models are viable, some models at the hot and cold end of the plausible temperature range can be ruled out. One class of previously unconsidered hot Jupiter atmospheric models that fit the data are those with C/O >~ 1 (as Jupiter may have), which have a significant paucity of H2O compared to solar abundance models with C/O = 0.5. The models indicate that HD209458b is in a situation intermediate between pure in situ reradiation and very efficient redistribution of heat; one which will require a careful treatment of atmospheric circulation. We discuss how future wavelength-dependent and phase-dependent observations will further constrain the atmospheric circulation regime. In the shorter term, additional planned measurements for HD209458b, especially Spitzer IRAC photometry, should lift many of the model degeneracies. Multiwavelength IR observations constrain the atmospheric structure and circulation properties of hot Jupiters and thus open a new chapter in quantitative extrasolar planetology.

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

We present numerical studies of 3-dimensional magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium in the extremely relativistic limit, in which the inertia of the charge carriers can be neglected. We have focused on strong Alfvenic turbulence in the limit. We have found the following results. First, the energy spectrum is consistent with a Kolmogorov spectrum: $E(k)\sim k^{-5/3}$. Second, turbulence shows a Goldreich-Sridhar type anisotropy: $k_{\|} \propto k_{\perp}^{2/3}$, where $k_{\|}$ and $k_{\perp}$ are wavenumbers along and perpendicular to the local mean magnetic field directions, respectively. These scalings are in agreement with earlier theoretical predictions by Thompson & Blaes.

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.

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

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.

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

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.

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.

Shallow-water numerical simulations show that the atmospheric circulation of the close-in extrasolar giant planet (EGP) HD 209458b is characterized by moving circumpolar vortices and few bands/jets (in contrast with ~10 bands/jets and absence of polar vortices on cloud-top Jupiter and Saturn). The large spatial scales of moving circulation structures on HD 209458b may generate detectable variability of the planet's atmospheric signatures. In this Letter, we generalize these results to other close-in EGPs, by noting that shallow-water dynamics is essentially specified by the values of the Rossby (Ro) and Burger (Bu) dimensionless numbers. The range of likely values of Ro (~ 0.01 - 10) and Bu (~ 1 - 200) for the atmospheric flow of known close-in EGPs indicates that their circulation should be qualitatively similar to that of HD 209458b. This results mostly from the slow rotation of these tidally-synchronized planets.

High-resolution atmospheric flow simulations of the tidally-locked extrasolar giant planet, HD 209458b, show large-scale spatio-temporal variability. This is in contrast to the simple, permanent day/night (i.e., hot/cold) picture. The planet's global circulation is characterized by a polar vortex in motion around each pole and a banded structure corresponding to ~3 broad zonal (east-west) jets. For very strong jets, the circulation-induced temperature difference between moving hot and cold regions can reach up to ~1000 K, suggesting that atmospheric variability could be observed in the planet's spectral and photometric signatures.

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.

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.

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.

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.

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.

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

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.

The study of incompressible magnetohydrodynamic (MHD) turbulence gives useful insights on many astrophysical problems. We describe a pseudo-spectral MHD code suitable for the study of incompressible turbulence. We review our recent works on direct three-dimensional numerical simulations for MHD turbulence in a periodic box. In those works, we use a pseudo-spectral code to solve the incompressible MHD equations. We first discuss the structure and properties of turbulence as functions of scale. The results are consistent with the scaling law recently proposed by Goldreich & Sridhar. The scaling law is based on the concept of scale-dependent isotropy: smaller eddies are more elongated than larger ones along magnetic field lines. This scaling law substantially changes our views on MHD turbulence. For example, as noted by Lazarian & Vishniac, the scaling law can provide a fast reconnection rate. We further discuss how the study of incompressible MHD turbulence can help us to understand physical processes in interstellar medium (ISM) by considering imbalanced cascade and viscous damped turbulence.

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 construct a magnetic helicity conserving dynamo theory which incorporates a calculated magnetic helicity current. In this model the fluid helicity plays a small role in large scale magnetic field generation. Instead, the dynamo process is dominated by a new quantity, derived from asymmetries in the second derivative of the velocity correlation function, closely related to the `twist and fold' dynamo model. The turbulent damping term is, as expected, almost unchanged. Numerical simulations with a spatially constant fluid helicity and vanishing resistivity are not expected to generate large scale fields in equipartition with the turbulent energy density. The prospects for driving a fast dynamo under these circumstances are uncertain, but if it is possible, then the field must be largely force-free. On the other hand, there is an efficient analog to the $\alpha-\Omega$ dynamo. Systems whose turbulence is driven by some anisotropic local instability in a shearing flow, like real stars and accretion disks, and some computer simulations, may successfully drive the generation of strong large scale magnetic fields, provided that $\partial_r\Omega< \partial_\theta v_z\omega_\theta>>0$. We show that this criterion is usually satisfied. Such dynamos will include a persistent, spatially coherent vertical magnetic helicity current with the same sign as $-\partial_r\Omega$, that is, positive for an accretion disk and negative for the Sun. We comment on the role of random magnetic helicity currents in storing turbulent energy in a disordered magnetic field, which will generate an equipartition, disordered field in a turbulent medium, and also a declining long wavelength tail to the power spectrum. As a result, calculations of the galactic `seed' field are largely irrelevant.

We perform direct 3-dimensional numerical simulations for magnetohydrodynamic (MHD) turbulence in a periodic box of size $2\pi$ threaded by strong uniform magnetic fields. We use a pseudo-spectral code with hyperviscosity and hyperdiffusivity to solve the incompressible MHD equations. We analyze the structure of the eddies as a function of scale. A straightforward calculation of anisotropy in wavevector space shows that the anisotropy is scale-{\it independent}. We discuss why this is {\it not} the true scaling law and how the curvature of large-scale magnetic fields affects the power spectrum and leads to the wrong conclusion. When we correct for this effect, we find that the anisotropy of eddies depends on their size: smaller eddies are more elongated than larger ones along {\it local} magnetic field lines. The results are consistent with the scaling law $\tilde{k}_{\parallel} \sim \tilde{k}_{\perp}^{2/3}$ proposed by Goldreich and Sridhar (1995, 1997). Here $\tilde{k}_{\|}$ (and $\tilde{k}_{\perp}$) are wavenumbers measured relative to the local magnetic field direction. However, we see some systematic deviations which may be a sign of limitations to the model, or our inability to fully resolve the inertial range of turbulence in our simulations.

We have tested the ability of driven turbulence to generate magnetic field structure from a weak uniform field using three dimensional numerical simulations of incompressible turbulence. We used a pseudo-spectral code with a numerical resolution of up to $144^3$ collocation points. We find that the magnetic fields are amplified through field line stretching at a rate proportional to the difference between the velocity and the magnetic field strength times a constant. Equipartition between the kinetic and magnetic energy densities occurs at a scale somewhat smaller than the kinetic energy peak. Above the equipartition scale the velocity structure is, as expected, nearly isotropic. The magnetic field structure at these scales is uncertain, but the field correlation function is very weak. At the equipartition scale the magnetic fields show only a moderate degree of anisotropy, so that the typical radius of curvature of field lines is comparable to the typical perpendicular scale for field reversal. In other words, there are few field reversals within eddies at the equipartition scale, and no fine-grained series of reversals at smaller scales. At scales below the equipartition scale, both velocity and magnetic structures are anisotropic; the eddies are stretched along the local magnetic field lines, and the magnetic energy dominates the kinetic energy on the same scale by a factor which increases at higher wavenumbers. We do not show a scale-free inertial range, but the power spectra are a function of resolution and/or the imposed viscosity and resistivity. Our results are consistent with the emergence of a scale-free inertial range at higher Reynolds numbers.