To explore the spatial variations of the regular (mean) magnetic field of the Andromeda galaxy (M31), we use Fourier analysis in azimuthal angle along four rings in the galaxy's plane. The Fourier coefficients give a quantitative measure of strength of the modes, enabling us to compare expectations from mean-field dynamo models of spiral galaxies. Earlier analyses indicated that the axisymmetric magnetic field (azimuthal Fourier mode m = 0) is sufficient to fit the observed polarization angles in a wide range of galactocentric distances. We apply a Bayesian inference approach to new, more sensitive radio continuum data at AA3.59, 6.18, and 11 . 33 cm and the earlier data at A20 . 46 cm to reveal subdominant contributions from the modes m = 1, 2, and 3 along with a dominant axisymmetric mode. Magnetic lines of the axisymmetric mode are close to trailing logarithmic spirals which are significantly more open than the spiral arms detectable in the interstellar dust and neutral hydrogen. The form of the m = 0 mode is consistent with galactic dynamo theory. Both the amplitudes and the pitch angles of the higher azimuthal modes ( m > 1) vary irregularly with r reflecting local variations in the magnetic field structure. The maximum strength of the mean magnetic field of 1 . 8-2 . 7 mu G (for the axisymmetric part of the field) occurs at 10-14 kpc but we find that its strength varies strongly along the azimuth; this variation gives rise to the m = 1 mode. We suggest a procedure of Bayesian inference which is independent of the specific nature of the depolarization and applies when the magneto-ionic layer observable in polarized emission is not symmetric along the line of sight because emission from its far side is completely depolarized.

Synchrotron radiation dominates the continuum emission of star-forming galaxies in the frequency range from a few MHz to about 30 GHz. We model the total synchrotron emission of a large population of evolving star-forming galaxies using the semi-analytic galaxy formation model GALFORM combined with the dynamo simulation code MAGNETIZER. Assuming local energy equipartition between cosmic rays and magnetic fields, we calculate the specific synchrotron luminosity L_ν for each simulated galaxy at various frequencies and find strong positive correlations between L_ν and both the star formation rate (SFR) and characteristic galaxy rotation speed V_ rot for redshifts up to z≃ 3. At low redshifts, the turbulent magnetic field is found to dominate in the synchrotron luminosity, but the contribution of the large-scale magnetic field increases with redshift and becomes important for z≳ 1. The correlation between L_ν and SFR arises from the tight correlation between the disc gas mass M_ gas and SFR, and the correlation between L_ν and V_ rot is additionally a consequence of the stellar mass Tully–Fisher relation for main-sequence galaxies. At low redshifts, the model predictions and observational data compiled for this work show remarkable agreement, but a discrepancy arises at higher redshifts, where modelled SFR values are systematically smaller than those previously inferred from observations. These theoretical models will aid the interpretation of next-generation radio surveys with the Square Kilometre Array and other telescopes.

The common envelope (CE) phase is a key stage in binary star evolution that is still not very well understood. Once engulfed by the giant star, the binary companion may accrete envelope material. For neutron star (NS) companions, such accretion may in principle occur at mass rates several orders of magnitude above the Eddington limit and may result in outflows dominated by powerful bi-polar jets with mass-loss rates similar to the accretion rates. Such jets would impact the morphology of the system and the rate of envelope unbinding, which affect the duration and outcome of the CE event. Employing 3D global hydrodynamic simulations, we study the role of such NS jets in a CE event involving a red giant branch star. The jets eventually drill through and break out of the envelope, producing prominent low-density bi-polar lobes. The jets cause about twice as much envelope mass to be unbound as compared to simulations of the same duration without NS jets. However, the rate of mass unbinding due to the jets decreases towards the ends of the simulations as the jets break out and energetically decouple from the envelope. Moreover, jet activity leads to slightly reduced drag on the binary, decreasing the rate of orbital energy transfer to the envelope. Hence, while such powerful jets can play an important role, negative feedback effects tend to prevent them from dominating envelope unbinding and dictating CE outcomes.

