Source author record

Kandaswamy Subramanian

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

49 published item(s)

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy.

preprint2022arXiv

Non-locality of the Turbulent Electromotive Force

The generation of large-scale magnetic fields ($\overline{\mathbf{B}}$) in astrophysical systems is driven by the mean turbulent electromotive force ($\overline{\cal E}$), the cross correlation between local fluctuations of velocity and magnetic fields. This can depend non-locally on $\overline{\mathbf{B}}$ through a convolution kernel $K_{ij}$. In a new approach to find $K_{ij}$, we directly fit the time series data of $\overline{\cal E}$ versus $\overline{\mathbf{B}}$ from a galactic dynamo simulation using singular value decomposition. We calculate the usual turbulent transport coefficients as moments of $K_{ij}$, show the importance of including non-locality over eddy length scales to fully capture their amplitudes and that higher order corrections to the standard transport coefficients are small in the present case.

preprint2021arXiv

Extragalactic magnetism with SOFIA (Legacy Program) -- II: A magnetically-driven flow in the starburst ring of NGC 1097

Galactic bars are frequent in disk galaxies and they may support the transfer of matter towards the central engine of active nuclei. The barred galaxy NGC 1097 has magnetic forces controlling the gas flow at several kpc scales, which suggest that magnetic fields (B-fields) are dynamically important along the bar and nuclear ring. However, the effect of the B-field on the gas flows in the central kpc scale has not been characterized. Using thermal polarized emission at $89$ $μ$m with HAWC+/SOFIA, here, we measure that the polarized flux is spatially located at the contact regions of the outer-bar with the starburst ring. The linear polarization decomposition analysis shows that the $89$ $μ$m and radio ($3.5$ and $6.2$ cm) polarization traces two different modes, $m$, of the B-field: a constant B-field orientation and dominated by $m=0$ at $89$ $μ$m, and a spiral B-field dominated by $m=2$ at radio. We show that the B-field at 89 $μ$m is concentrated in the warmest region of a shock driven by the galactic-bar dynamics in the contact regions between the outer-bar with the starburst ring. Radio polarization traces a superposition of the spiral B-field outside and within the starburst ring. According to Faraday rotation measures between $3.5$ and $6.2$ cm, the radial component of the B-field along the contact regions points toward the galaxy's center on both sides. We conclude that gas streams outside and within the starburst ring follow the B-field, which feeds the black hole with matter from the host galaxy.

preprint2021arXiv

Properties of polarized synchrotron emission from fluctuation-dynamo action -- I. Application to galaxy clusters

Using magnetohydrodynamic simulations of fluctuation dynamos, we perform broad-bandwidth synthetic observations to investigate the properties of polarized synchrotron emission and the role that Faraday rotation plays in inferring the polarized structures in the intracluster medium (ICM) of galaxy clusters. In the saturated state of the dynamo, we find a Faraday depth (FD) dispersion $σ_{\rm FD} \approx 100$ rad m$^{-2}$, in agreement with observed values in the ICM. Remarkably, the FD power spectrum is qualitatively similar to $M(k)/k$, where $M(k)$ is the magnetic spectrum and $k$ the wavenumber. However, this similarity is broken at high $k$ when FD is obtained by applying RM synthesis to polarized emission from the ICM due to poor resolution and complexities of spectrum in FD space. Unlike the Gaussian probability distribution function (PDF) obtained for FD, the PDF of the synchrotron intensity is lognormal. A relatively large $σ_{\rm FD}$ in the ICM gives rise to strong frequency-dependent variations of the pixel-wise mean and peak polarized intensities at low frequencies ($\lesssim 1.5\,{\rm GHz}$). The mean fractional polarization $\langle p \rangle$ obtained at the resolution of the simulations increases from $<0.1$ at 0.5 GHz to its intrinsic value of $\sim0.3$ at 6 GHz. Beam smoothing significantly affects the polarization properties below $\lesssim 1.5\,{\rm GHz}$, reducing $\langle p \rangle$ to $\lesssim 0.01$ at 0.5 GHz. At frequencies $\gtrsim 5\,{\rm GHz}$, polarization remains largely unaffected, even when recovered using RM synthesis. Thus, our results underline the need for high-frequency ($\gtrsim 5\,{\rm GHz}$) observations with future radio telescopes to effectively probe the properties of polarized emission in the ICM.

preprint2020arXiv

Characterising the Dynamo in a Radiatively Inefficient Accretion Flow

