The magnetic fields of spiral galaxies are so strong that they cannot be primordial. Their typical values are over one billion times higher than any value predicted for the early Universe. Explaining this immense growth and incorporating it in galaxy evolution theories is one of the long-standing challenges in astrophysics. So far, the most successful theory for the sustained growth of the galactic magnetic field is the alpha-omega dynamo. This theory predicts a characteristic dipolar or quadrupolar morphology for the galactic magnetic field, which has been observed in external galaxies. However, so far, there has been no direct demonstration of a mean-field dynamo operating in direct, multi-physics simulations of spiral galaxies. We do so in this work. We employ numerical models of isolated, star-forming spiral galaxies that include a magnetized gaseous disk, a dark matter halo, stars, and stellar feedback. Naturally, the resulting magnetic field has a complex morphology that includes a strong random component. Using a smoothing of the magnetic field on small scales, we are able to separate the mean from the turbulent component and analyze them individually. We find that a mean-field dynamo naturally occurs as a result of the dynamical evolution of the galaxy and amplifies the magnetic field by an order of magnitude over half a Gyr. Despite the highly dynamical nature of these models, the morphology of the mean component of the field is identical to analytical predictions. This result underlines the importance of the mean-field dynamo in galactic evolution. Moreover, by demonstrating the natural growth of the magnetic field in a complex galactic environment, it brings us a step closer to understanding the cosmic origin of magnetic fields.

We investigate the effect of magnetic helicity on the stability of buoyant magnetic cavities as found in the intergalactic medium. In these cavities, we insert helical magnetic fields and test whether or not helicity can increase their stability to shredding through the Kelvin-Helmholtz instability and, with that, their lifetime. This is compared to the case of an external vertical magnetic field which is known to reduce the growth rate of the Kelvin-Helmholtz instability. By comparing a low-helicity configuration with a high helicity one with the same magnetic energy we find that an internal helical magnetic field stabilizes the cavity. This effect increases as we increase the helicity content. Stabilizing the cavity with an external magnetic field requires instead a significantly stronger field at higher magnetic energy. We conclude that the presence of helical magnetic fields is a viable mechanism to explain the stability of intergalactic cavities on time scales longer than 100 Myr.

Proxima Centauri is known to host an earth-like planet in its habitable zone; very recently a second candidate planet was proposed based on radial velocities. At quadrature, the expected projected separation of this new candidate is larger than 1 arcsec, making it a potentially interesting target for direct imaging. While difficult, identification of the optical counterpart of this planet would allow detailed characterization of the closest planetary system. We searched for a counterpart in SPHERE images acquired during four years through the SHINE survey. In order to account for the large orbital motion of the planet, we used a method that assumes the circular orbit obtained from radial velocities and exploits the sequence of observations acquired close to quadrature in the orbit. We checked this with a more general approach that considers keplerian motion, K-stacker. We did not obtain a clear detection. The best candidate has S/N=6.1 in the combined image. A statistical test suggests that the probability that this detection is due to random fluctuation of noise is < 1% but this result depends on the assumption that distribution of noise is uniform over the image. The position of this candidate and the orientation of its orbital plane fit well with observations in the ALMA 12m array image. However, the astrometric signal expected from the orbit of the candidate we detected is 3-sigma away from the astrometric motion of Proxima as measured from early Gaia data. This, together with the unexpectedly high flux associated with our direct imaging detection, means we cannot confirm that our candidate is indeed Proxima c. On the other hand, if confirmed, this would be the first observation in imaging of a planet discovered from radial velocities and the second one (after Fomalhaut b) of reflecting circumplanetary material. Further confirmation observations should be done as soon as possible.

Proxima c, a candidate second planet orbiting Proxima Centauri, was detected with the radial velocity method. The announced long orbital period (5.21$^{+0.26}_{-0.22}$ years), and small semi-amplitude of the induced Doppler signal (1.2$\pm$0.4 m/s), make this detection challenging and the target worthy of a follow-up in the next years. We intend to evaluate the impact of future data on the statistical significance of the detection through realistic simulated radial velocities to be added to the published dataset, spanning up to one orbital period of Proxima c in the time range 2019-2023. We find that the detection significance of Proxima c increases depending not only on the amount of data collected, but also on the number of instruments used, and especially on the timespan covered by the observational campaign. However, on average we do not get strong statistical evidence and we predict that, in the best-case scenario, in the next 5 years the detection of Proxima c can become significant at 4$\sigma$ level. If instead Proxima c does not exist, the detected signal may lower its significance down to 2$\sigma$.

