Search

The origin and evolution of magnetic fields of neutron stars from birth has long been a source of debate. Here, motivated by recent simulations of the Hall cascade with magnetic helicity, we invoke a model where the large-scale magnetic field of neutron stars grows as a product of small-scale turbulence through an inverse cascade. We apply this model to a simulated population of neutron stars at birth and show how this model can account for the evolution of such objects across the $P\dot{P}$ diagram, explaining both pulsar and magnetar observations. Under the assumption that small-scale turbulence is responsible for large-scale magnetic fields, we place a lower limit on the spherical harmonic degree of the energy-carrying magnetic eddies of $\approx 40$. Our results favor the presence of a highly resistive pasta layer at the base of the neutron star crust. We further discuss the implications of this paradigm on direct observables, such as the nominal age and braking index of pulsars., Comment: Submitted. 6 pages, 3 figures

Mean-field dynamo theory has important applications in solar physics and galactic magnetism. We discuss some of the many turbulence effects relevant to the generation of large-scale magnetic fields in the solar convection zone. The mean-field description is then used to illustrate the physics of the $\alpha$ effect, turbulent pumping, turbulent magnetic diffusivity, and other effects on a modern solar dynamo model. We also discuss how turbulence transport coefficients are derived from local simulations of convection and then used in mean-field models., Comment: 64 pages, 20 figures, 2 tables, submitted to Space Science Reviews, special issue "Solar and stellar dynamos: a new era"

Turbulence is typically not in equilibrium, i.e. mean quantities such as the mean energy and helicity are typically time-dependent. The effect of non-stationarity on the turbulent hydromagnetic dynamo process is studied here with the use of the two-scale direct-interaction approximation (TSDIA), which allows to explicitly relate the mean turbulent Reynolds and Maxwell stresses and the mean electromotive force (EMF) to the spectral characteristics of turbulence, such as e.g. the mean energy, as well as kinetic and cross-helicity. It is demonstrated, that the non-equilibrium effects can enhance the dynamo process when the magnetohydrodynamic (MHD) turbulence is both helical and cross-helical. This effect is based on the turbulent infinitesimal-impulse cross-response functions, which do not affect turbulent flows in equilibrium. The evolution and sources of the cross-helicity in MHD turbulence is also discussed., 23 pages, 3 figures

We study the evolution of magnetic fields coupled with chiral fermion asymmetry in the framework of chiral magnetohydrodynamics with zero initial total chirality. The initial magnetic field has a turbulent spectrum peaking at a certain characteristic scale and is fully helical with positive helicity. The initial chiral chemical potential is spatially uniform and negative. We consider two opposite cases where the ratio of the length scale of the chiral plasma instability (CPI) to the characteristic scale of the turbulence is smaller and larger than unity. These initial conditions might be realized in cosmological models such as certain types of axion inflation. The magnetic field and chiral chemical potential evolve with inverse cascading in such a way that the magnetic helicity and chirality cancel each other at all times. The CPI time scale is found to determine mainly the time when the magnetic helicity spectrum attains negative values at high wave numbers. The turnover time of the energy-carrying eddies, on the other hand, determines the time when the peak of the spectrum starts to shift to smaller wave numbers via an inverse cascade. The onset of helicity decay is determined by the time when the chiral magnetic effect becomes efficient at the peak of the initial magnetic energy spectrum. When spin flipping is important, the chiral chemical potential vanishes and the magnetic helicity becomes constant, which leads to a faster increase of the correlation length, as expected from magnetic helicity conservation. This also happens when the initial total chirality is imbalanced. Our findings have important implications for baryogenesis after axion inflation., Comment: 21 pages, 20 figures, 4 tables

