ERC Adv. Grant Project 227952 AstroDyn

Second reporting period: 1 August 2010 - 31 January 2012 (months 19-36)

Project's full title: Astrophysical dynamos

Name of the PI: Axel Brandenburg

Name of PI's host institution for project: Nordita

Project website:

Management of relation between PI and HI:

Professor Larus Thorlacius has been the administrative officer since 2010. Since August 2011, Ms Marianne Persson Söderlind has been head of administration at Nordita and has managed Nordita's budget including that of the AstroDyn project.

Expenses incurred during reporting period 2:

Following the original plan, 4 PhD students, hired in 2009, are still being supported by the grant. All 4 received the licentiate degree in 2011 and are expected to finish in 2013. Two of the post-docs who started in 2009 finished in August 2011. One post-doc has been employed for 6 months (A.Hubbard). The assistant professor (Mitra) who started in 2010 continues until the end of the project, and the visiting professor (Rheinhardt) who started in 2009 had been prolonged for a second year and finished in August 2011. As already reported to you after sending you Report 1, this led to a shift of the planned expenses from Periods III and IV to Period II. According to the original plan, the expenses for post-docs where spread equally over the 5 year grant period, corresponding to a fraction of 1.5/5=0.3 of the total expense. This fraction now increased to 0.5. Likewise, the expenses for a visiting professor where spread equally over the 5 year grant period, corresponding here to a fraction of 0.2 of the total expense for this item. This fraction now increased to 0.35. This shift in expenses is justified by the enhanced scientific accomplishments in Period II. Indeed, as can be seen from Table 1 detailing the status of completion, we estimate that we have now completed the project to 70%. According to the original plan, a new post-doc should have been hired, but owing the reduced value of the Euro in Sweden (see figure below), this has not been done.

Personnel paid during reporting period 2:

S. Candelaresi 18 months PhD student F. Del Sordo 18 months PhD student K. Kemel 18 months PhD student J. Warnecke 18 months PhD student P. Chatterjee 12 months post-doc G. Guerrero 12 months post-doc A. Hubbard 6 months post-doc M. Rheinhardt 12.5 months visiting professor D. Mitra 18 months assistant professor

We have made some savings in equipment in consumables, and publication costs to cover the enhanced travel costs that occurred owing to an increased demand to attend meetings and to present our results.

**Objectives for the reporting period and corresponding achievements**

The following 14 items in italics are excerpts from the original proposal of 2008. The present status of the achievements is described for all items. All the papers that are quoted acknowledge the ERC grant. The results of Period I are summarized in small letters and documented by 26 peer-reviewed papers plus 8 additional papers that acknowledge the ERC, but whose scope falls outside the originally anticipated goals. The results of Period II are summarized in normal letters and documented by 34 peer-reviewed papers plus 10 additional papers addressing other goals.

*Code validation. Continue testing the spherical extension of the PENCIL CODE by comparing with other codes. Much of this has already been completed successfully, but there are some issues connected with the treatment of boundary conditions in large scale effect dynamos where the comparison is not yet satisfactory.*(*Phase 1*)**Results from reporting period I**. The implementation of spherical geometry in the PENCIL CODE has been developed and a number of additional tests and code enhancements have been carried out by Drs Mitra and Plasson as well as Mr. Svedin. Further tests have been performed by Dr. Babkovskaia in connection with applications to turbulent combustion. The following papers where mentioned in reporting period I:- Mitra, D., Tavakol, R., Brandenburg, A., & Moss, D.: 2009, ``Turbulent dynamos in spherical shell segments of varying geometrical extent,''
*Astrophys. J.***697**, 923-933 (arXiv:0812.3106, ADS, PDF) - Babkovskaia, N., Haugen, N. E. L., Brandenburg, A.: 2011, ``A high-order public domain code for direct numerical simulations of turbulent combustion,''
*J. Comp. Phys.***230**, 1-12 (arXiv:1005.5301, ADS, DOI, PDF)

**New achievements**. During the PENCIL CODE User Meeting in New York, an important development was undertaken to restructure the mean-field modelling as part of a substructure underneath the magnetic field modules. The main scientific outcomes of this are reported below in connection with the formation of magnetic structures in stratified turbulence. In addition, Drs Chatterjee and Mitra have continued working on the implementation of the anelastic solver.- Mitra, D., Tavakol, R., Brandenburg, A., & Moss, D.: 2009, ``Turbulent dynamos in spherical shell segments of varying geometrical extent,''
*Nonlinear testfield method. Determination of the quenching of the nonlinear effect and the turbulent diffusivity by large scale magnetic fields using the testfield method. The importance of the small scale current helicity for the effect is still not entirely settled, so it is important to extend work along those lines. The idea is to calculate not only the response of each test field on the small scale velocity, but also on the small scale magnetic field. For this work the Cartesian configuration of the PENCIL CODE will be used.*(*Phase 1*)**Results from reporting period I**. A new test-field method has been developed by Rheinhardt & Brandenburg that is able to account for the effects of MHD background turbulence. This method has been tested for the Roberts flow and the paper has now appeared. In addition, we have identified pitfalls in determining the correct effect using the more traditional imposed-field method. This work has led to new possibilities for determining magnetic helicity fluxes both in mean-field and in direct simulations, which in turn has led to the realization that diffusive magnetic helicity fluxes can be more important than previously thought.- Rheinhardt, M., & Brandenburg, A.: 2010, ``Test-field method for mean-field coefficients with MHD background,''
*Astron. Astrophys.***520**, A28 (arXiv:1004.0689, ADS, DOI, PDF) - Hubbard, A., Del Sordo, F., Käpylä, P. J., & Brandenburg, A.: 2009, ``The effect with imposed and dynamo-generated magnetic fields,''
*Monthly Notices Roy. Astron. Soc.***398**, 1891-1899 (arXiv:0904.2773, ADS, PDF) - Käpylä, P. J., Korpi, M. J., & Brandenburg, A.: 2010, ``The alpha effect in rotating convection with sinusoidal shear,''
*Monthly Notices Roy. Astron. Soc.***402**, 1458-1466 (arXiv:0908.2423, ADS, PDF)