Unveiling the evolution of toroidal field instability, known as Tayler instability, is essential to understand the strength and topology of the magnetic fields observed in early-type stars, in the core of the red giants, or in any stellar radiative zone. We want to study the non-linear evolution of the instability of a toroidal field stored in a stably stratified layer, in spherical symmetry and in the absence of rotation. In particular, we intend to quantify the suppression of the instability as a function of the Brunt-V\"ais\"ala ($\omega_{\rm BV}$) and the Alfv\'en ($\omega_{\rm A}$) frequencies. We use the MHD equations as implemented in the anelastic approximation in the EULAG-MHD code and perform a large series of numerical simulations of the instability exploring the parameter space for the $\omega_{\rm BV}$ and $\omega_{\rm A}$. We show that beyond a critical value gravity strongly suppress the instability, in agreement with the linear analysis. The intensity of the initial field also plays an important role: weaker fields show much slower growth rates. Moreover, in the case of very low gravity, the fastest growing modes have a large characteristic radial scale, at variance with the case of strong gravity, where the instability is characterized by horizontal displacements. Our results illustrate that the anelastic approximation can efficiently describe the evolution of toroidal field instability in stellar interiors. The suppression of the instability as a consequence of increasing values of $\omega_{\rm BV}$ might play a role to explain the magnetic desert in Ap/Bp stars since weak fields are only marginally unstable in the case of strong gravity.

We present further evidence of the presence of acoustic oscillations on the slowly-rotating, over-active G8 sub-giant EK Eri. This star was observed with the 1-m Hertzsprung SONG telescope, at the Observatorio del Teide for two different runs of 8 and 13 nights, respectively, and separated by about a year. We determined a significant excess of power around $\nu_\mathrm{max} = 253 \pm 3\,\mu$Hz in the first observing run and we were able to determine the large separation $\Delta\nu = 16.43 \pm 0.22\,\mu$Hz. No significant excess of power was instead detected in a subsequent SONG observing season, as also supported by our analysis of the simultaneous TESS photometric observations. We propose a new amplitude-luminosity relation in order to account for the missing power in the power spectrum. Based on the evolutionary stage of this object we argue that standard $\alpha^2\Omega$ dynamo cannot be excluded as a possible origin for the observed magnetic field.

The full-sky Planck polarization data at 850um revealed unexpected properties of the E and B mode power spectra of dust emission in the interstellar medium (ISM). The positive cross-correlation between the total dust intensity, T, with the B modes has raised new questions about the physical mechanisms that affect dust polarization, such as the Galactic magnetic-field structure. This is key both to better understanding ISM dynamics and to accurately describing Galactic foregrounds to the polarization of the Cosmic Microwave Background (CMB). In this theoretical paper we investigate the possibility that the observed cross-correlations in the dust polarization power spectra, and specifically between T and B, can be related to a parity-odd quantity in the ISM such as the magnetic helicity. We produce synthetic dust polarization data, derived from 3D analytical toy models of density structures and helical magnetic fields, to compare with the E and B modes of observations. Focusing on the observed T-B correlation, we propose a new line of interpretation of the Planck observations based on a large-scale helical component of the Galactic magnetic field in the solar neighborhood. Our analysis shows that: I) the sign of magnetic helicity does not affect E and B modes for isotropic magnetic-field configurations; II) helical magnetic fields threading interstellar filaments cannot reproduce the Planck results; III) a weak helical left-handed magnetic field structure in the solar neighborhood may explain the T-B correlation seen in the Planck data. This work suggests a new perspective for the interpretation of the dust polarization power spectra, which strongly supports the imprint of a large-scale structure of the Galactic magnetic field in the solar neighborhood.