The conversion of axionlike particles (ALPs) and photons in magnetised astrophysical environments provides a promising route to search for ALPs. The strongest limits to date on light ALPs use galaxy clusters as ALP--photon converters. However, such studies traditionally rely on simple models of the cluster magnetic fields, with the state-of-the-art being Gaussian random fields (GRFs). We present the first systematic study of ALP-photon conversion in more realistic, turbulent fields from dedicated magnetohydrodynamic (MHD) simulations, which we compare with GRF models. For GRFs, we analytically derive the distribution of conversion ratios at fixed energy and find that it follows an exponential law. We find that the MHD models agree with the exponential law for typical, small amplitude mixings but exhibit distinctly heavy tails for rare and large mixings. We explain how non-Gaussian, local spikes in the MHD magnetic field are mainly responsible for the heavy tail. Our results indicate that limits placed on ALPs using GRFs are conservative but that MHD models are necessary to reach the full potential of these searches., Comment: 8 pages, 3 figures

If the alpha effect plays a role in the generation of the Sun's magnetic field, the field should show evidence of magnetic helicity of opposite signs at large and small length scales. Measuring this faces two challenges: (i) in weak-field regions, horizontal field measurements are unreliable because of the pi ambiguity, and (ii) one needs a truly global approach to computing helicity spectra in the case where one expects a sign reversal across the equator at all wavenumbers. Here we develop such a method using spin-2 spherical harmonics to decompose the linear polarization in terms of the parity-even and parity-odd E and B polarizations, respectively. Using simple one- and two-dimensional models, we show that the product of the spectral decompositions of E and B, taken at spherical harmonic degrees that are shifted by one, can act as a proxy of the global magnetic helicity with a sign that represents that in the northern hemisphere. We then apply this method to the analysis of solar synoptic vector magnetograms, from which we extract a pseudo-polarization corresponding to a "pi-ambiguated" magnetic field, i.e., a magnetic field vector that has no arrow. We find a negative sign of the global EB helicity proxy at spherical harmonic degrees of around 6. This could indicate a positive magnetic helicity at large length scales, but the spectrum fails to capture clear evidence of the well-known negative magnetic helicity at smaller scales. This method might also be applicable to stellar and Galactic polarization data., 11 pages, 7 figures, 5 tables, ApJ, in press

We investigate the effect of turbulence on the combined condensational and collisional growth of cloud droplets by means of high resolution direct numerical simulations of turbulence and a superparticle approximation for droplet dynamics and collisions. The droplets are subject to turbulence as well as gravity, and their collision and coalescence efficiencies are taken to be unity. We solve the thermodynamic equations governing temperature, water-vapor mixing ratio, and the resulting supersaturation fields together with the Navier-Stokes equation. We find that the droplet-size distribution broadens with increasing Reynolds number and/or mean energy dissipation rate. Turbulence affects the condensational growth directly through supersaturation fluctuations, and it influences collisional growth indirectly through condensation. Our simulations show for the first time that, in the absence of the mean updraft cooling, supersaturation fluctuation-induced broadening of droplet-size distributions enhances the collisional growth. This is contrary to classical (non-turbulent) condensational growth, which leads to a growing mean droplet size, but a narrower droplet-size distribution. Our findings, instead, show that condensational growth facilitates collisional growth by broadening the size distribution in the tails at an early stage of rain formation. With increasing Reynolds numbers, evaporation becomes stronger. This counteracts the broadening effect due to condensation at late stages of rain formation. Our conclusions are consistent with results of laboratory experiments and field observations, and show that supersaturation fluctuations are important for precipitation., Comment: revised version, 17 pages, 1 table, 13 figures

In the presence of strong density stratification, hydromagnetic turbulence attains qualitatively new properties: the formation of magnetic flux concentrations. We review here the theoretical foundations of this mechanism in terms of what is now called the negative effective magnetic pressure instability. We also present direct numerical simulations of forced turbulence in strongly stratified layers and discuss the qualitative and quantitative similarities with corresponding mean-field simulations. Finally, the relevance to sunspot formation is discussed., Comment: 23 pages, 13 figures, special issue "Focus on Turbulence in Astrophysical and Laboratory Plasmas" of New Journal of Physics, published

