Source author record

Luke Chamandy

Luke Chamandy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2024arXiv

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

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.

preprint2022arXiv

Jets from main sequence and white dwarf companions during common envelope evolution

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.

preprint2021arXiv

Successive common envelope events from multiple planets

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.

preprint2021arXiv

The Formation of Discs in the Interior of AGB Stars from the Tidal Disruption of Planets and Brown Dwarfs

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.

preprint2020arXiv

Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution

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.

preprint2020arXiv

Common Envelope Evolution on the Asymptotic Giant Branch: Unbinding within a Decade?

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.

preprint2020arXiv

Parameters of the Supernova-Driven Interstellar Turbulence

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.

preprint2018arXiv

Accretion in Common Envelope Evolution

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.

preprint2016arXiv

An analytical dynamo solution for large-scale magnetic fields of galaxies

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.

preprint2016arXiv

Statistical tests of galactic dynamo theory

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.

preprint2015arXiv

Nonlinear Galactic Dynamos and the Magnetic Pitch Angle

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.

preprint2014arXiv

Magnetic Spiral Arms and Galactic Outflows

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.

preprint2014arXiv

Non-linear galactic dynamos: A toolbox

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.

preprint2013arXiv

Galactic Spiral Patterns and Dynamo Action II: Asymptotic Solutions

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.

preprint2013arXiv

Magnetic arms generated by multiple interfering galactic spiral patterns

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)

preprint2012arXiv

Galactic spiral patterns and dynamo action I: A new twist on magnetic arms

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.

preprint2001arXiv

Exact Solutions with w - modes

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.