We explore the MRI driven dynamo in a radiatively inefficient accretion flow (RIAF) using the mean field dynamo paradigm. Using singular value decomposition (SVD) we obtain the least squares fitting dynamo coefficients $α$ and $γ$ by comparing the time series of the turbulent electromotive force and the mean magnetic field. Our study is the first one to show the poloidal distribution of these dynamo coefficients in global accretion flow simulations. Surprisingly, we obtain a high value of the turbulent pumping coefficient $γ$ which transports the mean magnetic flux radially outward. This would have implications for the launching of magnetised jets which are produced efficiently in presence a large-scale poloidal magnetic field close to the compact object. We present a scenario of a truncated disc beyond the RIAF where a large scale dynamo-generated poloidal magnetic field can aid jet-launching close to the black hole. Magnitude of all the calculated coefficients decreases with radius. Meridional variations of $α_{ϕϕ}$, responsible for toroidal to poloidal field conversion, is very similar to that found in shearing box simulations using the `test field' (TF) method. By estimating the relative importance of $α$-effect and shear, we conclude that the MRI driven large-scale dynamo, which operates at high latitudes beyond a disc scale height, is essentially of the $α-Ω$ type.

preprint2019arXiv

Gravitational wave generation in a viable scenario of inflationary magnetogenesis

Generation of magnetic fields during inflation is a promising mechanism for the origin of the observed large scale magnetic fields in the universe. Among several attempts, a popular model is one where the inflaton and the electromagnetic field are coupled through a coupling function $f$ leading to a term in the Lagrangian density of the form, $f^2 F^{μν} F_{μν}$. A number of potential difficulties with such models have been raised in the literature. In our earlier work, we have suggested viable models of inflationary magnetogenesis which avoid these problems and at the same time can lead to either nonhelical or helical magnetic fields of astrophysical interest. Our models require a low energy scale for inflation and reheating (reheating temperature, $T_R < 10^4$ GeV) and generate a blue spectrum of electromagnetic (EM) field which peaks around the horizon scale of reheating. We show here that the anisotropic stress associated with these EM fields naturally source the production of a stochastic background of Gravitational waves (GW) with frequencies in the range of tens of nano Hertz to milli Hertz. These two extremes of the range can be probed respectively by pulsar timing arrays (PTA) experiments and the upcoming Laser Interferometric Space Array (LISA). The peak value of the GW spectrum energy represented by $d Ω_{GW}/d \ln k$ is $ 10^{-6} $ for the models which lead to nonhelical primordial fields and $2 \times 10^{-6} $ for the helical case for $T_R=100$ GeV. In this case the spectrum peaks at a frequency $30 μ$Hz for non helical case and at $40 μ$Hz for helical case. These values are obtained when the ratio of EM energy density to the cosmological density at reheating $ε\sim 1$ and decrease approximately as $ε^2$ for smaller values.

preprint2019arXiv

Turbulent transport coefficients in galactic dynamo simulations using singular value decomposition

Coherent magnetic fields in disc galaxies are thought to be generated by a large-scale (or mean-field) dynamo operating in their interstellar medium. A key driver of mean magnetic field growth is the turbulent electromotive force (EMF), which represents the influence of correlated small-scale (or fluctuating) velocity and magnetic fields on the mean field. The EMF is usually expressed as a linear expansion in the mean magnetic field and its derivatives, with the dynamo tensors as expansion coefficients. Here, we adopt the singular value decomposition (SVD) method to directly measure these turbulent transport coefficients in a simulation of the turbulent interstellar medium that realizes a large-scale dynamo. Specifically, the SVD is used to least-square fit the time series data of the EMF with that of the mean field and its derivatives, to determine these coefficients. We demonstrate that the spatial profiles of the EMF reconstructed from the SVD coefficients match well with that taken directly from the simulation. Also, as a direct test, we use the coefficients to simulate a 1-D mean-field dynamo model and find an overall similarity in the evolution of the mean magnetic field between the dynamo model and the direct simulation. We also compare the results with those which arise using simple regression and the ones obtained previously using the test-field (TF) method, to find reasonable qualitative agreement. Overall, the SVD method provides an effective post-processing tool to determine turbulent transport coefficients from simulations.

preprint2016arXiv

A unified large/small-scale dynamo in helical turbulence

We use high resolution direct numerical simulations (DNS) to show that helical turbulence can generate significant large-scale fields even in the presence of strong small-scale dynamo action. During the kinematic stage, the unified large/small-scale dynamo grows fields with a shape-invariant eigenfunction, with most power peaked at small scales or large $k$, as in \citet{SB14}. Nevertheless, the large-scale field can be clearly detected as an excess power at small $k$ in the negatively polarized component of the energy spectrum for a forcing with positively polarized waves. Its strength $bar{B}$, relative to the total rms field $B_{rms}$, decreases with increasing magnetic Reynolds number, $R_m$. However, as the Lorentz force becomes important, the field generated by the unified dynamo orders itself by saturating on successively larger scales. The magnetic integral scale for the positively polarized waves, characterizing the small-scale field, increases significantly from the kinematic stage to saturation. This implies that the small-scale field becomes as coherent as possible for a given forcing scale, which averts the $R_m$-dependent quenching of $\bar{B}/B_rms$. These results are obtained for $1024^3$ DNS with magnetic Prandtl numbers of $P_m=0.1$ and $10$. For $P_m=0.1$, $\bar{B}/B_{rms}$ grows from about $0.04$ to about $0.4$ at saturation, aided in the final stages by helicity dissipation. For $P_m=10$, $\bar{B}/B_{rms}$ grows from much less than 0.01 to values of the order the $0.2$. Our results confirm that there is a unified large/small-scale dynamo in helical turbulence.

preprint2016arXiv

Probing magnetic fields with Square Kilometre array and its precursors

Origin of magnetic fields, its structure and effects on dynamical processes in stars to galaxies are not well understood. Lack of a direct probe has hampered its study. The first phase of Square Kilometre Array (SKA-I), will have more than an order of magnitude higher sensitivity than existing radio telescopes. In this contribution, we discuss specific science cases that are of interest to the Indian community concerned with astrophysical turbulence and magnetic fields. The SKA-I will allow observations of a large number of background sources with detectable polarisation and measure their Faraday depths (FDs) through the Milky Way, other galaxies and their circum-galactic medium. This will probe line-of-sight magnetic fields in these objects well and provide field configurations. Detailed comparison of observational data with models which consider various processes giving rise to field amplification and maintenance will then be possible. Such observations will also provide the coherence scale of the fields and measure its random component. Measuring the random component is important to characterise turbulence in the medium. Observations of FDs with redshift will provide important information on magnetic field evolution as a function of redshift. The background sources could also be used to probe magnetic fields and its coherent scale in galaxy clusters and in bridges formed between interacting galaxies. Other than FDs, sensitive observations of synchrotron emission from galaxies will provide complimentary information on their magnetic field strengths in the sky plane. The core shift measurements of AGNs can provide more precise measurements of magnetic field very close (<pc) to the black hole and its evolution. The low band of SKA-I will also be useful to study circularly polarized emission from Sun and comparing various models of field configurations with observations.

preprint2015arXiv

On the resilience of helical magnetic fields to turbulent diffusion and the astrophysical implications

The extent to which large scale magnetic fields are susceptible to turbulent diffusion is important for interpreting the need for in situ large scale dynamos in astrophysics and for observationally inferring field strengths compared to kinetic energy. By solving coupled equations for magnetic energy and magnetic helicity in a system initiated with isotropic turbulence and an arbitrarily helical large scale field, we quantify the decay rate of the latter for a bounded or periodic system. The energy associated with the non-helical magnetic field rapidly decays by turbulent diffusion, but the decay rate of the helical component depends on whether the ratio of its magnetic energy to the turbulent kinetic energy exceeds a critical value given by M_{1,c} =(k_1/k_2)^2, where k_1 and k_2 are the wave numbers of the large and forcing scales. Turbulently diffusing helical fields to small scales while conserving magnetic helicity requires a rapid increase in total magnetic energy. As such, only when the helical fields are sub-critical can they so diffuse. When super-critical, the large scale helical field decays slowly, at a rate determined by microphysical dissipation even when macroscopic turbulence is present. Amplification of small scale magnetic helicity abates the turbulent diffusion. Two implications are that: (1) Standard arguments supporting the need for in situ large scale dynamos based on the otherwise rapid turbulent diffusion of large scale fields require re-thinking since only the non-helical field is so diffused in a closed system. Boundary terms could however provide potential pathways for rapid change of the large scale helical field. (2) Since M_{1,c} <<1 for k_1 << k_2, the presence of long-lived ordered large scale helical fields, as in extragalactic jets, does not guarantee that the magnetic field dominates the kinetic energy.

preprint2014arXiv

A physical model for the redshift evolution of high-z Lyman-Break Galaxies

We present a galaxy formation model to understand the evolution of stellar mass (M*) - UV luminosity relations, stellar mass functions and specific star formation rate (sSFR) of Lyman Break Galaxies (LBGs) along with their UV luminosity functions in the redshift range 3 < z < 8. Our models assume a physically motivated form for star formation in galaxies and model parameters are calibrated by fitting the observed UV luminosity functions (LFs) of LBGs. We find the fraction of baryons that gets converted into stars remains nearly constant for z < 4 but shows an increase for z < 4. However, the rate of converting baryons into stars does not evolve significantly in the redshift range 3 < z < 8. Our model further successfully explains the M* - UV luminosity (M_AB) correlations of LBGs. While our model predictions of stellar mass functions compare well with the inferred data from observations at the low mass end, we need to invoke the Eddington bias to fit the high mass end. At any given redshift, we find the sSFR to be constant over the stellar mass range 5 \times 10^8 -5 \times 10^9 M_\odot and the redshift evolution of sSFR is well approximated by a form (1+z)^2.4 for 3 < z < 8 which is consistent with observations. Thus we find that dark matter halo build up in the LCDM model is sufficient to explain the evolution of UV LFs of LBGs along with their M* - M_AB relations, the stellar mass functions and the sSFR for 3 < z < 8.

preprint2014arXiv

Fluctuation dynamo at finite correlation times and the Kazantsev spectrum

Fluctuation dynamos are generic to astrophysical systems. The only analytical model of the fluctuation dynamo is Kazantsev model which assumes a delta-correlated in time velocity field. We derive a generalized model of fluctuation dynamo with finite correlation time, $τ$, using renovating flows. For $τ\to 0$, we recover the standard Kazantsev equation for the evolution of longitudinal magnetic correlation, $M_L$. To the next order in $τ$, the generalized equation involves third and fourth spatial derivatives of $M_L$. It can be recast using the Landau-Lifschitz approach, to one with at most second derivatives of $M_L$. Remarkably, we then find that the magnetic power spectrum, remains the Kazantsev spectrum of $M(k) \propto k^{3/2}$, in the large $k$ limit, independent of $τ$.

preprint2014arXiv

Fluctuation dynamo at finite correlation times using renewing flows

Fluctuation dynamos are generic to turbulent astrophysical systems. The only analytical model of the fluctuation dynamo, due to Kazantsev, assumes the velocity to be delta-correlated in time. This assumption breaks down for any realistic turbulent flow. We generalize the analytic model of fluctuation dynamo to include the effects of a finite correlation time, $τ$, using renewing flows. The generalized evolution equation for the longitudinal correlation function $M_L$ leads to the standard Kazantsev equation in the $τ\to 0$ limit, and extends it to the next order in $τ$. We find that this evolution equation involves also third and fourth spatial derivatives of $M_L$, indicating that the evolution for finite $τ$ will be non-local in general. In the perturbative case of small-$τ$ (or small Strouhl number), it can be recast using the Landau-Lifschitz approach, to one with at most second derivatives of $M_L$. Using both a scaling solution and the WKBJ approximation, we show that the dynamo growth rate is reduced when the correlation time is finite. Interestingly, to leading order in $τ$, we show that the magnetic power spectrum, preserves the Kazantsev form, $M(k) \propto k^{3/2}$, in the large $k$ limit, independent of $τ$.

preprint2014arXiv

Helical cosmological magnetic fields from extra-dimensions

We study the inflationary generation of helical cosmological magnetic fields in a higher-dimensional generalization of the electromagnetic theory. For this purpose, we also include a parity breaking piece to the electromagnetic action. The evolution of extra-dimensional scale factor allows the breaking of conformal invariance of the effective electromagnetic action in $1+3$ dimensions required for such generation. Analytical solutions for the vector potential can be obtained in terms of Coulomb wave-functions for some special cases. We also present numerical solutions for the vector potential evolution in more general cases. In the presence of a higher-dimensional cosmological constant there exist solutions for the scale factors in which both normal and extra dimensional space either inflate or deflate simultaneously with the same rate. In such a scenario, with the number of extra dimensions $D=4$, a scale invariant spectrum of helical magnetic field is obtained. The net helicity arises, as one helical mode comes to dominate over the other at the superhorizon scales. A magnetic field strength of the order of $10^{-9}$ $G$ can be obtained for the inflationary scale $H\simeq 10^{-3}$ $M_{pl}$. Weaker fields will be generated for lower scales of inflation. Magnetic fields generated in this model respects the bounds on magnetic fields by Planck and $γ$-ray observations (i.e. $10^{-16}$ $G$ $<$ $B_{obs}<3.4\times 10^{-9}$ $G$).

preprint2014arXiv

Magnetic Spiral Arms and Galactic Outflows

Galactic magnetic arms have been observed between the gaseous arms of some spiral galaxies; their origin remains unclear. We suggest that magnetic spiral arms can be naturally generated in the interarm regions because the galactic fountain flow or wind is likely to be weaker there than in the arms. Galactic outflows lead to two countervailing effects: removal of small-scale magnetic helicity, which helps to avert catastrophic quenching of the dynamo, and advection of the large-scale magnetic field, which suppresses dynamo action. For realistic galactic parameters, the net consequence of outflows being stronger in the gaseous arms is higher saturation large-scale field strengths in the interarm regions as compared to in the arms. By incorporating rather realistic models of spiral structure and evolution into our dynamo models, an interlaced pattern of magnetic and gaseous arms can be produced.

preprint2014arXiv

Non-linear galactic dynamos: A toolbox

We compare various models and approximations for non-linear mean-field dynamos in disc galaxies to assess their applicability and accuracy, and thus to suggest a set of simple solutions suitable to model the large-scale galactic magnetic fields in various contexts. The dynamo saturation mechanisms considered are the magnetic helicity balance involving helicity fluxes (the dynamical $α$-quenching) and an algebraic $α$-quenching. The non-linear solutions are then compared with the marginal kinematic and asymptotic solutions. We also discuss the accuracy of the no-$z$ approximation. Although these tools are very different in the degree of approximation and hence complexity, they all lead to remarkably similar solutions for the mean magnetic field. In particular, we show that the algebraic $α$-quenching non-linearity can be obtained from a more physical dynamical $α$-quenching model in the limit of nearly azimuthal magnetic field. This suggests, for instance, that earlier results on galactic disc dynamos based on the simple algebraic non-linearity are likely to be reliable, and that estimates based on simple, even linear models are often a good starting point. We suggest improved no-$z$ and algebraic $α$-quenching models, and also incorporate galactic outflows into a simple analytical dynamo model to show that the outflow can produce leading magnetic spirals near the disc surface. The simple dynamo models developed are applied to estimate the magnetic pitch angle and the arm-interarm contrast in the saturated magnetic field strength for realistic parameter values.

preprint2014arXiv

Traces of large-scale dynamo action in the kinematic stage

Using direct numerical simulations (DNS) we verify that in the kinematic regime, a turbulent helical dynamo grows in such a way that the magnetic energy spectrum remains to high precision shape-invariant, i.e., at each wavenumber $k$ the spectrum grows with the same growth rate. Signatures of large-scale dynamo action can be identified through the excess of magnetic energy at small $k$, of one of the two oppositely polarized constituents. Also a suitably defined planar average of the magnetic field can be chosen such that its rms value isolates the strength of the mean field. However, these different means of analysis suggest that the strength of the large-scale field diminishes with increasing magnetic Reynolds number ${\rm Re}_{\rm M}$ like ${\rm Re}_{\rm M}^{-1/2}$ for intermediate values and like ${\rm Re}_{\rm M}^{-3/4}$ for larger ones. Both an analysis from the Kazantsev model including helicity and the DNS show that this arises due to the magnetic energy spectrum still peaking at resistive scales, even when helicity is present. As expected, the amplitude of the large-scale field increases with increasing fractional helicity, enabling us to determine the onset of large scale dynamo action and distinguishing it from that of the small-scale dynamo. Our DNS show that, contrary to earlier results for smaller scale separation (only 1.5 instead of now 4), the small-scale dynamo can still be excited at magnetic Prandtl numbers of 0.1 and only moderate values of the magnetic Reynolds numbers ($\sim 160$).

preprint2013arXiv

A study of large scale dynamo growth rates from numerical simulations and implications for mean field theories

Understanding large scale magnetic field growth in turbulent plasmas in the magnetohydrodynamic limit is a goal of magnetic dynamo theory. In particular, assessing how well large scale helical field growth and saturation in simulations matches that predicted by existing theories is important for progress. Using numerical simulations of isotropically forced turbulence without large scale shear with the implications, we focus on several aspects of this comparison that have not been previously tested: (1) Leading mean field dynamo theories which break the field into large and small scales predict that large scale helical field growth rates are determined by the difference between kinetic helicity and current helicity with no dependence on the non-helical energy in small scale magnetic fields. Our simulations show that the growth rate of the large scale field from fully helical forcing is indeed unaffected by the presence or absence of small scale magnetic fields amplified in a precursor non-helical dynamo. However, because the precursor non helical dynamo in our simulations produced fields that were strongly sub-equipartition with respect to the kinetic energy, we cannot yet rule out the potential influence of stronger non- helical small scale fields. (2) We have identified two features in our simulations which cannot be explained by the most minimalist versions of two-scale mean field theory: (i) fully helical small scale forcing produces significant non-helical large scale magnetic energy and (ii) the saturation of the large scale field growth is time-delayed with respect to what minimalist theory predicts. We comment on desirable generalizations to the theory in this context and future desired work.

preprint2013arXiv

Clustering at high redshift: The connection between Lyman Alpha emitters and Lyman break galaxies

We present a physically motivated semi-analytic model to understand the clustering of high redshift Lyman Alpha Emitters (LAEs). We show that the model parameters constrained by the observed luminosity functions, can be used to predict large scale bias and angular correlation function of LAEs. These predictions are shown to reproduce the observations remarkably well. We find that average masses of dark matter halos hosting LAEs brighter than threshold narrow band magnitude ~ 25 are ~ 10^11 M_\odot. These are smaller than that of typical Lyman Break Galaxies (LBGs) brighter than similar threshold continuum magnitude by a factor ~ 10. This results in a smaller clustering strength of LAEs compared to LBGs. However, using the observed relationship between UV continuum and Lyman-alpha luminosity of LAEs, we show that both LAEs and LBGs belong to the same parent galaxy population with narrow band techniques having greater efficiency in picking up galaxies with low UV luminosity. We also show that the lack of evidence for the presence of the one halo term in the observed LAE angular correlation functions can be attributed to sub-Poisson distribution of LAEs in dark matter halos as a result of their low halo occupations.

preprint2013arXiv

Cosmological Magnetogenesis From Extra-dimensional Gauss Bonnet Gravity

Generation of primordial magnetic fields during inflation typically requires the breaking of conformal invariance of the Electromagnetic action. In this paper this has been achieved naturally in a higher dimensional cosmological model with a Gauss-Bonnet term in the action. The evolution of the scale factor of the extra dimension (whose dynamics is influenced by the Gauss-Bonnet term) acts as the cause for the breaking of conformal invariance. Different cases have been investigated, each of which is characterized by the number of higher dimensions, the value of the Gauss-Bonnet parameter, and the cosmological constant. Many of the scenarios considered are highly constrained by the requirements that the cosmic evolution is stable, that the normal dimensions expand and that there is no back reaction due to growing electric fields. However there do exist scenarios which satisfy the above requirements and are well suited for magnetogenesis. In particular, a scenario where the number of extra dimensions $D=4$ and the cosmological constant is non-zero, turns out to be best suited for generating primordial magnetic fields. It is shown that for these values of parameters, a scale invariant magnetic field of the order of $10^{-10}-10^{-9}$ $G$ can be produced. Even in these most favorable scenarios, the higher dimensional space expands during inflation at the same rate as the normal dimension. Hence if a mechanism could freeze the evolution of the higher dimension, this seems to be a viable mechanism to produce acceptable primordial magnetic fields.

preprint2013arXiv

Galactic Spiral Patterns and Dynamo Action II: Asymptotic Solutions

The exploration of mean-field galactic dynamos affected by a galactic spiral pattern, begun in Chamandy et al. (2013, hereafter Paper I) with numerical simulations, is continued here with an asymptotic solution. The mean-field dynamo model used generalizes the standard theory to include the delayed response of the mean electromotive force to variations of the mean magnetic field and turbulence parameters (the temporal non-locality, or τ effect). The effect of the spiral pattern on the dynamo considered is the enhancement of the α-effect in spiral-shaped regions (which may overlap the gaseous spiral arms or be located in the interarm regions). The axisymmetric and enslaved non-axisymmetric modes of the mean magnetic field are studied semi-analytically to clarify and strengthen the numerical results. Good qualitative agreement is obtained between the asymptotic solution and numerical solutions of Paper I for a global, rigidly rotating material spiral (density wave). At all galactocentric distances except for the co-rotation radius, we find magnetic arms displaced in azimuth from the α-arms, so that the ridges of magnetic field strength are more tightly wound than the α-arms. Moreover, the effect of a finite dynamo relaxation time τ (related to the turbulence correlation time) is to phase-shift the magnetic arms in the direction opposite to the galactic rotation even at the co-rotation radius. This mechanism can be used to explain the phase shifts between magnetic and material arms observed in some spiral galaxies.

preprint2013arXiv

Magnetic arms generated by multiple interfering galactic spiral patterns

Interfering two- and three-arm spiral patterns have previously been inferred to exist in many galaxies and also in numerical simulations, and invoked to explain important dynamical properties, such as lack of symmetry, kinks in spiral arms, and star formation in armlets. The non-axisymmetric galactic mean-field dynamo model of Chamandy et al. is generalized to allow for such multiple co-existing spiral patterns in the kinetic alpha_k effect, leading to the existence of magnetic spiral arms in the large-scale magnetic field with several new properties. The large-scale magnetic field produced by an evolving superposition of two- and three-arm (or two- and four-arm) patterns evolves with time along with the superposition. Magnetic arms can be stronger and more extended in radius and in azimuth when produced by two interfering patterns rather than by one pattern acting alone. Transient morphological features arise in the magnetic arms, including bifurcations, disconnected armlets, and temporal and spatial variation in arm strength and pitch angles. Pitch angles of the large-scale magnetic field and magnetic arm structures (ridges) are smaller than those typically inferred from observations of spiral galaxies for model parameters of Chamandy et al., but can become comparable to typically inferred values for certain (still realistic) parameters. The magnetic field is sometimes strongest in between the alpha_k-arms, unlike in standard models with a single pattern, where it is strongest within the alpha_k-arms. Moreover, for models with a two- and three-arm pattern, some amount of m=1 azimuthal symmetry is found to be present in the magnetic field, which is generally not the case for forcing by single two- or three-arm patterns. (abridged)

preprint2013arXiv

Primordial Magnetic Field Limits from CMB Trispectrum - Scalar Modes and Planck Constraints

Cosmic magnetic fields are observed to be coherent on large scales and could have a primordial origin. Non-Gaussian signals in the cosmic microwave background (CMB) are generated by primordial magnetic fields as the magnetic stresses and temperature anisotropy they induce depend quadratically on the magnetic field. We compute the CMB scalar trispectrum on large angular scales, for nearly scale-invariant magnetic fields, sourced via the Sachs-Wolfe effect. The trispectra induced by magnetic energy density and by magnetic scalar anisotropic stress are found to have typical magnitudes of approximately $10^{-29}$ and $10^{-19}$, respectively. The scalar anisotropic stress trispectrum is also calculated in the flat-sky approximation and yields a similar result. Observational limits on CMB non-Gaussianity from the Planck mission data allow us to set upper limits of $B_0 \lesssim 0.6 $ nG on the present value of the primordial cosmic magnetic field. Considering the inflationary magnetic curvature mode in the trispectrum can further tighten the magnetic field upper limit to $B_0 \lesssim 0.05 $ nG. These sub-nanoGauss constraints from the magnetic trispectrum are the most stringent limits so far on the strength of primordial magnetic fields, on megaparsec scales, significantly better than the limits obtained from the CMB bispectrum and the CMB power spectrum.

preprint2013arXiv

Resilience of helical fields to turbulent diffusion II: direct numerical simulations

The recent study of Blackman and Subramanian (Paper I) indicates that large scale helical magnetic fields are resilient to turbulent diffusion in the sense that helical fields stronger than a critical value, decay on slow (~resistively mediated), rather than fast ($\sim$ turbulent) time scales. This gives more credence to potential fossil field origin models of the magnetic fields in stars, galaxies and compact objects. We analyze a suite of direct numerical simulations (DNS) of decaying large scale helical magnetic fields in the presence of non-helical turbulence to further study the physics of helical field decay. We study two separate cases: (1) the initial field is large enough to decay resistively, is tracked until it transitions to decay fast, and the critical large scale helical field at that transition is sought; (2) the case of Paper I, wherein there is a critical initial helical field strength below which the field undergoes fast decay right from the beginning. In case (1), both DNS and solutions of the two scale model (from Paper 1), reveal that the transition energy, $E_{c1}$, is independent of the scale of the turbulent forcing, within a small range of $\Rm$. We also find that the kinetic alpha, $α_K$, is subdominant to magnetic alpha, $α_M$, in the DNS, justifying an assumption in the two scale model. For case (2), we show exact solutions of two scale model in the limit of $η\rightarrow 0$ in fully helical case, leading to the transition energy, $E_{c2} = (k_1/\kf)^2 M_{eq}$, where $k_1$ and $\kf$ are the large scale and small turbulent forcing scale respectively and $M_{eq}$ is the equipartition magnetic energy. The DNS in this case agree qualitatively with the two scale model but the $R_M$ currently achievable, is too small to satisfy a condition $3/R_M << (k_1/k_f)^2$, necessary to robustly reveal the transition, $E_{c2}$ (Abridged).

preprint2012arXiv

Cosmic Microwave Background Trispectrum and Primordial Magnetic Field Limits

Primordial magnetic fields will generate non-Gaussian signals in the cosmic microwave background (CMB) as magnetic stresses and the temperature anisotropy they induce depend quadratically on the magnetic field. We compute a new measure of magnetic non-Gaussianity, the CMB trispectrum, on large angular scales, sourced via the Sachs-Wolfe effect. The trispectra induced by magnetic energy density and by magnetic scalar anisotropic stress are found to have typical magnitudes of approximately a few times 10^{-29} and 10^{-19}, respectively. Observational limits on CMB non-Gaussianity from WMAP data allow us to conservatively set upper limits of a nG, and plausibly sub-nG, on the present value of the primordial cosmic magnetic field. This represents the tightest limit so far on the strength of primordial magnetic fields, on Mpc scales, and is better than limits from the CMB bispectrum and all modes in the CMB power spectrum. Thus, the CMB trispectrum is a new and more sensitive probe of primordial magnetic fields on large scales.

preprint2012arXiv

Current status of turbulent dynamo theory: From large-scale to small-scale dynamos

Several recent advances in turbulent dynamo theory are reviewed. High resolution simulations of small-scale and large-scale dynamo action in periodic domains are compared with each other and contrasted with similar results at low magnetic Prandtl numbers. It is argued that all the different cases show similarities at intermediate length scales. On the other hand, in the presence of helicity of the turbulence, power develops on large scales, which is not present in non-helical small-scale turbulent dynamos. At small length scales, differences occur in connection with the dissipation cutoff scales associated with the respective value of the magnetic Prandtl number. These differences are found to be independent of whether or not there is large-scale dynamo action. However, large-scale dynamos in homogeneous systems are shown to suffer from resistive slow-down even at intermediate length scales. The results from simulations are connected to mean field theory and its applications. Recent work on helicity fluxes to alleviate large-scale dynamo quenching, shear dynamos, nonlocal effects and magnetic structures from strong density stratification are highlighted. Several insights which arise from analytic considerations of small-scale dynamos are discussed.

preprint2012arXiv

Fluctuation dynamos and their Faraday rotation signatures

Turbulence is ubiquitous in many astrophysical systems like galaxies, galaxy clusters and possibly even the IGM filaments. We study fluctuation dynamo action in turbulent systems focusing on one observational signature; the Faraday rotation measure (RM) from background radio sources seen through the magnetic field generated by such a dynamo. We simulate the fluctuation dynamo (FD) in periodic boxes up to resolutions of 512^3, with varying fluid and magnetic Reynolds numbers, and measure the resulting random RMs. We show that, even though the magnetic field generated is intermittent, it still allows for contributions to the RM to be significant. When the dynamo saturates, it is of order 40%-50% of the value expected in a model where fields of strength B_rms uniformly fill cells of the largest turbulent eddy but are randomly oriented from one cell to another. This level of RM dispersion obtains across different values of magnetic Reynolds number and Prandtl number explored. We also use the random RMs to probe the structure of the generated fields to distinguish the contribution from intense and diffuse field regions. We find that the strong field regions (say with B > 2B_rms) contribute only of order 15%-20% to the RM. Thus rare structures do not dominate the RM; rather the general 'sea' of volume filling fluctuating fields are the dominant contributors. We also show that the magnetic integral scale, L_{int}, which is directly related to the RM dispersion, increases in all the runs, as Lorentz forces become important to saturate the dynamo. It appears that due to the ordering effect of the Lorentz forces, L_{int} of the saturated field tends to a modest fraction, 1/2-1/3 of the integral scale of the velocity field, for all our runs. These results are then applied to discuss the RM signatures of FD generated fields in young galaxies, galaxy clusters and intergalactic filaments.

preprint2012arXiv

Galactic spiral patterns and dynamo action I: A new twist on magnetic arms

We generalise the theory of mean-field galactic dynamos by allowing for temporal non-locality in the mean electromotive force (emf). This arises in random flows due to a finite response time of the mean emf to changes in the mean magnetic field, and leads to the telegraph equation for the mean field. The resulting dynamo model also includes the nonlinear dynamo effects arising from magnetic helicity balance. Within this framework, coherent large-scale magnetic spiral arms superimposed on the dominant axially symmetric magnetic structure are considered. A non-axisymmetric forcing of the mean-field dynamo by a spiral pattern (either stationary or transient) is invoked, with the aim of explaining the phenomenon of magnetic arms. For a stationary dynamo forcing by a rigidly rotating material spiral, we find corotating non-axisymmetric magnetic modes enslaved to the axisymmetric modes and strongly peaked around the corotation radius. For a forcing by transient material arms wound up by the galactic differential rotation, the magnetic spiral is able to adjust to the winding so that it resembles the material spiral at all times. There are profound effects associated with the temporal non-locality, i.e. finite `dynamo relaxation time'. For the case of a rigidly rotating spiral, a finite relaxation time causes each magnetic arm to mostly lag the corresponding material arm with respect to the rotation. For a transient material spiral that winds up, the finite dynamo relaxation time leads to a large, negative (in the sense of the rotation) phase shift between the magnetic and material arms, similar to that observed in NGC 6946 and other galaxies. We confirm that sufficiently strong random seed fields can lead to global reversals of the regular field along the radius whose long-term survival depends on specific features of a given galaxy.

