Dynamo saturation in direct simulations of the multi‐phase turbulent interstellar medium

The ordered magnetic field observed via polarised synchrotron emission in nearby disc galaxies can be explained by a mean‐field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical‐flux initial conditions are potentially able to influence this dynamo via the occurrence of the magnetorotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean‐field dynamo. This is motivated by the observation that different saturation behaviour was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three‐dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing‐periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean‐field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched α) dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxy's star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub‐grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Minimal Spanning Tree algorithm for source detection in γ-ray images

We developed a source detection algorithm based on the Minimal Spanning Tree (MST), that is a graph-theoretical method useful for finding clusters in a given set of points. This algorithm is applied to γ-ray bi-dimensional images where the points correspond to the arrival direction of photons, and the possible sources are associated with the regions where they clusterize. Some filters to select these clusters and to reduce the spurious detections are introduced. An empirical study of the statistical properties of MST on random fields is carried out in order to derive some criteria to estimate the best filter values. We also introduce two parameters useful to verify the goodness of candidate sources. To show how the MST algorithm works in practice, we present an application to an EGRET observation of the Virgo field, at high Galactic latitude and with a low and rather uniform background, in which several sources are detected.     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

The dependence of galaxy formation on cosmological parameters: can we distinguish between the WMAP1 and WMAP3 parameter sets?

We combine N -body simulations of structure growth with physical modelling of galaxy evolution to investigate whether the shift in cosmological parameters between the first- and third-year results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) affects predictions for the galaxy population. Structure formation is significantly delayed in the WMAP3 cosmology, because the initial matter fluctuation amplitude is lower on the relevant scales. The decrease in dark matter clustering strength is, however, almost entirely offset by an increase in halo bias, so predictions for galaxy clustering are barely altered. In both cosmologies, several combinations of physical parameters can reproduce observed, low-redshift galaxy properties; the star formation, supernova feedback and active galactic nucleus feedback efficiencies can be played off against each other to give similar results. Models which fit observed luminosity functions predict projected two-point correlation functions which scatter by about 10–20 per cent on large scale and by larger factors on small scale, depending both on cosmology and on details of galaxy formation. Measurements of the pairwise velocity distribution prefer the WMAP1 cosmology, but careful treatment of the systematics is needed. Given present modelling uncertainties, it is not easy to distinguish between the WMAP1 and WMAP3 cosmologies on the basis of low-redshift galaxy properties. Model predictions diverge more dramatically at high redshift. Better observational data at   z > 2  will better constrain galaxy formation and perhaps also cosmological parameters.