Many spiral galaxies host magnetic fields with energy densities comparable to those of the turbulent and thermal motions of their interstellar gas. However, quantitative comparison between magnetic field properties inferred from observation and those obtained from theoretical modeling has been lacking. In L. Chamandy et al. (Paper I ), we developed a simple, axisymmetric galactic dynamo model that uses various observational data as input. Here, we apply our model to calculate radial profiles of azimuthally and vertically averaged magnetic field strength and pitch angle, gas velocity dispersion and scale height, turbulent correlation time and length, and the sizes of supernova remnants for the galaxies M31, M33, M51, and NGC 6946, using input data collected from the literature. Scaling factors are introduced to account for a lack of precision in both theory and observation. Despite the simplicity of our model, its outputs agree fairly well with galaxy properties inferred from observation. Additionally, we find that most of the parameter values are similar between galaxies. We extend the model to predict the magnetic field pitch angles arising from a combination of mean-field dynamo action and the winding up of the random small-scale field owing to the large-scale radial shear. We find their magnitudes to be much smaller than those of the pitch angles measured in polarized radio and far-infrared emission. This suggests that effects not included in our model, such as effects associated with spiral arms, are needed to explain the pitch angle values.

Common envelope (CE) evolution is largely governed by the drag torque applied on the in-spiralling stellar components by the envelope. Previous work has shown that idealized models of the torque based on a single body moving in rectilinear motion through an unperturbed atmosphere can be highly inaccurate. Progress requires new models for the torque that account for binarity. Towards this end, we perform a new 3D global hydrodynamic CE simulation with the mass of the companion point particle set equal to the mass of the asymptotic giant branch star core particle to maximize symmetry and facilitate interpretation. First, we find that a region around the particles of a scale comparable to their separation contributes essentially all of the torque. Second, the density pattern of the torque-dominating gas and, to an extent, this gas itself, is roughly in corotation with the binary. Third, approximating the spatial distribution of the torquing gas as a uniform-density prolate spheroid whose major axis resides in the orbital plane and lags the line joining the binary components by a constant phase angle reproduces the torque evolution remarkably well, analogous to studies of binary supermassive black holes. Fourth, we compare the torque measured in the simulation with the predictions of a model that assumes two weak point-mass perturbers undergoing circular motion in a uniform background without gas self-gravity and find remarkable agreement with our results if the background density is taken to be equal to a fixed fraction ( approximate to 0.44 ) of the density at the spheroid surface. Overall, this work makes progress towards developing simple time-dependent models of the CE phase, for example by informing the development of drag force prescriptions for 1D spherically symmetric CE simulations, which could be used to explore the parameter space of luminous red novae or in binary population synthesis studies.

We propose that certain white dwarf (WD) planets, such as WD 1856+534 b, may form out of material from a stellar companion that tidally disrupts from common envelope evolution with the WD progenitor star. The disrupted companion shreds into an accretion disc, out of which a gas giant protoplanet forms due to gravitational instability. To explore this scenario, we make use of detailed stellar evolution models consistent with WD 1856+534. The minimum mass companion that produces a gravitationally-unstable disc after tidal disruption is $\sim0.15\,\mathrm{M}_\odot$. In this scenario, WD 1856+534 b might have formed at or close to its present separation, in contrast to other proposed scenarios where it would have migrated in from a much larger separation. Planet formation from tidal disruption is a new channel for producing second-generation planets around WDs.

ABSTRACT We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\, \mathrm{M_\odot }$ red giant branch primary and a $1\, \mathrm{M_\odot }$ companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of ${\sim}2\, {\rm yr}$, after which the envelope would be unbound.

Galactic dynamo models have generally relied on input parameters that are very challenging to constrain. We address this problem by developing a model that uses observable quantities as input: the galaxy rotation curve, the surface densities of the gas, stars and star formation rate, and the gas temperature. The model can be used to estimate parameters of the random and mean components of the magnetic field, as well as the gas scale height, root-mean-square velocity and the correlation length and time of the interstellar turbulence, in terms of the observables. We use our model to derive theoretical scaling relations for the quantities of interest, finding reasonable agreement with empirical scaling relations inferred from observation. We assess the dependence of the results on different assumptions about turbulence driving, finding that agreement with observations is improved by explicitly modeling the expansion and energetics of supernova remnants. The model is flexible enough to include alternative prescriptions for the physical processes involved, and we provide links to two open-source PYTHON programs that implement it.