Stellar dynamos are affected by boundary conditions imposed by stellar coronae. Under some approximations it is possible to find analytical solutions. Interior dynamo models often consider a current-free coronae without taking into account the constraints imposed by the presence of currents in the corona. We aim to analytically evaluate the effect of coronal currents and of an outer boundary condition on the efficiency of an $\alpha^2$ dynamo. We intend to estimate the change in geometry and dinamo excitation numbers with respect to the current-free case. We analytically solve the turbulent dynamo induction equation for a homogeneous, non-mirror symmetric turbulence, in a spherical domain surrounded by a linear force-free corona with magnetic field B satisfying curl B = \beta B. The main result is that the dynamo number is a decreasing function of \beta. Moreover, if the current is parallel to the field (\beta > 0) the dynamo number is smaller than in the force-free case. On the contrary (\beta < 0) the dynamo number is greater than the in force-free case. We conclude that the presence of currents in the corona needs to be taken into account because it affects the condition for excitation of a dynamo.

The detection and characterisation of Earth-like planets with Doppler signals of the order of 1 m/s currently represent one of the greatest challenge for extrasolar-planet hunters. As results for such findings are often controversial, it is desirable to provide independent confirmations of the discoveries. Testing different models for the suppression of non-Keplerian stellar signals usually plaguing radial velocity data is essential to ensuring findings are robust and reproducible. Using an alternative treatment of the stellar noise to that discussed in the discovery paper, we re-analyse the radial velocity dataset that led to the detection of a candidate terrestrial planet orbiting the star Proxima Centauri. We aim to confirm the existence of this outstanding planet, and test the existence of a second planetary signal. Our technique jointly modelled Keplerian signals and residual correlated signals in radial velocities using Gaussian processes. We analysed only radial velocity measurements without including other ancillary data in the fitting procedure. In a second step, we have compared our outputs with results coming from photometry, to provide a consistent physical interpretation. Our analysis was performed in a Bayesian framework to quantify the robustness of our findings. We show that the correlated noise can be successfully modelled as a Gaussian process regression, and contains a periodic term modulated on the stellar rotation period and characterised by an evolutionary timescale of the order of one year. Both findings appear to be robust when compared with results obtained from archival photometry, thus providing a reliable description of the noise properties. We confirm the existence of a coherent signal described by a Keplerian orbit equation that can be attributed to the planet Proxima b [...]

Owing to technological advances, the number of exoplanets discovered has risen dramatically in the last few years. However, when trying to observe Earth analogs, it is often difficult to test the veracity of detection. We have developed a new approach to the analysis of exoplanetary spectral observations based on temporal multifractality, which identifies time scales that characterize planetary orbital motion around the host star, and those that arise from stellar features such as spots. Without fitting stellar models to spectral data, we show how the planetary signal can be robustly detected from noisy data using noise amplitude as a source of information. For observation of transiting planets, combining this method with simple geometry allows us to relate the time scales obtained to primary and secondary eclipse of the exoplanets. Making use of data obtained with ground-based and space-based observations we have tested our approach on HD 189733b. Moreover, we have investigated the use of this technique in measuring planetary orbital motion via Doppler shift detection. Finally, we have analyzed synthetic spectra obtained using the SOAP 2.0 tool, which simulates a stellar spectrum and the influence of the presence of a planet or a spot on that spectrum over one orbital period. We have demonstrated that, so long as the signal-to-noise-ratio $\ge$ 75, our approach reconstructs the planetary orbital period, as well as the rotation period of a spot on the stellar surface.