We review our current understanding of galactic dynamo theory, paying particular attention to numerical simulations both of the mean-field equations and the original three-dimensional equations relevant to describing the magnetic field evolution for a turbulent flow. We emphasize the theoretical difficulties in explaining non-axisymmetric magnetic fields in galaxies and discuss the observational basis for such results in terms of rotation measure analysis. Next, we discuss nonlinear theory, the role of magnetic helicity conservation and magnetic helicity fluxes. This leads to the possibility that galactic magnetic fields may be bi-helical, with opposite signs of helicity and large and small length scales. We discuss their observational signatures and close by discussing the possibilities of explaining the origin of primordial magnetic fields., Comment: 28 pages, 15 figure, to appear in Lecture Notes in Physics "Magnetic fields in diffuse media", Eds. E. de Gouveia Dal Pino and A. Lazarian

In the presence of strong density stratification, turbulence can lead to a large-scale instability of a horizontal magnetic field if its strength is in a suitable range (within a few percent of the turbulent equipartition value). This instability is related to a suppression of the turbulent pressure so that the turbulence contribution to the mean magnetic pressure becomes negative. This results in the excitation of a negative effective magnetic pressure instability (NEMPI). This instability has so far only been studied for an imposed magnetic field. We want to know how NEMPI works when the mean magnetic field is generated self-consistently by an $\alpha^2$ dynamo, whether it is affected by global spherical geometry, and whether it can influence the properties of the dynamo itself. We adopt the mean-field approach which has previously been shown to provide a realistic description of NEMPI in direct numerical simulations. We assume axisymmetry and solve the mean-field equations with the Pencil-Code for an adiabatic stratification at a total density contrast in the radial direction of approximately 4 orders of magnitude. NEMPI is found to work when the dynamo-generated field is about 4% of the equipartition value, which is achieved through strong $\alpha$ quenching. This instability is excited in the top 5% of the outer radius provided the density contrast across this top layer is at least 10. NEMPI is found to occur at lower latitudes when the mean magnetic field is stronger. For weaker fields, NEMPI can make the dynamo oscillatory with poleward migration. NEMPI is a viable mechanism for producing magnetic flux concentrations in a strongly stratified spherical shell in which a magnetic field is generated by a strongly quenched $\alpha$ effect dynamo., Comment: published in A&A

Studies of solar and stellar convection often employ simple polytropic setups using the diffusion approximation instead of solving the proper radiative transfer equation. This allows one to control separately the polytropic index of the hydrostatic reference solution, the temperature contrast between top and bottom, and the Rayleigh and Peclet numbers. We extend such studies by including radiative transfer in the gray approximation using a Kramers-like opacity with freely adjustable coefficients. We study the properties of such models and compare them with results from the diffusion approximation. We use the Pencil Code, which is a high-order finite difference code where radiation is treated using the method of long characteristics. The source function is given by the Planck function. The opacity is written as kappa=kappa_0 rho^a T^b, where b is varied from -3.5 to +5, and kappa_0 is varied by four orders of magnitude. We consider sets of one dimensional models and perform a comparison with the diffusion approximation. Except for the case where b=5, we find one-dimensional hydrostatic equilibria with a nearly polytropic stratification and a polytropic index close to n=(3-b)/(1+a), covering both convectively stable (n>3/2) and unstable (n, 16 pages, 17 figures, submitted to A&A

We report on turbulent dynamo simulations in a spherical wedge with an outer coronal layer. We apply a two-layer model where the lower layer represents the convection zone and the upper layer the solar corona. This setup is used to study the coronal influence on the dynamo action beneath the surface. Increasing the radial coronal extent gradually to three times the solar radius and changing the magnetic Reynolds number, we find that dynamo action benefits from the additional coronal extent in terms of higher magnetic energy in the saturated stage. The flux of magnetic helicity can play an important role in this context., 4 pages, 1 figure, To appear in proceedings of IAUS 302: Magnetic fields throughout stellar evolution (August 2013, Biarritz, France)

We consider mean-field dynamo models with fluctuating ��effect, both with and without shear. The ��effect is chosen to be Gaussian white noise with zero mean and given covariance. We show analytically that the mean magnetic field does not grow, but, in an infinitely large domain, the mean-squared magnetic field shows exponential growth of the fastest growing mode at a rate proportional to the shear rate, which agrees with earlier numerical results of Yousef et al (2008) and recent analytical treatment by Heinemann et al (2011) who use a method different from ours. In the absence of shear, an incoherent ��^2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third and fourth order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian., Mysteriously missing a line in introduction, ("Heinemann et al. (2011) have recently proposed a simple analytically tractable model ..") which has now been added

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNS) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in formation of large-scale inhomogeneous magnetic structures from initial uniform magnetic field. Earlier investigations of this instability in DNS of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L^2/eta_t, where eta_t is the turbulent diffusivity and L is the scale of the domain. DNS are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number., Comment: 14 pages, 8 figures, submitted to special issue "Advances of European Solar Physics" in Solar Physics

