MNRAS 000 ,

1-7 (2019)

## Hydrodynamic Simulations of Disrupted Planetary Accretion Discs Inside the Core of an AGB Star

- G. Guidarelli 1 ? , J. Nordhaus 1 2 ; y , L. Chamandy 3 , Z. Chen , E. G. Blackman , 4 3
- A. Frank , J. Carroll-Nellenback 3 3 5 6 ; ; , B. Liu 3 5 6 ; ; ,
- 1 Center for Computational Relativity and Gravitation, Rochester Institute of Technology, Rochester, NY, 14623, USA
- 2 National Technical Institute for the Deaf, Rochester Institute of Technology, Rochester, NY, 14623, USA
- 3 Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA

4

Department of Physics, University of Alberta, Edmonton, AB, T6G 2E1, Canada

- 5 Center for Integrated Research Computing, University of Rochester, Rochester, NY 14627, USA
- 6 Laboratory for Laser Energetics, University of Rochester, Rochester, NY 14623, USA

Accepted XXX. Received YYY; in original form ZZZ

## 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 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.

Key words: hydrodynamics - planet-star interactions - stars: AGB and post-AGB - binaries: general - stars: interiors - stars: magnetic field

## 1 INTRODUCTION

White dwarfs (WDs) are the terminal end products of stellar evolution for main-sequence stars with initial masses glyph[lessorsimilar] 8 M GLYPH&lt;12&gt; . Zeeman spectroscopy has revealed a peculiar subset of WDs with abnormally high surface-averaged magnetic fields; aptly termed high-field magnetic white dwarfs (HFMWD) (Preston 1970). Constituting GLYPH&lt;24&gt; 10 % of all WDs, HFMWDs have magnetic fields of GLYPH&lt;24&gt; 10 4 GLYPH&lt;0&gt; 10 9 Gauss while typical WDs have a magnetic field of GLYPH&lt;24&gt; 10 1 GLYPH&lt;0&gt; 10 3 Gauss (Schmidt et al. 2003; Liebert et al. 2003; Kepler et al. 2013; Kepler et al. 2017; Ferrario et al. 2015). They are more massive than generic WDs (Kawka et al. 2007) and none have been found in detached binary systems with M dwarf companions (Liebert et al. 2005; Silvestri et al. 2006, 2007).

field is a remnant from the progenitor main-sequence star (Angel et al. 1981; Wickramasinghe &amp; Ferrario 2005) and more recently the plausibility of magnetic field generation via an internal dynamo as the white dwarf crystallizes has been explored (Isern et al. 2017). However, neither of these theories easily explain the mass disparity or the lack of HFMWDs in non-interacting binaries. An alternative is that the magnetic field is correlated with close binary interactions such as mergers or common envelope evolution (CEE).

The origin of the strong magnetic fields present in HFMWDs is debated. There are arguments the magnetic

- ? E-mail: gcg3642@g.rit.edu
- y E-mail: nordhaus@astro.rit.edu

CEE creates a pathway for reducing binary separations and leads to some of the most interesting and energetic events in the observable universe (Paczynski 1976; Nordhaus &amp; Blackman 2006; Staff et al. 2016; Belczynski et al. 2016; Wilson &amp; Nordhaus 2019). CEE may be a key process in forming a HFMWD as it can lead to WD-WD mergers and the engulfment and accretion of lower-mass companions. A hot corona formed after a WD-WD merger might last long enough to generate a strong magnetic field consistent with observations of HFMWDs in the solar neighborhood (Garc´ ıa-Berro et al. 2012). This scenario requires two

common envelope interactions as both main-sequence stars evolve to become white dwarfs. The consistency of the mass distribution of the merged WD remnants with the mass distribution of isolated HFMWDs remains a subject of debate.

Alternatively, the magnetic field might also be the result of CEE between the pre-WD AGB star and a low-mass companion (Nordhaus et al. 2011). As the AGB star expands to form the CE, the relatively slow moving circumbinary material would dynamically drag on the companion thereby reducing its orbit. If the companion mass is sufficiently low, the liberated orbital energy during inspiral will not exceed the envelope's binding energy, resulting in continued orbital decay until the companion reaches the tidal shredding radius. At this radius, the self gravity of the companion is overwhelmed by the differential potential across its surface resulting in extreme deformation and eventual disruption. The subsequent accretion disc formed by the tidally shredded material could amplify (via a dynamo), advect and anchor a magnetic field to the surface of the white dwarf core. A long Ohmic decay timescale ensures that a HFMWD emerges at a later stage of evolution after the star has shed its outer layers. We focus on this scenario, as it is consistent with the stringent observational constraints provided by the Sloan Digital Sky Survey (Nordhaus et al. 2011).