We explore the redshift evolution of the radio luminosity function (RLF) of star-forming galaxies using GALFORM, a semi-analytic model of galaxy formation (SAMGF) and MAGNETIZER, a dynamo model of the magnetic field evolving in a galaxy through its formation history. Assuming energy equipartition between the magnetic field and cosmic rays, we derive the synchrotron luminosity of each sample galaxy. In a model where the turbulent speed is correlated with the star formation rate, the RLF is in fair agreement with observations in the redshift range $0 \leq z \leq 2$. At larger redshifts, the structure of galaxies, their interstellar matter and turbulence appear to be rather different from those at $z\lesssim2$, so that the turbulence and magnetic field models applicable at low redshifts become inadequate. The strong redshift evolution of the RLF at $0 \leq z \leq 2$ can be attributed to an increased number, at high redshift, of galaxies with large disc volumes and strong magnetic fields. On the other hand, in models where the turbulent speed in galaxies is a constant or assumed to be an explicit function of the redshift, the redshift evolution of the RLF is poorly captured. The evolution of the interstellar turbulence and outflow parameters appear to be major (but not the only) drivers of the RLF changes. We find that both the small- and large-scale magnetic fields contribute significantly to the RLF but the small-scale field dominates at high redshifts. Polarisation observations will therefore be important to distinguish these two components and understand better the evolution of galaxies and their nonthermal constituents.

ABSTRACT It has long been speculated that jet feedback from accretion on to the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model without a jet. Our binary consists of a 2-M⊙ red giant branch primary and a 1- or 0.5-M⊙ main sequence (MS) or white dwarf (WD) secondary companion modelled as a point particle. We run the simulations for 10 orbits (40 d). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of the order of the Eddington rate, with maximum velocity vj of the order of the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on the jet mass-loss rate, launch speed, companion mass, opening angle, and accretion rate. In line with our theoretical estimates, jets are choked around the time of first periastron passage and remain choked thereafter. Subsequent to choking, but not before, jets efficiently transfer energy to bound envelope material. This leads to increases in unbound mass of up to $\sim 10{{\ \rm per\ cent}}$, as compared to the simulations without jets. We also estimate the cumulative effects of jets over a full CE phase, finding that jets launched by MS and WD companions are unlikely to dominate envelope unbinding.

Many stars harbour multi-planet systems. As these stars expand late in their evolutions, the innermost planet may be engulfed, leading to a common envelope (CE) event. Even if this is insufficient to eject the envelope, it may expand the star further, causing additional CE events, with the last one unbinding what remains of the envelope. This multi-planet CE scenario may have broad implications for stellar and planetary evolution across a range of systems. We develop a simplified version and show that it may be able to explain the recently observed planet WD 1856 b.

Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primaries. But CE evolution with AGB progenitors is desirable to simulate because such events are the likely progenitors of most bi-polar planetary nebulae (PNe), and prominent observational testing grounds for CE physics. Here we present a high-resolution global simulation of CE evolution involving an AGB primary and 1-M-circle dot secondary, evolved for 20 orbital revolutions. During the last 16 of these orbits, the envelope unbinds at an almost constant rate of about 0.1-0.2M(circle dot) yr(-1). If this rate were maintained, the envelope would be unbound in less than 10 yr. The dominant source of this unbinding is consistent with inspiral; we assess the influence of the ambient medium to be subdominant. We compare this run with a previous run that used an RGB phase primary evolved from the same 2-M-circle dot main-sequence star to assess the influence of the evolutionary state of the primary. When scaled appropriately, the two runs are quite similar, but with some important differences.

Galactic dynamo models take as input certain parameters of the interstellar turbulence, most essentially the correlation time τ, root-mean-square turbulent speed u, and correlation scale l. However, these quantities are difficult, or, in the case of τ, impossible, to directly observe, and theorists have mostly relied on order of magnitude estimates. Here we present an analytic model to derive these quantities in terms of a small set of more accessible parameters. in our model, turbulence is assumed to be driven concurrently by isolated supernovae (SNe) and superbubbles (SBs), but clustering of SNe to form SBs can be turned off if desired, which reduces the number of model parameters by about half. In general, we find that isolated SNe and SBs can inject comparable amounts of turbulent energy into the interstellar medium, but SBs do so less efficiently. this results in rather low overall conversion rates of SN energy into turbulent energy of ∼1–3%. the results obtained for l, u and τ for model parameter values representative of the Solar neighbourhood are consistent with those determined from direct numerical simulations. Our analytic model can be combined with existing dynamo models to predict more directly the magnetic field properties for nearby galaxies or for statistical populations of galaxies in cosmological models.