**New achievements**. We have now extended the determination of turbulent transport coefficients to contributions that do not depend on the magnetic field, but depend just on rotation (the Yoshizawa effect). In this connection we have for the first time determined numerically its quenching due to self-consistently generated magnetic fields. We have also extended the test-field method to irrotational flows, to helical shear flows, and to passive scalar transport. The following publications have recently emerged:- Brandenburg, A., & Rädler, K.-H.: 2012, ``Yoshizawa's cross-helicity effect and its quenching,''
*Geophys. Astrophys. Fluid Dyn.*, submitted (arXiv:1112.1237, HTML, PDF) - Rädler, K.-H., Brandenburg, A., Del Sordo, F., & Rheinhardt, M.: 2011, ``Mean-field diffusivities in passive scalar and magnetic transport
in irrotational flows,''
*Phys. Rev. E***84**, 4 (arXiv:1104.1613, ADS, DOI, PDF) - Snellman, J. E., Rheinhardt, M., Käpylä, P. J., Mantere, M. J., & Brandenburg, A.: 2012, ``Mean-field closure parameters for passive scalar turbulence,''
*Phys. Scr.*, submitted (arXiv:1112.4777, HTML, PDF) - Snellman, J. E., Brandenburg, A., Käpylä, P. J., & Mantere, M. J.: 2012, ``Verification of Reynolds stress parameterizations from simulations,''
*Astron. Nachr.***333**, 78-83 (arXiv:1109.4857, ADS, DOI, PDF) - Rogachevskii, I., Kleeorin, N., Käpylä, P. J., & Brandenburg, A.: 2011, ``Pumping velocity in homogeneous helical turbulence with shear,''
*Phys. Rev. E***84**, 056314 (arXiv:1105.5785, ADS, DOI, PDF)

- Rheinhardt, M., & Brandenburg, A.: 2010, ``Test-field method for mean-field coefficients with MHD background,''
*Catastrophic quenching in a spherical shell. Reproduce the catastrophic quenching behavior in a closed sphere or spherical shell sector using perfectly conducting boundary conditions and forced turbulence. Some work in this direction has already been done, but the results are not yet well understood nor entirely conclusive.*(*Phase 1*)**Results from reporting period I**. In pursuit of this problem, Dr. Mitra has come across a new type of solution that yields equatorward migration even without shear. This result is quite surprising and has now been published. Additional work is in progress and has been combined with dynamos in spherical shells (item 7 in this list).- Mitra, M., Tavakol, R., Käpylä, P. J., & Brandenburg, A.: 2010, ``Oscillatory migrating magnetic fields in helical turbulence
in spherical domains,''
*Astrophys. J. Lett.***719**, L1-L4 (arXiv:0901.2364, PDF)

**New achievements**. To understand catastrophic quenching, we need to understand the magnetically generated effect, whose value is characterized by magnetic helicity, which is a conserved quantity at large magnetic Reynolds numbers. There is the possibility that there might be higher order invariants that could also matter. This possibility has been tested in a recent paper where some evidence for the need of higher order invariants has been found.- Candelaresi, S., & Brandenburg, A.: 2011, ``Decay of helical and non-helical magnetic knots,''
*Phys. Rev. E***84**, 016406 (arXiv:1103.3518, ADS, DOI, PDF)