Analytic estimates for the magnetic fields generated via a dynamo operating in discs formed from GLYPH&lt;24&gt; 1-500 M J companions 1 were sufficient to explain the full range of observed HFMWDs (Nordhaus et al. 2011). However, this scenario requires the disc to survive at least several rotation times for the fields to amplify. Generation of strong magnetic fields is a necessary but not necessarily sufficient condition to produce a HFMWD. The disc must also persist and sustain steady-state magnetic fields long enough to accrete onto the proto-WD. In this purely hydrodynamic study, we focus on stability of the disc in the presence of the harsh interior of the AGB star and leave details of magnetic field amplification and advection to future work.

The disc is initially cold ( GLYPH&lt;24&gt; 10 4 GLYPH&lt;0&gt; 5 K) and dense ( GLYPH&lt;24&gt; 1 g cm GLYPH&lt;0&gt; 3 ), compared to the hot ( GLYPH&lt;24&gt; 10 7 GLYPH&lt;0&gt; 8 K) and less dense ( GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 2 g cm GLYPH&lt;0&gt; 3 ) interior of the AGB star. Vertical shear and steep temperature gradients may threaten disc stability as entrainment of hot material will eventually fully mix the disc into the stellar envelope. The central question that motivates this work is to determine if the discs survive the AGB interior long enough so that a strong magnetic field could develop.

In this paper, we present hydrodynamic simulations of accretion discs formed from the disruption of planetary companions in the interior of a 2 M GLYPH&lt;12&gt; AGB star. Section 2 details some additional aspects of the formation scenario while Section 3 outlines our numerical approach, setup and simulation parameters. Section 4 presents our simulation results while we conclude and discuss future directions in Section 5.

## 2 ADDITIONAL ASPECTS OF THE FORMATION SCENARIO

Before detailing our numerical approach, we note that this scenario requires the companion to survive until the point

Figure 1. Stellar density and pressure profiles before and after modification.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-d768b34e-d3f4-4631-baf7-344407a37366.jpg)

at which it is tidally disrupted and forms a disc. While it is likely that low-mass companions survive inspiral until they are tidally disrupted, there are a few effects that might lead to the planet's destruction before this point which we comment on here.

Villaver &amp; Livio (2007) previously suggested that thermal evaporation during a common envelope phase destroys companions of less than 15-Jupiter masses. However, this scenario requires that the radiative flux penetrate the large optical depth of the planet on a timescale shorter than orbital decay. At 10 11 cm from the center of an extreme-AGB star, the evaporation timescale is on the order of GLYPH&lt;24&gt; 10 3 GLYPH&lt;0&gt; 4 years compared to an orbital decay timescale of a few years, suggesting that the planet likely shreds before significant evaporation occurs (Watson et al. 1981; Murray-Clay et al. 2009; Lopez et al. 2012). Furthermore, brown dwarfs around white dwarfs in post-CE orbits show no evidence that the common envelope phase affected their masses (Maxted et al. 2006; Casewell et al. 2018; Longstaff et al. 2019). In contrast to brown dwarfs, all but the most massive planets are unlikely to survive the common envelope phase based on theoretical arguments and observational searches for planets in white dwarf habitable zones (Nordhaus &amp; Spiegel 2013; Sandhaus et al. 2016; Xu et al. 2015; van Sluijs &amp; Van Eylen 2017; Vanderburg et al. 2015).

Another possibility is that the orbital energy dissipated by drag might be communicated to the interior of the planet, unbinding it before it tidally disrupts. This is unlikely to occur via mechanical means in light of the large density contrast between the AGB envelope and the planetary interior. Aerodynamic drag exerts a steady pressure on the leading face of the planet, and turbulent mixing will be inhibited by the steep entropy gradient at the planet's surface. While it appears likely that even low-mass planets inspiral until they are tidally disrupted, high-resolution CE simulations with self-gravitating planets in hydrostatic equilibrium (HSE) are needed to ultimately verify these assumptions.