preprint2012arXiv

Mean field dynamo action in renovating shearing flows

We study mean field dynamo action in renovating flows with finite and non zero correlation time ($τ$) in the presence of shear. Previous results obtained when shear was absent are generalized to the case with shear. The question of whether the mean magnetic field can grow in the presence of shear and non helical turbulence, as seen in numerical simulations, is examined. We show in a general manner that, if the motions are strictly non helical, then such mean field dynamo action is not possible. This result is not limited to low (fluid or magnetic) Reynolds numbers nor does it use any closure approximation; it only assumes that the flow renovates itself after each time interval $τ$. Specifying to a particular form of the renovating flow with helicity, we recover the standard dispersion relation of the $α^2 Ω$ dynamo, in the small $τ$ or large wavelength limit. Thus mean fields grow even in the presence of rapidly growing fluctuations, surprisingly, in a manner predicted by the standard quasilinear closure, even though such a closure is not strictly justified. Our work also suggests the possibility of obtaining mean field dynamo growth in the presence of helicity fluctuations, although having a coherent helicity will be more efficient.

preprint2012arXiv

Spatial Clustering of High Redshift Lyman Break Galaxies

We present a physically motivated semi-analytic model to understand the clustering of high redshift LBGs. We show that the model parameters constrained by the observed luminosity function, can be used to predict large scale (θ> 80 arcsec) bias and angular correlation function of galaxies. These predictions are shown to reproduce the observations remarkably well. We then adopt these model parameters to calculate the halo occupation distribution (HOD) using the conditional mass function. The halo model using this HOD is shown to provide a reasonably good fit to the observed clustering of LBGs at both large (θ>80") and small (θ< 10") angular scales for z=3-5 and several limiting magnitudes. However, our models underpredict the clustering amplitude at intermediate angular scales, where quasi-linear effects are important. The average mass of halos contributing to the observed clustering is found to be 6.2 x 10^{11} M_\odot and the characteristic mass of a parent halo hosting satellite galaxies is 1.2 \times 10^{12} M_\odot for a limiting absolute magnitude of -20.5 at z=4. For a given threshold luminosity these masses decrease with increasing z and at any given z these are found to increase with increasing value of threshold luminosity. We find that approximately 40 % of the halos above a minimum mass M_{min}, can host detectable central galaxies and about 5-10 % of these halos are likely to also host a detectable satellite. The satellites form typically a dynamical timescale prior to the formation of the parent halo. The small angular scale clustering is due to central-satellite pairs and is quite sensitive to changes in the duration of star formation in a halo. The present data favor star formation in a halo lasting typically for a few dynamical time-scales. Our models also reproduce different known trends between parameters related to star formation.