- Mitra, M., Tavakol, R., Käpylä, P. J., & Brandenburg, A.: 2010, ``Oscillatory migrating magnetic fields in helical turbulence
in spherical domains,''
*Dynamo effect from the MRI. Calculate the nonlinear effect and the turbulent diffusivity for turbulence driven by the magneto-rotational instability (MRI). Some work in this direction has already been done, but only a few representative test cases at relatively low resolution were done. This work is primarily relevant to accretion discs. However, understanding this case may also teach us general aspects of magnetically driven dynamos that may in some form also work in the Sun.*(*Phase 1*)**Results from reporting period I**. This work has been started with the help of a student from the ENS in Paris, Emeric Bron, who visited us for 1/2 year on an internship. Preparatory work on this topic has already been published with another student who also came on an internship. Both works have led to new issues concerning the importance of using open boundary conditions for the magnetic field. This has also led to new work that helped resolving the question of the dependence of the onset of MRI on the value of the magnetic Prandtl number. In now turns out that with open boundary conditions the onset is independent of the magnetic Prandtl number. This following work appeared in the previous reporting period.- Vermersch, V., & Brandenburg, A.: 2009, ``Shear-driven magnetic buoyancy oscillations,''
*Astron. Nachr.***330**, 797-806 (arXiv:0909.0324, ADS, PDF) - Käpylä, P. J., & Korpi, M. J.: 2011, ``Magnetorotational instability driven dynamos at low magnetic Prandtl numbers,''
*Monthly Notices Roy. Astron. Soc.***413**, 901-907 (arXiv:1004.2417, ADS

**New achievements**. The possibility of explaing MRI dynamo action as the result of an incoherent -shear dynamo has now been explored further and has been shown to obey the observed linear scaling of growth rate with shear rate. The MRI has now been explored further in systems where the Hall effect is important.- Mitra, D., & Brandenburg, A.: 2012, ``(,''
*Monthly Notices Roy. Astron. Soc.***420**, 2170-Scaling and intermittency in incoherent -shear dynamoarXiv:1107.2419, DOI, HTML, PDF) - Bejarano, C., Gomez, D. O., & Brandenburg, A.: 2011, ``Shear-driven instabilities in Hall-MHD plasmas,''
*Astrophys. J.***737**, 62 (arXiv:1012.5284, ADS, DOI, PDF)

**Main objectives originally scheduled for the second reporting period**:- Vermersch, V., & Brandenburg, A.: 2009, ``Shear-driven magnetic buoyancy oscillations,''
*Testfield method in spherical geometry. Adapt the testfield method to spherical coordinates. Originally the testfield method was developed in connection with full spheres, and then the testfields consisted of field components of constant value or constant slope. However, only afterwards it became clear that the scale (or wavenumber) of the field components must be the same for one set of all tensor components, and so it is necessary to work with spherical harmonic functions as testfields. In other words, constant and linearly varying field components are insufficient.*(*Phase 2*)**Results from reporting period I**. In preparation of this task, Dr. Mitra has implemented a helical flow in spherical geometry that will allow us to validate the testfield method in spherical geometry. We plan to implement this method soon.**New achievements**. This topic had to be postponed because of recent advances in the field of magnetic flux concentrations (item 8 below).*Alpha effect from convection. The calculation of the effect in convective turbulence is at the moment rather unclear. There are some results suggesting that goes to zero in the limit of large magnetic Reynolds numbers even for kinematically weak magnetic fields. There remain however several open questions regarding the amount of stratification ( should be proportional to the local stratification gradient and should hence be absent in Boussinesq convection) and regarding the degree of scale separation.*(*Phase 2*)**Results from reporting period I**. In the mean time the situation has changed dramatically. Significant progress in this direction has been made by Dr. Käpylä and collaborators in a series of papers:- Käpylä, P. J., Korpi, M. J., & Brandenburg, A.: 2009, ``Alpha effect and turbulent diffusion from convection,''
*Astron. Astrophys.***500**, 633-646 (arXiv:0812.1792, ADS, PDF)

**New achievements**. We have now found a surprising occurrence of oscillatory large-scale dynamo action from Cartesian convection simulations. This work has now been submitted.- Käpylä, P. J., Mantere, M. J., & Brandenburg, A.: 2012, ``Oscillatory large-scale dynamos from Cartesian convection simulations,''
*Geophys. Astrophys. Fluid Dyn.*, submitted (arXiv:1111.6894, PDF)

- Käpylä, P. J., Korpi, M. J., & Brandenburg, A.: 2009, ``Alpha effect and turbulent diffusion from convection,''
*Dynamo in open shells with and without shear. Calculate the saturation of the magnetic field and the underlying dynamo effects with open boundary conditions in a spherical shell sector with and without shear. One expects low saturation amplitude with magnetic energy of the mean field being inversely proportional to the magnetic Reynolds number in the absence of shear, but of order unity in the presence of shear. The shear is here critical, because it is responsible for the local driving of small scale magnetic helicity fluxes.*(*Phase 2*)**Results from reporting period I**. Significant progress has also been made in understanding the nature of magnetic helicity and its fluxes. Additional items connected with understanding magnetic helicity fluxes.- Mitra, D., Käpylä, P. J., Tavakol, R., & Brandenburg, A.: 2009, ``Alpha effect and diffusivity in helical turbulence with shear,''
*Astron. Astrophys.***495**, 1-8 (arXiv:0806.1608, ADS, PDF) - Mitra, D., Candelaresi, S., Chatterjee, P., Tavakol, R., &
Brandenburg, A.: 2010, ``Equatorial magnetic helicity flux in simulations with different gauges,''
*Astron. Nachr.***331**, 130-135 (arXiv:0911.0969, ADS, PDF) - Del Sordo, F., Candelaresi, S., & Brandenburg, A.: 2010, ``Magnetic field decay of three interlocked flux rings with zero linking number,''
*Phys. Rev. E***81**, 036401 (arXiv:0910.3948, ADS, PDF) - Käpylä, P. J., Korpi, M. J., & Brandenburg, A.: 2010, ``Open and closed boundaries in large-scale convective dynamos,''
*Astron. Astrophys.***518**, A22 (arXiv:0911.4120, ADS, DOI, PDF) - Hubbard, A., & Brandenburg, A.: 2010, ``Magnetic helicity fluxes in an dynamo embedded in a halo,''
*Geophys. Astrophys. Fluid Dyn.***104**, 577-590 (arXiv:1004.4591, ADS, DOI, PDF) - Hubbard, A., & Brandenburg, A.: 2011, ``Magnetic helicity flux in the presence of shear,''
*Astrophys. J.***727**, 11 (arXiv:1006.3549, ADS, DOI, PDF)

**New achievements**. While the progress has been significant, our work also showed that the magnetic helicity fluxes are still small compared with microscopic diffusion unless the magnetic Reynolds number exceeds values around to . Most surprisingly, there is now evidence that the Vishniac-Cho flux may not exist. We have also established that in a stationary but open system the divergence of magnetic helicity fluxes (even separately for those of the small-scale field) is gauge-invariant if the magnetic helicity density is found to be statistically steady. The following papers have now be written:- Hubbard, A., Rheinhardt, M. & Brandenburg, A.: 2011, ``The fratricide of dynamos by their siblings,''
*Astron. Astrophys.***535**, A48 (arXiv:1102.2617, ADS, DOI, PDF) - Mitra, D., Moss, D., Tavakol, R., & Brandenburg, A.: 2011, ``Alleviating alpha quenching by solar wind and meridional flow,''
*Astron. Astrophys.***526**, A138 (arXiv:1008.4226, ADS, DOI, PDF) - Hubbard, A., & Brandenburg, A.: 2012, ``Catastrophic quenching in dynamos revisited,''
*Astrophys. J.*, in press (arXiv:1107.0238, HTML) - Brandenburg, A.: 2011, ``Chandrasekhar-Kendall functions in astrophysical dynamos,''
*Pramana J. Phys.***77**, 67-76 (arXiv:1103.4976, ADS, DOI, PDF) - Candelaresi, S., Hubbard, A., Brandenburg, A., & Mitra, D.: 2011, ``Magnetic helicity transport in the advective gauge family,''
*Phys. Plasmas***18**, 012903 (arXiv:1010.6177, ADS, DOI, PDF) - Kemel, K., Brandenburg, A., & Ji, H.: 2011, ``A model of driven and decaying magnetic turbulence in a cylinder,''
*Phys. Rev. E***84**, 056407 (arXiv:1106.1129, ADS, DOI, PDF)

**Objectives for the remainder of the grant period and status report, with emphasis on already obtained results**:- Mitra, D., Käpylä, P. J., Tavakol, R., & Brandenburg, A.: 2009, ``Alpha effect and diffusivity in helical turbulence with shear,''
*Magnetic flux concentrations near the surface. Test the scenario that the emergence of active regions and sunspots can be explained as the result of flux concentrations from local dynamo action via negative turbulent magnetic pressure effects or turbulent flux collapse.*(*Phase 2*)**Results from reporting period I**. This work constitutes one of the corner stones of our project in that we must explore scenarios for being able to explain the formation of magnetic flux concentrations in the absence of deep-rooted hypothetical flux loops at the bottom of the convection zone. This work has been started with Professors Kleeorin and Rogachevskii, as well as Mr. Kemel (one of our PhD students). One paper has now been published and one has been submitted.- Brandenburg, A., Kleeorin, N., & Rogachevskii, I.: 2010, ``Large-scale magnetic flux concentrations from turbulent stresses,''
*Astron. Nachr.***331**, 5-13 (arXiv:0910.1835, ADS, PDF) - Brandenburg, A., Kemel, K., Kleeorin, N., & Rogachevskii, I.: 2012, ``The negative effective magnetic pressure in stratified forced turbulence,''
*Astrophys. J.*, in press (arXiv:1005.5700, HTML, PDF)

**New achievements**. The field has now seen a major transformation with the detection of the negative effective magnetic pressure instability in turbulence simulations. This has been a major milestone for this project that contributed to putting this effect on the map. Furthermore, Dr. Käpylä has started look at the possibility of producing magnetic flux concentration in stratified convection.- Brandenburg, A., Kemel, K., Kleeorin, N., Mitra, D., & Rogachevskii, I.: 2011, ``Detection of negative effective magnetic pressure instability in turbulence simulations,''
*Astrophys. J. Lett.***740**, L50 (arXiv:1109.1270, ADS, DOI, HTML, PDF) - Käpylä, P. J., Brandenburg, A., Kleeorin, N., Mantere, M. J., & Rogachevskii, I.: 2012, ``Negative effective magnetic pressure in turbulent convection,''
*Monthly Notices Roy. Astron. Soc.*, in press (arXiv:1104.4541, HTML, PDF) - Kemel, K., Brandenburg, A., Kleeorin, N., Mitra, D., & Rogachevskii, I.: 2012, ``Spontaneous formation of magnetic flux concentrations in stratified turbulence,''
*Solar Phys.*, in press (arXiv:1112.0279, HTML, PDF) - Kemel, K., Brandenburg, A., Kleeorin, N., & Rogachevskii, I.: 2012, ``Properties of the negative effective magnetic pressure instability,''
*Astron. Nachr.***333**, 95-100 (arXiv:1107.2752, DOI, PDF)

- Brandenburg, A., Kleeorin, N., & Rogachevskii, I.: 2010, ``Large-scale magnetic flux concentrations from turbulent stresses,''
*CME-like features above the surface. Analyze the nature of the expelled magnetic field in simulations that couple to a simplified version of the lower solar wind. It is possible that the magnetic field above the surface might resemble coronal mass ejections (CMEs), in which case more detailed comparisons with actual coronal mass ejections would be beneficial.*(*Phase 3*)**Results from reporting period I**. This project has been started with Mr. Jörn Warnecke, one of our PhD students who arrived in August 2009. Our first steps in this direction include a simple Cartesian model with a force-free outer layer above the turbulence zone. A revised version of the paper

- Warnecke, J., & Brandenburg, A.: 2010, ``Surface appearance of dynamo-generated large-scale fields,''
*Astron. Astrophys.***523**, A19 (arXiv:1002.3620, ADS, DOI, PDF)

**New achievements**. As planned, this work has now been extended to include spherical geometry. We have also now included convection. Using observations from the Ulysses space craft we have for the first time determined the magnetic helicity spectrum in the solar wind. Similar results have now also been verified using the simulations of Warnecke et al.- Warnecke, J., Käpylä, P. J., Mantere, M. J., & Brandenburg, A.: 2012, ``Coronal ejections driven by a convective dynamo with differential rotation in a sphere,''
*Solar Phys.*, submitted (arXiv:1112.0505, HTML, PDF) - Warnecke, J., Brandenburg, A., & Mitra, D.: 2011, ``Dynamo-driven plasmoid ejections above a spherical surface,''
*Astron. Astrophys.***534**, A11 (arXiv:1104.0664, ADS, DOI, HTML, PDF) - Brandenburg, A., Subramanian, K., Balogh, A., & Goldstein, M. L.: 2011, ``Scale-dependence of magnetic helicity in the solar wind,''
*Astrophys. J.***734**, 9 (arXiv:1101.1709, ADS, DOI, PDF)

- Warnecke, J., & Brandenburg, A.: 2010, ``Surface appearance of dynamo-generated large-scale fields,''
*Convective dynamo in spherical shell. Set up convection in the spherical shell. If the resulting scale of the flow is small enough and there is scale separation it would be useful to simulate the resulting magnetic field, compare with forced turbulence simulations in spherical shells and see whether contact can be made both with the Sun and with improved mean field models.*(*Phase 3*)**Results from reporting period I**. this work has been started under the initiative of Dr. Käpylä and first results have been published.- Käpylä, P. J., Korpi, M. J., & Brandenburg, A., Mitra, D.,
& Tavakol, R.: 2010, ``Convective dynamos in spherical wedge geometry,''
*Astron. Nachr.***331**, 73-81 (arXiv:0909.1330, ADS, PDF)

**New achievements**. We have now obtained convection-driven dynamo in spherical wedges. Selfconsistently driven differential rotation has been studied for different degree of stratification and rotation rates.- Käpylä, P. J., Mantere, M. J., & Brandenburg, A.: 2011, ``Effects of stratification in spherical shell convection,''
*Astron. Nachr.***332**, 883-890 (arXiv:1109.4625, ADS, DOI, PDF) - Käpylä, P. J., Mantere, M. J., Guerrero, G., Brandenburg, A., & Chatterjee, P.: 2011, ``Reynolds stress and heat flux in spherical shell convection,''
*Astron. Astrophys.***531**, A162 (arXiv:1010.1250, ADS, DOI, PDF)

- Käpylä, P. J., Korpi, M. J., & Brandenburg, A., Mitra, D.,
& Tavakol, R.: 2010, ``Convective dynamos in spherical wedge geometry,''
*Buoyancy-driven dynamo. The turbulence in accretion discs is believed to be driven by the magnetorotational instability. It was one of the first examples showing cyclic dynamo action somewhat reminiscent of the solar dynamo. It was believed to be a prototype of magnetically driven dynamos. In the mean time, another example of a magnetically driven dynamo has emerged, where magnetic buoyancy works in the presence of shear and stratification alone. This phenomenon is superficially similar to a magnetically dominated version of the shear-current effect. We are now in a good position to identify the governing mechanism by using the recently developed testfield method.*(*Phase 3*)**Results from reporting period I**. Dr. Chatterjee has started with this project. First results have been reported in the following review paper:- Brandenburg, A., Chatterjee, P., Del Sordo, F., Hubbard, A., Käpylä, P. J., & Rheinhardt, M.: 2010, ``Turbulent transport in hydromagnetic flows,''
*Phys. Scr.***T142**, 014028 (arXiv:1004.5380, ADS, DOI, PDF)

**New achievements**. Surprisingly, we have found that in a completely mirror-symmetric system, an effect can still emerge as a result of spontaneous symmetry breaking. Furthermore, Drs Guerrero and Käpylä have modeled convection with a strong shear layer (tachocline) at the bottom. They find the emergence of flux tubes at the top of the domain, but the field has become rather weak by the time it reaches the surface.- Chatterjee, P., Mitra, D., Rheinhardt, & M. Brandenburg, A.: 2011, ``Alpha effect due to buoyancy instability of a magnetic layer,''
*Astron. Astrophys.***534**, A46 (arXiv:1011.1218, ADS, DOI, PDF) - Chatterjee, P., Mitra, D., Brandenburg, A., & Rheinhardt, M.: 2011, ``Spontaneous chiral symmetry breaking by hydromagnetic buoyancy,''
*Phys. Rev. E***84**, 025403R (arXiv:1011.1251, ADS, DOI, PDF) - Guerrero, G., & Kpylä, P. J.: 2011, ``Dynamo action and magnetic buoyancy in convection simulations with vertical shear,''
*Astron. Astrophys.***533**, A40

- Brandenburg, A., Chatterjee, P., Del Sordo, F., Hubbard, A., Käpylä, P. J., & Rheinhardt, M.: 2010, ``Turbulent transport in hydromagnetic flows,''
*Deep convection dynamo. The deeper layers of the Sun are characterized by rather low values of the energy flux relative to the natural units given by , where is the density and is the speed of sound. In the Sun this is accomplished by nearly perfectly adiabatic conditions, which implies low Mach numbers on the order of . Such conditions cannot be economically simulated with compressible codes, so it is necessary to turn to an anelastic configuration of the PENCIL CODE. This should not be so hard to do because a Poisson solver has already been implemented in connection with solving for self-gravitating flows. Another possibility would be multigrid solvers. One such multigrid solver is also already present in the PENCIL CODE, but this subroutine still need to be parallelized. In discussions with Professor J Toomre from Boulder concerning near-future peta-flop computing it became clear that there is great interest in mesh-based codes that are able to solve anelastic flows in spherical shells.*(*Phase 4*)**Results from reporting period I**. Dr. Chatterjee has started implementing an anelastic solver into the PENCIL CODE. This work has been presented at the last PENCIL CODE User Meeting in New York (http://www.nordita.org/software/pencil-code/UserMeetings/2010/).**New achievements**. Much of this topic had to be postponed because of recent advances in the field of magnetic flux concentrations (item 8 below). However, we did make some progress in understanding the effect of strong stratification in the stably stratified tachocline and the layers beneath.- Kitchatinov, L. L., & Brandenburg, A.: 2012, ``Transport of angular momentum and chemical species by anisotropic mixing in stellar radiative interiors,''
*Astron. Nachr.*, in press (arXiv:1201.2484, HTML, PDF) - Brandenburg, A., Rädler, K.-H., & Kemel, K.: 2012, ``Mean-field transport in stratified and/or rotating turbulence,''
*Astron. Astrophys.***539**, A35 (arXiv:1108.2264, HTML)

- Kitchatinov, L. L., & Brandenburg, A.: 2012, ``Transport of angular momentum and chemical species by anisotropic mixing in stellar radiative interiors,''
*Solar dynamo models and solar cycle forecast. Among the popular applications of solar dynamo theory and solar magnetohydrodynamics are solar cycle predictions, solar subsurface weather, and space weather. Also of interest are predictions of solar activity during its first 500 thousand years. This has great relevance for predicting the loss of volatile elements from the Earth's atmosphere, for example, and for understanding the conditions on Earth during the time when life began colonizing the planet. In this connection it is also of interest to calculate the deflection of cosmic ray particles by the Sun's magnetic field and on the scale of the galaxy which is relevant for galactic cosmic rays.*(*Phase 4*)**Results from reporting period I**. In a preparatory step of this work, Mr. Svedin has started developing a data assimilation package for the PENCIL CODE. The first steps of this work are currently being written up. For future models of the solar dynamo, the effects of magnetic helicity fluxes have now been studied in more detail both in Cartesian as well as on spherical mean-field models:- Brandenburg, A., Candelaresi, S., & Chatterjee, P.: 2009, ``Small-scale magnetic helicity losses from a mean-field dynamo,''
*Monthly Notices Roy. Astron. Soc.***398**, 1414-1422 (arXiv:0905.0242, ADS, PDF) - Guerrero, G., Chatterjee, P., & Brandenburg, A.,: 2010, ``Shear-driven and diffusive helicity fluxes in dynamos,''
*Monthly Notices Roy. Astron. Soc.***409**, 1619-1630 (arXiv:1005.4818, ADS, DOI, PDF) - Chatterjee, P., Brandenburg, A., & Guerrero, G.: 2010, ``Can catastrophic quenching be alleviated by separating shear and effect?''
*Geophys. Astrophys. Fluid Dyn.***104**, 591-599 (arXiv:1005.5708, ADS, DOI, PDF) - Chatterjee, P., Guerrero, G., & Brandenburg, A.: 2011, ``Magnetic helicity fluxes in interface and flux transport dynamos,''
*Astron. Astrophys.***525**, A5 (arXiv:1005.5335, ADS, DOI, PDF)

**New achievements**. While the effects of the near-surface shear layer have still not been taken into account into much of this modeling, new insights have been gained by modelling the meridional advection of magnetic structures on the solar surface.- Guerrero, G., Rheinhardt, M., Brandenburg, A., & Dikpati, M.: 2012, ``Plasma flow vs. magnetic feature-tracking speeds in the Sun,''
*Monthly Notices Roy. Astron. Soc.***420**, L1-L5 (arXiv:1107.4801, ADS, DOI, HTML, PDF) - Rheinhardt, M., & Brandenburg, A.: 2012, ``Modeling spatio-temporal nonlocality in mean-field dynamos,''
*Astron. Nachr.***333**, 71-77 (arXiv:1110.2891, ADS, DOI, HTML, PDF)

- Brandenburg, A., Candelaresi, S., & Chatterjee, P.: 2009, ``Small-scale magnetic helicity losses from a mean-field dynamo,''
*Applications to laboratory liquid sodium dynamos. Unexpected beneficial insights have come from recent laboratory dynamo experiments. Unlike numerical dynamos, experimental liquid metal dynamos are able to address the regime of rather low values of the magnetic Prandtl number of the order of , which is of interest for solar and stellar conditions. At the same time the magnetic Reynolds number can be large enough (above 100) to allow for dynamo action. (This magnetic Prandtl number is not to be confused with the*(*turbulent*magnetic Prandtl number that was mentioned before in connection with flux transport dynamos!) The Cadarache experiment is particularly interesting to us. Simulations of this flow have been attempted by various groups using the Taylor-Green flow as a model. Again, the nature of the resulting dynamo effect has not yet been elucidated. It would be useful to analyze the resulting flows using the testfield method. It is hoped that such work can teach us important aspects about small-scale dynamos at low magnetic Prandtl number, which is relevant to the Sun, but hard to address numerically. Another relevant application is precession-driven dynamos. Preparatory simulations in Cartesian simulations have been carried out in collaboration with Agris Gailitis from Riga/Latvia.*Phase 4*)**New achievements**. There is significant hope to be able to determine for the first time the alpha effect in a turbulent liquid-metal plane Couette flow. Preparations have been performed through analytical and numerical calculations. We have also persued turbulent dynamo calculations at low magnetic Prandtl numbers that are relevant to liquid metal dynamos.- Rüdiger, G., & Brandenburg, A.: 2012, ``The alpha-effect in a turbulent liquid-metal plane Couette flow,''
*Phys. Rev. E*, submitted (arXiv:1201.0652, PDF) - Brandenburg, A.: 2011, ``Nonlinear small-scale dynamos at low magnetic Prandtl numbers,''
*Astrophys. J.***741**, 92 (arXiv:1106.5777, ADS, DOI, HTML, PDF) - Brandenburg, A.: 2011, ``Dissipation in dynamos at low and high magnetic Prandtl numbers,''
*Astron. Nachr.***332**, 51-56 (arXiv:1010.4805, ADS, DOI, PDF)

- Rüdiger, G., & Brandenburg, A.: 2012, ``The alpha-effect in a turbulent liquid-metal plane Couette flow,''

The ultimate goal of the project is of course to establish the cause of the
equatorward migration of magnetic activity belts at low solar latitudes.
Is it the rather feeble meridional circulation, as assumed in the
now rather popular flux transport models, even though one
has to assume unrealistic values of the turbulent magnetic Prandtl number,
or is it perhaps the near-surface shear layer, which would have
indeed the right sign?

Two reviews have been published that outline our current thinking:

- Brandenburg, A.: 2009, ``The critical role of magnetic helicity in astrophysical dynamos,''
*Plasma Phys. Control. Fusion***51**, 124043 (arXiv:0909.4377, PDF) - Brandenburg, A., & Nordlund, Å.: 2011, ``Astrophysical turbulence modeling,''
*Rep. Prog. Phys.***74**, 046901 (arXiv:0912.1340, ADS, DOI, PDF)

The success of our project is further evidenced by a number of
publications on other timely aspects of dynamo theory.
**Papers from reporting period I**.

- Rüdiger, G., Kitchatinov, L. L., & Brandenburg, A.: 2011, ``Cross helicity and turbulent magnetic diffusivity in the solar convection zone,''
*Solar Phys.***269**, 3-12 (arXiv:1004.4881, ADS, DOI, PDF) - Käpylä, P. J., Brandenburg, A., Korpi, M. J., Snellman, J. E., &
Narayan, R.: 2010, ``Angular momentum transport in convectively unstable shear flows,''
*Astrophys. J.***719**, 67-76 (arXiv:1003.0900, PDF) - Madarassy, E. J. M., & Brandenburg, A.: 2010, ``Calibrating passive scalar transport in shear-flow turbulence,''
*Phys. Rev. E***82**, 016304 (arXiv:0906.3314, ADS, PDF) - Kahniashvili, T., Brandenburg, A., Tevzadze, A. G., & Ratra, B.: 2010, ``Numerical simulations of the decay of primordial magnetic turbulence,''
*Phys. Rev. D***81**, 123002 (arXiv:1004.3084, ADS, PDF) - Brandenburg, A.: 2010, ``Magnetic field evolution in simulations with Euler potentials,''
*Monthly Notices Roy. Astron. Soc.***401**, 347-354 (arXiv:0907.1906, ADS, PDF) - Rädler, K.-H., & Brandenburg, A.: 2010, ``Mean electromotive force proportional to mean flow in mhd turbulence,''
*Astron. Nachr.***331**, 14-21 (arXiv:0910.0071, ADS, PDF) - Hubbard, A., & Brandenburg, A.: 2009, ``Memory effects in turbulent transport,''
*Astrophys. J.***706**, 712-726 (arXiv:0811.2561, ADS, PDF) - Sur, S., & Brandenburg, A.: 2009, ``The role of the Yoshizawa effect in the Archontis dynamo,''
*Monthly Notices Roy. Astron. Soc.***399**, 273-280 (arXiv:0902.2394, ADS, PDF)

- Chan, C. K., Mitra, D., & Brandenburg, A.: 2012, ``Dynamics of saturated energy condensation in two-dimensional turbulence,''
*Phys. Rev. E*, in press (arXiv:1109.6937, PDF) - Haugen, N. E. L., Kleeorin, N., Rogachevskii, I., & Brandenburg, A.: 2012, ``Detection of the phenomenon of turbulent thermal diffusion in numerical simulations,''
*Phys. Fluids*, submitted (arXiv:1101.4188, HTML, PDF) - Brandenburg, A., & Petrosyan, A.: 2012, ``Reynolds number dependence of kinetic helicity decay in linearly forced turbulence,''
*Astron. Nachr.*, submitted (arXiv:1012.1464, HTML) - Dosopoulou, F., Del Sordo, F., Tsagas, C. G., Brandenburg A.: 2012, ``Vorticity production and survival in viscous and magnetized cosmologies,''
*Phys. Rev. D*, in press (arXiv:1112.6164, HTML, PDF) - Rempel, E. L., Chian, A. C.-L., & Brandenburg, A.: 2012, ``Lagrangian chaos in an ABC-forced Nonlinear Dynamo,''
*Phys. Scr.*, in press (arXiv:1201.4324, HTML) - Rempel, E. L., Chian, A. C.-L., & Brandenburg, A.: 2011, ``Lagrangian coherent structures in a nonlinear dynamo,''
*Astrophys. J.***735**, L9 (arXiv:1011.6327, ADS, DOI, PDF) - Del Sordo, F., & Brandenburg, A.: 2011, ``Vorticity production through rotation, shear, and baroclinicity,''
*Astron. Astrophys.***528**, A145 (arXiv:1008.5281, ADS, DOI, PDF) - Brandenburg, A., Haugen, N. E. L., & Babkovskaia, N.: 2011, ``Turbulent front speed in the Fisher equation: dependence on Damköhler number,''
*Phys. Rev. E***83**, 016304 (arXiv:1008.5145, ADS, DOI, PDF) - Ray, S. S., Mitra, D. Perlekar, P. & Pandit, R.: 2011, ``Dynamic multiscaling in two-dimensional fluid turbulence,''
*Phys. Rev. Lett.***107**, 184503 (arXiv:1105.5160) - Perlekar,P., Ray, S. S., Mitra, D. & Pandit, R.: 2011, ``Persistence Problem in Two-Dimensional Fluid Turbulence,''
*Phys. Rev. Lett.***106**, 054501 (arXiv:1009.1494)

**3. Explanation of the use of resources**

A detailed working plan is given in the extended synopsis of Section 2.
We summarize here the intermediate goals as described in detail
in that section, where the goals were ordered by the phase within
the grant period.

- Code validation, nonlinear testfield method,
catastrophic quenching in a spherical shell,
dynamo effect from the MRI
(
**items 1-4 in Sect. 2**).*Estimated completion*: 65%, resources consumed in Period II: 225,000.00 EUR (246,000.00 in Period I). Completion date from Annex I: July 2010.*Comments*: other topics such as near-surface magnetic flux concentrations and CME-like features have become extremely timely, so some of the resources have been diverted to those topics. - Testfield method in spherical geometry, alpha effect from convection,
dynamo in open shells with and without shear
(
**items 5-7 in Sect. 2**).*Estimated completion*: 45%, resources consumed in Period II: 102,000.00 EUR (164,000.00 in Period I). Completion date from Annex I: February 2012.*Comments*: Studies of the alpha effect from convection have been nearly completed. - Magnetic flux concentrations near the surface,
CME-like features above the surface, convective dynamo in spherical shell,
buoyancy-driven dynamo
(
**items 8-11 in Sect. 2**).*Estimated completion*: 14%, resources consumed in Period II: 245,000.00 EUR (131,000.00 in Period I). Completion date from Annex I: July 2013.*Comments*: studies of near-surface magnetic flux concentrations and CME-like features have already been 40% completed. - Deep convection dynamo, solar cycle forecast,
applications to laboratory liquid sodium dynamos
(
**items 12-14 in Sect. 2**).*Estimated completion*: 6%, resources consumed in Period II: 143,000.00 EUR (33,000.00 in Period I). Completion date from Annex I: February 2014.*Comments*: a number of preparatory steps toward solar cycle forecast have been completed.