Additionally, we note that the disc masses considered in this work are less than their typical progenitor companions. According to StarSmasher simulations of tidal interactions of Jupiter-mass companions around low-mass/intermediatemass black holes, GLYPH&lt;24&gt; 50% of the companion mass remains bound after disruption (Perets et al. 2016). We expect a similar retention ratio for our systems.

The fraction of MS stars with one or more companions has been estimated to be GLYPH&lt;24&gt; 0 44 : (Raghavan et al. 2010). The mass fraction distribution of intermediate-mass mainsequence stars follows f ' q ' / q GLYPH&lt;27&gt; where q is the mass ratio and GLYPH&lt;27&gt; = GLYPH&lt;0&gt; 0 5 : GLYPH&lt;6&gt; 0 2 : (Duchˆ ene &amp; Kraus 2013). This implies that GLYPH&lt;24&gt; 25% of the main sequence stars with companions have q glyph[greaterorsimilar] 0 125 : . This is on the order of mass ratios predicted to be involved with this common envelope scenario. Therefore, the fraction of main-sequence stars with low-mass companions is of order the fraction of HFMWDs.

## 3 NUMERICAL METHODS

We performed three-dimensional hydrodynamic simulations of accretion discs ( 1 GLYPH&lt;0&gt; 10MJ ) around the core of a 2M GLYPH&lt;12&gt; AGB star. The discs were given an initial Keplerian rotation profile with a radius of 2 GLYPH&lt;2&gt; 10 10 cm and height of 5 GLYPH&lt;2&gt; 10 9 cm . The radial extent of the accretion disc is chosen to be the tidal shredding radius, which is the location where the differential gravitational potential of the proto-WD core exceeds the self gravity of a constant density sphere, and is given by:

$$R _ { s } & \simeq R _ { c } \left ( \frac { 2 M } { M _ { c } } \right ) ^ { 1 / 3 } & ( 1 ) & & \underset { \text{mod} } { \text{smoc} }$$

where R c and M c are the radius and mass of the companion, M is the mass of the primary interior to R s (Nordhaus &amp; Blackman 2006).

The initial disc masses were selected according to previous analytic estimates of systems that produce HFMWDs Nordhaus et al. (2011). The ambient AGB density distribution was determined with the one-dimensional stellar evolution code MESA. With MESA, we evolved a 2 M GLYPH&lt;12&gt; zero-agemain-sequence star with solar metalicity ( z = 0 02 : ) through all phases of its stellar lifetime until it began cooling as a white dwarf. The density and pressure profiles were extracted at the time when the star's radius was at its maximum. This is a reasonable time at which a companion would be engulfed as the physical volume of the star is maximized and tidal torques are the strongest (Nordhaus et al. 2010; Nordhaus &amp; Spiegel 2013).

In order to ensure that boundary conditions do not affect the disc dynamics it is necessary to resolve beyond the disc by at least an order of magnitude ( GLYPH&lt;24&gt; 10 11 cm ). The simulation box size combined with limited computation power means that the innermost region of the AGB density profile cannot be properly resolved. We employ a point particle at the center of the grid to accurately match the gravitational potential. This requires a modification to the original profile to conserve total mass and maintain hydro-static equilibrium. We use the modified Lane-Emden equation to adjust the profile:

$$\frac { 1 } { \xi } \frac { d } { d \xi } \left ( \xi ^ { 2 } \frac { d \theta } { d \xi } + \xi ^ { 2 } \frac { g _ { c } ( \alpha \xi ) } { 4 \pi G \rho _ { c } \alpha } \right ) + \theta ^ { n } = 0. & & ( 2 ) & & \text{(} 2 ) & & \text{app}$$

The Lane-Emden equation is modified by the second term on the left-hand side in Equation 2, representing the gravitational attraction of the point particle (Ohlmann et al. 2017). The polytropic index is given as n , with the re-scaled radial coordinate, GLYPH&lt;24&gt; , and density, GLYPH&lt;18&gt; , defined by:

$$\xi & \equiv \frac { r } { \alpha } & ( 3 ) & \quad \text{$\text{$\text{$\text{$no\quad}$} } \\ & 2 \times 1$$

GLYPH&lt;11&gt;

and

$$\underset { \text{ratio} } { 1 0 }. \underset { \text{ratio} } { \begin{array} { c } \theta \end{array} \end{array} } \begin{pmatrix} \theta \end{pmatrix} \begin{pmatrix} \rho \\ \overline { \rho } _ { c } \end{pmatrix} ^ { 2 / 3 }$$

where

$$\text{have} \\ \text{to be} \\ \text{+, the}$$

Lastly, gc is the smoothed gravitational acceleration from the point particle and y GLYPH&lt;17&gt; r GLYPH&lt;157&gt; h :

$$\text{ups} \\ \text{tions} \\ \text{OR} - \text{R} \quad. \quad \begin{cases} \quad \frac { 1 } { 2 } ^ { 2 } } { \gamma ( \frac { 6 4 } { 3 } + y ( - 4 8 + y ( \frac { 1 9 2 } { 5 } + y \frac { - 3 2 } { 3 } ) ) ) - \frac { 2 } { 3 0 y ^ { 2 } } } & \quad r > = h \\ \frac { y ( \frac { 6 4 } { 3 } + y ^ { 2 } ( \frac { - 1 9 2 } { 5 } + 3 2 y ) ) } { h ^ { 2 } } & \quad h / 2 < = r < h \\ \frac { y ( \frac { 3 2 } { 3 } + y ^ { 2 } ( \frac { - 1 9 2 } { 5 } + 3 2 y ) ) } { h ^ { 2 } } & \quad r < h / 2 \\ \text{OR} - \text{R} \quad. \quad \end{cases} \quad. \quad. \quad. \quad. \quad \begin{cases} ( 6 ) \\ \gamma ( \frac { 3 2 } { 3 } + y ^ { 2 } ( \frac { - 1 9 2 } { 5 } + 3 2 y ) ) \\ h ^ { 2 } } & \quad r < h / 2 \\ \gamma ( \frac { 6 ) } { 3 } \end{cases}.$$

Equation 2 was solved using a third order Runge-Kutta integrator such that solutions were accepted when the slope of the density profile matched at the smoothing radius. The resultant profiles in Figure 1 have a reduced constant central density where the profile was previously underresolved. The added point particle ensures the total mass inside the smoothing radius remains the same as the initial model.

A point particle also requires the choice of sub-grid model for losing and gaining mass. The current sub-grid models for accretion available in AstroBEAR remove thermal energy from the grid, a fact that is only appropriate in contexts where the gas exterior to the point particle is optically thin. In our environment this would lead to unphysical accretion rates (Krumholz et al. 2004; Federrath et al. 2010; Chamandy et al. 2018). Therefore, we choose a non-accreting point particle which allows the gas to accumulate naturally in the center. As the pressure builds, the accretion rate decreases until it is eventually halted. Advection of the field and formation of a HFMWD requires that the central pressure be removed, most likely through a bipolar outflow or jet (Nordhaus et al. 2011) and that the quasi-steady state field sustain in the disc until that time.

The fluid in our simulation abides by

$$\int _ { t \text{ a f-} } \frac { \varpi } { \partial t } + \nabla \cdot ( \rho \mathfrak { x } ) = 0,$$

$$\underset { \text{new} } { \overset { \partial \rho \mathbf v } { \underset { \text{$\pi$} } } } \quad \frac { \partial \rho \mathbf v } { \partial t } + \nabla \cdot ( \rho \mathbf v \mathbf v ) = - \nabla P - \rho g _ { c } ( r ) \hat { r },$$

and

$$\frac { \infty } { \text{onal} } \quad P = \tilde { n } k _ { b } T,$$

where GLYPH&lt;26&gt; is the density, t is time, v is the velocity, p is the pressure, r is the radial coordinate, ˜ n is the number density of particles, k b is Boltzmann constant and T is the temperature. With monatomic the GLYPH&lt;13&gt; -law gas equation of state where GLYPH&lt;13&gt; = 5 3 GLYPH&lt;157&gt; .

The total simulated time for each disc mass is GLYPH&lt;24&gt; 22000 s , approximately 10 orbits. The simulations were completed with the three-dimensional, multi-physics, AMR code, AstroBEAR AstroBEAR . utilizes a Riemann solver to solve the fluid equations.

We use 128 3 computational cells with 5 levels of refinement for an effective resolution of 4096 3 with extrapolated boundary conditions. The simulated box side lengths are 2 GLYPH&lt;2&gt; 10 11 cm which means the smallest distance resolved is