preprint2011arXiv

Cosmic ray transport in galaxy clusters: implications for radio halos and gamma-rays

Observations of giant radio halos provide unambiguous evidence for the existence of cosmic ray (CR) electrons and magnetic fields in galaxy clusters. The physical mechanism generating radio halos is still heavily debated. We critically discuss the proposed models for the radio halo emission and highlight the weaknesses underlying each explanation. We present an idea how the interplay of CR propagation and turbulent advection selects a bimodal spatial CR distribution that is characteristic for the dynamical state of a cluster. As a result, strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma ray emission of clusters. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons.

preprint2011arXiv

Scale-dependence of magnetic helicity in the solar wind

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. (1982, Phys. Rev. Lett. 48, 1256) by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k~2 AU^{-1} (or frequency nu~2 uHz) at distances below 2.8 AU and at k~30 AU^{-1} (or nu~25 uHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10^{45} Mx^2/cycle at large scales and to a 3 times lower value at smaller scales.

preprint2011arXiv

The First Magnetic Fields

We review current ideas on the origin of galactic and extragalactic magnetic fields. We begin by summarizing observations of magnetic fields at cosmological redshifts and on cosmological scales. These observations translate into constraints on the strength and scale magnetic fields must have during the early stages of galaxy formation in order to seed the galactic dynamo. We examine mechanisms for the generation of magnetic fields that operate prior during inflation and during subsequent phase transitions such as electroweak symmetry breaking and the quark-hadron phase transition. The implications of strong primordial magnetic fields for the reionization epoch as well as the first generation of stars is discussed in detail. The exotic, early-Universe mechanisms are contrasted with astrophysical processes that generate fields after recombination. For example, a Biermann-type battery can operate in a proto-galaxy during the early stages of structure formation. Moreover, magnetic fields in either an early generation of stars or active galactic nuclei can be dispersed into the intergalactic medium.

preprint2011arXiv

Weighing neutrinos using high redshift galaxy luminosity functions

Laboratory experiments measuring neutrino oscillations, indicate small mass differences between different mass eigenstates of neutrinos. The absolute mass scale is however not determined, with at present the strongest upper limits coming from astronomical observations rather than terrestrial experiments. The presence of massive neutrinos suppresses the growth of perturbations below a characteristic mass scale, thereby leading to a decreased abundance of collapsed dark matter halos. Here we show that this effect can significantly alter the predicted luminosity function (LF) of high redshift galaxies. In particular we demonstrate that a stringent constraint on the neutrino mass can be obtained using the well measured galaxy LF and our semi-analytic structure formation models. Combining the constraints from the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) data with the LF data at z = 4, we get a limit on the sum of the masses of 3 degenerate neutrinos Σm_ν< 0.52 eV at the 95 % CL. The additional constraints using the prior on Hubble constant strengthens this limit to Σm_ν< 0.29 eV at the 95 % CL. This neutrino mass limit is a factor of order 4 improvement compared to the constraint based on the WMAP7 data alone, and as stringent as known limits based on other astronomical observations. As different astronomical measurements may suffer from different set of biases, the method presented here provides a complementary probe of Σm_ν. We suggest that repeating this exercise on well measured luminosity functions over different redshift ranges can provide independent and tighter constraints on Σm_ν.