At large magnetic Reynolds numbers, magnetic helicity evolution plays an important role in astrophysical large-scale dynamos. The recognition of this fact led to the development of the dynamical alpha quenching formalism, which predicts catastrophically low mean fields in open systems. Here we show that in oscillatory alpha-Omega dynamos this formalism predicts an unphysical magnetic helicity transfer between scales. An alternative technique is proposed where this artifact is removed by using the evolution equation for the magnetic helicity of the total field in the shearing--advective gauge. In the traditional dynamical alpha quenching formalism, this can be described as an additional magnetic helicity flux of small-scale fields that does not appear in homogeneous alpha-squared dynamos. In alpha-Omega dynamos, the alternative formalism is shown to lead to larger saturation fields than previously obtained with the traditional formalism., Comment: 9 pages, 7 figures, Published, ApJ

The decay of kinetic helicity is studied in numerical models of forced turbulence using either an externally imposed forcing function as an inhomogeneous term in the equations or, alternatively, a term linear in the velocity giving rise to a linear instability. The externally imposed forcing function injects energy at the largest scales, giving rise to a turbulent inertial range with nearly constant energy flux while for linearly forced turbulence the spectral energy is maximum near the dissipation wavenumber. Kinetic helicity is injected once a statistically steady state is reached, but it is found to decay on a turbulent time scale regardless of the nature of the forcing and the value of the Reynolds number., Comment: 7 pages, 8 figures, Astron. Nachr. (accepted)

Autocatalysis is a fundamental concept, used in a wide range of domains. From the most general definition of autocatalysis, that is a process in which a chemical compound is able to catalyze its own formation, several different systems can be described. We detail the different categories of autocatalyses, and compare them on the basis of their mechanistic, kinetic, and dynamic properties. It is shown how autocatalytic patterns can be generated by different systems of chemical reactions. The notion of autocatalysis covering a large variety of mechanistic realisations with very similar behaviors, it is proposed that the key signature of autocatalysis is its kinetic pattern expressed in a mathematical form., 10 pages, 5 figures, Submitted to J Phys Chem A

We investigate the emergence of a large-scale magnetic field. This field is dynamo-generated by turbulence driven with a helical forcing function. Twisted arcade-like field structures are found to emerge in the exterior above the turbulence zone. Time series of the magnetic field structure show recurrent plasmoid ejections., 2 pages, 2 figures, conference proceeding

Magnetic buoyancy is believed to drive the transport of magnetic flux tubes from the convection zone to the surface of the Sun. The magnetic fields form twisted loop-like structures in the solar atmosphere. In this paper we use helical forcing to produce a large-scale dynamo-generated magnetic field, which rises even without magnetic buoyancy. A two layer system is used as computational domain where the upper part represents the solar atmosphere. Here, the evolution of the magnetic field is solved with the stress--and--relax method. Below this region a magnetic field is produced by a helical forcing function in the momentum equation, which leads to dynamo action. We find twisted magnetic fields emerging frequently to the outer layer, forming arch-like structures. In addition, recurrent plasmoid ejections can be found by looking at space--time diagrams of the magnetic field. Recent simulations in spherical coordinates show similar results., 4 pages, 8 figures, To appear in the proceedings of the IAU273 "Physics of Sun and Star Spots"