Many astrophysical bodies harbor magnetic fields that are thought to be sustained by a dynamo process. However, it has been argued that the production of large-scale magnetic fields by mean-field dynamo action is strongly suppressed at large magnetic Reynolds numbers owing to the conservation of magnetic helicity. This phenomenon is known as {\it catastrophic quenching}. Advection of magnetic fields by stellar and galactic winds toward the outer boundaries and away from the dynamo is expected to alleviate such quenching. Here we explore the relative roles played by advective and turbulent--diffusive fluxes of magnetic helicity in the dynamo. In particular, we study how the dynamo is affected by advection. We do this by performing direct numerical simulations of a turbulent dynamo of $\alpha^2$ type driven by forced turbulence in a Cartesian domain in the presence of a flow away from the equator where helicity changes sign. Our results indicate that in the presence of advection, the dynamo, otherwise stationary, becomes oscillatory. We confirm an earlier result for turbulent--diffusive magnetic helicity fluxes that for small magnetic Reynolds numbers ($\Rm\lesssim 100...200$, based on the wavenumber of the energy-carrying eddies) the magnetic helicity flux scales less strongly with magnetic Reynolds number ($\Rm^{-1/2}$) than the term describing magnetic helicity destruction by resistivity ($\Rm^{-1}$). Our new results now suggest that for larger $\Rm$ the former becomes approximately independent of $\Rm$, while the latter falls off more slowly. We show for the first time that both for weak and stronger winds, the magnetic helicity flux term becomes comparable to the resistive term for $\Rm\gtrsim 1000$, which is necessary for alleviating catastrophic quenching.

We study spontaneous breakdown of chiral symmetry during the nonlinear evolution of the Tayler instability. We start with an initial steady state of zero helicity. Within linearized perturbation calculations, helical perturbations of this initial state have the same growth rate for either sign of helicity. Direct numerical simulations (DNS) of the fully nonlinear equations, however, show that an infinitesimal excess of one sign of helicity in the initial perturbation gives rise to a saturated helical state. We further show that this symmetry breaking can be described by weakly nonlinear finite--amplitude equations with undetermined coefficients which can be deduced solely from symmetry consideration. By fitting solutions of the amplitude equations to data from DNS we further determine the coefficients of the amplitude equations.

We study the role of viscosity and the effects of a magnetic field on a rotating, self-gravitating fluid, using Newtonian theory and adopting the ideal magnetohydrodynamic approximation. Our results confirm that viscosity can generate vorticity in inhomogeneous environments, while the magnetic tension can produce vorticity even in the absence of fluid pressure and density gradients. Linearizing our equations around an Einstein-de Sitter cosmology, we find that viscosity adds to the diluting effect of the universal expansion. Typically, however, the dissipative viscous effects are confined to relatively small scales. We also identify the characteristic length bellow which the viscous dissipation is strong and beyond which viscosity is essentially negligible. In contrast, magnetism seems to favor cosmic rotation. The magnetic presence is found to slow down the standard decay-rate of linear vortices, thus leading to universes with more residual rotation than generally anticipated.

The chiral symmetry breaking properties of the Tayler instability are discussed. Effective amplitude equations are determined in one case. This model has three free parameters that are determined numerically. Comparison with chiral symmetry breaking in biochemistry is made.