preprint2010arXiv

Cosmic ray transport in galaxy clusters: implications for radio halos, gamma-ray signatures, and cool core heating

We investigate the interplay of cosmic ray (CR) propagation and advection in galaxy clusters. Propagation in form of CR diffusion and streaming tends to drive the CR radial profiles towards being flat, with equal CR number density everywhere. Advection of CR by the turbulent gas motions tends to produce centrally enhanced profiles. We assume that the CR streaming velocity is of the order of the sound velocity. This is motivated by plasma physical arguments. The CR streaming is then usually larger than typical advection velocities and becomes comparable or lower than this only for periods with trans- and super-sonic cluster turbulence. As a consequence a bimodality of the CR spatial distribution results. Strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma-ray emission of clusters, since more centrally concentrated CR will find higher target densities for hadronic CR proton interactions, higher plasma wave energy densities for CR electron and proton re-acceleration, and stronger magnetic fields. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons. Energy dependence of the CR propagation should lead to spectral steepening of dying radio halos. Furthermore, we show that the interplay of CR diffusion with advection implies first order CR re-acceleration in the pressure-stratified atmospheres of galaxy clusters. Finally, we argue that CR streaming could be important in turbulent cool cores of galaxy clusters since it heats preferentially the central gas with highest cooling rate.