Figure 2. Evolution of the one Jupiter mass disc at 2 and 10 orbits. The left-most panel is a mid-plane slice of the density of the disc from a top down view. The other panels show a slice of the temperature of the disc from the side with the right-most panel a magnified view indicated by the black box. Velocity vectors are shown on all plots with magnitude indicated by colour and the projected magnitude by the length.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-4d604e66-311a-4d69-9758-52ab3affc388.jpg)

GLYPH&lt;24&gt; 10 8 cm. This resolution accurately resolves the modified ambient profile transition which is important for maintaining HSE and is typically two orders of magnitude smaller than the smallest grid cells used in CE simulations that resolve the full primary star (Ohlmann et al. 2016; Chamandy et al. 2018; Chamandy et al. 2019; Passy et al. 2012). The computations were run on Stampede 2 at TACC (Towns et al. 2014).

The ambient velocities initially appear stochastic at 2 orbits and develop more structure at 10 orbits as the gas is spun up and circulates with the disc.

## 4 RESULTS

Figures 2-5 present density and temperature snapshots from each simulation at 2 and 10 orbits 2 . In each figure, the left panel shows a face-on view of the density of the mid-plane of the disc. The middle panel presents an edge-on view of temperature while the right panel details the zoomed-in region of the middle panel. The top three panels show the disc structure with corresponding velocity vectors at 2 orbits while the bottom three panels present the same information at 10 orbits.

The temperature plots for each of the discs show pronounced Kelvin-Helmholtz instabilities at 2 orbits. Comparing the profile of each disc at 2 and 10 orbits, it is clear that the discs expand vertically-upward and radiallyoutward with the more massive discs expanding the most.

2 An orbit is defined as the rotation period of the disc's outer radius at the start of the simulation.

Each disc is initially rotating at the Keplerian speed with the ambient stellar interior stationary. As such, the first few orbits are dominated by shear at the interface of the disc edges with the AGB interior. Because the disc is not in perfect HSE with the ambient at the start of the computation, the disc experiences some compression during the first quarter-orbit as seen in Figure 6 after which it reaches a quasi-steady state. As the point particle cannot accrete, the central pressure increases, effectively halting accretion. The result is mass outflow from the disc radially.

Figure 6 shows the disc mass inside of hollow cylinders with outer radius given by the colour key in units of 10 9 cm, and with a thickness of 5 GLYPH&lt;2&gt; 10 8 cm. In Figure 7, we present the mass outflow rate of the disc, i.e. the rate of change of the total disc mass inside a cylinder of radius 5 GLYPH&lt;2&gt; 10 9 cm. In aggregate, these simulations show that the discs become stable after a brief relaxation phase even in an environment where the initial shear is maximized. Note that global 3D simulations of CE phases (the precursor of our systems) demonstrate that as the companion inspirals, the ambient stellar gas spins up via angular momentum transfer and can approach co-rotation with the disc (Ricker &amp; Taam 2008; Ricker &amp; Taam 2012; Iaconi et al. 2018). The faster the rotation of the surrounding AGB star, the less shear there would be between the disc and ambient material, and therefore more likely that the discs survive on long

Figure 3. Evolution of the three Jupiter mass disc at 2 and 10 orbits. (Similar to 2)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-95d726e8-77fe-464f-b28d-811f05f23464.jpg)

timescales. Thus, we conclude that the stability of the discs formed is not threatened by the entrainment of hot material in a strongly sheared, realistic AGB environment.

Our selected resolution for the simulations was the highest allowable given our computational budget and significantly higher than the highest resolution seen in global CE simulations (Ohlmann et al. 2016; Chamandy et al. 2018; Chamandy et al. 2019; Passy et al. 2012). We check the robustness of our results by lowering the resolution of the 3 M J simulation by two levels of refinement, i.e. a factor of 4. The mass outflow rate for the low resolution simulation is shown in Figure 7. Because the mass outflow rate is larger with lower resolution, it is likely that increasing our resolution would further decrease the mass outflow rate and thus increase the disc survival timescales. Note that the lower resolution simulation is far less steady compared to the higher resolution simulations which also suggests that the discs may further stabilize at higher resolution. In summary, it is apparent that the discs persist for at least 10 orbits and perhaps glyph[greaterorsimilar] 100 orbits in all cases if the outflow rates remain constant. Thus any magnetic field generated in the disc would likely have sufficient time to reach the surface of the white dwarf if the central pressure can be relieved via an outflow as was suggested in Nordhaus et al. (2011).