The degree of intermittency of the magnetic field of a large-scale dynamo is considered. Based on simulations it is argued that there is a tendency for the field to become more diffuse and non-intermittent as the dynamo saturates. The simulations are idealized in that the turbulence is strongly helical and shear is strong, so the tendency for the field to become more diffuse is somewhat exaggerated. Earlier results concerning the effects of magnetic buoyancy are discussed. It is emphasized that the resulting magnetic buoyancy is weak compared with the stronger effects of simultaneous downward pumping. These findings are used to support the notion that the solar dynamo might operate in a distributed fashion where the near-surface shear layer could play an important role., Comment: 10 pages, 7 figures, to be published in Solar-stellar dynamos as revealed by helio- and asteroseismology, ed. M. Dikpati et al

Several aspects of mathematical astrobiology are discussed. It is argued that around the time of the origin of life the handedness of biomolecules must have established itself through an instability. Possible pathways of producing a certain handedness include mechanisms involving either autocatalysis or, alternatively, epimerization as governing effects. Concepts for establishing hereditary information are discussed in terms of the theory of hypercycles. Instabilities toward parasites and possible remedies by invoking spatial extent are reviewed. Finally, some effects of early life are discussed that contributed to modifying and regulating atmosphere and climate of the Earth, and that could have contributed to the highly oxidized state of its crust., Comment: 32 pages, 16 figures, to appear in Encyclopedia of Complexity and System Science

Star-disc coupling is considered in numerical models where the stellar field is not an imposed perfect dipole, but instead a more irregular self-adjusting dynamo-generated field. Using axisymmetric simulations of the hydromagnetic mean-field equations, it is shown that the resulting stellar field configuration is more complex, but significantly better suited for driving a stellar wind. In agreement with recent findings by a number of people, star-disc coupling is less efficient in braking the star than previously thought. Moreover, stellar wind braking becomes equally important. In contrast to a perfect stellar dipole field, dynamo-generated stellar fields favor field-aligned accretion with considerably higher velocity at low latitudes, where the field is weaker and originating in the disc. Accretion is no longer nearly periodic (as it is in the case of a stellar dipole), but it is more irregular and episodic., Comment: 19 pages, 15 figures, some errors corrected, Astron. Nach. (submitted). For higher quality images, see http://www.nordita.dk/~brandenb/papers/stellardyn/

The location of the solar dynamo is discussed in the context of new insights into the theory of nonlinear turbulent dynamos. It is argued that, from a dynamo-theoretic point of view, the bottom of the convection zone is not a likely location and that the solar dynamo may be distributed over the convection zone. The near surface shear layer produces not only east-west field alignment, but it also helps the dynamo disposing of its excess small scale magnetic helicity., Comment: 6 pages, 3 figures, in Solar MHD: Theory and Observations, eds: J. W. Leibacher, H. Uitenbroek, & R. F. Stein, Astron. Soc. Pac. Conf. Ser

The excitation of gravity waves by penetrative convective plumes is investigated using 2D direct simulations of compressible convection. The oscillation field is measured by a new technique based on the projection of our simulation data onto the theoretical g-modes solutions of the associated linear eigenvalue problem. This allows us to determine both the excited modes and their corresponding amplitudes accurately., Comment: 4 pages, 2 figures, proceeding of the Porto conference "Asteroseismology Across the HR Diagram", edited by M.J. Thompson, M.S Cunha & M.J.P.G. Monteiro, Kluwer Academic Publishers (also in Ap&SS, 284, 237)

The dimensionless kinetic energy dissipation rate C_epsilon is estimated from numerical simulations of statistically stationary isotropic box turbulence that is slightly compressible. The Taylor microscale Reynolds number Re_lambda range is 20 < Re_lambda < 220 and the statistical stationarity is achieved with a random phase forcing method. The strong Re_lambda dependence of C_epsilon abates when Re_lambda approx. 100 after which C_epsilon slowly approaches approx 0.5 a value slightly different to previously reported simulations but in good agreement with experimental results. If C_epsilon is estimated at a specific time step from the time series of the quantities involved it is necessary to account for the time lag between energy injection and energy dissipation. Also, the resulting value can differ from the ensemble averaged value by up to +-30%. This may explain the spread in results from previously published estimates of C_epsilon., 7 pages, 7 figures. Submitted to Phys. Rev. E