preprint2010arXiv

Magnetic fields in the early universe

We give a pedagogical introduction to two aspects of magnetic fields in the early universe. We first focus on how to formulate electrodynamics in curved space time, defining appropriate magnetic and electric fields and writing Maxwell equations in terms of these fields. We then specialize to the case of magnetohydrodynamics in the expanding universe. We emphasize the usefulness of tetrads in this context. We then review the generation of magnetic fields during the inflationary era, deriving in detail the predicted magnetic and electric spectra for some models. We discuss potential problems arising from back reaction effects and from the large variation of the coupling constants required for such field generation.

preprint2010arXiv

Nonperturbative quasilinear approach to the shear dynamo problem

We study large-scale dynamo action due to turbulence in the presence of a linear shear flow. Our treatment is quasilinear and equivalent to the standard `first order smoothing approximation'. However it is non perturbative in the shear strength. We first derive an integro-differential equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invariance of the linear shear flow. We show that, for non helical turbulence, the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves; this is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence, to all orders in the shear parameter, there is no shear-current type effect for non helical turbulence in a linear shear flow, in quasilinear theory in the limit of zero resistivity. We then develop a systematic approximation of the integro-differential equation for the case when the mean magnetic field varies slowly compared to the turbulence correlation time. For non-helical turbulence, the resulting partial differential equations can again be solved by making a shearing coordinate transformation in Fourier space. The resulting solutions are in the form of shearing waves, labeled by the wavenumber in the sheared coordinates. These shearing waves can grow at early and intermediate times but are expected to decay in the long time limit.

preprint2010arXiv

Primordial Magnetic Field Limits from Cosmic Microwave Background Bispectrum of Magnetic Passive Scalar Modes

Primordial magnetic fields lead to non-Gaussian signals in the cosmic microwave background (CMB) even at the lowest order, as magnetic stresses and the temperature anisotropy they induce depend quadratically on the magnetic field. In contrast, CMB non-Gaussianity due to inflationary scalar perturbations arises only as a higher order effect. Apart from a compensated scalar mode, stochastic primordial magnetic fields also produce scalar anisotropic stress that remains uncompensated till neutrino decoupling. This gives rise to an adiabatic-like scalar perturbation mode that evolves passively thereafter (called the passive mode). We compute the CMB reduced bispectrum ($b_{l_{_1}l_{_2}l_{_3}}$) induced by this passive mode, sourced via the Sachs-Wolfe effect, on large angular scales. For any configuration of bispectrum, taking a partial sum over mode-coupling terms, we find a typical value of $l_1(l_1+1)l_3(l_3+1) b_{l_{_1}l_{_2}l_{_3}} \sim 6-9 \times 10^{-16}$, for a magnetic field of $B_0 \sim 3$ nG, assuming a nearly scale-invariant magnetic spectrum . We also evaluate, in full, the bispectrum for the squeezed collinear configuration over all angular mode-coupling terms and find $l_1(l_1+1)l_3(l_3+1) b_{l_{_1}l_{_2}l_{_3}} \approx -1.4 \times 10^{-16}$. These values are more than $\sim 10^6$ times larger than the previously calculated magnetic compensated scalar mode CMB bispectrum. Observational limits on the bispectrum from WMAP7 data allow us to set upper limits of $B_0 \sim 2$ nG on the present value of the cosmic magnetic field of primordial origin. This is over 10 times more stringent than earlier limits on $B_0$ based on the compensated mode bispectrum.

preprint2009arXiv

A Reconnecting Flux Rope Dynamo

We develop a new model of the fluctuation dynamo in which the magnetic field is confined in thin flux ropes advected by a multi-scale flow modeling turbulence. Magnetic dissipation occurs only via reconnection of the flux ropes. We investigate the kinetic energy release into heat, mediated by the dynamo action, both in our model and by solving the induction equation with the same flow. We find that a flux rope dynamo is an order of magnitude more efficient at converting mechanical energy into heat. The probability density of the magnetic energy release in reconnections has a power-law form with the slope -3, consistent with the Solar corona heating by nanoflares.

preprint2009arXiv

Cosmic ray driven outflows from high redshift galaxies

We study winds in high redshift galaxies driven by a relativistic cosmic ray (proton) component in addition to the hot thermal gas component. Cosmic rays (CRs) are likely to be efficiently generated in supernova shocks inside galaxies. We obtain solutions of such CR driven free winds in a gravitational potential of the NFW form, relevant to galaxies. Cosmic rays naturally provide the extra energy and/or momentum input to the system, needed for a transonic wind solution in a gas with adiabatic index $γ=5/3$. We show that CRs can effectively drive winds even when the thermal energy of the gas is lost due to radiative cooling. These wind solutions predict an asymptotic wind speed closely related to the circular velocity of the galaxy. Furthermore, the mass outflow rate per unit star formation rate (eta_w) is predicted to be ~ 0.2-0.5 for massive galaxies, with masses $M \sim 10^{11}-10^{12} M_\odot$. We show eta_w to be inversely proportional to the square of the circular velocity. Magnetic fields at the $μ$G levels are also required in these galaxies to have a significant mass loss. A large eta_w for small mass galaxies implies that CR driven outflows could provide a strong negative feedback to the star formation in dwarf galaxies. Further, our results will also have important implications to the metal enrichment of the IGM. These conclusions are applicable to the class of free wind models where the source region is confined to be within the sonic point.

preprint2009arXiv

Fluctuation dynamo based on magnetic reconnections

We develop a new model of the fluctuation dynamo in which the magnetic field is confined to thin flux ropes advected by a multi-scale flow which models turbulence. Magnetic dissipation occurs only via reconnections of flux ropes. The model is particularly suitable for rarefied plasma, such as the Solar corona or galactic halos. We investigate the kinetic energy release into heat, mediated by dynamo action, both in our model and by solving the induction equation with the same flow. We find that the flux rope dynamo is more than an order of magnitude more efficient at converting mechanical energy into heat. The probability density of the magnetic energy released during reconnections has a power-law form with the slope -3, consistent with the Solar corona heating by nanoflares. We also present a nonlinear extension of the model. This shows that a plausible saturation mechanism of the fluctuation dynamo is the suppression of turbulent magnetic diffusivity, due to suppression of random stretching at the location of the flux ropes. We confirm that the probability distribution function of the magnetic line curvature has a power-law form suggested by (Schekochihin et al., 2002). We argue, however, using our results that this does not imply a persistent folded structure of magnetic field, at least in the nonlinear stage.

preprint2009arXiv

Models of high redshift luminosity functions and galactic outflows: The dependence on halo mass function

The form of the halo mass function is a basic ingredient in any semi-analytical galaxy formation model. We study the existing forms of the mass functions in the literature and compare their predictions for semi-analytical galaxy formation models. Two methods are used in the literature to compute the net formation rate of halos, one by simply taking the derivative of the halo mass function and the other using the prescription due to Sasaki (1994). For the Press-Schechter (PS) mass function, we compare various model predictions, using these two methods. However, as the Sasaki formalism cannot be easily generalized for other mass functions, we use the derivative while comparing model predictions of different mass functions. We show that the reionization history and UV luminosity function of Lyman break galaxies (LBGs) predicted by the PS mass function differs from those using any other existing mass function, like Sheth-Tormen (ST) mass function.In particular the reionization efficiency of molecular cooled halos has to be substantially reduced when one uses the ST and other mass functions obtained from the simulation instead of the PS mass function. Using $χ^2$-minimization, we find that the observed UV luminosity functions of LBGs at $3.0\le z\le 7.4$ are better reproduced by models using the ST mass function compared to models that use the PS mass function. On the other hand, the volume filling factor of the metals expelled from the galaxies through supernovae driven outflows differs very little between models with different mass functions. It depends on the way we treat merging outflows. We also show that the porosity weighted average quantities related to the outflow are not very sensitive to the differences in the halo mass function.

preprint2009arXiv

Primordial magnetic fields and the HI signal from the epoch of reionization

The implication of primordial magnetic-field-induced structure formation for the HI signal from the epoch of reionization is studied. Using semi-analytic models, we compute both the density and ionization inhomogeneities in this scenario. We show that: (a) The global HI signal can only be seen in emission, unlike in the standard $Λ$CDM models, (b) the density perturbations induced by primordial fields, leave distinctive signatures of the magnetic field Jeans' length on the HI two-point correlation function, (c) the length scale of ionization inhomogeneities is $\la 1 \rm Mpc$. We find that the peak expected signal (two-point correlation function) is $\simeq 10^{-4} \rm K^2$ in the range of scales $0.5\hbox{-}3 \rm Mpc$ for magnetic field strength in the range $5 \times 10^{-10} \hbox{-}3 \times 10^{-9} \rm G$. We also discuss the detectability of the HI signal. The angular resolution of the on-going and planned radio interferometers allows one to probe only the largest magnetic field strengths that we consider. They have the sensitivity to detect the magnetic field-induced features. We show that thefuture SKA has both the angular resolution and the sensitivity to detect the magnetic field-induced signal in the entire range of magnetic field values we consider, in an integration time of one week.

preprint2009arXiv

Shear dynamo problem: Quasilinear kinematic theory

Large-scale dynamo action due to turbulence in the presence of a linear shear flow is studied. Our treatment is quasilinear and kinematic but is non perturbative in the shear strength. We derive the integro-differential equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invariance of the linear shear flow. For non helical turbulence the time evolution of the cross-shear components of the mean field do not depend on any other components excepting themselves. This is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence the shear-current assisted dynamo is essentially absent, although large-scale non helical dynamo action is not ruled out.

preprint2009arXiv

Understanding the redshift evolution of the luminosity functions of Lyman-alpha emitters

We present a semi-analytical model of star formation which explains simultaneously the observed UV luminosity function of high redshift Lyman break galaxies (LBGs) and luminosity functions of Lyman-alpha emitters. We consider both models that use the Press-Schechter (PS) and Sheth-Tormen (ST) halo mass functions to calculate the abundances of dark matter halos. The Lyman-alpha luminosity functions at z < 4 are well reproduced with only <10% of the LBGs emitting Lyman-alpha lines with rest equivalent width greater than the limiting equivalent width of the narrow band surveys. However, the observed luminosity function at z > 5 can be reproduced only when we assume that nearly all LBGs are Lyman-alpha emitters. Thus it appears that 4 < z < 5 marks the epoch when a clear change occurs in the physical properties of the high redshift galaxies. As Lyman-alpha escape depends on dust and gas kinematics of the inter stellar medium (ISM), this could mean that on an average the ISM at z > 5 could be less dusty, more clumpy and having more complex velocity field. All of these will enable easier escape of the Lyman-alpha photons. At z > 5 the observed Lyman-alpha luminosity function are well reproduced with the evolution in the halo mass function along with very minor evolution in the physical properties of high redshift galaxies. In particular, upto z = 6.5, we do not see the effect of evolving inter galactic medium (IGM) opacity on the Lyman-alpha escape from these galaxies.

preprint2006arXiv

Magnetic helicity density and its flux in weakly inhomogeneous turbulence

A gauge invariant and hence physically meaningful definition of magnetic helicity density for random fields is proposed, using the Gauss linking formula, as the density of correlated field line linkages. This definition is applied to the random small scale field in weakly inhomogeneous turbulence, whose correlation length is small compared with the scale on which the turbulence varies. For inhomogeneous systems, with or without boundaries, our technique then allows one to study the local magnetic helicity density evolution in a gauge independent fashion, which was not possible earlier. This evolution equation is governed by local sources (owing to the mean field) and by the divergence of a magnetic helicity flux density. The role of magnetic helicity fluxes in alleviating catastrophic quenching of mean field dynamos is discussed.

preprint2005arXiv

CMB Anisotropy Due to Tangled magnetic fields in re-ionized models

Primordial tangled cosmological Magnetic Fields source rotational velocity perturbations of the baryon fluid, even in the post-recombination universe. These vortical modes inturn leave a characteristic imprint on the temperature anisotropy of the Cosmic Microwave Background (CMB), if the CMB photons can be re-scatterred after recombination. Observations from WMAP indicate that the Universe underwent a relatively early re-ionization ($z_{ri} \sim 15$), which does indeed lead to a significant optical depth for re-scattering of CMB photons after the re-ionization epoch. We compute the resulting additional temperature anisotropies, induced by primordial magnetic fields in the post-recombination universe. We show that in models with early re-ionization, a nearly scale-invariant spectrum of tangled magnetic fields which redshift to a present value of $B_0 \sim 3 \times 10^{-9}$ Gauss, produce vector mode perturbations which in turn induce additional temperature anisotropy of about 0.3 to 0.4 $ μ$K over very small angular scales, with $ l$ upto $\sim 10000 $ or so.

preprint1998arXiv

Microwave Background Signals from Tangled Magnetic Fields

An inhomogeneous cosmological magnetic field will create Alfven-wave modes that induce a small rotational velocity perturbation on the last scattering surface of the microwave background radiation. The Alfven-wave mode survives Silk damping on much smaller scales than the compressional modes. This, in combination with its rotational nature, ensures that there will be no sharp cut-off in anisotropy on arc-minute scales. We estimate that a magnetic field which redshifts to a present value of $3\times 10^{-9}$ Gauss produces temperature anisotropies at the 10 micro Kelvin level at and below 10 arc-min scales. A tangled magnetic field, which is large enough to influence the formation of large scale structure is therefore potentially detectable by future observations.

preprint1997arXiv

Magnetohydrodynamics in the Early Universe and the Damping of Non-linear Alfven Waves

The evolution and viscous damping of cosmic magnetic fields in the early universe, is analysed. Using the fact that the fluid, electromagnetic, and shear viscous energy-momentum tensors are all conformally invariant, the evolution is transformed from the expanding universe setting into that in flat spacetime. Particular attention is paid to the evolution of nonlinear Alfven modes. For a small enough magnetic field, which satisfies our observational constraints, these wave modes either oscillate negligibly or, when they do oscillate, become overdamped. Hence they do not suffer Silk damping on galactic and subgalactic scales. The smallest scale which survives damping depends on the field strength and is of order a dimensionless Alfven velocity times the usual baryon-photon Silk damping scale. After recombination, nonlinear effects can convert the Alfven mode into compressional, gravitationally unstable waves and seed cosmic structures if the cosmic magnetic field is sufficiently strong.