Certain aspects of the mean-field theory of turbulent passive scalar transport and of mean-field electrodynamics are considered with particular emphasis on aspects of compressible fluids. It is demonstrated that the total mean-field diffusivity for passive scalar transport in a compressible flow may well be smaller than the molecular diffusivity. This is in full analogy to an old finding regarding the magnetic mean-field diffusivity in an electrically conducting turbulently moving compressible fluid. These phenomena occur if the irrotational part of the motion dominates the vortical part, the P\`eclet or magnetic Reynolds number is not too large, and, in addition, the variation of the flow pattern is slow. For both the passive scalar and the magnetic cases several further analytical results on mean-field diffusivities and related quantities found within the second-order correlation approximation are presented, as well as numerical results obtained by the test-field method, which applies independently of this approximation. Particular attention is paid to non-local and non-instantaneous connections between the turbulence-caused terms and the mean fields. Two examples of irrotational flows, in which interesting phenomena in the above sense occur, are investigated in detail. In particular, it is demonstrated that the decay of a mean scalar in a compressible fluid under the influence of these flows can be much slower than without any flow, and can be strongly influenced by the so-called memory effect, that is, the fact that the relevant mean-field coefficients depend on the decay rates themselves.

Two setups with interlocked magnetic flux tubes are used to study the evolution of magnetic energy and helicity on magnetohydrodynamical (MHD) systems like plasmas. In one setup the initial helicity is zero while in the other it is finite. To see if it is the actual linking or merely the helicity content that influences the dynamics of the system we also consider a setup with unlinked field lines as well as a field configuration in the shape of a trefoil knot. For helical systems the decay of magnetic energy is slowed down by the helicity which decays slowly. It turns out that it is the helicity content, rather than the actual linking, that is significant for the dynamics.

Observations have shown that the Sun's magnetic field has helical structures. The helicity content in magnetic field configurations is a crucial constraint on the dynamical evolution of the system. Since helicity is connected with the number of links we investigate configurations with interlocked magnetic flux rings and one with unlinked rings. It turns out that it is not the linking of the tubes which affects the magnetic field decay, but the content of magnetic helicity.

In the interstellar medium the turbulence is believed to be forced mostly through supernova explosions. In a first approximation these flows can be written as a gradient of a potential being thus devoid of vorticity. There are several mechanisms that could lead to vorticity generation, like viscosity and baroclinic terms, rotation, shear and magnetic fields, but it is not clear how effective they are, neither is it clear whether the vorticity is essential in determining the turbulent diffusion acting in the ISM. Here we present a study of the role of rotation, shear and baroclinicity in the generation of vorticity in the ISM.

In the absence of rotation and shear, and under the assumption of constant temperature or specific entropy, purely potential forcing by localized expansion waves is known to produce irrotational flows that have no vorticity. Here we study the production of vorticity under idealized conditions when there is rotation, shear, or baroclinicity, to address the problem of vorticity generation in the interstellar medium in a systematic fashion. We use three-dimensional periodic box numerical simulations to investigate the various effects in isolation. We find that for slow rotation, vorticity production in an isothermal gas is small in the sense that the ratio of the root-mean-square values of vorticity and velocity is small compared with the wavenumber of the energy-carrying motions. For Coriolis numbers above a certain level, vorticity production saturates at a value where the aforementioned ratio becomes comparable with the wavenumber of the energy-carrying motions. Shear also raises the vorticity production, but no saturation is found. When the assumption of isothermality is dropped, there is significant vorticity production by the baroclinic term once the turbulence becomes supersonic. In galaxies, shear and rotation are estimated to be insufficient to produce significant amounts of vorticity, leaving therefore only the baroclinic term as the most favorable candidate. We also demonstrate vorticity production visually as a result of colliding shock fronts.

A spherical hydrodynamical expansion flow can be described as the gradient of a potential. In that case no vorticity should be produced, but several additional mechanisms can drive its production. Here we analyze the effects of baroclinicity, rotation and shear in the case of a viscous fluid. Those flows resemble what happens in the interstellar medium. In fact in this astrophysical environment supernovae explosion are the dominant flows and, in a first approximation, they can be seen as spherical. One of the main difference is that in our numerical study we examine only weakly supersonic flows, while supernovae explosions are strongly supersonic.

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

The resistive decay of chains of three interlocked magnetic flux rings is considered. Depending on the relative orientation of the magnetic field in the three rings, the late-time decay can be either fast or slow. Thus, the qualitative degree of tangledness is less important than the actual value of the linking number or, equivalently, the net magnetic helicity. Our results do not suggest that invariants of higher order than that of the magnetic helicity need to be considered to characterize the decay of the field.

Estimates for the nonlinear alpha effect in helical turbulence with an applied magnetic field are presented using two different approaches: the imposed-field method where the electromotive force owing to the applied field is used, and the test-field method where separate evolution equations are solved for a set of different test fields. Both approaches agree for stronger fields, but there are apparent discrepancies for weaker fields that can be explained by the influence of dynamo-generated magnetic fields on the scale of the domain that are referred to as meso-scale magnetic fields. Examples are discussed where these meso-scale fields can lead to both drastically overestimated and underestimated values of alpha compared with the kinematic case. It is demonstrated that the kinematic value can be recovered by resetting the fluctuating magnetic field to zero in regular time intervals. It is concluded that this is the preferred technique both for the imposed-field and the test-field methods.

Using the test-field method for nearly irrotational turbulence driven by spherical expansion waves it is shown that the turbulent magnetic diffusivity increases with magnetic Reynolds numbers. Its value levels off at several times the rms velocity of the turbulence multiplied by the typical radius of the expansion waves. This result is discussed in the context of the galactic mean-field dynamo.