The presence of magnetic fields in many astrophysical objects is due to dynamo action, whereby a part of the kinetic energy is converted into magnetic energy. A turbulent dynamo that produces magnetic field structures on the same scale as the turbulent flow is known as the fluctuation dynamo. We use numerical simulations to explore the nonlinear, statistically steady state of the fluctuation dynamo in driven turbulence. We demonstrate that as the magnetic field growth saturates, its amplification and diffusion are both affected by the back-reaction of the Lorentz force upon the flow. The amplification of the magnetic field is reduced due to stronger alignment between the velocity field, magnetic field, and electric current density. Furthermore, we confirm that the amplification decreases due to a weaker stretching of the magnetic field lines. The enhancement in diffusion relative to the field line stretching is quantified by a decrease in the computed local value of the magnetic Reynolds number. Using the Minkowski functionals, we quantify the shape of the magnetic structures produced by the dynamo as magnetic filaments and ribbons in both kinematic and saturated dynamos and derive the scalings of the typical length, width, and thickness of the magnetic structures with the magnetic Reynolds number. We show that all three of these magnetic length scales increase as the dynamo saturates. The magnetic intermittency, strong in the kinematic dynamo (where the magnetic field strength grows exponentially) persists in the statistically steady state, but intense magnetic filaments and ribbons are more volume-filling.

Energy equipartition between cosmic rays and magnetic fields is often assumed to infer magnetic field properties from the synchrotron observations of star-forming galaxies. However, there is no compelling physical reason to expect the same. We aim to explore the validity of the energy equipartition assumption. After describing popular arguments in favour of the assumption, we first discuss observational results which support it at large scales and how certain observations show significant deviations from equipartition at scales smaller than $\approx 1 \, {\rm kpc}$, probably related to the propagation length of the cosmic rays. Then we test the energy equipartition assumption using test-particle and MHD simulations. From the results of the simulations, we find that the energy equipartition assumption is not valid at scales smaller than the driving scale of the ISM turbulence ($\approx 100 \, {\rm pc}$ in spiral galaxies), which can be regarded as the lower limit for the scale beyond which equipartition is valid. We suggest that one must be aware of the dynamical scales in the system before assuming energy equipartition to extract magnetic field information from synchrotron observations. Finally, we present ideas for future observations and simulations to investigate in more detail under which conditions the equipartition assumption is valid or not.

Interpretations of synchrotron observations often assume a tight correlation between magnetic and cosmic ray energy densities. We examine this assumption using both test-particle simulations of cosmic rays and MHD simulations which include cosmic rays as a diffusive fluid. We find no spatial correlation between the cosmic rays and magnetic field energy densities at turbulent scales. Moreover, the cosmic ray number density and magnetic field energy density are statistically independent. Nevertheless, the cosmic ray spatial distribution is highly inhomogeneous, especially at low energies because the particles are trapped between random magnetic mirrors. These results can significantly change the interpretation of synchrotron observations and thus our understanding of the strength and structure of magnetic fields in the Milky Way and nearby spiral galaxies.

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)

Synchrotron radiation from cosmic rays is a key observational probe of the galactic magnetic field. Interpreting synchrotron emission data requires knowledge of the cosmic ray number density, which is often assumed to be in energy equipartition (or otherwise tightly correlated) with the magnetic field energy. However, there is no compelling observational or theoretical reason to expect such tight correlation to hold across all scales. We use test particle simulations, tracing the propagation of charged particles (protons) through a random magnetic field, to study the cosmic ray distribution at scales comparable to the correlation scale of the turbulent flow in the interstellar medium ($\simeq 100\,{\rm pc}$ in spiral galaxies). In these simulations, we find that there is no spatial correlation between the cosmic ray number density and the magnetic field energy density. In fact, their distributions are approximately statistically independent. We find that low-energy cosmic rays can become trapped between magnetic mirrors, whose location depends more on the structure of the field lines than on the field strength.

Zeldovich's stretch-twist fold (STF) dynamo provided a breakthrough in conceptual understanding of fast dynamos, including fluctuation or small scale dynamos. We study the evolution and saturation behaviour of two types of Baker's map dynamos, which have been used to model Zeldovich's STF dynamo process. Using such maps allows one to analyze dynamos at much higher magnetic Reynolds numbers $R_M$ as compared to direct numerical simulations. In the 2-strip map dynamo there is constant constructive folding while the 4-strip map dynamo also allows the possibility of field reversal. Incorporating a diffusive step parameterised by $R_M$, we find that the magnetic field $B(x)$ is amplified only above a critical $R_M=R_{crit} \sim 4$ for both types of dynamos. We explore the saturation of these dynamos in 3 ways; by a renormalized decrease of the effective $R_M$ (Case I) or due to a decrease in the efficiency of field amplification by stretching (Case II), or a combination of both effects (Case III). For Case I, we show that $B(x)$ in the saturated state, for both types of maps, goes back to the marginal eigenfunction, which is obtained for the critical $R_M=R_{crit}$. This is independent of the initial $R_M=R_{M0}$. On the other hand in Case II, for the 2-strip map, we show that $B(x)$ now saturates preserving the structure of the kinematic eigenfunction. Thus the energy is transferred to larger scales in Case I but remains at the smallest resistive scales in Case II. For the 4-strip map, the $B(x)$ oscillates with time, although with a structure similar to the kinematic eigenfunction. Interestingly, the saturated state for Case III shows an intermediate behaviour, with $B(x)$ now similar to the kinematic eigenfunction for an intermediate $R_M=R_{sat}$, with $R_{M0}>R_{sat}>R_{crit}$. These saturation properties are akin to the ones discussed in the context of fluctuation dynamos.

The propagation of cosmic rays in turbulent magnetic fields is a diffusive process driven by the scattering of the charged particles by random magnetic fluctuations. Such fields are usually highly intermittent, consisting of intense magnetic filaments and ribbons surrounded by weaker, unstructured fluctuations. Studies of cosmic ray propagation have largely overlooked intermittency, instead relying on Gaussian random magnetic fields. Using test particle simulations, we investigate cosmic ray diffusivity in intermittent, dynamo-generated magnetic fields. The results are compared with those obtained from non-intermittent magnetic fields having identical power spectra. The presence of magnetic intermittency significantly enhances cosmic ray diffusion over a wide range of particle energies. We demonstrate that the results can be interpreted in terms of a correlated random walk.

We investigate the anomalies in the Earth - Moon system using ancient eclipse data. We identify nine groups of anomalous eclipses between 400 and 1800 AD recorded in parts of India that should have completely missed the subcontinent as per NASA simulations (Espenak and Meeus, 2011). We show that the typical correction to the lunar location required to reconcile the anomalous eclipses is relatively small and consistent with the fluctuations in the length of day that are observed in recent periods. We then investigate the change in Earth's moment of inertia due to differential acceleration of land and water that can account for this discrepancy. We show that 80 percent of these discrepancies occur when the Moon is at declinations greater than 10 deg and closer to its major standstill of 28 deg while it spends 46 percent of the time in this region. We simulate the differential interaction of the Moon's gravity with landmass and water using finite element method to account for landmass and water mass. We show that the results of eclipse error are consistent with the estimate of a small differential acceleration when the Moon is over land at high latitudes. However, we encounter some examples where the results from simulations studies cannot explain the phenomena. Hence we propose that dT corrections have to be coupled with some other mechanism possibly a small vertical oscillation in the Moon's rotational plane with a period of the order of a few hundred years to achieve the required adjustment in the Eclipse maps.