-
Understanding the radio luminosity function of star-forming galaxies and its cosmological evolution
Authors:
Charles Jose,
Luke Chamandy,
Anvar Shukurov,
Kandaswamy Subramanian,
Luiz Felippe S. Rodrigues,
Carlton M. Baugh
Abstract:
We explore the redshift evolution of the radio luminosity function (RLF) of star-forming galaxies using GALFORM, a semi-analytic model of galaxy formation and a dynamo model of the magnetic field evolving in a galaxy. 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 correlat…
▽ More
We explore the redshift evolution of the radio luminosity function (RLF) of star-forming galaxies using GALFORM, a semi-analytic model of galaxy formation and a dynamo model of the magnetic field evolving in a galaxy. 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 is a constant or an explicit function of $z$, the observed 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 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.
△ Less
Submitted 11 June, 2024; v1 submitted 23 February, 2024;
originally announced February 2024.
-
Galactic magnetic fields I. Theoretical model and scaling relations
Authors:
Luke Chamandy,
Rion Glenn Nazareth,
Gayathri Santhosh
Abstract:
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 magn…
▽ More
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.
△ Less
Submitted 3 April, 2024; v1 submitted 27 November, 2023;
originally announced November 2023.
-
How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution
Authors:
Luke Chamandy,
Jonathan Carroll-Nellenback,
Eric G. Blackman,
Adam Frank,
Yisheng Tu,
Baowei Liu,
Yangyuxin Zou,
Jason Nordhaus
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\,M_\odot$ red giant branch primary and a $1\,M_\odot$ companion modeled as a particle. We compare a run employing a tabulated equation of state (EOS) tha…
▽ More
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\,M_\odot$ red giant branch primary and a $1\,M_\odot$ companion modeled 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, $\sim15$ 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 $\sim2\,\mathrm{yr}$, after which the envelope would be unbound.
△ Less
Submitted 4 January, 2024; v1 submitted 28 April, 2023;
originally announced April 2023.
-
Jets from main sequence and white dwarf companions during common envelope evolution
Authors:
Yangyuxin Zou,
Luke Chamandy,
Jonathan Carroll-Nellenback,
Eric G. Blackman,
Adam Frank
Abstract:
It has long been speculated that jet feedback from accretion onto the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process, but this conjecture has heretofore remained largely untested. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model w…
▽ More
It has long been speculated that jet feedback from accretion onto the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process, but this conjecture has heretofore remained largely untested. 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 $2M_\odot$ red giant branch primary and a $1M_\odot$ or $0.5M_\odot$ main sequence or white dwarf secondary companion modeled as a point particle. We run the simulations for 10 orbits (40 days). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of order the Eddington rate, with maximum velocity $v_\mathrm{j}$ of order the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane, and supplies kinetic energy at the rate $\sim\dot{M}_\mathrm{j} v_\mathrm{j}^2/40$. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on jet mass-loss rate, launch speed, companion mass, opening angle, and whether or not subgrid accretion is turned on. In line with our theoretical estimates, we find that in all cases the jet becomes choked around the time of first periastron passage. We also find that jets lead to increases in unbound mass of up to $\sim10\%$, as compared to simulations which do not include a jet.
△ Less
Submitted 7 September, 2022; v1 submitted 11 February, 2022;
originally announced February 2022.
-
The Formation of Discs in the Interior of AGB Stars from the Tidal Disruption of Planets and Brown Dwarfs
Authors:
Gabriel Guidarelli,
Jason Nordhaus,
Jonathan Carroll-Nellenback,
Luke Chamandy,
Eric G. Blackman,
Adam Frank
Abstract:
A significant fraction of isolated white dwarfs host magnetic fields in excess of a MegaGauss. Observations suggest that these fields originate in interacting binary systems where the companion is destroyed thus leaving a singular, highly-magnetized white dwarf. In post-main-sequence evolution, radial expansion of the parent star may cause orbiting companions to become engulfed. During the common…
▽ More
A significant fraction of isolated white dwarfs host magnetic fields in excess of a MegaGauss. Observations suggest that these fields originate in interacting binary systems where the companion is destroyed thus leaving a singular, highly-magnetized white dwarf. In post-main-sequence evolution, radial expansion of the parent star may cause orbiting companions to become engulfed. During the common envelope phase, as the orbital separation rapidly decreases, low-mass companions will tidally disrupt as they approach the giant's core. We hydrodynamically simulate the tidal disruption of planets and brown dwarfs, and the subsequent accretion disc formation, in the interior of an asymptotic giant branch star. These dynamically formed discs are commensurate with previous estimates, suggesting strong magnetic fields may originate from these tidal disruption events.
△ Less
Submitted 30 April, 2021;
originally announced May 2021.
-
Successive common envelope events from multiple planets
Authors:
Luke Chamandy,
Eric G. Blackman,
Jason Nordhaus,
Emily Wilson
Abstract:
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 fo…
▽ More
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.
△ Less
Submitted 18 February, 2021; v1 submitted 22 November, 2020;
originally announced November 2020.
-
Parameters of the Supernova-Driven Interstellar Turbulence
Authors:
Luke Chamandy,
Anvar Shukurov
Abstract:
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 der…
▽ More
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 $\sim1$-$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.
△ Less
Submitted 15 August, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Common Envelope Evolution on the Asymptotic Giant Branch: Unbinding within a Decade?
Authors:
Luke Chamandy,
Eric G. Blackman,
Adam Frank,
Jonathan Carroll-Nellenback,
Yisheng Tu
Abstract:
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) primar…
▽ More
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\,\mathrm{M}_\odot$ 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.2\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$. If this rate were maintained, the envelope would be unbound in less than $10\,\mathrm{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\,\mathrm{M}_\odot$ 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.
△ Less
Submitted 27 May, 2020; v1 submitted 14 April, 2020;
originally announced April 2020.
-
Synthesizing Observations and Theory to Understand Galactic Magnetic Fields: Progress and Challenges
Authors:
Rainer Beck,
Luke Chamandy,
Ed Elson,
Eric G. Blackman
Abstract:
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 fo…
▽ More
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. Further, we measure the velocity dispersion of HI gas in arm and inter-arm regions in three galaxies, M 51, M 74, and NGC 6496, 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.
△ Less
Submitted 21 December, 2019; v1 submitted 18 December, 2019;
originally announced December 2019.
-
Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution
Authors:
Yangyuxin Zou,
Adam Frank,
Zhuo Chen,
Thomas Reichardt,
Orsola De Marco,
Eric G. Blackman,
Jason Nordhaus,
Bruce Balick,
Jonathan Carroll-Nellenback,
Luke Chamandy,
Baowei Liu
Abstract:
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 ev…
▽ More
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.
△ Less
Submitted 3 June, 2020; v1 submitted 3 December, 2019;
originally announced December 2019.
-
How Drag Force Evolves in Global Common Envelope Simulations
Authors:
Luke Chamandy,
Eric G. Blackman,
Adam Frank,
Jonathan Carroll-Nellenback,
Yangyuxin Zou,
Yisheng Tu
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 October, 2019; v1 submitted 16 August, 2019;
originally announced August 2019.
-
Hydrodynamic Simulations of Disrupted Planetary Accretion Discs Inside the Core of an AGB Star
Authors:
Gabriel Guidarelli,
Jason Nordhaus,
Luke Chamandy,
Zhuo Chen,
Eric G. Blackman,
Adam Frank,
Jonathan Carroll-Nellenback,
Baowei Liu
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 31 July, 2019;
originally announced August 2019.
-
Energy Budget and Core-Envelope Motion in Common Envelope Evolution
Authors:
Luke Chamandy,
Yisheng Tu,
Eric G. Blackman,
Jonathan Carroll-Nellenback,
Adam Frank,
Baowei Liu,
Jason Nordhaus
Abstract:
We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total e…
▽ More
We analyze a 3D hydrodynamic simulation of common envelope evolution to understand how energy is transferred between various forms and whether theory and simulation are mutually consistent given the setup. Virtually all of the envelope unbinding in the simulation occurs before the end of the rapid plunge-in phase, here defined to coincide with the first periastron passage. In contrast, the total envelope energy is nearly constant during this time because positive energy transferred to the gas from the core particles is counterbalanced by the negative binding energy from the closer proximity of the inner layers to the plunged-in secondary. During the subsequent slow spiral-in phase, energy continues to transfer to the envelope from the red giant core and secondary core particles. We also propose that relative motion between the centre of mass of the envelope and the centre of mass of the particles could account for the offsets of planetary nebula central stars from the nebula's geometric centre.
△ Less
Submitted 25 March, 2019; v1 submitted 28 December, 2018;
originally announced December 2018.
-
New Insights on Galactic Dynamos
Authors:
Luke Chamandy,
Anvar Shukurov,
A. Russ Taylor
Abstract:
We argue that pitch angles of the azimuthally averaged large-scale or mean magnetic fields in nearby spiral galaxies inferred from observations can tentatively be explained with simple galactic dynamo models. Agreement is not perfect, but is reasonable considering the uncertainty in dynamo parameters.
We argue that pitch angles of the azimuthally averaged large-scale or mean magnetic fields in nearby spiral galaxies inferred from observations can tentatively be explained with simple galactic dynamo models. Agreement is not perfect, but is reasonable considering the uncertainty in dynamo parameters.
△ Less
Submitted 16 October, 2018;
originally announced October 2018.
-
Accretion in Common Envelope Evolution
Authors:
Luke Chamandy,
Adam Frank,
Eric G. Blackman,
Jonathan Carroll-Nellenback,
Baowei Liu,
Yisheng Tu,
Jason Nordhaus,
Zhuo Chen,
Bo Peng
Abstract:
Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simula…
▽ More
Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.
△ Less
Submitted 10 October, 2018;
originally announced October 2018.
-
Evolution of galactic magnetic fields
Authors:
Luiz Felippe S. Rodrigues,
Luke Chamandy,
Anvar Shukurov,
Carlton M. Baugh,
A. Russ Taylor
Abstract:
We study the cosmic evolution of the magnetic fields of a large sample of spiral galaxies in a cosmologically representative volume by employing a semi-analytic galaxy formation model and numerical dynamo solver in tandem. We start by deriving time- and radius-dependent galaxy properties using the galform galaxy formation model, which are then fed into the nonlinear mean-field dynamo equations. Th…
▽ More
We study the cosmic evolution of the magnetic fields of a large sample of spiral galaxies in a cosmologically representative volume by employing a semi-analytic galaxy formation model and numerical dynamo solver in tandem. We start by deriving time- and radius-dependent galaxy properties using the galform galaxy formation model, which are then fed into the nonlinear mean-field dynamo equations. These are solved to give the large-scale (mean) field as a function of time and galactocentric radius for a thin disc, assuming axial symmetry. A simple prescription for the evolution of the small-scale (random) magnetic field component is also adopted. We find that, while most massive galaxies are predicted to have large-scale magnetic fields at redshift z=0, a significant fraction of them are expected to contain negligible large-scale field. Our model indicates that, for most of the galaxies containing large-scale magnetic fields today, the mean-field dynamo becomes active at z<3. We compute the radial profiles of pitch angle, and find broad agreement with observational data for nearby galaxies.
△ Less
Submitted 5 December, 2018; v1 submitted 27 September, 2018;
originally announced September 2018.
-
Non-linear Galactic Dynamos and the Magnetic Rädler Effect
Authors:
Luke Chamandy,
Nishant K. Singh
Abstract:
We show that the magnetic analogue of the Rädler effect of mean-field dynamo theory leads to a non-linear backreaction that quenches a large-scale galactic dynamo, and can result in saturation of the large-scale magnetic field at near-equipartition with turbulent kinetic energy density. In a rotating fluid containing small-scale magnetic fluctuations, anisotropic terms in the mean electromotive fo…
▽ More
We show that the magnetic analogue of the Rädler effect of mean-field dynamo theory leads to a non-linear backreaction that quenches a large-scale galactic dynamo, and can result in saturation of the large-scale magnetic field at near-equipartition with turbulent kinetic energy density. In a rotating fluid containing small-scale magnetic fluctuations, anisotropic terms in the mean electromotive force are induced via the Coriolis effect and these terms lead to a reduction of the growth rate in a predominantly $αΩ$-type galactic dynamo (Chamandy & Singh 2017). By including the generation of small-scale magnetic fluctuations by turbulent tangling of the large-scale magnetic field, one obtains a negative feedback effect that quenches the dynamo and leads to the saturation of the large-scale field. This saturation mechanism is found to be competitive with the dynamical $α$-quenching mechanism for realistic galactic parameter values. Furthermore, in the context of the dynamical $α$-quenching model, a separate non-linear term is obtained which has the same form as the helicity flux term of Vishniac & Cho (2001), but which depends on the strength of small-scale magnetic fluctuations. We briefly discuss the observational implications of the magnetic Rädler effect for galaxies.
△ Less
Submitted 24 August, 2018; v1 submitted 14 May, 2018;
originally announced May 2018.
-
Accretion in common envelope evolution
Authors:
Luke Chamandy,
Adam Frank,
Eric G. Blackman,
Jonathan Carroll-Nellenback,
Baowei Liu,
Yisheng Tu,
Jason Nordhaus,
Zhuo Chen,
Bo Peng
Abstract:
Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the…
▽ More
Common envelope evolution (CEE) is presently a poorly understood, yet critical, process in binary stellar evolution. Characterizing the full 3D dynamics of CEE is difficult in part because simulating CEE is so computationally demanding. Numerical studies have yet to conclusively determine how the envelope ejects and a tight binary results, if only the binary potential energy is used to propel the envelope. Additional power sources might be necessary and accretion onto the inspiraling companion is one such source. Accretion is likely common in post-asymptotic giant branch (AGB) binary interactions but how it operates and how its consequences depend on binary separation remain open questions. Here we use high resolution global 3D hydrodynamic simulations of CEE with the adaptive mesh refinement (AMR) code AstroBEAR, to bracket the range of CEE companion accretion rates by comparing runs that remove mass and pressure via a subgrid accretion model with those that do not. The results show that if a pressure release valve is available, super-Eddington accretion may be common. Jets are a plausible release valve in these environments, and they could also help unbind and shape the envelopes.
△ Less
Submitted 9 May, 2018;
originally announced May 2018.
-
Derivation and precision of mean field electrodynamics with mesoscale fluctuations
Authors:
Hongzhe Zhou,
Eric G. Blackman,
Luke Chamandy
Abstract:
Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two ques…
▽ More
Mean field electrodynamics (MFE) facilitates practical modeling of secular, large scale properties of astrophysical or laboratory systems with fluctuations.Practitioners commonly assume wide scale separation between mean and fluctuating quantities, to justify equality of ensemble and spatial or temporal averages.Often however, real systems do not exhibit such scale separation. This raises two questions: (I) what are the appropriate generalized equations of MFE in the presence of mesoscale fluctuations? (II) how precise are theoretical predictions from MFE? We address both by first deriving the equations of MFE for different types of averaging, along with mesoscale correction terms that depend on the ratio of averaging scale to variation scale of the mean. We then show that even if these terms are small, predictions of MFE can still have a significant precision error. This error has an intrinsic contribution from the dynamo input parameters and a filtering contribution from differences in the way observations and theory are projected through the measurement kernel.Minimizing the sum of these contributions can produce an optimal scale of averaging that makes the theory maximally precise.The precision error is important to quantify when comparing to observations because it quantifies the resolution of predictive power. We exemplify these principles for galactic dynamos, comment on broader implications, and identify possibilities for further work.
△ Less
Submitted 21 March, 2019; v1 submitted 6 October, 2017;
originally announced October 2017.
-
A new constraint on mean-field galactic dynamo theory
Authors:
Luke Chamandy,
Nishant K. Singh
Abstract:
Appealing to an analytical result from mean-field theory, we show, using a generic galaxy model, that galactic dynamo action can be suppressed by small-scale magnetic fluctuations. This is caused by the magnetic analogue of the Rädler or $Ω\times J$ effect, where rotation-induced corrections to the mean-field turbulent transport result in what we interpret to be an effective reduction of the stand…
▽ More
Appealing to an analytical result from mean-field theory, we show, using a generic galaxy model, that galactic dynamo action can be suppressed by small-scale magnetic fluctuations. This is caused by the magnetic analogue of the Rädler or $Ω\times J$ effect, where rotation-induced corrections to the mean-field turbulent transport result in what we interpret to be an effective reduction of the standard $α$ effect in the presence of small-scale magnetic fields.
△ Less
Submitted 10 April, 2017; v1 submitted 15 December, 2016;
originally announced December 2016.
-
Statistical tests of galactic dynamo theory
Authors:
Luke Chamandy,
Anvar Shukurov,
A. Russ Taylor
Abstract:
Mean-field galactic dynamo theory is the leading theory to explain the prevalence of regular magnetic fields in spiral galaxies, but its systematic comparison with observations is still incomplete and fragmentary. Here we compare predictions of mean-field dynamo models to observational data on magnetic pitch angle and the strength of the mean magnetic field. We demonstrate that a standard $α^2Ω$ d…
▽ More
Mean-field galactic dynamo theory is the leading theory to explain the prevalence of regular magnetic fields in spiral galaxies, but its systematic comparison with observations is still incomplete and fragmentary. Here we compare predictions of mean-field dynamo models to observational data on magnetic pitch angle and the strength of the mean magnetic field. We demonstrate that a standard $α^2Ω$ dynamo model produces pitch angles of the regular magnetic fields of nearby galaxies that are reasonably consistent with available data. The dynamo estimates of the magnetic field strength are generally within a factor of a few of the observational values. Reasonable agreement between theoretical and observed pitch angles generally requires the turbulent correlation time $τ$ to be in the range 10-20 Myr, in agreement with standard estimates. Moreover, good agreement also requires that the ratio of the ionized gas scale height to root-mean-square turbulent velocity increases with radius. Our results thus widen the possibilities to constrain interstellar medium (ISM) parameters using observations of magnetic fields. This work is a step toward systematic statistical tests of galactic dynamo theory. Such studies are becoming more and more feasible as larger datasets are acquired using current and up-and-coming instruments.
△ Less
Submitted 19 September, 2016;
originally announced September 2016.
-
An analytical dynamo solution for large-scale magnetic fields of galaxies
Authors:
Luke Chamandy
Abstract:
We present an effectively global analytical asymptotic galactic dynamo solution for the regular magnetic field of an axisymmetric thin disc in the saturated state. This solution is constructed by combining two well-known types of local galactic dynamo solution, parameterized by the disc radius. Namely, the critical (zero growth) solution obtained by treating the dynamo equation as a perturbed diff…
▽ More
We present an effectively global analytical asymptotic galactic dynamo solution for the regular magnetic field of an axisymmetric thin disc in the saturated state. This solution is constructed by combining two well-known types of local galactic dynamo solution, parameterized by the disc radius. Namely, the critical (zero growth) solution obtained by treating the dynamo equation as a perturbed diffusion equation is normalized using a non-linear solution that makes use of the `no-$z$' approximation and the dynamical $α$-quenching non-linearity. This overall solution is found to be reasonably accurate when compared with detailed numerical solutions. It is thus potentially useful as a tool for predicting observational signatures of magnetic fields of galaxies. In particular, such solutions could be painted onto galaxies in cosmological simulations to enable the construction of synthetic polarized synchrotron and Faraday rotation measure (RM) datasets. Further, we explore the properties of our numerical solutions, and their dependence on certain parameter values. We illustrate and assess the degree to which numerical solutions based on various levels of approximation, common in the dynamo literature, agree with one another.
△ Less
Submitted 1 September, 2016; v1 submitted 3 August, 2016;
originally announced August 2016.
-
Nonlinear Galactic Dynamos and the Magnetic Pitch Angle
Authors:
Luke Chamandy,
A. Russ Taylor
Abstract:
Pitch angles $p$ of the large-scale magnetic fields $\overline{\bf{\it{B}}}$ of spiral galaxies have previously been inferred from observations to be systematically larger in magnitude than predicted by standard mean-field dynamo theory. This discrepancy is more pronounced if dynamo growth has saturated, which is reasonable to assume given that such fields are generally inferred to be close to ene…
▽ More
Pitch angles $p$ of the large-scale magnetic fields $\overline{\bf{\it{B}}}$ of spiral galaxies have previously been inferred from observations to be systematically larger in magnitude than predicted by standard mean-field dynamo theory. This discrepancy is more pronounced if dynamo growth has saturated, which is reasonable to assume given that such fields are generally inferred to be close to energy equipartition with the interstellar turbulence. This 'pitch angle problem' is explored using local numerical mean-field dynamo solutions as well as asymptotic analytical solutions. It is first shown that solutions in the saturated or kinematic regimes depend on only five dynamo parameters, two of which are tightly constrained by observations of galaxy rotation curves. The remaining 3-dimensional (dimensionless) parameter space can be constrained to some extent using theoretical arguments. Predicted values of $|p|$ can be as large as $\sim40^\circ$, which is similar to the largest values inferred from observations, but only for a small and non-standard region of parameter space. We argue, based on independent evidence, that such non-standard parameter values are plausible. However, these values are located toward the boundary of the allowed parameter space, suggesting that additional physical effects may need to be incorporated. We therefore suggest possible directions for extending the basic model considered.
△ Less
Submitted 17 July, 2015; v1 submitted 10 June, 2015;
originally announced June 2015.
-
Magnetic Spiral Arms and Galactic Outflows
Authors:
Luke Chamandy,
Anvar Shukurov,
Kandaswamy Subramanian
Abstract:
Galactic magnetic arms have been observed between the gaseous arms of some spiral galaxies; their origin remains unclear. We suggest that magnetic spiral arms can be naturally generated in the interarm regions because the galactic fountain flow or wind is likely to be weaker there than in the arms. Galactic outflows lead to two countervailing effects: removal of small-scale magnetic helicity, whic…
▽ More
Galactic magnetic arms have been observed between the gaseous arms of some spiral galaxies; their origin remains unclear. We suggest that magnetic spiral arms can be naturally generated in the interarm regions because the galactic fountain flow or wind is likely to be weaker there than in the arms. Galactic outflows lead to two countervailing effects: removal of small-scale magnetic helicity, which helps to avert catastrophic quenching of the dynamo, and advection of the large-scale magnetic field, which suppresses dynamo action. For realistic galactic parameters, the net consequence of outflows being stronger in the gaseous arms is higher saturation large-scale field strengths in the interarm regions as compared to in the arms. By incorporating rather realistic models of spiral structure and evolution into our dynamo models, an interlaced pattern of magnetic and gaseous arms can be produced.
△ Less
Submitted 17 September, 2014; v1 submitted 18 August, 2014;
originally announced August 2014.
-
Non-linear galactic dynamos: A toolbox
Authors:
Luke Chamandy,
Anvar Shukurov,
Kandaswamy Subramanian,
Katherine Stoker
Abstract:
We compare various models and approximations for non-linear mean-field dynamos in disc galaxies to assess their applicability and accuracy, and thus to suggest a set of simple solutions suitable to model the large-scale galactic magnetic fields in various contexts. The dynamo saturation mechanisms considered are the magnetic helicity balance involving helicity fluxes (the dynamical $α$-quenching)…
▽ More
We compare various models and approximations for non-linear mean-field dynamos in disc galaxies to assess their applicability and accuracy, and thus to suggest a set of simple solutions suitable to model the large-scale galactic magnetic fields in various contexts. The dynamo saturation mechanisms considered are the magnetic helicity balance involving helicity fluxes (the dynamical $α$-quenching) and an algebraic $α$-quenching. The non-linear solutions are then compared with the marginal kinematic and asymptotic solutions. We also discuss the accuracy of the no-$z$ approximation. Although these tools are very different in the degree of approximation and hence complexity, they all lead to remarkably similar solutions for the mean magnetic field. In particular, we show that the algebraic $α$-quenching non-linearity can be obtained from a more physical dynamical $α$-quenching model in the limit of nearly azimuthal magnetic field. This suggests, for instance, that earlier results on galactic disc dynamos based on the simple algebraic non-linearity are likely to be reliable, and that estimates based on simple, even linear models are often a good starting point. We suggest improved no-$z$ and algebraic $α$-quenching models, and also incorporate galactic outflows into a simple analytical dynamo model to show that the outflow can produce leading magnetic spirals near the disc surface. The simple dynamo models developed are applied to estimate the magnetic pitch angle and the arm-interarm contrast in the saturated magnetic field strength for realistic parameter values.
△ Less
Submitted 24 June, 2014; v1 submitted 11 March, 2014;
originally announced March 2014.
-
Magnetic arms generated by multiple interfering galactic spiral patterns
Authors:
Luke Chamandy,
Kandaswamy Subramanian,
Alice Quillen
Abstract:
Interfering two- and three-arm spiral patterns have previously been inferred to exist in many galaxies and also in numerical simulations, and invoked to explain important dynamical properties, such as lack of symmetry, kinks in spiral arms, and star formation in armlets. The non-axisymmetric galactic mean-field dynamo model of Chamandy et al. is generalized to allow for such multiple co-existing s…
▽ More
Interfering two- and three-arm spiral patterns have previously been inferred to exist in many galaxies and also in numerical simulations, and invoked to explain important dynamical properties, such as lack of symmetry, kinks in spiral arms, and star formation in armlets. The non-axisymmetric galactic mean-field dynamo model of Chamandy et al. is generalized to allow for such multiple co-existing spiral patterns in the kinetic alpha_k effect, leading to the existence of magnetic spiral arms in the large-scale magnetic field with several new properties. The large-scale magnetic field produced by an evolving superposition of two- and three-arm (or two- and four-arm) patterns evolves with time along with the superposition. Magnetic arms can be stronger and more extended in radius and in azimuth when produced by two interfering patterns rather than by one pattern acting alone. Transient morphological features arise in the magnetic arms, including bifurcations, disconnected armlets, and temporal and spatial variation in arm strength and pitch angles. Pitch angles of the large-scale magnetic field and magnetic arm structures (ridges) are smaller than those typically inferred from observations of spiral galaxies for model parameters of Chamandy et al., but can become comparable to typically inferred values for certain (still realistic) parameters. The magnetic field is sometimes strongest in between the alpha_k-arms, unlike in standard models with a single pattern, where it is strongest within the alpha_k-arms. Moreover, for models with a two- and three-arm pattern, some amount of m=1 azimuthal symmetry is found to be present in the magnetic field, which is generally not the case for forcing by single two- or three-arm patterns. (abridged)
△ Less
Submitted 28 October, 2013; v1 submitted 2 August, 2013;
originally announced August 2013.
-
Galactic Spiral Patterns and Dynamo Action II: Asymptotic Solutions
Authors:
Luke Chamandy,
Kandaswamy Subramanian,
Anvar Shukurov
Abstract:
The exploration of mean-field galactic dynamos affected by a galactic spiral pattern, begun in Chamandy et al. (2013, hereafter Paper I) with numerical simulations, is continued here with an asymptotic solution. The mean-field dynamo model used generalizes the standard theory to include the delayed response of the mean electromotive force to variations of the mean magnetic field and turbulence par…
▽ More
The exploration of mean-field galactic dynamos affected by a galactic spiral pattern, begun in Chamandy et al. (2013, hereafter Paper I) with numerical simulations, is continued here with an asymptotic solution. The mean-field dynamo model used generalizes the standard theory to include the delayed response of the mean electromotive force to variations of the mean magnetic field and turbulence parameters (the temporal non-locality, or τ effect). The effect of the spiral pattern on the dynamo considered is the enhancement of the α-effect in spiral-shaped regions (which may overlap the gaseous spiral arms or be located in the interarm regions). The axisymmetric and enslaved non-axisymmetric modes of the mean magnetic field are studied semi-analytically to clarify and strengthen the numerical results. Good qualitative agreement is obtained between the asymptotic solution and numerical solutions of Paper I for a global, rigidly rotating material spiral (density wave). At all galactocentric distances except for the co-rotation radius, we find magnetic arms displaced in azimuth from the α-arms, so that the ridges of magnetic field strength are more tightly wound than the α-arms. Moreover, the effect of a finite dynamo relaxation time τ (related to the turbulence correlation time) is to phase-shift the magnetic arms in the direction opposite to the galactic rotation even at the co-rotation radius. This mechanism can be used to explain the phase shifts between magnetic and material arms observed in some spiral galaxies.
△ Less
Submitted 11 July, 2013; v1 submitted 21 January, 2013;
originally announced January 2013.
-
Galactic spiral patterns and dynamo action I: A new twist on magnetic arms
Authors:
Luke Chamandy,
Kandaswamy Subramanian,
Anvar Shukurov
Abstract:
We generalise the theory of mean-field galactic dynamos by allowing for temporal non-locality in the mean electromotive force (emf). This arises in random flows due to a finite response time of the mean emf to changes in the mean magnetic field, and leads to the telegraph equation for the mean field. The resulting dynamo model also includes the nonlinear dynamo effects arising from magnetic helici…
▽ More
We generalise the theory of mean-field galactic dynamos by allowing for temporal non-locality in the mean electromotive force (emf). This arises in random flows due to a finite response time of the mean emf to changes in the mean magnetic field, and leads to the telegraph equation for the mean field. The resulting dynamo model also includes the nonlinear dynamo effects arising from magnetic helicity balance. Within this framework, coherent large-scale magnetic spiral arms superimposed on the dominant axially symmetric magnetic structure are considered. A non-axisymmetric forcing of the mean-field dynamo by a spiral pattern (either stationary or transient) is invoked, with the aim of explaining the phenomenon of magnetic arms. For a stationary dynamo forcing by a rigidly rotating material spiral, we find corotating non-axisymmetric magnetic modes enslaved to the axisymmetric modes and strongly peaked around the corotation radius. For a forcing by transient material arms wound up by the galactic differential rotation, the magnetic spiral is able to adjust to the winding so that it resembles the material spiral at all times. There are profound effects associated with the temporal non-locality, i.e. finite `dynamo relaxation time'. For the case of a rigidly rotating spiral, a finite relaxation time causes each magnetic arm to mostly lag the corresponding material arm with respect to the rotation. For a transient material spiral that winds up, the finite dynamo relaxation time leads to a large, negative (in the sense of the rotation) phase shift between the magnetic and material arms, similar to that observed in NGC 6946 and other galaxies. We confirm that sufficiently strong random seed fields can lead to global reversals of the regular field along the radius whose long-term survival depends on specific features of a given galaxy.
△ Less
Submitted 30 October, 2012; v1 submitted 26 July, 2012;
originally announced July 2012.
-
Exact Solutions with w - modes
Authors:
Mustapha Ishak,
Luke Chamandy,
Nicholas Neary,
Kayll Lake
Abstract:
An explicit necessary condition for the occurrence of resonance scattering of axial gravitational waves, along with the internal trapping of null geodesics, is proposed for static spherically symmetric perfect fluid solutions to Einstein's equations. Some exact inhomogeneous solutions which exhibit this trapping are given with special attention to boundary conditions and the physical acceptabili…
▽ More
An explicit necessary condition for the occurrence of resonance scattering of axial gravitational waves, along with the internal trapping of null geodesics, is proposed for static spherically symmetric perfect fluid solutions to Einstein's equations. Some exact inhomogeneous solutions which exhibit this trapping are given with special attention to boundary conditions and the physical acceptability of the space times. In terms of the tenuity ($α= R/M$ at the boundary) all the examples given lie in the narrow range $2.8 < α< 2.9$. The tenuity can be raised to more interesting values by the addition of an envelope without altering the trapping.
△ Less
Submitted 30 March, 2001; v1 submitted 27 July, 2000;
originally announced July 2000.