## 5 CONCLUSIONS

In this work, we present 3D adaptive-mesh-refinement hydrodynamic simulations of accretion discs around the core of an AGB star. These discs are expected to form from the tidal disruption of a low-mass companion inside a post-main- sequence star during a common envelope interaction. Such discs are initially cold and dense compared to the hot stellar ambient meaning entrainment of hot gas could dissolve the discs. Our simulations show that despite significant shear and temperature gradients, planetary mass discs could survive and operate on timescales long enough to amplify strong magnetic fields.

Transport of the magnetic fields to the white dwarf surface requires a valve that can relieve central pressure. As opposed to neutron stars and black holes who can remove pressure via neutrino cooling and advection through an event horizon, white dwarfs can in principle relieve pressure through strong outflows or jets. Developing a sub-grid point-particle model that conserves thermal energy and appropriately couples to outflows would allow one to investigate whether the strong magnetic fields generated in the disc can anchor to the proto-white dwarf.

Future studies could also improve upon this work in several aspects. Instead of starting with a well-formed disc, simulations that follow the inspiral of a self-gravitating planet as it tidally disrupts and settles into a disc would improve estimates of the initial disc mass and structure. Furthermore, utilizing the full-MHD capability of AstroBEAR in a tidal disruption simulation, or in the simulation setup described in this work, would allow one to study the amplification and dynamics of the magnetic field.

## ACKNOWLEDGEMENTS

GG and JN acknowledge support from the following grants: NASAHST-AR-15044, NASA HST-AR-14563, NTID SPDI-

## 6 Guidarelli et al.

Figure 4. Evolution of the six Jupiter mass disc at 2 and 10 orbits. (Similar to 2)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-4c1ae5af-09b6-4daa-a765-c5852748b25a.jpg)

Figure 5. Evolution of the ten Jupiter mass disc at 2 and 10 orbits. (Similar to 2)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-2314342c-7e63-46df-8f1b-09fed6b59308.jpg)

Figure 6. disc mass inside a cylindrical shell of thickness 0 5 : GLYPH&lt;2&gt; 10 9 cm for the 3 Jupiter mass disc. The colours indicate the outer radius in units of 10 9 cm . The values are independent of the cylindrical shell height of the cylinder as long as it is greater than the disc height

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-14156ca2-88ec-4178-b6b3-07699950f6d8.jpg)

Figure 7. disc mass outflow rate inside a cylinder of radius 5 GLYPH&lt;2&gt; 10 9 cm centered on the AGB core after the disc reaches a steady state. The violet line is a lower resolution version of the 3 Jupiter mass disc and shows that with the increased resolution we experience slower and more stable mass outflow.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00157v1-figure-0fdc9e0e-9f1c-49e5-adb4-cfe1ccad63c8.jpg)

15992. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562. The Center for Integrated Research Computing (CIRC) at the University of Rochester provided additional computational and visualization resources.

VisIt is supported by the Department of Energy with funding from the Advanced Simulation and Computing Program and the Scientific Discovery through Advanced Computing Program.

## REFERENCES

Angel J. R. P., Borra E. F., Landstreet J. D., 1981, ApJS, 45, 457 Belczynski K., Holz D. E., Bulik T., O'Shaughnessy R., 2016, Nature, 534, 512

Casewell S. L., et al., 2018, MNRAS, 476, 1405

Chamandy L., et al., 2018, MNRAS, 480, 1898

Chamandy L., Tu Y., Blackman E. G., Carroll-Nellenback J., Frank A., Liu B., Nordhaus J., 2019, Monthly Notices of the Royal Astronomical Society, 486, 1070

Duchˆ ene G., Kraus A., 2013, ARA&amp;A, 51, 269

Federrath C., Banerjee R., Clark P. C., Klessen R. S., 2010, ApJ,

713, 269

Ferrario L., de Martino D., G¨nsicke B. T., 2015, Space Science a Reviews, 191, 111

Garc´ ıa-Berro E., et al., 2012, The Astrophysical Journal, 749, 25 Iaconi R., De Marco O., Passy J.-C., Staff J., 2018, MNRAS, 477, 2349