Constraining dynamo theories of magnetic field origin by observation is indispensable but challenging, in part because the basic quantities measured by observers and predicted by modelers are different. We clarify these differences and sketch out ways to bridge the divide. Based on archival and previously unpublished data, we then compile various important properties of galactic magnetic fields for nearby spiral galaxies. We consistently compute strengths of total, ordered, and regular fields, pitch angles of ordered and regular fields, and we summarize the present knowledge on azimuthal modes, field parities, and the properties of non-axisymmetric spiral features called magnetic arms. We review related aspects of dynamo theory, with a focus on mean-field models and their predictions for large-scale magnetic fields in galactic discs and halos. Furthermore, we measure the velocity dispersion of H i gas in arm and inter-arm regions in three galaxies, M 51, M 74, and NGC 6946, since spiral modulation of the root-mean-square turbulent speed has been proposed as a driver of non-axisymmetry in large-scale dynamos. We find no evidence for such a modulation and place upper limits on its strength, helping to narrow down the list of mechanisms to explain magnetic arms. Successes and remaining challenges of dynamo models with respect to explaining observations are briefly summarized, and possible strategies are suggested. With new instruments like the Square Kilometre Array (SKA), large data sets of magnetic and non-magnetic properties from thousands of galaxies will become available, to be compared with theory.

The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes.

Many stars harbour multi-planet systems. As these stars expand late in their evolution, the innermost planet may be engulfed, leading to a common envelope (CE) event. Even if this CE interaction is insufficient to eject the envelope, it may cause the star to expand further, leading to additional CE events, with the last one unbinding what remains of the envelope. This multi-planet CE scenario may have broad implications for stellar and planetary evolution across a range of systems. We develop a simplified version and show that it may be able to explain the recently observed planet WD 1856 b.

During the common envelope binary interaction, the expanding layers of the gaseous common envelope recombine and the resulting recombination energy has been suggested as a contributing factor to the ejection of the envelope. In this paper we perform a comparative study between simulations with and without the inclusion of recombination energy. We use two distinct setups, comprising 0.88-M$_{\odot}$ and 1.8-M$_{\odot}$ giants, that have been studied before and can serve as benchmarks. In so doing we conclude that (i) the final orbital separation is not affected by the choice of equation of state. In other words, simulations that unbind but a small fraction of the envelope result in similar final separations to those that, thanks to recombination energy, unbind a far larger fractions. (ii) The adoption of a tabulated equation of state results in a much greater fraction of unbound envelope and we demonstrate the cause of this to be the release of recombination energy. (iii) The fraction of hydrogen recombination energy that is allowed to do work should be about half of that which our adiabatic simulations use. (iv) However, for the heavier star simulation we conclude that it is helium and not hydrogen recombination energy that unbinds the gas and we determine that all helium recombination energy is thermalised in the envelope and does work. (v)~The outer regions of the expanding common envelope are likely to see the formation of dust. This dust would promote additional unbinding and shaping of the ejected envelope into axisymmetric morphologies.

We compute the forces, torque, and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably approximated by conventional analytic theory and peaks at values proportional to the companion mass. Good agreement between global and local 3D `wind tunnel' simulations, including similar net drag force and flow pattern, is obtained for comparable regions of parameter space. However, subsequent to the first periastron passage, the drag force is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and depends only weakly on companion mass. The discrepancy is exacerbated for larger companion mass and when the inter-particle separation reduces to the Bondi-Hoyle-Lyttleton accretion radius, creating a turbulent thermalized region. Greater flow symmetry during this phase leads to near balance of opposing gravitational forces in front of and behind the companion, hence a small net drag. The reduced drag force at late times helps explain why companion-core separations necessary for envelope ejection are not reached by the end of limited duration CE simulations.

Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system, suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc. Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.