The advantages of high-order finite difference scheme for astrophysical MHD and turbulence simulations are highlighted. A number of one-dimensional test cases are presented ranging from various shock tests to Parker-type wind solutions. Applications to magnetized accretion discs and their associated outflows are discussed. Particular emphasis is placed on the possibility of dynamo action in three-dimensional turbulent convection and shear flows, which is relevant to stars and astrophysical discs. The generation of large scale fields is discussed in terms of an inverse magnetic cascade and the consequences imposed by magnetic helicity conservation are reviewed with particular emphasis on the issue of alpha-quenching., Comment: 68 pages, 36 figures, to be published in "The Fluid Mechanics of Astrophysics and Geophysics", Vol. 8, ed. Antonio Ferriz-Mas

The emergence of a large scale magnetic field from randomly forced isotropic strongly helical flows is discussed in terms of the inverse cascade of magnetic helicity and the alpha-effect. In simulations of such flows the maximum field strength exceeds the equipartition field strength for large scale separation. However, helicity conservation controls the speed at which this final state is reached. In the presence of open boundaries magnetic helicity fluxes out of the domain are possible. This reduces the timescales of the field growth, but it also tends to reduce the maximum attainable field strength., Comment: 8 pages, 4 figs, to be published in "Dynamo and dynamics, a mathematical challenge", Ed: P. Chossat, Nato ASI Series

We investigate the behaviour of $\alpha\Omega$ dynamos with a dynamic $\alpha$, whose evolution is governed by the imbalance between a driving and a damping term. We focus on truncated versions of such dynamo models which are often studied in connection with solar and stellar variability. Given the approximate nature of such models, it is important to study how robust they are with respect to reasonable changes in the formulation of the driving and damping terms. For each case, we also study the effects of changes of the dynamo number and its sign, the truncation order and initial conditions. Our results show that changes in the formulation of the driving term have important consequences for the dynamical behaviour of such systems, with the detailed nature of these effects depending crucially on the form of the driving term assumed, the value and the sign of the dynamo number and the initial conditions. On the other hand, the change in the damping term considered here seems to produce little qualitative effect., Comment: Figures included. Also in http://www.maths.qmw.ac.uk/~eoc

We study axisymmetric mean field spherical and spherical shell dynamo models, with both dynamic and algebraic $��$--quenchings. Our results show that there are qualitative as well as quantitative differences and similarities between these models. Regarding similarities, both groups of models exhibit symmetric, antisymmetric and mixed modes of behaviour. As regards differences, the important feature in the full sphere models is the occurrence of chaotic behaviour in the algebraic $��$--quenching models. For the spherical shell models with dynamic $��$ the main features include the possibility of multi-attractor regimes with final state sensitivity with respect to small changes in the magnitude of $��$ and the initial parity. We find the effect of introducing a dynamic $��$ is likely to be complicated and depend on the region of the parameter space considered, rather than a uniform change towards simplicity or complexity., Latex with full figures, also available on web page http://www.maths.qmw.ac.uk/~eoc; Astronomy and Astrophysics, accepted

We use the mean-field dynamo equations to show that an incoherent alpha effect in mirror-symmetric turbulence in a shearing flow can generate a large scale, coherent magnetic field. We illustrate this effect with simulations of a few simple systems. In accretion disks, this process can lead to axisymmetric magnetic domains whose radial and vertical dimensions will be comparable to the disk height. This process may be responsible for observations of dynamo activity seen in simulations of dynamo-generated turbulence involving, for example, the Balbus-Hawley instability. In this case the magnetic field strength will saturate at $\sim (h/r)^2$ times the ambient pressure in real accretion disks. The resultant dimensionless viscosity will be of the same order. In numerical simulations the azimuthal extent of the simulated annulus should be substituted for $r$. We compare the predictions of this model to numerical simulations previously reported by Brandenburg et al. (1995). In a radiation pressure dominated environment this estimate for viscosity should be reduced by a factor of $(P_{gas}/P_{radiation})^6$ due to magnetic buoyancy., 23 pages, uses aastex