- Isern J., Garc´ ıa-Berro E., K ˜ ijlebi A B., Lor´ en-Aguilar P., 2017, The Astrophysical Journal, 836, L28
- Kawka A., Vennes S., Schmidt G. D., Wickramasinghe D. T., Koch R., 2007, The Astrophysical Journal, 654, 499

Kepler S. O., et al., 2013, Monthly Notices of the Royal Astro- nomical Society, 429, 2934

Kepler S. O., et al., 2017, VizieR Online Data Catalog, 742 Krumholz M. R., McKee C. F., Klein R. I., 2004, ApJ, 611, 399 Liebert J., Bergeron P., Holberg J. B., 2003, The Astronomical Journal, 125, 348

- Liebert J., et al., 2005, The Astronomical Journal, 129, 2376 Longstaff E. S., Casewell S. L., Wynn G. A., Page K. L., Williams

P. K. G., Braker I., Maxted P. F. L., 2019, MNRAS, 484, 2566

Lopez E. D., Fortney J. J., Miller N., 2012, ApJ, 761, 59

- Maxted P. F. L., Napiwotzki R., Dobbie P. D., Burleigh M. R., 2006, Nature, 442, 543
- Murray-Clay R. A., Chiang E. I., Murray N., 2009, ApJ, 693, 23 Nordhaus J., Blackman E. G., 2006, Monthly Notices of the Royal Astronomical Society, 370, 2004
- Nordhaus J., Spiegel D. S., 2013, MNRAS, 432, 500
- Nordhaus J., Spiegel D. S., Ibgui L., Goodman J., Burrows A., 2010, MNRAS, 408, 631

Nordhaus J., Wellons S., Spiegel D. S., Metzger B. D., Blackman E. G., 2011, Proceedings of the National Academy of Science, 108, 3135

- Ohlmann S. T., R¨pke F. K., Pakmor R., Springel V., 2016, The o Astrophysical Journal, 816, L9
- Ohlmann S. T., R¨pke F. K., Pakmor R., Springel V., 2017, A&amp;A, o 599, A5

Paczynski B., 1976, in Eggleton P., Mitton S., Whelan J., eds, , Vol. 73, Structure and Evolution of Close Binary Systems, iau sympos edn, Dordrecht: Reidel, pp 75-80

Passy J. C., et al., 2012, Astrophysical Journal, 744, 52

- Perets H. B., Li Z., Lombardi James C. J., Milcarek Stephen R. J., 2016, ApJ, 823, 113

Preston G. W., 1970, ApJ, 160, L143

- Raghavan D., et al., 2010, ApJS, 190, 1
- Ricker P. M., Taam R. E., 2008, ApJ, 672, L41
- Ricker P. M., Taam R. E., 2012, The Astrophysical Journal, 746, 74
- Sandhaus P. H., Debes J. H., Ely J., Hines D. C., Bourque M., 2016, ApJ, 823, 49

Schmidt G. D., et al., 2003, The Astrophysical Journal, 595, 1101 Silvestri N. M., et al., 2006, The Astronomical Journal, 131, 1674 Silvestri N. M., et al., 2007, The Astronomical Journal, 134, 741 Staff J. E., De Marco O., Wood P., Galaviz P., Passy J.-C., 2016, Monthly Notices of the Royal Astronomical Society, 458, 832

- Towns J., et al., 2014, Computing in Science &amp; Engineering, 16, 62

Vanderburg A., et al., 2015, Nature, 526, 546

Villaver E., Livio M., 2007, ApJ, 661, 1192

Watson A. J., Donahue T. M., Walker J. C. G., 1981, Icarus, 48, 150

Wickramasinghe D. T., Ferrario L., 2005, MNRAS, 356, 1576

Wilson E. C., Nordhaus J., 2019, MNRAS, 485, 4492

- Xu S., Ertel S., Wahhaj Z., Milli J., Scicluna P., Bertrang G. H. M., 2015, A&amp;A, 579, L8
- van Sluijs L., Van Eylen V., 2017, Monthly Notices of the Royal Astronomical Society, 474, 4603

This paper has been typeset from a T E X/LT E X file prepared by A the author.
Gabriel Guidarelli Jason Nordhaus Luke Chamandy Zhuo Chen Eric G. Blackman Adam Frank Jonathan Carroll-Nellenback Baowei Liu
Hydrodynamic Simulations of Disrupted Planetary Accretion Discs Inside the Core of an AGB Star
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.