# Magnetically driven winds from accretion disks in post-asymptotic giant branch binaries

Olivier Verhamme, Jacques Kluska, Jonathan Ferreira, Dylan Bollen, Toon de Prins, Devika Kamath, Hans van Winckel

► To cite this version:

Olivier Verhamme, Jacques Kluska, Jonathan Ferreira, Dylan Bollen, Toon de Prins, et al.. Magnetically driven winds from accretion disks in post-asymptotic giant branch binaries. Astron.Astrophys., 2024, 684, pp.A79. 10.1051/0004-6361/202347708 . hal-04432345

HAL Id: hal-04432345

https://hal.science/hal-04432345

Submitted on 5 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Magnetically driven winds from accretion disks in post-asymptotic giant branch binaries

Olivier Verhamme, Jacques Kluska, Jonathan Ferreira, Dylan Bollen, Toon de Prins, Devika Kamath, Hans van Winckel

► To cite this version:

Olivier Verhamme, Jacques Kluska, Jonathan Ferreira, Dylan Bollen, Toon de Prins, et al. Magnetically driven winds from accretion disks in post-asymptotic giant branch binaries. Astronomy & Astrophysics - A&A, 2024, 684, pp.A79. 10.1051/0004-6361/202347708. hal-04432345

HAL Id: hal-04432345

https://hal.science/hal-04432345v1

Submitted on 5 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Magnetically driven winds from accretion disks in post-asymptotic giant branch binaries

Olivier Verhamme<fnref n="1" />, Jacques Kluska<fnref n="1" />, Jonathan Ferreira<fnref n="2" />, Dylan Bollen<fnref n="1" />,<fnref n="3" />, Toon De Prins<fnref n="1" />, Devika Kamath<fnref n="3" />,<fnref n="4" />, and Hans Van Winckel<fnref n="1" />

Received 11 August 2023 / Accepted 8 December 2023

## ABSTRACT

**Context.** Jets are commonly detected in post-asymptotic giant branch (post-AGB) binaries and originate from an accretion process onto the companion of the post-AGB primary. These jets are revealed by high-resolution spectral time series.

**Aims.** This paper is part of a series. In this work, we move away from our previous parametric modelling and include a self-similar wind model that allows the physical properties of post-AGB binaries to be characterised. This model describes magnetically driven jets from a thin accretion disk threaded by a large-scale, near equipartition vertical field.

**Methods.** We expanded our methodology in order to simulate the high-resolution dynamic spectra coming from the obscuration of the primary by the jets launched by the companion. We present the framework to exploit the self-similar jet models for post-AGB binaries. We performed a parameter study to investigate the impact of different parameters (inclination, accretion rate, inner and outer launching radius) on the synthetic spectra.

**Results.** We successfully included the physical jet models into our framework. The synthetic spectra have a very similar orbital phase coverage and absorption strengths as the observational data. The magnetohydrodynamic (MHD) jet models provide a good representation of the actual jet creation process in these evolved binaries. Challenges remain, however, as the needed high-accretion rate would induce accretion disks that are too hot in comparison to the data. Moreover, the rotational signature of the models is not detected in the observations. In future research, we will explore models with a higher disk ejection efficiency and even lower magnetisation in order to solve some of the remaining discrepancies between the observed and synthetic dynamic spectra.

**Key words.** accretion, accretion disks – stars: AGB and post-AGB – binaries: spectroscopic – stars: jets

## 1. Introduction

In this paper, we discuss disk winds launched by main-sequence companions in a binary with a post-asymptotic giant branch (post-AGB) primary (Van Winckel 2003). We use the term “post-AGB star” to refer to a star that has just evolved off the AGB branch, meaning the star is contracting and heating at a constant luminosity up to a value of $10^5 L_\odot$, with mass-loss rate estimates between $10^{-9}$ and $10^{-7} M_\odot \text{yr}^{-1}$ (Oomen et al. 2019; Van Winckel 2003, 2018). The most important building blocks (see Fig. 1) are the evolved primary and its main-sequence companion, which are in an orbit ranging from a couple hundred to a few thousand days and are surrounded by a stable circumbinary disk of gas and dust. A catalogue of Galactic binary post-AGB systems was recently published by Kluska et al. (2022).

The luminosity of the primary dominates the optical spectrum, but the thermal contribution of the accretion disk around the companion can be resolved with near-infrared interferometry (Hillen et al. 2016; Kluska et al. 2018). The observation of an accretion disk strengthens the belief that the bipolar outflow is launched from the circum-companion accretion disk. This outflow is not resolved, but its presence can be deduced indirectly by time-resolved high-resolution optical spectra. The outflow is

seen as a high-velocity and phase-dependent absorption feature in the Balmer lines, which also follow the movement of the secondary companion. This absorption component appears around a superior conjunction (Van Winckel 2007; Gorlova et al. 2012; Bollen et al. 2017). When the companion moves in front of the primary, one of the cones of the jet gradually enters the line of sight (LOS) towards the primary. This creates an absorption component as continuum photons from the background post-AGB primary are selectively removed from the LOS. We illustrate this in Fig. 2. In the figure, we show our time series of high-resolution optical spectra as obtained by our HERMES spectrograph (Raskin 2011; Raskin et al. 2011) mounted on the Mercator telescope at the Roque de los Muchachos observatory in Spain. We folded the spectra around $H\alpha$ on the orbital period and display the spectra in velocity space. At the superior conjunction, one can see the outflow in absorption when the bicone is in the LOS towards the primary.

The observed collimated outflows can be described as extended disk winds launched from this accretion disk (Blandford & Payne 1982; Ferreira & Casse 2013). This type of magneto-centrifugal outflow is formed when the magnetic field lines are threading the disk at an angle acute enough for the material to be transported by the centrifugal force over the field lines

<footnote n="1">¹ Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
e-mail: olivier.verhamme@kuleuven.be</footnote>

<footnote n="2">² Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="3">³ Department of Physics & Astronomy, School of Mathematical and Physical Sciences, Macquarie University, Sydney, NSW 2109, Australia</footnote>

<footnote n="4">⁴ Astronomy, Astrophysics and Astrophotonics Research Centre, Macquarie University, Sydney, NSW 2109, Australia</footnote>
# Magnetically driven winds from accretion disks in post-asymptotic giant branch binaries

Olivier Verhamme<fnref n="1" />, Jacques Kluska<fnref n="1" />, Jonathan Ferreira<fnref n="2" />, Dylan Bollen<fnref n="1" />,<fnref n="3" />, Toon De Prins<fnref n="1" />, Devika Kamath<fnref n="3" />,<fnref n="4" />, and Hans Van Winckel<fnref n="1" />

Received 11 August 2023 / Accepted 8 December 2023

## ABSTRACT

**Context.** Jets are commonly detected in post-asymptotic giant branch (post-AGB) binaries and originate from an accretion process onto the companion of the post-AGB primary. These jets are revealed by high-resolution spectral time series.

**Aims.** This paper is part of a series. In this work, we move away from our previous parametric modelling and include a self-similar wind model that allows the physical properties of post-AGB binaries to be characterised. This model describes magnetically driven jets from a thin accretion disk threaded by a large-scale, near equipartition vertical field.

**Methods.** We expanded our methodology in order to simulate the high-resolution dynamic spectra coming from the obscuration of the primary by the jets launched by the companion. We present the framework to exploit the self-similar jet models for post-AGB binaries. We performed a parameter study to investigate the impact of different parameters (inclination, accretion rate, inner and outer launching radius) on the synthetic spectra.

**Results.** We successfully included the physical jet models into our framework. The synthetic spectra have a very similar orbital phase coverage and absorption strengths as the observational data. The magnetohydrodynamic (MHD) jet models provide a good representation of the actual jet creation process in these evolved binaries. Challenges remain, however, as the needed high-accretion rate would induce accretion disks that are too hot in comparison to the data. Moreover, the rotational signature of the models is not detected in the observations. In future research, we will explore models with a higher disk ejection efficiency and even lower magnetisation in order to solve some of the remaining discrepancies between the observed and synthetic dynamic spectra.

**Key words.** accretion, accretion disks – stars: AGB and post-AGB – binaries: spectroscopic – stars: jets

## 1. Introduction

In this paper, we discuss disk winds launched by main-sequence companions in a binary with a post-asymptotic giant branch (post-AGB) primary (Van Winckel 2003). We use the term “post-AGB star” to refer to a star that has just evolved off the AGB branch, meaning the star is contracting and heating at a constant luminosity up to a value of $10^5 L_\odot$, with mass-loss rate estimates between $10^{-9}$ and $10^{-7} M_\odot \text{yr}^{-1}$ (Oomen et al. 2019; Van Winckel 2003, 2018). The most important building blocks (see Fig. 1) are the evolved primary and its main-sequence companion, which are in an orbit ranging from a couple hundred to a few thousand days and are surrounded by a stable circumbinary disk of gas and dust. A catalogue of Galactic binary post-AGB systems was recently published by Kluska et al. (2022).

The luminosity of the primary dominates the optical spectrum, but the thermal contribution of the accretion disk around the companion can be resolved with near-infrared interferometry (Hillen et al. 2016; Kluska et al. 2018). The observation of an accretion disk strengthens the belief that the bipolar outflow is launched from the circum-companion accretion disk. This outflow is not resolved, but its presence can be deduced indirectly by time-resolved high-resolution optical spectra. The outflow is

seen as a high-velocity and phase-dependent absorption feature in the Balmer lines, which also follow the movement of the secondary companion. This absorption component appears around a superior conjunction (Van Winckel 2007; Gorlova et al. 2012; Bollen et al. 2017). When the companion moves in front of the primary, one of the cones of the jet gradually enters the line of sight (LOS) towards the primary. This creates an absorption component as continuum photons from the background post-AGB primary are selectively removed from the LOS. We illustrate this in Fig. 2. In the figure, we show our time series of high-resolution optical spectra as obtained by our HERMES spectrograph (Raskin 2011; Raskin et al. 2011) mounted on the Mercator telescope at the Roque de los Muchachos observatory in Spain. We folded the spectra around $H\alpha$ on the orbital period and display the spectra in velocity space. At the superior conjunction, one can see the outflow in absorption when the bicone is in the LOS towards the primary.

The observed collimated outflows can be described as extended disk winds launched from this accretion disk (Blandford & Payne 1982; Ferreira & Casse 2013). This type of magneto-centrifugal outflow is formed when the magnetic field lines are threading the disk at an angle acute enough for the material to be transported by the centrifugal force over the field lines

<footnote n="1">¹ Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
e-mail: olivier.verhamme@kuleuven.be</footnote>

<footnote n="2">² Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="3">³ Department of Physics & Astronomy, School of Mathematical and Physical Sciences, Macquarie University, Sydney, NSW 2109, Australia</footnote>

<footnote n="4">⁴ Astronomy, Astrophysics and Astrophotonics Research Centre, Macquarie University, Sydney, NSW 2109, Australia</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Conceptual figure of post-AGB binaries. Figure is from Bollen et al. (2022).</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Observable influence of the orbit and jet on the spectra. The top panel gives the dynamic observed spectra of the $H\alpha$ profile of the post-AGB star HD 52961 folded onto the orbital period. The initial phase is random. The lower panel gives the same but for the C I multiplet around 7115 Å. The orbital phase is shown twice to guide the eye.</figcaption></figure>

away from the accretion disk. This type of disk wind can occur even at low magnetisation (Jacquemin-Ide et al. 2019, 2021).

This paper is a continuation of a series of papers (Bollen et al. 2017, 2019, 2020, 2022) describing both the observational effort as well as the modelling efforts to exploit the spectral time series in order to model the jets and to use these systems to study the launching physics of jets. In our previous papers,

we focussed on developing a geometric parametric model to fit the observed time series. This was followed by a line-tracing radiative transfer module for quantifying the total mass loss in the jets. The methodology is outlined in Bollen et al. (2022) and was applied to a total of 16 objects. The jet model used in these papers is purely geometric, with a jet axis and an angle-dependent density and velocity law. The results of these studies have also shown that high mass-ejection rates are needed (reaching values up to $10^{-4} M_{\odot} \text{yr}^{-1}$) and hence so are high accretion rates onto the accretion disk from the circumbinary disk. To further study these surprisingly large mass accretion rates, further investigation is necessary.

To that end, we use the same setup in this work but move away from the geometric models and include self-similar physical jet models. These were developed in the framework jets in YSOs (Jacquemin-Ide et al. 2019). While the long-term goal is to reproduce the data quantitatively, in this work, we present the framework to use these self-similar jet models for post-AGB binary systems. We qualitatively compare the results with the observables of one target and observe the influence of the different model parameters without fine-tuning the model to reproduce the data. We organised the paper as follows: In Sect. 2, we give a short overview of the methodology employed up to now using geometric biconal models. Next, we focus on the properties of the self-similar jet models and how these are implemented (Sect. 3). We perform a parameter study of the physical parameters and focus on their impact on the observables, namely, the time series of high-resolution $H\alpha$ profiles. We analyse our parameter study in Sect. 5. We conclude our paper by highlighting our most important findings and by isolating some avenues for future research.

## 2. Initial framework: The parametric model

We give a brief description of the methodology used in previous studies (Bollen et al. 2022, and references therein). The parametric model is a purely geometric model of the jet, and we use it to generate synthetic dynamical spectra that can be compared to the observed spectra. The first step in this approach is to define the initial spectrum that will serve as the basis on which the absorption from the jet model is added (Sect. 2.1). This initial spectrum was also used for the self-similar models in this work. Then, the parameters of the jet model are constrained as described in Sect. 2.2. Finally, the absorption profile is obtained by performing a radiative transfer step (Sect. 2.4).

### 2.1. Initial spectra

To estimate the impact of an obscuring jet bicone, it is necessary to determine the unobscured spectrum at all orbital phases. This unobscured spectrum is assumed to be the spectrum observed with no circumsecondary jet and serves as a background spectrum from which the jet obscuration is computed. One of the issues in determining the background spectra is that the unabsorbed base spectra of the $H\alpha$ profile does not correspond to a purely photospheric spectrum. Outside the conjunction, the $H\alpha$ profile shows a double-peaked profile, the origin of which is still unclear. It is necessary to create a background spectra that takes this into account. In this process, we assume that this background spectrum is orbit independent.

An example of how the dynamic background spectra appear without the absorption of the jet is given in Fig. 3. To determine this background, we used the average of observations in phases where no jet absorption was seen and assumed that this is a valid
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Conceptual figure of post-AGB binaries. Figure is from Bollen et al. (2022).</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Observable influence of the orbit and jet on the spectra. The top panel gives the dynamic observed spectra of the $H\alpha$ profile of the post-AGB star HD 52961 folded onto the orbital period. The initial phase is random. The lower panel gives the same but for the C I multiplet around 7115 Å. The orbital phase is shown twice to guide the eye.</figcaption></figure>

away from the accretion disk. This type of disk wind can occur even at low magnetisation (Jacquemin-Ide et al. 2019, 2021).

This paper is a continuation of a series of papers (Bollen et al. 2017, 2019, 2020, 2022) describing both the observational effort as well as the modelling efforts to exploit the spectral time series in order to model the jets and to use these systems to study the launching physics of jets. In our previous papers,

we focussed on developing a geometric parametric model to fit the observed time series. This was followed by a line-tracing radiative transfer module for quantifying the total mass loss in the jets. The methodology is outlined in Bollen et al. (2022) and was applied to a total of 16 objects. The jet model used in these papers is purely geometric, with a jet axis and an angle-dependent density and velocity law. The results of these studies have also shown that high mass-ejection rates are needed (reaching values up to $10^{-4} M_{\odot} \text{yr}^{-1}$) and hence so are high accretion rates onto the accretion disk from the circumbinary disk. To further study these surprisingly large mass accretion rates, further investigation is necessary.

To that end, we use the same setup in this work but move away from the geometric models and include self-similar physical jet models. These were developed in the framework jets in YSOs (Jacquemin-Ide et al. 2019). While the long-term goal is to reproduce the data quantitatively, in this work, we present the framework to use these self-similar jet models for post-AGB binary systems. We qualitatively compare the results with the observables of one target and observe the influence of the different model parameters without fine-tuning the model to reproduce the data. We organised the paper as follows: In Sect. 2, we give a short overview of the methodology employed up to now using geometric biconal models. Next, we focus on the properties of the self-similar jet models and how these are implemented (Sect. 3). We perform a parameter study of the physical parameters and focus on their impact on the observables, namely, the time series of high-resolution $H\alpha$ profiles. We analyse our parameter study in Sect. 5. We conclude our paper by highlighting our most important findings and by isolating some avenues for future research.

## 2. Initial framework: The parametric model

We give a brief description of the methodology used in previous studies (Bollen et al. 2022, and references therein). The parametric model is a purely geometric model of the jet, and we use it to generate synthetic dynamical spectra that can be compared to the observed spectra. The first step in this approach is to define the initial spectrum that will serve as the basis on which the absorption from the jet model is added (Sect. 2.1). This initial spectrum was also used for the self-similar models in this work. Then, the parameters of the jet model are constrained as described in Sect. 2.2. Finally, the absorption profile is obtained by performing a radiative transfer step (Sect. 2.4).

### 2.1. Initial spectra

To estimate the impact of an obscuring jet bicone, it is necessary to determine the unobscured spectrum at all orbital phases. This unobscured spectrum is assumed to be the spectrum observed with no circumsecondary jet and serves as a background spectrum from which the jet obscuration is computed. One of the issues in determining the background spectra is that the unabsorbed base spectra of the $H\alpha$ profile does not correspond to a purely photospheric spectrum. Outside the conjunction, the $H\alpha$ profile shows a double-peaked profile, the origin of which is still unclear. It is necessary to create a background spectra that takes this into account. In this process, we assume that this background spectrum is orbit independent.

An example of how the dynamic background spectra appear without the absorption of the jet is given in Fig. 3. To determine this background, we used the average of observations in phases where no jet absorption was seen and assumed that this is a valid
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>Fig. 3. Modelling effect on the dynamic spectra. The top panel shows the initial Hα line-profile of HD 52961 which serves as background spectrum. The lower panel shows the synthetic dynamic spectra of the Hα line using the stellar disk wind parametric model.</figcaption></figure>

background spectrum also at phases of conjunction. To use these spectra correctly, it is important to understand the orbital motion of the system. The wavelength shift of the spectra reveals the radial velocity of the post-AGB star around the centre of mass. We used the orbital parameters of Oomen et al. (2018) given in Table 1. Three of the spectra outside the conjunction contribute to a mean spectrum that serves as a new base spectrum. This spectrum is then used to create a synthetic dynamic spectrum by reintroducing the correct wavelength shift for the specific orbital phase of the post-AGB star.

## 2.2. Modelled components

The model consists of two main components: the post-AGB star itself and the biconical jet centred on the companion. The post-AGB star is orders of magnitude brighter than its typical main-sequence low-mass companion, so we assumed that this binary partner will not contribute to the resulting spectra of the system. In optical wavelengths, the same holds true for the circumbinary disk and the secondary accretion disk. The post-AGB primary and the biconical jet orbit each other with certain binary parameters. At some orbital phases, a LOS from the surface of the post-AGB star passes through the bicone. This ray, characterised by the background spectrum, interacts with the material in the jet at the given Doppler velocity, and this forms the basis of this model.

The post-AGB star is represented by a uniform disk perpendicular to the LOS. This disk is divided into small regions via a Fibonacci grid in which every grid point covers an equal surface area. We then determined the LOS of every grid point in the disk.

The disk wind itself is represented by a cone, which is shown in Fig. 4. This cone has its vertex not located at the secondary star but somewhere under the secondary star, depending on the outer launching radius of the disk wind and the opening angle. The disk wind has two parts: one central fast part and an outer slower part with a higher density. A fit angle determines where the fast region stops and the slow region begins. This angle, together with the offset of the cone, corresponds to a certain inner launching radius.

The inner and outer launching radii cannot be chosen freely. To start, the inner launching radius has to be larger than the stellar radius. The outer disk radius must be smaller than the Roche lobe radius of our main-sequence star, as material outside of this lobe will not be bound to the star (Eggleton 1983):

$$ \frac{r_L}{a} = \frac{0.49q^{2/3}}{0.69q^{2/3} + \ln(1 + q^{1/3})} \quad (1) $$

Here, q is the mass ratio between the two stars, and a is the semi-major axis of the primary. To estimate the radius of the main-sequence companion, we used $R_2 = 1.01(M_2/M_\odot)^{0.724}$ (Demircan & Kahraman 1991).

## 2.3. Fitting procedure

For each phase, the global Hα profiles were determined using the jet parameters of inclination, opening angle, velocity structure, and scaled density distribution as well as the radius of the post-AGB star and its orbit. A Markov chain Monte Carlo (MCMC) fitting routine was used to fit these synthetic spectral time series to the observations. We illustrate this in Fig. 3 panel b, where we show the simulated dynamic spectrum compared with the top panel of Fig. 2, which are the observables themselves. This best fit is based on a disk wind model. The parameters of this best fit are given in the right side of Table 1.

## 2.4. Radiative transfer

To determine the actual mass-loss rates, we needed to compute the absolute densities and assumed temperatures instead of using scaling parameters. To do so, the parametric model was followed by a radiative transfer module using a temperature and density fit while all other parameters from the parametric fitting were kept constant. The jet was assumed to be isothermal. The observables to be reproduced were then the impact of the jet on the Balmer series (Hα to Hδ). The equivalent width of the absorption was used as a tracer of the total opacity.

The source function $S_v$ is described by the Planck function $B_v$. As we only needed to find the best fits for two parameters, a grid of jet temperatures and densities were computed, and each of these models was compared to the data. This procedure is described in detail in Bollen et al. (2020) and allowed us to determine the total mass loss from the jet while assuming the bicone is symmetric in both lobes. The found best-fit disk wind mass-ejection rate for the model described in Table 1 is $2.6 \times 10^{-7} M_\odot \text{yr}^{-1}$, with a $1-\sigma$ interval of $[1.3, 21] \times 10^{-7} M_\odot \text{yr}^{-1}$.
<figure><img src="image_4.png" /><figcaption>Fig. 3. Modelling effect on the dynamic spectra. The top panel shows the initial Hα line-profile of HD 52961 which serves as background spectrum. The lower panel shows the synthetic dynamic spectra of the Hα line using the stellar disk wind parametric model.</figcaption></figure>

background spectrum also at phases of conjunction. To use these spectra correctly, it is important to understand the orbital motion of the system. The wavelength shift of the spectra reveals the radial velocity of the post-AGB star around the centre of mass. We used the orbital parameters of Oomen et al. (2018) given in Table 1. Three of the spectra outside the conjunction contribute to a mean spectrum that serves as a new base spectrum. This spectrum is then used to create a synthetic dynamic spectrum by reintroducing the correct wavelength shift for the specific orbital phase of the post-AGB star.

## 2.2. Modelled components

The model consists of two main components: the post-AGB star itself and the biconical jet centred on the companion. The post-AGB star is orders of magnitude brighter than its typical main-sequence low-mass companion, so we assumed that this binary partner will not contribute to the resulting spectra of the system. In optical wavelengths, the same holds true for the circumbinary disk and the secondary accretion disk. The post-AGB primary and the biconical jet orbit each other with certain binary parameters. At some orbital phases, a LOS from the surface of the post-AGB star passes through the bicone. This ray, characterised by the background spectrum, interacts with the material in the jet at the given Doppler velocity, and this forms the basis of this model.

The post-AGB star is represented by a uniform disk perpendicular to the LOS. This disk is divided into small regions via a Fibonacci grid in which every grid point covers an equal surface area. We then determined the LOS of every grid point in the disk.

The disk wind itself is represented by a cone, which is shown in Fig. 4. This cone has its vertex not located at the secondary star but somewhere under the secondary star, depending on the outer launching radius of the disk wind and the opening angle. The disk wind has two parts: one central fast part and an outer slower part with a higher density. A fit angle determines where the fast region stops and the slow region begins. This angle, together with the offset of the cone, corresponds to a certain inner launching radius.

The inner and outer launching radii cannot be chosen freely. To start, the inner launching radius has to be larger than the stellar radius. The outer disk radius must be smaller than the Roche lobe radius of our main-sequence star, as material outside of this lobe will not be bound to the star (Eggleton 1983):

$$ \frac{r_L}{a} = \frac{0.49q^{2/3}}{0.69q^{2/3} + \ln(1 + q^{1/3})} \quad (1) $$

Here, q is the mass ratio between the two stars, and a is the semi-major axis of the primary. To estimate the radius of the main-sequence companion, we used $R_2 = 1.01(M_2/M_\odot)^{0.724}$ (Demircan & Kahraman 1991).

## 2.3. Fitting procedure

For each phase, the global Hα profiles were determined using the jet parameters of inclination, opening angle, velocity structure, and scaled density distribution as well as the radius of the post-AGB star and its orbit. A Markov chain Monte Carlo (MCMC) fitting routine was used to fit these synthetic spectral time series to the observations. We illustrate this in Fig. 3 panel b, where we show the simulated dynamic spectrum compared with the top panel of Fig. 2, which are the observables themselves. This best fit is based on a disk wind model. The parameters of this best fit are given in the right side of Table 1.

## 2.4. Radiative transfer

To determine the actual mass-loss rates, we needed to compute the absolute densities and assumed temperatures instead of using scaling parameters. To do so, the parametric model was followed by a radiative transfer module using a temperature and density fit while all other parameters from the parametric fitting were kept constant. The jet was assumed to be isothermal. The observables to be reproduced were then the impact of the jet on the Balmer series (Hα to Hδ). The equivalent width of the absorption was used as a tracer of the total opacity.

The source function $S_v$ is described by the Planck function $B_v$. As we only needed to find the best fits for two parameters, a grid of jet temperatures and densities were computed, and each of these models was compared to the data. This procedure is described in detail in Bollen et al. (2020) and allowed us to determine the total mass loss from the jet while assuming the bicone is symmetric in both lobes. The found best-fit disk wind mass-ejection rate for the model described in Table 1 is $2.6 \times 10^{-7} M_\odot \text{yr}^{-1}$, with a $1-\sigma$ interval of $[1.3, 21] \times 10^{-7} M_\odot \text{yr}^{-1}$.
<table><caption>**Table 1.** Orbital parameter and best spatiokinematic fit.</caption><thead><tr><th></th><th>Orbital parameters</th><th></th><th>Parametric fit</th></tr></thead><tbody><tr><td>Period [days]</td><td>1279.4 ± 0.23</td><td><i>i</i> [°]</td><td>69.88</td></tr><tr><td>Eccentricity</td><td>0.22 ± 0.00</td><td><i>θ</i><sub>out</sub> [°]</td><td>49.56</td></tr><tr><td><i>T</i><sub>0</sub> [days]</td><td>2456309 ± 112</td><td><i>θ</i><sub>in</sub> [°]</td><td>49.31</td></tr><tr><td><i>ω</i> [°]</td><td>11 ± 73</td><td><i>θ</i><sub>cav</sub> [°]</td><td>25.14</td></tr><tr><td><i>K</i><sub>1</sub> [km s<sup>-1</sup>]</td><td>13.264 ± 0.005</td><td><i>φ</i><sub>tilt</sub> [°]</td><td>-4.65</td></tr><tr><td><i>γ</i> [km s<sup>-1</sup>]</td><td>6.1 ± 0.1</td><td><i>v</i><sub>in</sub> [km s<sup>-1</sup>]</td><td>158</td></tr><tr><td>Semi-major axis (a) [R<sub>☉</sub>]</td><td>586</td><td><i>v</i><sub>out</sub> [km s<sup>-1</sup>]</td><td>37</td></tr><tr><td>Roche-lobe radius (<i>R</i><sub>L</sub>) [R<sub>☉</sub>]</td><td>395</td><td><i>c</i><sub>v</sub></td><td>0.04</td></tr><tr><td></td><td></td><td><i>p</i><sub>v</sub></td><td>2.32</td></tr><tr><td></td><td></td><td><i>p</i><sub>ρ,in</sub></td><td>14.92</td></tr><tr><td></td><td></td><td><i>p</i><sub>ρ,out</sub></td><td>-9.89</td></tr><tr><td></td><td></td><td><i>c</i><sub>τ</sub></td><td>2.29</td></tr><tr><td></td><td></td><td><i>R</i><sub>1</sub> [R<sub>☉</sub>]</td><td>105.4</td></tr><tr><td></td><td></td><td><i>R</i><sub>2</sub> [R<sub>☉</sub>]</td><td>1.04</td></tr><tr><td></td><td></td><td><i>R</i><sub>in</sub> [R<sub>2</sub>]</td><td>8.6</td></tr><tr><td></td><td></td><td><i>R</i><sub>out</sub> [R<sub>L</sub>]</td><td>1.00</td></tr></tbody></table>

**Notes.** Listed in the table are all orbital parameters of the object HD 52961. Here, “Period” is the time of one revolution in days; “Eccentricity” defines the elongation of the orbit; $T_0$ sets the position of the orbit at phase zero; $\omega$ is the argument of periapsis, which defines the orientation of the ellipse; $K_1$ is the highest radial velocity offset from the median; $\gamma$ is the radial velocity of the mass centre. The semi-major axis and the Roche-lobe radius was calculated assuming an inclination of 60°. The second column shows the best-fit parameters of HD 52961 found using the parametric fit described in Bollen et al. (2022). The tabulated parameters are: (i) inclination angle of the binary system; ($\theta_{out}$) jet outer angle; ($\theta_{in}$) jet inner angle; ($\theta_{cav}$) jet cavity angle; ($\phi_{tilt}$) jet tilt; ($v_{in}$) inner jet velocity; ($v_{out}$) jet velocity at the jet edges; ($c_v$) velocity scaling parameter; ($p_v$) exponent for the velocity profile; ($p_{\rho,in}$, $p_{\rho,out}$) exponent for the density profile for the outer and inner region; ($c_\tau$) optical depth scaling parameter; ($R_1$) the radius of the post-AGB star. The term $R_2$ is the radius of the secondary. Underneath the line, the inner and outer launching radius is shown.

## 2.5. Overview results

A full overview of the results obtained for 16 objects can be found in Bollen et al. (2022, and references therein). The main results to highlight are as follows:

- The geometric models are versatile and are able to fit a wide variety of time series.

- The jets are wide (>30°), with a slow and dense outflow component along the sides of the jet cone (1–80 km s⁻¹) and a faster inner jet component (>150 km s⁻¹). The best fits are found when the centre has a very low density.

- The maximal outflow velocity is similar to the escape velocities of main-sequence stars and not similar to those of compact objects, such as white dwarfs. This can be seen directly from the broadness of the Hα absorption feature once the inclination has been determined. The resulting maximum outflow velocities correspond to a main-sequence star and not to a compact object, such as a white dwarf.

- The best-fit mass accretion rates are higher than expected for all of the systems. The current mass loss of the primary is too low to account for this accretion, so it is assumed that the circumbinary disk is the main source of accreted material. The estimated mass accretion rates are high, and using the estimated total mass of the circumbinary disk, we can estimate the lifetime of the circumbinary disk to be on the order of 100–1000 years (Bollen et al. 2020). The jet phenomenon should, thus, be transient. This is in tension with the high number of jets detected in post-AGB binaries (Kluska et al. 2022; Kamath et al. 2015; Manick et al. 2017).

The models used so far are a purely geometric description of the jets. In this work, we report on our implementation of physical, self-similar jet models (Jacquemin-Ide et al. 2019) in the same

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Disk wind model. In the parametric disk wind model, the disk wind is divided into two parts: the inner wind, up to an opening angle of $\alpha_{in}$, and a disk wind, from $\alpha_{in}$ up to $\alpha_{out}$. The position of the point of the cone is set by the combination of the outer launching radius of the disk wind ($d_{out}$) and $\alpha_{out}$.</figcaption></figure>

framework presented above. We evaluate the models on the basis of the synthetic observables that imprint onto the spectroscopic time series.
<table><caption>**Table 1.** Orbital parameter and best spatiokinematic fit.</caption><thead><tr><th></th><th>Orbital parameters</th><th></th><th>Parametric fit</th></tr></thead><tbody><tr><td>Period [days]</td><td>1279.4 ± 0.23</td><td><i>i</i> [°]</td><td>69.88</td></tr><tr><td>Eccentricity</td><td>0.22 ± 0.00</td><td><i>θ</i><sub>out</sub> [°]</td><td>49.56</td></tr><tr><td><i>T</i><sub>0</sub> [days]</td><td>2456309 ± 112</td><td><i>θ</i><sub>in</sub> [°]</td><td>49.31</td></tr><tr><td><i>ω</i> [°]</td><td>11 ± 73</td><td><i>θ</i><sub>cav</sub> [°]</td><td>25.14</td></tr><tr><td><i>K</i><sub>1</sub> [km s<sup>-1</sup>]</td><td>13.264 ± 0.005</td><td><i>φ</i><sub>tilt</sub> [°]</td><td>-4.65</td></tr><tr><td><i>γ</i> [km s<sup>-1</sup>]</td><td>6.1 ± 0.1</td><td><i>v</i><sub>in</sub> [km s<sup>-1</sup>]</td><td>158</td></tr><tr><td>Semi-major axis (a) [<i>R</i><sub>☉</sub>]</td><td>586</td><td><i>v</i><sub>out</sub> [km s<sup>-1</sup>]</td><td>37</td></tr><tr><td>Roche-lobe radius (<i>R</i><sub>L</sub>) [<i>R</i><sub>☉</sub>]</td><td>395</td><td><i>c</i><sub>v</sub></td><td>0.04</td></tr><tr><td></td><td></td><td><i>p</i><sub>v</sub></td><td>2.32</td></tr><tr><td></td><td></td><td><i>p</i><sub>ρ,in</sub></td><td>14.92</td></tr><tr><td></td><td></td><td><i>p</i><sub>ρ,out</sub></td><td>-9.89</td></tr><tr><td></td><td></td><td><i>c</i><sub>τ</sub></td><td>2.29</td></tr><tr><td></td><td></td><td><i>R</i><sub>1</sub> [<i>R</i><sub>☉</sub>]</td><td>105.4</td></tr><tr><td></td><td></td><td><i>R</i><sub>2</sub> [<i>R</i><sub>☉</sub>]</td><td>1.04</td></tr><tr><td></td><td></td><td><i>R</i><sub>in</sub> [<i>R</i><sub>2</sub>]</td><td>8.6</td></tr><tr><td></td><td></td><td><i>R</i><sub>out</sub> [<i>R</i><sub>L</sub>]</td><td>1.00</td></tr></tbody></table>

**Notes.** Listed in the table are all orbital parameters of the object HD 52961. Here, “Period” is the time of one revolution in days; “Eccentricity” defines the elongation of the orbit; $T_0$ sets the position of the orbit at phase zero; $\omega$ is the argument of periapsis, which defines the orientation of the ellipse; $K_1$ is the highest radial velocity offset from the median; $\gamma$ is the radial velocity of the mass centre. The semi-major axis and the Roche-lobe radius was calculated assuming an inclination of 60°. The second column shows the best-fit parameters of HD 52961 found using the parametric fit described in Bollen et al. (2022). The tabulated parameters are: (i) inclination angle of the binary system; ($\theta_{out}$) jet outer angle; ($\theta_{in}$) jet inner angle; ($\theta_{cav}$) jet cavity angle; ($\phi_{tilt}$) jet tilt; ($v_{in}$) inner jet velocity; ($v_{out}$) jet velocity at the jet edges; ($c_v$) velocity scaling parameter; ($p_v$) exponent for the velocity profile; ($p_{\rho,in}$, $p_{\rho,out}$) exponent for the density profile for the outer and inner region; ($c_\tau$) optical depth scaling parameter; ($R_1$) the radius of the post-AGB star. The term $R_2$ is the radius of the secondary. Underneath the line, the inner and outer launching radius is shown.

## 2.5. Overview results

A full overview of the results obtained for 16 objects can be found in Bollen et al. (2022, and references therein). The main results to highlight are as follows:

- The geometric models are versatile and are able to fit a wide variety of time series.

- The jets are wide (>30°), with a slow and dense outflow component along the sides of the jet cone (1–80 km s⁻¹) and a faster inner jet component (>150 km s⁻¹). The best fits are found when the centre has a very low density.

- The maximal outflow velocity is similar to the escape velocities of main-sequence stars and not similar to those of compact objects, such as white dwarfs. This can be seen directly from the broadness of the Hα absorption feature once the inclination has been determined. The resulting maximum outflow velocities correspond to a main-sequence star and not to a compact object, such as a white dwarf.

- The best-fit mass accretion rates are higher than expected for all of the systems. The current mass loss of the primary is too low to account for this accretion, so it is assumed that the circumbinary disk is the main source of accreted material. The estimated mass accretion rates are high, and using the estimated total mass of the circumbinary disk, we can estimate the lifetime of the circumbinary disk to be on the order of 100–1000 years (Bollen et al. 2020). The jet phenomenon should, thus, be transient. This is in tension with the high number of jets detected in post-AGB binaries (Kluska et al. 2022; Kamath et al. 2015; Manick et al. 2017).

The models used so far are a purely geometric description of the jets. In this work, we report on our implementation of physical, self-similar jet models (Jacquemin-Ide et al. 2019) in the same

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Disk wind model. In the parametric disk wind model, the disk wind is divided into two parts: the inner wind, up to an opening angle of $\alpha_{in}$, and a disk wind, from $\alpha_{in}$ up to $\alpha_{out}$. The position of the point of the cone is set by the combination of the outer launching radius of the disk wind ($d_{out}$) and $\alpha_{out}$.</figcaption></figure>

framework presented above. We evaluate the models on the basis of the synthetic observables that imprint onto the spectroscopic time series.
# 3. Inclusion of a physical jet model

## 3.1. Accretion-ejection model

We wish to use a magnetically driven jet model that is self-consistent with the underlying accretion disk physics. In this accretion-ejection model, a large-scale vertical magnetic field $B_z$ is assumed to thread the disk. An MHD turbulence is then triggered within the disk, which provides an anomalous (turbulent) viscosity but also magnetic diffusivity, thereby allowing accretion to proceed (see e.g. Balbus 2003, and references therein). However, the existence of a vertical field also leads to the launching of jets, which provide another means of disk angular momentum removal (Blandford & Payne 1982). Semi-analytical models are possible thanks to the use of a self-similar Ansatz where all axisymmetric and steady-state quantities are assumed to be power laws of the cylindrical radius (see e.g. Ferreira & Pelletier 1995; Ferreira 1997). This mathematical framework allows for an exact variable separation and transforms the full set of MHD steady-state partial differential equations into ordinary equations on a self-similar variable $x = z/r$. The turbulent disk of isothermal scale height $h(r) = \epsilon r$, where the disk aspect ratio $\epsilon$ is used as a free constant, is then described within a resistive MHD framework. The disk is in a magneto-hydrostatic vertical balance that is slightly sub-Keplerian accretes mostly due to the dominant jet torque. The steady-state is then possible only if the assumed turbulence allows plasma to slip through the magnetic field within the disk. Above the disk, turbulence is expected to be quenched, and the ejected plasma enters an ideal MHD regime. Since bipolar jets are emitted, the disk accretion rate is a function of the radius, $M_a(r) \propto r^{\xi}$, where $\xi$ is the local disk ejection efficiency. The integration starts from the disk equatorial plane at $z = 0$ in the midst of the resistive MHD zone, and it proceeds in the vertical direction until the flow becomes an ideal MHD, achieves a super slow magnetosonic speed (SM critical point) above the disk surface then a super-Alfvénic poloidal speed (A critical point), and finally obtains a (conventional) super-fast magnetosonic speed. We refer the interested reader to the work of Jacquemin-Ide et al. (2019) and references therein, which is the last published paper of a series on these exact MHD accretion-ejection solutions.

In this paper, we do not explore the parameter space of such MHD solutions and instead use one particular solution that has been already used in the context of Young Stellar Objects (YSOs) by Panoglou et al. (2012). Although it has been shown to provide slightly too fast outflows for YSOs, it has been well documented and explains most YSO jet features. Moreover, it provides a suitable reference for our post-AGB binary system since the companion object in our case is also a solar mass star. Such an MHD solution is described by the following set of parameters: a disk aspect ratio $\epsilon = 0.03$, a disk ejection efficiency $\xi = 0.04$, MHD turbulence level $\alpha_m = 2$, magnetic Prandtl number $P_m = 1$, MHD turbulence anisotropy $\chi_m = 1.42$. We also assumed some additional heating at the disk surface so that outflows are actually magneto-thermal with a parameter $f = 0.0066$, corresponding to 0.66% of the released power that goes into heating Casse & Ferreira (2000a). As a consequence of a large disk magnetisation $\mu = 0.351$, accretion is super-sonic within the disk (with a sonic Mach number $m_s = 1.3$), jets carry away nearly all the angular momentum of the underlying disk, and the total jet power is about 88% of the released accretion power, with only ~12% of that power being radiated by the disk. We note that this energy budget applies only to the portion of the disk (located between radii $r_{in}$ and $r_{out}$) that launches the MHD jets. For any radius $r_o$ inside this range, the rotation

velocity, density, and magnetic field at the disk midplane are respectively:

$$ \Omega_o r_o = \delta \Omega_{K_o} r_o = 29.8 \delta m^{1/2} \left(\frac{r_o}{1 \text{ AU}}\right)^{-1/2} \text{km s}^{-1}, \quad (2) $$

$$ n_o = \frac{4.7}{m_s} 10^{-14} m^{-1/2} \dot{m}_{-5} \left(\frac{\epsilon}{0.01}\right)^{-2} \left(\frac{r_o}{1 \text{ AU}}\right)^{\xi-3/2} \left(\frac{r_{in}}{1 \text{ AU}}\right)^{-\xi} \text{cm}^{-3}, \quad (3) $$

$$ B_z = 2.9 m^{1/4} \dot{m}_{-5}^{1/2} \left(\frac{\mu}{m_s}\right)^{1/2} \left(\frac{r_o}{1 \text{ AU}}\right)^{-5/4+\xi/2} \left(\frac{r_{in}}{1 \text{ AU}}\right)^{-\xi/2} \text{G}, \quad (4) $$

where $m = M_2/M_\odot$, $\dot{m}_{-5} = \dot{M}_a(r_{in})/(10^{-5} M_\odot \text{yr}^{-1})$. For the solution used here, the Keplerian deviation was found to be $\delta = 0.9$. At the Alfvén critical point (located at $r_A/r_o = 3.8$, $z_A/r_o = 3.1$), the flow has already achieved a poloidal velocity that is 2.1 times the Keplerian velocity ($\Omega_{K_o} r_o$) at the anchoring radius $r_o$ and 1.9 times the local rotational speed. Asymptotically, the jet speed almost reaches its maximum velocity $v_{jet} = \sqrt{2\lambda - 3\Omega_{K_o} r_o} = 4.8\Omega_{K_o} r_o$, with a magnetic lever arm $\lambda = 1 + \frac{1}{2\xi} \approx 13.1$.

## 3.2. From line to 3D grid

Originally these models have mostly been used in single star objects. As a result, the models have been stored in such a way that allows for the direct calculation of the required values along a LOS towards the central object of the accretion. However, it is possible to combine the results from all possible angles and compute the density and velocity field in a grid as long as the points in the grid fall within the magnetic field lines piercing the edges of the accretion disk. This allowed us to compute along a new LOS that is not towards the central accretor. This is required, as the disk wind is backlit by the actual object we observe, namely, the luminous post-AGB star. The shape, velocity, and density that are typical for these extended disk winds are shown in Figs. 5–7. By converting this grid to a Cartesian grid, the determination of the position of each point of the jet becomes less computationally intensive. Figures 6 and 7 show the Cartesian integrated velocity in the three dimensions over the jet, while Fig. 5 shows the density of the jet.

From the 2D jet model described above, we produced a 3D grid by axial symmetry with the z-axis. As the jet follows the secondary, the motion of this secondary influences the velocity field of the jet with respect to the observer. The velocity in the models is given in three directions (radial velocity from the centre of the secondary star $v_R$, angular velocity $v_\phi$, and vertical velocity $v_z$). We converted the cylindrical velocities to Cartesian, as the LOS is already set in Cartesian coordinates, using the following conversion:

$$ v_x = v_R \cdot \cos\phi + v_\phi \cdot \cos(90^\circ + \phi) \quad (5) $$

$$ v_y = v_R \cdot \sin\phi + v_\phi \cdot \sin(90^\circ + \phi) \quad (6) $$

$$ v_z = v_z. \quad (7) $$

Figure 5 shows the new geometry that is centred on the post-AGB star instead of the main-sequence star. Here, a typical LOS is also shown with a black line. Outside our computed jet, the density is just set to zero, thereby neglecting any other wind source (post-AGB wind as well as any wind from the companion star).

# 4. Parameter study

In our parameter study, we studied the impact of the physical parameters of the jet model on our synthetic observables.
# 3. Inclusion of a physical jet model

## 3.1. Accretion-ejection model

We wish to use a magnetically driven jet model that is self-consistent with the underlying accretion disk physics. In this accretion-ejection model, a large-scale vertical magnetic field $B_z$ is assumed to thread the disk. An MHD turbulence is then triggered within the disk, which provides an anomalous (turbulent) viscosity but also magnetic diffusivity, thereby allowing accretion to proceed (see e.g. Balbus 2003, and references therein). However, the existence of a vertical field also leads to the launching of jets, which provide another means of disk angular momentum removal (Blandford & Payne 1982). Semi-analytical models are possible thanks to the use of a self-similar Ansatz where all axisymmetric and steady-state quantities are assumed to be power laws of the cylindrical radius (see e.g. Ferreira & Pelletier 1995; Ferreira 1997). This mathematical framework allows for an exact variable separation and transforms the full set of MHD steady-state partial differential equations into ordinary equations on a self-similar variable $x = z/r$. The turbulent disk of isothermal scale height $h(r) = \epsilon r$, where the disk aspect ratio $\epsilon$ is used as a free constant, is then described within a resistive MHD framework. The disk is in a magneto-hydrostatic vertical balance that is slightly sub-Keplerian accretes mostly due to the dominant jet torque. The steady-state is then possible only if the assumed turbulence allows plasma to slip through the magnetic field within the disk. Above the disk, turbulence is expected to be quenched, and the ejected plasma enters an ideal MHD regime. Since bipolar jets are emitted, the disk accretion rate is a function of the radius, $M_a(r) \propto r^{\xi}$, where $\xi$ is the local disk ejection efficiency. The integration starts from the disk equatorial plane at $z = 0$ in the midst of the resistive MHD zone, and it proceeds in the vertical direction until the flow becomes an ideal MHD, achieves a super slow magnetosonic speed (SM critical point) above the disk surface then a super-Alfvénic poloidal speed (A critical point), and finally obtains a (conventional) super-fast magnetosonic speed. We refer the interested reader to the work of Jacquemin-Ide et al. (2019) and references therein, which is the last published paper of a series on these exact MHD accretion-ejection solutions.

In this paper, we do not explore the parameter space of such MHD solutions and instead use one particular solution that has been already used in the context of Young Stellar Objects (YSOs) by Panoglou et al. (2012). Although it has been shown to provide slightly too fast outflows for YSOs, it has been well documented and explains most YSO jet features. Moreover, it provides a suitable reference for our post-AGB binary system since the companion object in our case is also a solar mass star. Such an MHD solution is described by the following set of parameters: a disk aspect ratio $\epsilon = 0.03$, a disk ejection efficiency $\xi = 0.04$, MHD turbulence level $\alpha_m = 2$, magnetic Prandtl number $P_m = 1$, MHD turbulence anisotropy $\chi_m = 1.42$. We also assumed some additional heating at the disk surface so that outflows are actually magneto-thermal with a parameter $f = 0.0066$, corresponding to 0.66% of the released power that goes into heating Casse & Ferreira (2000a). As a consequence of a large disk magnetisation $\mu = 0.351$, accretion is super-sonic within the disk (with a sonic Mach number $m_s = 1.3$), jets carry away nearly all the angular momentum of the underlying disk, and the total jet power is about 88% of the released accretion power, with only ~12% of that power being radiated by the disk. We note that this energy budget applies only to the portion of the disk (located between radii $r_{in}$ and $r_{out}$) that launches the MHD jets. For any radius $r_o$ inside this range, the rotation

velocity, density, and magnetic field at the disk midplane are respectively:

$$ \Omega_o r_o = \delta \Omega_{K_o} r_o = 29.8 \delta m^{1/2} \left(\frac{r_o}{1 \text{ AU}}\right)^{-1/2} \text{km s}^{-1}, \quad (2) $$

$$ n_o = \frac{4.7}{m_s} 10^{-14} m^{-1/2} \dot{m}_{-5} \left(\frac{\epsilon}{0.01}\right)^{-2} \left(\frac{r_o}{1 \text{ AU}}\right)^{\xi-3/2} \left(\frac{r_{\text{in}}}{1 \text{ AU}}\right)^{-\xi} \text{cm}^{-3}, \quad (3) $$

$$ B_z = 2.9 m^{1/4} \dot{m}_{-5}^{1/2} \left(\frac{\mu}{m_s}\right)^{1/2} \left(\frac{r_o}{1 \text{ AU}}\right)^{-5/4+\xi/2} \left(\frac{r_{\text{in}}}{1 \text{ AU}}\right)^{-\xi/2} \text{G}, \quad (4) $$

where $m = M_2/M_\odot$, $\dot{m}_{-5} = \dot{M}_a(r_{\text{in}})/(10^{-5} M_\odot \text{yr}^{-1})$. For the solution used here, the Keplerian deviation was found to be $\delta = 0.9$. At the Alfvén critical point (located at $r_A/r_o = 3.8$, $z_A/r_o = 3.1$), the flow has already achieved a poloidal velocity that is 2.1 times the Keplerian velocity ($\Omega_{K_o} r_o$) at the anchoring radius $r_o$ and 1.9 times the local rotational speed. Asymptotically, the jet speed almost reaches its maximum velocity $v_{\text{jet}} = \sqrt{2\lambda - 3\Omega_{K_o} r_o} = 4.8\Omega_{K_o} r_o$, with a magnetic lever arm $\lambda = 1 + \frac{1}{2\xi} \approx 13.1$.

## 3.2. From line to 3D grid

Originally these models have mostly been used in single star objects. As a result, the models have been stored in such a way that allows for the direct calculation of the required values along a LOS towards the central object of the accretion. However, it is possible to combine the results from all possible angles and compute the density and velocity field in a grid as long as the points in the grid fall within the magnetic field lines piercing the edges of the accretion disk. This allowed us to compute along a new LOS that is not towards the central accretor. This is required, as the disk wind is backlit by the actual object we observe, namely, the luminous post-AGB star. The shape, velocity, and density that are typical for these extended disk winds are shown in Figs. 5–7. By converting this grid to a Cartesian grid, the determination of the position of each point of the jet becomes less computationally intensive. Figures 6 and 7 show the Cartesian integrated velocity in the three dimensions over the jet, while Fig. 5 shows the density of the jet.

From the 2D jet model described above, we produced a 3D grid by axial symmetry with the z-axis. As the jet follows the secondary, the motion of this secondary influences the velocity field of the jet with respect to the observer. The velocity in the models is given in three directions (radial velocity from the centre of the secondary star $v_R$, angular velocity $v_\phi$, and vertical velocity $v_z$). We converted the cylindrical velocities to Cartesian, as the LOS is already set in Cartesian coordinates, using the following conversion:

$$ v_x = v_R \cdot \cos \phi + v_\phi \cdot \cos(90^\circ + \phi) \quad (5) $$

$$ v_y = v_R \cdot \sin \phi + v_\phi \cdot \sin(90^\circ + \phi) \quad (6) $$

$$ v_z = v_z. \quad (7) $$

Figure 5 shows the new geometry that is centred on the post-AGB star instead of the main-sequence star. Here, a typical LOS is also shown with a black line. Outside our computed jet, the density is just set to zero, thereby neglecting any other wind source (post-AGB wind as well as any wind from the companion star).

# 4. Parameter study

In our parameter study, we studied the impact of the physical parameters of the jet model on our synthetic observables.
<figure><img src="image_7.png" /></figure>

**Fig. 5.** MHD disk wind density with LOS. Schematic of the central region of a post-AGB binary showing the post-AGB star in yellow. The jet is shown as starting from the accretion disk with the densities plotted in the colour plot. The black line shows how a hypothetical LOS would pass through the jet. For this figure as well as Figs. 6 and 7, the wind values are $\dot{M}(r_{\text{in}}) = 10^{-4} M_{\odot} \text{yr}^{-1}$, $r_{\text{in}} = 20 R_{2}$, and $r_{\text{out}} = 0.5 R_{\text{L}}$ with an inclination of 60°.

<figure><img src="image_8.png" /></figure>

**Fig. 6.** Rotational velocity of MHD disk wind. The figure shows an example of the rotational velocity in the jet ($v_{\phi}$); the direction of this velocity is into the page. Up to the Alfvén surface (green line), the magnetic field almost enforces isorotation to the plasma despite the opening of the flow. The decrease of the rotation velocity occurs only beyond the Alfvén surface. The wind values for this plot are the same as in Fig. 5.

We used the synthetic dynamic $H\alpha$ profiles of the binary post-AGB star HD 52961 for illustration. We utilised the orbital parameters from Oomen et al. (2018, Table 1). We limited ourselves to the four parameters that impact the observables of a given model, which are the inner and outer launching radius, the inclination, and the mass accretion rate.

We first defined the base model that would be modified by changing the parameters mentioned above. The dynamic synthetic spectrum of this base model is shown in Fig. 8. Since our goal was only to assess the first implications of using a physical disk wind model instead of a parametric one, we used only one MHD model, the parameters of which have been specified in Sect. 2.2. The quantities that needed to be specified were (1) the innermost jet launching radius $r_{\text{in}}$, (2) the outermost jet launching radius $r_{\text{out}}$, (3) the disk accretion rate $\dot{M}_{\text{a}}(r_{\text{in}})$ at $r_{\text{in}}$, and (4) the LOS inclination angle $i$ (0° is for pole-on objects, 90° for edge-on). The mass of the secondary was computed from

<figure><img src="image_9.png" /></figure>

**Fig. 7.** Velocity distribution of MHD disk wind. The figure shows an example of the velocity of a jet in the plane ($v_{\text{R}}$ and $v_{\text{z}}$). The wind values for this plot are the same as in Fig. 5.

<figure><img src="image_10.png" /></figure>

**Fig. 8.** Base model synthetic spectra. The figure shows the synthetic dynamic spectrum using the base model. This model has an ejection efficiency of $\xi = 0.04$, a magnetisation of $\mu = 0.351$, a scale height of $\epsilon = 0.03$, and the temperature of 5000 K. All the other relevant parameters can be found in Tables 1 and 2.

the mass function and the inclination assuming a post-AGB star mass of $0.6 M_{\odot}$. Additionally, the jet was set at a temperature of 5000 K based on the fact that the jet feature is in absorption (see the fit from Bollen 2021). The number of free parameters as a result has been reduced greatly from 12 in the parametric models to four in these new models, with all four of these models having direct links to the physical behaviour of the jet or how we observed it.

We defined a base model with an inclination of 60°, a mass accretion rate of $\dot{M}_{\text{a}}(r_{\text{in}}) = 10^{-4} M_{\odot} \text{yr}^{-1}$, an inner launching radius of $20 R_{2} = r_{\text{in}} = 0.09$ AU, and an outer launching radius of $r_{\text{out}} = 0.9 R_{\text{L}} = 1.6$ AU. As the inclination sets the mass of the secondary, it also influences the stellar radius and the Roche-lobe radius. The choice of the orbital parameters as well as the disk parameters of the base model are given in Table 1. The choice of these parameters is based on the results obtained from the study of this object using the parametric jet model used in Bollen et al. (2022) and recalled in Table 1. In the spatiokinematic study, a high inclination as well as a very wide jet were found, which we induced by choosing a large outer launching radius. For the inner jet launching radius $r_{\text{in}}$, we used instead a constraint on the maximum observable velocity, which we obtained along the
<figure><img src="image_6.png" /></figure>

**Fig. 5.** MHD disk wind density with LOS. Schematic of the central region of a post-AGB binary showing the post-AGB star in yellow. The jet is shown as starting from the accretion disk with the densities plotted in the colour plot. The black line shows how a hypothetical LOS would pass through the jet. For this figure as well as Figs. 6 and 7, the wind values are $\dot{M}(r_{\text{in}}) = 10^{-4} M_{\odot} \text{yr}^{-1}$, $r_{\text{in}} = 20 R_{2}$, and $r_{\text{out}} = 0.5 R_{\text{L}}$ with an inclination of 60°.

<figure><img src="image_7.png" /></figure>

**Fig. 6.** Rotational velocity of MHD disk wind. The figure shows an example of the rotational velocity in the jet ($v_{\phi}$); the direction of this velocity is into the page. Up to the Alfvén surface (green line), the magnetic field almost enforces isorotation to the plasma despite the opening of the flow. The decrease of the rotation velocity occurs only beyond the Alfvén surface. The wind values for this plot are the same as in Fig. 5.

We used the synthetic dynamic $H\alpha$ profiles of the binary post-AGB star HD 52961 for illustration. We utilised the orbital parameters from Oomen et al. (2018, Table 1). We limited ourselves to the four parameters that impact the observables of a given model, which are the inner and outer launching radius, the inclination, and the mass accretion rate.

We first defined the base model that would be modified by changing the parameters mentioned above. The dynamic synthetic spectrum of this base model is shown in Fig. 8. Since our goal was only to assess the first implications of using a physical disk wind model instead of a parametric one, we used only one MHD model, the parameters of which have been specified in Sect. 2.2. The quantities that needed to be specified were (1) the innermost jet launching radius $r_{\text{in}}$, (2) the outermost jet launching radius $r_{\text{out}}$, (3) the disk accretion rate $\dot{M}_{\text{a}}(r_{\text{in}})$ at $r_{\text{in}}$, and (4) the LOS inclination angle $i$ (0° is for pole-on objects, 90° for edge-on). The mass of the secondary was computed from

<figure><img src="image_8.png" /></figure>

**Fig. 7.** Velocity distribution of MHD disk wind. The figure shows an example of the velocity of a jet in the plane ($v_{\text{R}}$ and $v_{\text{z}}$). The wind values for this plot are the same as in Fig. 5.

<figure><img src="image_9.png" /></figure>

**Fig. 8.** Base model synthetic spectra. The figure shows the synthetic dynamic spectrum using the base model. This model has an ejection efficiency of $\xi = 0.04$, a magnetisation of $\mu = 0.351$, a scale height of $\epsilon = 0.03$, and the temperature of 5000 K. All the other relevant parameters can be found in Tables 1 and 2.

the mass function and the inclination assuming a post-AGB star mass of $0.6 M_{\odot}$. Additionally, the jet was set at a temperature of 5000 K based on the fact that the jet feature is in absorption (see the fit from Bollen 2021). The number of free parameters as a result has been reduced greatly from 12 in the parametric models to four in these new models, with all four of these models having direct links to the physical behaviour of the jet or how we observed it.

We defined a base model with an inclination of 60°, a mass accretion rate of $\dot{M}_{\text{a}}(r_{\text{in}}) = 10^{-4} M_{\odot} \text{yr}^{-1}$, an inner launching radius of $20 R_{2} = r_{\text{in}} = 0.09$ AU, and an outer launching radius of $r_{\text{out}} = 0.9 R_{\text{L}} = 1.6$ AU. As the inclination sets the mass of the secondary, it also influences the stellar radius and the Roche-lobe radius. The choice of the orbital parameters as well as the disk parameters of the base model are given in Table 1. The choice of these parameters is based on the results obtained from the study of this object using the parametric jet model used in Bollen et al. (2022) and recalled in Table 1. In the spatiokinematic study, a high inclination as well as a very wide jet were found, which we induced by choosing a large outer launching radius. For the inner jet launching radius $r_{\text{in}}$, we used instead a constraint on the maximum observable velocity, which we obtained along the
<table><caption>**Table 2.** All the different parameters that deviate from the base model.</caption><thead><tr><th></th><th colspan="6">Base model</th></tr></thead><tbody><tr><td>Inclination [°]</td><td>60</td><td>18</td><td>36</td><td>54</td><td>60</td><td>72</td></tr><tr><td>Mass accretion rate [M<sub>☉</sub> yr<sup>-1</sup>]</td><td>10<sup>-4</sup></td><td>10<sup>-3</sup></td><td>10<sup>-4</sup></td><td>10<sup>-5</sup></td><td>10<sup>-6</sup></td><td>10<sup>-7</sup></td></tr><tr><td>Outer launching radius [R<sub>L</sub>]</td><td>0.9</td><td>0.5</td><td>0.6</td><td>0.7</td><td>0.8</td><td>0.9</td></tr><tr><td>Inner launching radius [R<sub>2</sub>]</td><td>20</td><td>1</td><td>11</td><td>21</td><td>31</td><td>41</td></tr></tbody></table>

**Notes.** The base model column shows the base model from where parameters where changed. The orbital parameters of this model are shown in Table 1.

innermost streamline. Since it is observed to be on the order of 100 km s⁻¹ (Fig. 2), we chose $r_{in} = 20 R_2$ where the rotation velocity is comparable. Nonetheless, even if the maximal jet velocity reaches $v_{jet} = 468$ km s⁻¹, the densest zone close to the disk probably dominates absorption and is associated only with poloidal projected velocities of 130–150 km s⁻¹, in agreement with observations.

We subsequently undertook a parametric study of 18 more models by varying the inclination angle $i$, disk accretion rate, and the inner and outer jet launching radii around the mentioned fiducial values. The inner and outer launching radii cannot be chosen completely freely. As before, $r_{in}$ has to be larger than the stellar radius, and $r_{out}$ must be smaller than the Roche lobe radius ($R_L$). The ranges of parameters we explored are shown in Table 2.

# 5. Results

In this section, we first describe the results obtained with our base model. After that, we analyse the effects of the different parameters on the synthetic dynamic spectra.

## 5.1. Base model

The base model (see Fig. 8) shows a broad absorption feature similar to what is observed in Fig. 2. The absorption feature is present over about 35% of the phase starting at a phase of zero. The velocity range of the base model reaches values up to 300 km s⁻¹, which is a lot higher than the observations show. The most remarkable difference is, however, that the synthetic spectra also show a smaller but redshifted feature at an orbital phase of about 0.3. This feature is less wide, only visible over 10% of the phase, and reaches velocities of only 100 km s⁻¹. This redshifted absorption is a result of the rotation being included in these models. At one point in the orbit, the LOS goes through one side of the disk wind that is rotating towards the post-AGB star. This causes the absorption feature to be asymmetric as well. There is also a lack of absorption in the middle of the broad absorption feature, which is also not observed. This is a result of not adding a minimum density of the jet and instead choosing the region outside of the jet to be void of any material.

## 5.2. Inclination

We modelled five different inclinations ranging from almost pole-on ($i = 18°$) to almost edge-on ($i = 72°$). The inclination is naturally limited in these systems, as at very high inclinations ($i > 70°$), the circumbinary disk obscures the post-AGB star, at which point, no direct light from the post-AGB primary will be seen anymore.

The resulting dynamic spectra are shown in Fig. 9. The inclination at which we observed the system strongly influences the resulting absorption feature, as it determines through which part of the jet the LOS passes. At lower inclinations, the LOS passes through a higher part of the jet, further from the launching region. This results in a weaker absorption feature, as the density decreases with altitude (z-direction). As the jet gets wider at higher altitudes, the phase-coverage of the jet becomes broader.

Additionally, at low inclinations, the rotational component of the velocity in the jet is closer to the perpendicular direction of the LOS, and the impact of this rotation to the observables disappears. Finally, the rotational velocity is also strongest when closest to the accreting object, and by looking through higher angles, the rotational velocity decreases (Fig. 6).

## 5.3. Disk accretion rate

The dynamic spectra in Fig. 10 show that as the mass accretion is lowered, the absorption feature becomes less strong, with no difference in the other characteristics of the dynamic spectrum. The accretion rate directly influences the density of the jet, as can be seen in Eq. (3). This is to be expected, as the disk accretion rate introduces only a scaling factor in the disk density, thereby on the ejected material. Since the disk ejection efficiency, $\xi = 0.04$, is small, the (one-sided) mass loss in the wind is $\dot{M}_j = (\xi/2) \ln(r_{out}/r_{in}) \dot{M}_a(r_{in}) = 5.8 \times 10^{-6} M_{\odot} \text{yr}^{-1}$ for the reference model. For the absorption feature to be noticeable, the mass accretion rate needs to be around $10^{-4} M_{\odot} \text{yr}^{-1}$ in the base-model setup assumed temperature.

## 5.4. Launching region

The inner and outer launching radius of the disk was changed independently for this parameter study (Figs. 11 and 12). Within our self-similar framework, the jet becomes wider as the launching region within the disk gets larger. As a consequence, changing the outer jet launching radius $r_{out}$ (while keeping $r_{in}$ constant) has a strong impact on the width of the jet and therefore on the phase coverage of the absorption feature (see Fig. 11). Increasing $r_{out}$ appears somewhat degenerate with lowering the inclination angle in the synthetic dynamic spectra. However, and contrary to the effect of inclination, not only does the maximum velocity remain untouched when increasing $r_{out}$, but the absorption feature also becomes slightly stronger (especially at low velocities). Figure 12 shows the influence of the innermost jet launching radius $r_{in}$ on the dynamic spectra while keeping $r_{out} = 0.9 R_L$. With no surprise, the most salient effect of $r_{in}$ is to change the maximum velocity seen in the absorption feature. A low inner launching region allows the fastest part of the jet closest to the poles to be included in the jet model, which results in a larger
<table><caption>**Table 2.** All the different parameters that deviate from the base model.</caption><thead><tr><th colspan="7">Base model</th></tr></thead><tbody><tr><td>Inclination [°]</td><td>60</td><td>18</td><td>36</td><td>54</td><td>60</td><td>72</td></tr><tr><td>Mass accretion rate [M<sub>☉</sub> yr<sup>-1</sup>]</td><td>10<sup>-4</sup></td><td>10<sup>-3</sup></td><td>10<sup>-4</sup></td><td>10<sup>-5</sup></td><td>10<sup>-6</sup></td><td>10<sup>-7</sup></td></tr><tr><td>Outer launching radius [R<sub>L</sub>]</td><td>0.9</td><td>0.5</td><td>0.6</td><td>0.7</td><td>0.8</td><td>0.9</td><td>1.0</td></tr><tr><td>Inner launching radius [R<sub>2</sub>]</td><td>20</td><td>1</td><td>11</td><td>21</td><td>31</td><td>41</td></tr></tbody></table>

**Notes.** The base model column shows the base model from where parameters where changed. The orbital parameters of this model are shown in Table 1.

innermost streamline. Since it is observed to be on the order of 100 km s⁻¹ (Fig. 2), we chose $r_{in} = 20 R_2$ where the rotation velocity is comparable. Nonetheless, even if the maximal jet velocity reaches $v_{jet} = 468$ km s⁻¹, the densest zone close to the disk probably dominates absorption and is associated only with poloidal projected velocities of 130–150 km s⁻¹, in agreement with observations.

We subsequently undertook a parametric study of 18 more models by varying the inclination angle $i$, disk accretion rate, and the inner and outer jet launching radii around the mentioned fiducial values. The inner and outer launching radii cannot be chosen completely freely. As before, $r_{in}$ has to be larger than the stellar radius, and $r_{out}$ must be smaller than the Roche lobe radius ($R_L$). The ranges of parameters we explored are shown in Table 2.

# 5. Results

In this section, we first describe the results obtained with our base model. After that, we analyse the effects of the different parameters on the synthetic dynamic spectra.

## 5.1. Base model

The base model (see Fig. 8) shows a broad absorption feature similar to what is observed in Fig. 2. The absorption feature is present over about 35% of the phase starting at a phase of zero. The velocity range of the base model reaches values up to 300 km s⁻¹, which is a lot higher than the observations show. The most remarkable difference is, however, that the synthetic spectra also show a smaller but redshifted feature at an orbital phase of about 0.3. This feature is less wide, only visible over 10% of the phase, and reaches velocities of only 100 km s⁻¹. This redshifted absorption is a result of the rotation being included in these models. At one point in the orbit, the LOS goes through one side of the disk wind that is rotating towards the post-AGB star. This causes the absorption feature to be asymmetric as well. There is also a lack of absorption in the middle of the broad absorption feature, which is also not observed. This is a result of not adding a minimum density of the jet and instead choosing the region outside of the jet to be void of any material.

## 5.2. Inclination

We modelled five different inclinations ranging from almost pole-on ($i = 18°$) to almost edge-on ($i = 72°$). The inclination is naturally limited in these systems, as at very high inclinations ($i > 70°$), the circumbinary disk obscures the post-AGB star, at which point, no direct light from the post-AGB primary will be seen anymore.

The resulting dynamic spectra are shown in Fig. 9. The inclination at which we observed the system strongly influences the resulting absorption feature, as it determines through which part of the jet the LOS passes. At lower inclinations, the LOS passes through a higher part of the jet, further from the launching region. This results in a weaker absorption feature, as the density decreases with altitude (z-direction). As the jet gets wider at higher altitudes, the phase-coverage of the jet becomes broader.

Additionally, at low inclinations, the rotational component of the velocity in the jet is closer to the perpendicular direction of the LOS, and the impact of this rotation to the observables disappears. Finally, the rotational velocity is also strongest when closest to the accreting object, and by looking through higher angles, the rotational velocity decreases (Fig. 6).

## 5.3. Disk accretion rate

The dynamic spectra in Fig. 10 show that as the mass accretion is lowered, the absorption feature becomes less strong, with no difference in the other characteristics of the dynamic spectrum. The accretion rate directly influences the density of the jet, as can be seen in Eq. (3). This is to be expected, as the disk accretion rate introduces only a scaling factor in the disk density, thereby on the ejected material. Since the disk ejection efficiency, $\xi = 0.04$, is small, the (one-sided) mass loss in the wind is $\dot{M}_j = (\xi/2) \ln(r_{out}/r_{in}) \dot{M}_a(r_{in}) = 5.8 \times 10^{-6} M_{\odot} \text{yr}^{-1}$ for the reference model. For the absorption feature to be noticeable, the mass accretion rate needs to be around $10^{-4} M_{\odot} \text{yr}^{-1}$ in the base-model setup assumed temperature.

## 5.4. Launching region

The inner and outer launching radius of the disk was changed independently for this parameter study (Figs. 11 and 12). Within our self-similar framework, the jet becomes wider as the launching region within the disk gets larger. As a consequence, changing the outer jet launching radius $r_{out}$ (while keeping $r_{in}$ constant) has a strong impact on the width of the jet and therefore on the phase coverage of the absorption feature (see Fig. 11). Increasing $r_{out}$ appears somewhat degenerate with lowering the inclination angle in the synthetic dynamic spectra. However, and contrary to the effect of inclination, not only does the maximum velocity remain untouched when increasing $r_{out}$, but the absorption feature also becomes slightly stronger (especially at low velocities). Figure 12 shows the influence of the innermost jet launching radius $r_{in}$ on the dynamic spectra while keeping $r_{out} = 0.9 R_L$. With no surprise, the most salient effect of $r_{in}$ is to change the maximum velocity seen in the absorption feature. A low inner launching region allows the fastest part of the jet closest to the poles to be included in the jet model, which results in a larger
<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Inclination effect. The dynamical spectra in this list all utilise the orbital parameters and accretion disk parameters in Table 1 while changing the inclination according to Table 2. From top left to bottom right, the inclinations used are 18°, 36°, 54°, and 72°.</figcaption></figure>

velocity coverage at the higher end of the velocities. On the contrary, if the inner jet launching region becomes large, the central empty core of the jet increases, and the lack of absorption present at the middle of the phase coverage becomes larger. This occurs because our model keeps the medium void where there is no wind. By adding a different and additional outflow component on the symmetry axis of the bipolar jet (e.g. a stellar wind), this gap in velocity space in the synthetic spectra can be removed, but this is not included in our current model.

## 5.5. Disk versus jet rotation

To confirm that the redshifted absorption is due to rotation, we include Fig. 13. Here, a dynamic spectra based on a self-similar jet model for which the rotation has been set to zero and all other parameters are the same as in the base model (Fig. 8) is shown. The resulting dynamic spectra indeed appear very different compared to the base model: the redshifted absorption is no longer present, and the blueshifted part is more symmetrical. This clearly shows the importance of the rotation velocity in the production of the redshifted absorption.

All models with rotation produce a redshifted absorption in the spectra, something that is not observed in the data nor in the parametric models (as they do not include rotation). To understand what part of the wind contributes to the redshifted absorption, we divided the wind in two regions: the sub-Alfvénic region and the super-Alfvénic region. The sub-Alfvénic region includes the disk region and the starting wind, while the super-Alfvénic

region is the launched wind. We created two synthetic spectra where only the sub-Alfvénic region is included or only the super-Alfvénic region is included. This is illustrated in Fig. 14.

Below the Alfvén point, there is no broad blueshifted absorption present. The redshifted absorption is comparable to the base model. Above the Alfvén point, in the jet itself, there is an extended blueshifted absorption. The redshifted absorption is still present but appears less pronounced. Just like the base model, there is a gap in the blueshifted part. The region below the Alfvén point has the largest contribution to the redshifted absorption component.

We made these models using an inclination of 60°, which is close to the highest possible inclination. At a lower inclination, the effect of the disk rotation is less severe, which can be seen in the inclination parameter study (Fig. 9), as the lower inclinations have less redshifted absorption. This strengthens the idea that the regions close to the disk are the main source of the redshifted rotational feature.

# 6. Discussion

The use of a physical steady-state jet model combined with our methodology to create synthetic observables has been proven successful. The steady-state MHD disk wind model produces, for our base parameter set, velocities and phase coverage of the absorption component in the synthetic dynamic spectra in the $H\alpha$ region that are very similar to our previous parametric model. The base model, shown in Fig. 8, shows a reasonable
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Inclination effect. The dynamical spectra in this list all utilise the orbital parameters and accretion disk parameters in Table 1 while changing the inclination according to Table 2. From top left to bottom right, the inclinations used are 18°, 36°, 54°, and 72°.</figcaption></figure>

velocity coverage at the higher end of the velocities. On the contrary, if the inner jet launching region becomes large, the central empty core of the jet increases, and the lack of absorption present at the middle of the phase coverage becomes larger. This occurs because our model keeps the medium void where there is no wind. By adding a different and additional outflow component on the symmetry axis of the bipolar jet (e.g. a stellar wind), this gap in velocity space in the synthetic spectra can be removed, but this is not included in our current model.

## 5.5. Disk versus jet rotation

To confirm that the redshifted absorption is due to rotation, we include Fig. 13. Here, a dynamic spectra based on a self-similar jet model for which the rotation has been set to zero and all other parameters are the same as in the base model (Fig. 8) is shown. The resulting dynamic spectra indeed appear very different compared to the base model: the redshifted absorption is no longer present, and the blueshifted part is more symmetrical. This clearly shows the importance of the rotation velocity in the production of the redshifted absorption.

All models with rotation produce a redshifted absorption in the spectra, something that is not observed in the data nor in the parametric models (as they do not include rotation). To understand what part of the wind contributes to the redshifted absorption, we divided the wind in two regions: the sub-Alfvénic region and the super-Alfvénic region. The sub-Alfvénic region includes the disk region and the starting wind, while the super-Alfvénic

region is the launched wind. We created two synthetic spectra where only the sub-Alfvénic region is included or only the super-Alfvénic region is included. This is illustrated in Fig. 14.

Below the Alfvén point, there is no broad blueshifted absorption present. The redshifted absorption is comparable to the base model. Above the Alfvén point, in the jet itself, there is an extended blueshifted absorption. The redshifted absorption is still present but appears less pronounced. Just like the base model, there is a gap in the blueshifted part. The region below the Alfvén point has the largest contribution to the redshifted absorption component.

We made these models using an inclination of 60°, which is close to the highest possible inclination. At a lower inclination, the effect of the disk rotation is less severe, which can be seen in the inclination parameter study (Fig. 9), as the lower inclinations have less redshifted absorption. This strengthens the idea that the regions close to the disk are the main source of the redshifted rotational feature.

# 6. Discussion

The use of a physical steady-state jet model combined with our methodology to create synthetic observables has been proven successful. The steady-state MHD disk wind model produces, for our base parameter set, velocities and phase coverage of the absorption component in the synthetic dynamic spectra in the $H\alpha$ region that are very similar to our previous parametric model. The base model, shown in Fig. 8, shows a reasonable
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 10. Mass accretion rate effect. These dynamic spectra all have the same input parameters (Tables 1 and 2) except for the mass accretion rate. These rates are, from top to bottom: $10^{-3} M_{\odot}yr^{-1}$, $10^{-4} M_{\odot}yr^{-1}$, and $10^{-5} M_{\odot}yr^{-1}$.</figcaption></figure>

accordance with the observations of HD 52961 when keeping in mind that we are not fitting the observations and used only one MHD model. However, since we are now using a physically consistent model that links the jet to the circum-companion accretion disk, several new questions arise.

An inner jet launching radius $r_{in}$ can be no smaller than ~20 R$_2$. Otherwise, the blueshifted absorption feature would extend to velocities far larger than those observed. This is of

course related to the reason why there should be no detectable jet arising from smaller disk regions.

As with the parametric wind models, we found a rather large disk accretion rate $M_a(r_{in})$ of about $10^{-4} M_{\odot}yr^{-1}$. If there is no way to keep wind density high while decreasing the disk accretion rate, then it raises the question as to what the lifetime of such an accretion-ejection event is.

Our MHD model predicts a redshifted absorption profile due to rotation, reaching velocities of about 100 km s$^{-1}$. While rotation was simply ignored in the parametric jet model, it is of course to be expected from any wind emitted from a rotating body. However, in none of the 16 objects published by Bollen et al. (2022), and references therein, has a redshifted absorption feature been detected.

Finally, a natural property of our MHD disk wind model is that it no longer exists inside $r_{in}$, leading to the creation of an empty core. This is indeed one of the common features in the disk winds creating systems (e.g. Bollen et al. 2022). However, this is also not detected in the various systems of interest. These four items are discussed further in the following sections.

## 6.1. Inner launching radius

Within our physical framework, the companion star is surrounded by an accretion disk extending up to ~R$_L$ for simplicity and fed by the circumbinary disk. Our main assumption is that there is a large-scale vertical field that threads the accretion disk and allows for the production of winds and jets. Such a magnetic field is probably advected along with the ionised disk material, as in the case of X-ray binaries or cataclysmic variables. What determines the extent of the jet emitting disk (hereafter, JED; Ferreira et al. 2006) is probably a matter of disk magnetisation, namely, whether the magnetic field is strong enough to launch jets. It sounds quite dubious to assume that such a magnetisation would decrease towards the star so that there is no clear reason why $r_{in}$ should be larger than the innermost disk radius $r_d$ itself<fnref n="1" />.

If the companion star has a dipolar magnetic field, such a field would be coupled to the accretion disk and truncate it, enforcing the disk material to accrete along magnetospheric columns or funnels. The position of this truncation radius thereby determines $r_d$ as is expressed in (Bessolaz et al. 2008; Zanni & Ferreira 2009, and references therein)

$$ \frac{r_d}{R_*} \approx 2m_s^{2/7} \left(\frac{B_*}{140 \text{ G}}\right)^{4/7} \left(\frac{\dot{M}_{a,2}}{10^{-8} M_{\odot} \text{yr}^{-1}}\right)^{-2/7} \left(\frac{M_2}{0.8 M_{\odot}}\right)^{-1/7} \left(\frac{R_2}{2 R_{\odot}}\right)^{5/7} \quad (8) $$

$$ \approx 0.005 m_s^{2/7} \left(\frac{B_*}{1 \text{ G}}\right)^{4/7} \left(\frac{\dot{M}_{a,2}}{10^{-4} M_{\odot} \text{yr}^{-1}}\right)^{-2/7} \left(\frac{M_2}{1 M_{\odot}}\right)^{-1/7} \left(\frac{R_2}{1 R_{\odot}}\right)^{5/7}, \quad (9) $$

where $\dot{M}_{a,2}$ is the disk accretion rate onto the star and $B_*$ is the stellar dipole field strength (see also Ireland et al. 2022 for numerical fits). Since the companion star is believed to be an old main-sequence solar-mass star, we do not expect field strengths of kilo Gauss and used instead the typical ~1 G value of the Sun. The above expression shows that for the parameters derived from our modelling, the disk would be expected to continue down to the stellar surface, namely, $r_{in} = r_d \sim R_2 = 4.6 \times 10^{-3} \text{ AU}$.

<footnote n="1">¹ We note that Bollen (2021) have already shown that the observed absorption features require wind velocities that are inconsistent with radiatively or thermally driven winds.</footnote>
<figure><img src="image_11.png" /><figcaption>Fig. 10. Mass accretion rate effect. These dynamic spectra all have the same input parameters (Tables 1 and 2) except for the mass accretion rate. These rates are, from top to bottom: $10^{-3} M_{\odot}yr^{-1}$, $10^{-4} M_{\odot}yr^{-1}$, and $10^{-5} M_{\odot}yr^{-1}$.</figcaption></figure>

accordance with the observations of HD 52961 when keeping in mind that we are not fitting the observations and used only one MHD model. However, since we are now using a physically consistent model that links the jet to the circum-companion accretion disk, several new questions arise.

An inner jet launching radius $r_{in}$ can be no smaller than ~20 R<sub>2</sub>. Otherwise, the blueshifted absorption feature would extend to velocities far larger than those observed. This is of

course related to the reason why there should be no detectable jet arising from smaller disk regions.

As with the parametric wind models, we found a rather large disk accretion rate $M_a(r_{in})$ of about $10^{-4} M_{\odot}yr^{-1}$. If there is no way to keep wind density high while decreasing the disk accretion rate, then it raises the question as to what the lifetime of such an accretion-ejection event is.

Our MHD model predicts a redshifted absorption profile due to rotation, reaching velocities of about 100 km s<sup>-1</sup>. While rotation was simply ignored in the parametric jet model, it is of course to be expected from any wind emitted from a rotating body. However, in none of the 16 objects published by Bollen et al. (2022), and references therein, has a redshifted absorption feature been detected.

Finally, a natural property of our MHD disk wind model is that it no longer exists inside $r_{in}$, leading to the creation of an empty core. This is indeed one of the common features in the disk winds creating systems (e.g. Bollen et al. 2022). However, this is also not detected in the various systems of interest. These four items are discussed further in the following sections.

## 6.1. Inner launching radius

Within our physical framework, the companion star is surrounded by an accretion disk extending up to ~R<sub>L</sub> for simplicity and fed by the circumbinary disk. Our main assumption is that there is a large-scale vertical field that threads the accretion disk and allows for the production of winds and jets. Such a magnetic field is probably advected along with the ionised disk material, as in the case of X-ray binaries or cataclysmic variables. What determines the extent of the jet emitting disk (hereafter, JED; Ferreira et al. 2006) is probably a matter of disk magnetisation, namely, whether the magnetic field is strong enough to launch jets. It sounds quite dubious to assume that such a magnetisation would decrease towards the star so that there is no clear reason why $r_{in}$ should be larger than the innermost disk radius $r_d$ itself<fnref n="1" />.

If the companion star has a dipolar magnetic field, such a field would be coupled to the accretion disk and truncate it, enforcing the disk material to accrete along magnetospheric columns or funnels. The position of this truncation radius thereby determines $r_d$ as is expressed in (Bessolaz et al. 2008; Zanni & Ferreira 2009, and references therein)

$$ \frac{r_d}{R_*} \approx 2m_s^{2/7} \left(\frac{B_*}{140 \text{ G}}\right)^{4/7} \left(\frac{\dot{M}_{a,2}}{10^{-8} M_{\odot}yr^{-1}}\right)^{-2/7} \left(\frac{M_2}{0.8 M_{\odot}}\right)^{-1/7} \left(\frac{R_2}{2 R_{\odot}}\right)^{5/7} \quad (8) $$

$$ \approx 0.005 m_s^{2/7} \left(\frac{B_*}{1 \text{ G}}\right)^{4/7} \left(\frac{\dot{M}_{a,2}}{10^{-4} M_{\odot}yr^{-1}}\right)^{-2/7} \left(\frac{M_2}{1 M_{\odot}}\right)^{-1/7} \left(\frac{R_2}{1 R_{\odot}}\right)^{5/7}, \quad (9) $$

where $\dot{M}_{a,2}$ is the disk accretion rate onto the star and $B_*$ is the stellar dipole field strength (see also Ireland et al. 2022 for numerical fits). Since the companion star is believed to be an old main-sequence solar-mass star, we do not expect field strengths of kilo Gauss and used instead the typical ~1 G value of the Sun. The above expression shows that for the parameters derived from our modelling, the disk would be expected to continue down to the stellar surface, namely, $r_{in} = r_d \sim R_2 = 4.6 \times 10^{-3}$ AU.

<footnote n="1">¹ We note that Bollen (2021) have already shown that the observed absorption features require wind velocities that are inconsistent with radiatively or thermally driven winds.</footnote>
<figure><img src="image_15.png" /><figcaption>Fig. 11. Outer launching radius effect. These two dynamic spectra share all input parameters (Tables 1 and 2) expect for the outer launching radius, which are, from left to right: 0.5 R<sub>L</sub>, and 1.0 R<sub>L</sub>.</figcaption></figure>

At such a small radius, the velocity space coverage would be too large for the model used, resulting in an extremely strong redshifted absorption. This is a strong hint that this MHD model, which has been used for YSOs, does not reproduce the physical conditions around the post-AGB companion star. We note that the parametric fit also provided an r<sub>in</sub> much larger than R<sub>2</sub> (see Table 1), so we believe this is quite generic.

## 6.2. Decreasing the disk accretion rate

The energy budget of an accretion disk undergoing mass loss through winds and jets can be written as

$$P_{\text{acc}} = 2P_{\text{jet}} + P_{\text{rad}} + P_{\text{adv}}, \qquad (10)$$

where P<sub>jet</sub> is the power feeding one jet, P<sub>rad</sub> is the disk luminosity, P<sub>adv</sub> is the power that is advected by the flow and released at the accretion shock onto the hard surface of the star, and

$$P_{\text{acc}} = \left[ \frac{GM_* \dot{M}_a}{2r} \right]_{\text{out}}^{\text{in}} = \frac{GM_2 \dot{M}_{a,in}}{2r_{\text{in}}} - \frac{GM_2 \dot{M}_{a,out}}{2r_{\text{out}}} \qquad (11) \\ \approx \frac{GM_2 \dot{M}_{a,in}}{2r_{\text{in}}} = 80 \left( \frac{\dot{M}_{a,in}}{10^{-4} M_{\odot} \text{yr}^{-1}} \right) \left( \frac{M_2}{1 M_{\odot}} \right) \left( \frac{r_{\text{in}}}{20 R_2} \right)^{-1} L_{\odot} \qquad (12)$$

is the released accretion power. Now P<sub>adv</sub>/P<sub>acc</sub> ≈ ε<sup>2</sup> and is therefore neglected, as we assumed the inner disk to be geometrically thin.

According to the jet model we used, the accretion power released up to r<sub>in</sub> = 20 R<sub>2</sub> amounts to ~80 L<sub>☉</sub>, 88% of which is fed into the two jets. However, mass within the disk must proceed down to the star with r<sub>d</sub> = R<sub>2</sub> and at the same rate. This will lead to an extra ~1500 L<sub>☉</sub> power, which is mostly released as radiation in the disk. Defining the disk radiative efficiency as η<sub>rad</sub> = P<sub>rad</sub>/P<sub>acc</sub> and assuming that the disk is optically thick, one gets the following rough estimate of its innermost effective temperature:

$$T_{\text{eff}} = \left( \eta_{\text{rad}} \frac{GM_2 \dot{M}_a}{8\pi\sigma r^3} \right)^{1/4} \qquad (13) \\ = 3.64 \times 10^4 \left( \frac{\eta_{\text{rad}} \dot{M}_a(r)}{10^{-4} M_{\odot} \text{yr}^{-1}} \right)^{1/4} \left( \frac{M_2}{M_{\odot}} \right)^{1/4} \left( \frac{r}{R_{\odot}} \right)^{-3/4} \text{ K.} \qquad (14)$$

This estimate shows that an accretion disk extending down to r<sub>d</sub> around a solar-like star should have a temperature of η<sub>rad</sub><sup>1/4</sup>

36 400 K. Observationally, the contribution of the accretion disk is only a minor fraction of the total emissivity in all wavelengths. Using high spatial resolution near-IR interferometric imaging, the contribution of the accretion disk was spatially resolved in one of the jet-producing systems (IRAS 08544–4431). This yielded a temperature of ~4000 K only for the accretion disk in that system (Hillen et al. 2014; Kluska et al. 2018). Thus, unless η<sub>rad</sub> ~ 10<sup>-4</sup>, accretion proceeding down to the star at such high accretion rates is inconsistent with known observational constraints. Even in the case of a JED established down to the star, and for the solution used, η<sub>rad</sub> = 0.12 is a value orders of magnitude larger than required.

The other obvious solution, namely, a much smaller disk accretion rate, seems more plausible. Inverting Eq. (14) and assuming T<sub>eff</sub> = 4000 K leads to a maximum accretion rate of

$$\dot{M}_a(r_d) = 2 \times 10^{-8} \eta_{\text{rad}}^{-1} \left( \frac{T_{\text{eff}}}{4000 \text{ K}} \right)^4 \left( \frac{M_2}{M_{\odot}} \right)^{-1} \left( \frac{r_d}{R_{\odot}} \right)^3 M_{\odot} \text{yr}^{-1}, \qquad (15)$$

namely, around 10<sup>-7</sup> M yr<sup>-1</sup> for a typical radiation efficiency, and a disk touching the star<sup>2</sup>. This is however three orders of magnitude smaller than the value used in our base model, which is somewhat problematic.

Although they are both model dependent, our two independent constraints (wind absorption features and disk temperature) are in strong tension. The circumbinary disk, which has a mass of approximately 10<sup>-2</sup> M (Bujarrabal et al. 2015), is believed to be the main fuel source for the circum-companion accretion disk (Bollen et al. 2022). With the derived high disk-accretion rate, the lifetime would be only ~10<sup>2</sup> years. If the lifetime of the circumbinary disk is this short, we would expect to find very few of these systems. In contradiction, a large number of the post-AGB binaries have a detectable disk (Kluska et al. 2022).

However, as seen above, a much smaller accretion rate than what used in our wind modelling would better fit the classical picture of a continuous accretion-ejection flow down to the star. But in that case, one would need to recover approximately the same wind density as the one found in our base model. Assuming ejection starts from r<sub>in</sub> = R<sub>2</sub> ensures that the midplane density will already be ~70–90 times larger than the midplane density at 20 R<sub>2</sub> for the same accretion rate. But the required decrease of the disk accretion rate needs to be compensated for by a much

<footnote n="2">2 We note that inserting this expression into Eq. (8) indeed shows that unless the stellar dipole field is much larger than 1 G, it is unable to truncate the disk, and r<sub>d</sub> = R<sub>2</sub> is a good approximation.</footnote>
<figure><img src="image_12.png" /><figcaption>Fig. 11. Outer launching radius effect. These two dynamic spectra share all input parameters (Tables 1 and 2) expect for the outer launching radius, which are, from left to right: 0.5 R<sub>L</sub>, and 1.0 R<sub>L</sub>.</figcaption></figure>

At such a small radius, the velocity space coverage would be too large for the model used, resulting in an extremely strong redshifted absorption. This is a strong hint that this MHD model, which has been used for YSOs, does not reproduce the physical conditions around the post-AGB companion star. We note that the parametric fit also provided an r<sub>in</sub> much larger than R<sub>2</sub> (see Table 1), so we believe this is quite generic.

## 6.2. Decreasing the disk accretion rate

The energy budget of an accretion disk undergoing mass loss through winds and jets can be written as

$$P_{\text{acc}} = 2P_{\text{jet}} + P_{\text{rad}} + P_{\text{adv}}, \qquad (10)$$

where P<sub>jet</sub> is the power feeding one jet, P<sub>rad</sub> is the disk luminosity, P<sub>adv</sub> is the power that is advected by the flow and released at the accretion shock onto the hard surface of the star, and

$$P_{\text{acc}} = \left[ \frac{GM_* \dot{M}_a}{2r} \right]_{\text{out}}^{\text{in}} = \frac{GM_2 \dot{M}_{a,in}}{2r_{\text{in}}} - \frac{GM_2 \dot{M}_{a,out}}{2r_{\text{out}}} \qquad (11) \\ \approx \frac{GM_2 \dot{M}_{a,in}}{2r_{\text{in}}} = 80 \left( \frac{\dot{M}_{a,in}}{10^{-4} M_{\odot} \text{yr}^{-1}} \right) \left( \frac{M_2}{1 M_{\odot}} \right) \left( \frac{r_{\text{in}}}{20 R_2} \right)^{-1} L_{\odot} \qquad (12)$$

is the released accretion power. Now P<sub>adv</sub>/P<sub>acc</sub> ≈ ε<sup>2</sup> and is therefore neglected, as we assumed the inner disk to be geometrically thin.

According to the jet model we used, the accretion power released up to r<sub>in</sub> = 20 R<sub>2</sub> amounts to ~80 L<sub>☉</sub>, 88% of which is fed into the two jets. However, mass within the disk must proceed down to the star with r<sub>d</sub> = R<sub>2</sub> and at the same rate. This will lead to an extra ~1500 L<sub>☉</sub> power, which is mostly released as radiation in the disk. Defining the disk radiative efficiency as η<sub>rad</sub> = P<sub>rad</sub>/P<sub>acc</sub> and assuming that the disk is optically thick, one gets the following rough estimate of its innermost effective temperature:

$$T_{\text{eff}} = \left( \eta_{\text{rad}} \frac{GM_2 \dot{M}_a}{8\pi\sigma r^3} \right)^{1/4} \qquad (13) \\ = 3.64 \times 10^4 \left( \frac{\eta_{\text{rad}} \dot{M}_a(r)}{10^{-4} M_{\odot} \text{yr}^{-1}} \right)^{1/4} \left( \frac{M_2}{M_{\odot}} \right)^{1/4} \left( \frac{r}{R_{\odot}} \right)^{-3/4} \text{ K.} \qquad (14)$$

This estimate shows that an accretion disk extending down to r<sub>d</sub> around a solar-like star should have a temperature of η<sub>rad</sub><sup>1/4</sup>

36 400 K. Observationally, the contribution of the accretion disk is only a minor fraction of the total emissivity in all wavelengths. Using high spatial resolution near-IR interferometric imaging, the contribution of the accretion disk was spatially resolved in one of the jet-producing systems (IRAS 08544–4431). This yielded a temperature of ~4000 K only for the accretion disk in that system (Hillen et al. 2014; Kluska et al. 2018). Thus, unless η<sub>rad</sub> ~ 10<sup>-4</sup>, accretion proceeding down to the star at such high accretion rates is inconsistent with known observational constraints. Even in the case of a JED established down to the star, and for the solution used, η<sub>rad</sub> = 0.12 is a value orders of magnitude larger than required.

The other obvious solution, namely, a much smaller disk accretion rate, seems more plausible. Inverting Eq. (14) and assuming T<sub>eff</sub> = 4000 K leads to a maximum accretion rate of

$$\dot{M}_a(r_d) = 2 \times 10^{-8} \eta_{\text{rad}}^{-1} \left( \frac{T_{\text{eff}}}{4000 \text{ K}} \right)^4 \left( \frac{M_2}{M_{\odot}} \right)^{-1} \left( \frac{r_d}{R_{\odot}} \right)^3 M_{\odot} \text{yr}^{-1}, \quad (15)$$

namely, around 10<sup>-7</sup> M yr<sup>-1</sup> for a typical radiation efficiency, and a disk touching the star<sup>2</sup>. This is however three orders of magnitude smaller than the value used in our base model, which is somewhat problematic.

Although they are both model dependent, our two independent constraints (wind absorption features and disk temperature) are in strong tension. The circumbinary disk, which has a mass of approximately 10<sup>-2</sup> M (Bujarrabal et al. 2015), is believed to be the main fuel source for the circum-companion accretion disk (Bollen et al. 2022). With the derived high disk-accretion rate, the lifetime would be only ~10<sup>2</sup> years. If the lifetime of the circumbinary disk is this short, we would expect to find very few of these systems. In contradiction, a large number of the post-AGB binaries have a detectable disk (Kluska et al. 2022).

However, as seen above, a much smaller accretion rate than what used in our wind modelling would better fit the classical picture of a continuous accretion-ejection flow down to the star. But in that case, one would need to recover approximately the same wind density as the one found in our base model. Assuming ejection starts from r<sub>in</sub> = R<sub>2</sub> ensures that the midplane density will already be ~70–90 times larger than the midplane density at 20 R<sub>2</sub> for the same accretion rate. But the required decrease of the disk accretion rate needs to be compensated for by a much

<footnote n="2">2 We note that inserting this expression into Eq. (8) indeed shows that unless the stellar dipole field is much larger than 1 G, it is unable to truncate the disk, and r<sub>d</sub> = R<sub>2</sub> is a good approximation.</footnote>
<figure><img src="image_16.png" /><figcaption>**Fig. 12.** Inner launching radius effect. All input parameters of these plots are the same (Tables 1 and 2) excluding the inner launching radius. The values are, from top to bottom: 1 $R_*$, 11 $R_*$, and 41 $R_*$.</figcaption></figure>

larger disk ejection efficiency $\xi$, namely, by using another MHD model.

This necessary modification of the MHD jet model may also
compensate for another consequence of putting $r_{\text{in}}$ at the stellar
surface. As shown in Fig. 11, a smaller $r_{\text{in}}$ leads to an impor-
tant modification of the velocity range since at $r_{\text{in}} = R_2$, the
Keplerian rotation is already 440 km s$^{-1}$. The MHD wind there-
fore needs to provide only a mild acceleration, which is possible
only with $\xi > 0.3$ or so (since $v_{\text{jet}} \approx \sqrt{\frac{1}{\xi} - 1\Omega_{K_0}r_0}$). But this

<figure><img src="image_17.png" /><figcaption>**Fig. 13.** Synthetic dynamic spectra of a self-similar model with no rotation.</figcaption></figure>

requires absorption signatures to mostly arise from regions close
to the accretion disk, where the poloidal velocity is still small.
Only a proper calculation of dynamic spectra will allow for con-
firmation of whether new MHD models launched close to the
stellar surface can indeed reproduce the observations.

Finally, another possibility would be to better compute the jet
temperature. For simplicity, we assumed an isothermal jet, but a
temperature distribution is actually expected, as a consequence
of thermal and ionisation equilibrium within an expanding out-
flow. By changing the jet temperature, one could, for instance,
increase the occupation of the first excited state of hydrogen,
which would create a stronger absorption feature at a given
jet density.

However, the temperature of the jet is limited to the tem-
perature of the post-AGB star (which is well determined) since
the jet is always seen in absorption. Indeed, the jet would be
seen in emission instead of absorption if a higher jet temperature
was assumed compared to the effective temperature of the post-
AGB star, which provides the background radiation. Incorporat-
ing non-local thermal equilibrium (NLTE) effects in the deter-
mination of the occupation numbers of the first excited state is
therefore another viable complexity increase. This has the bene-
fit of being able to pump hydrogen to the excitation stage with-
out needing extra density, and by using departure coefficients,
this could be done in a computationally inexpensive manner. The
inclusion of NLTE effects is postponed for future work.

6.3. *The absence of rotation signatures*

The main discrepancy between our physical jet model and obser-
vations is the presence of a redshifted absorption in the syn-
thetic spectra. This is present in nearly all our 19 models but
is missing in all dynamic spectra observed in post-AGB binaries
(Bollen 2021). As shown before, this is a natural consequence of
a wind being emitted from a rotating body. This feature was not
present in our previous investigations simply because our para-
metric wind model assumed zero rotation.

However, when artificially taking rotation out of our phys-
ical jet model, this redshifted feature does indeed disappear,
but the absorption pattern also becomes too symmetric (see
Fig. 13). This may be a sign that some rotation is indeed present
in real systems but that it is never as large as in the models.
Therefore, this puts interesting stringent constraints on wind
models. This is especially critical since, as discussed before,
<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>**Fig. 12.** Inner launching radius effect. All input parameters of these plots are the same (Tables 1 and 2) excluding the inner launching radius. The values are, from top to bottom: 1 $R_*$, 11 $R_*$, and 41 $R_*$.</figcaption></figure>

larger disk ejection efficiency $\xi$, namely, by using another MHD model.

This necessary modification of the MHD jet model may also
compensate for another consequence of putting $r_{\text{in}}$ at the stellar
surface. As shown in Fig. 11, a smaller $r_{\text{in}}$ leads to an impor-
tant modification of the velocity range since at $r_{\text{in}} = R_2$, the
Keplerian rotation is already 440 km s$^{-1}$. The MHD wind there-
fore needs to provide only a mild acceleration, which is possible
only with $\xi > 0.3$ or so (since $v_{\text{jet}} \approx \sqrt{\frac{1}{\xi} - 1\Omega_{K_0}r_o}$). But this

<figure><img src="image_16.png" /><figcaption>**Fig. 13.** Synthetic dynamic spectra of a self-similar model with no rotation.</figcaption></figure>

requires absorption signatures to mostly arise from regions close
to the accretion disk, where the poloidal velocity is still small.
Only a proper calculation of dynamic spectra will allow for con-
firmation of whether new MHD models launched close to the
stellar surface can indeed reproduce the observations.

Finally, another possibility would be to better compute the jet
temperature. For simplicity, we assumed an isothermal jet, but a
temperature distribution is actually expected, as a consequence
of thermal and ionisation equilibrium within an expanding out-
flow. By changing the jet temperature, one could, for instance,
increase the occupation of the first excited state of hydrogen,
which would create a stronger absorption feature at a given
jet density.

However, the temperature of the jet is limited to the tem-
perature of the post-AGB star (which is well determined) since
the jet is always seen in absorption. Indeed, the jet would be
seen in emission instead of absorption if a higher jet temperature
was assumed compared to the effective temperature of the post-
AGB star, which provides the background radiation. Incorporat-
ing non-local thermal equilibrium (NLTE) effects in the deter-
mination of the occupation numbers of the first excited state is
therefore another viable complexity increase. This has the bene-
fit of being able to pump hydrogen to the excitation stage with-
out needing extra density, and by using departure coefficients,
this could be done in a computationally inexpensive manner. The
inclusion of NLTE effects is postponed for future work.

6.3. *The absence of rotation signatures*

The main discrepancy between our physical jet model and obser-
vations is the presence of a redshifted absorption in the syn-
thetic spectra. This is present in nearly all our 19 models but
is missing in all dynamic spectra observed in post-AGB binaries
(Bollen 2021). As shown before, this is a natural consequence of
a wind being emitted from a rotating body. This feature was not
present in our previous investigations simply because our para-
metric wind model assumed zero rotation.

However, when artificially taking rotation out of our phys-
ical jet model, this redshifted feature does indeed disappear,
but the absorption pattern also becomes too symmetric (see
Fig. 13). This may be a sign that some rotation is indeed present
in real systems but that it is never as large as in the models.
Therefore, this puts interesting stringent constraints on wind
models. This is especially critical since, as discussed before,
<figure><img src="image_18.png" /><figcaption>Fig. 14. Disk winds separated into sub- or super-Alfvénic. The synthetic dynamic spectra of the models only includes either sub- or super-Alfvénic flow. In the left panel, all the material that is above the Alfvénic point has been removed. In the right panel, all the material below the Alfvén point has been removed.</figcaption></figure>

we expect ejection to take place up to the companion surface, where Keplerian speed will be 440 km s⁻¹.

There are two different means of decreasing the rotation rate of an MHD wind. The first one is to simply decrease the rotation of the disk itself so that it becomes significantly sub-Keplerian. The deviation to the Keplerian rotation law at the disk midplane is

$$ \delta^2 = \frac{\Omega_o^2}{\Omega_{K_o}^2} = 1 - \left(\frac{5}{2} - \xi + \frac{m_s^2}{2}\right) \epsilon^2 - p\mu\epsilon, \quad (16) $$

where $p \sim B_r^+/B_z$ is a measure of the magnetic field curvature at the disk surface and is on the order of unity (Casse & Ferreira 2000a). As can be seen from this expression, rotation depends mostly on the disk aspect ratio $\epsilon$ and the disk magnetisation $\mu$. The thicker (hotter) the disk becomes, the slower the rotation gets. In order to reach values of $\delta \sim 0.5$ or less, the disk must become geometrically thick with $\epsilon \sim 0.3$ or so. This might, however, be problematic for purely centrifugally driven winds. Indeed, for such a driving mechanism, the magnetic field lines need to be bent by more than 30° at the disk surface (i.e. $p \ge 1$, Blandford & Payne 1982), which results in a decrease of the rotation speed of the plasma at the disk surface. But if that surface rotation becomes too small, the centrifugal push becomes inefficient to drive ejection. This is the reason why purely magnetically driven winds cannot be launched from thick accretion disks (Casse & Ferreira 2000a). However, the accretion disk in our case is irradiated by the post-AGB primary star. This certainly provides additional heating at the disk surface so that the wind is probably also thermally driven. The MHD solution used here already includes such a thermal effect, but it remained negligible (less than 1%, parameter $f$). Therefore, to capture the influence of the post-AGB star on the secondary accretion disk, one needs to seek thermomagnetic MHD winds with a large thermal input ($f$). We note that this condition is also consistent with the constraint of having a denser wind, as thermal effects are known to enhance the disk ejection efficiency $\xi$ (Casse & Ferreira 2000b).

The second way to decrease the wind rotation is by lowering the MHD efficiency of the wind acceleration. This is done by increasing the jet mass loading parameter $\kappa$ (Blandford & Payne 1982), which is related to the disk by Ferreira (1997), Jacquemin-Ide et al. (2019):

$$ \kappa \approx \frac{\xi m_s}{\mu} = \alpha_m p \xi \mu^{-1/2}. \quad (17) $$

Thus, the larger $\kappa$ is, the more mass that is loaded per magnetic flux unit, implying that the magnetic field becomes less able to shape the outflow. It can be seen right away that the more massive the outflow (i.e. $\xi$ large) and/or the smaller the disk magnetisation (i.e. $\mu$ small), the more the outflow becomes matter dominated. In the large $\kappa$ limit, the magnetic lever arm $\lambda$ decreases close to unity, and the field lines become initially more conical, leading to a drastic and fast decrease of the wind rotation velocity (Jacquemin-Ide et al. 2019).

To summarise all the above constraints, a new wind model setup would require ejection down to the stellar surface $R_2$ with a disk accretion rate two or three orders of magnitude smaller than that used here, requiring an MHD wind model that is deeply modified. These modifications include: a much thicker accretion disk with $\epsilon = 0.3$ or more; a much larger disk ejection efficiency, $\xi = 0.3$ or more, probably with a significant thermal content as well (larger f value); and a lower disk magnetisation $\mu$. This exploration will be done in a subsequent work.

### 6.4. The inner axial wind

Our disk wind model has, by construction, a central empty space between the vertical axis and the innermost field line anchored at $r_{in}$. In order to assess the effect of only the disk wind, we arbitrarily set the density to zero in this inner region. This void can actually be seen in the synthetic dynamical spectra as a slightly redshifted zone of no absorption (see Fig. 8). However, observations do not show any evidence of such a feature. Moreover, our exploration of the effect of $r_{in}$ shows that even for a wind with $r_{in} = R_2$ launched near the companion star, our MHD model still displays this lack of absorption. Again, only an approach using a different MHD model will firmly assess whether this is systematic to disk winds or dependant on the specific model used. According to the exploration of many objects done by Bollen et al. (2022), this is systematic to disk winds. If this is so, then it means that an extra stellar wind component must also be included in the model. There would be no special difficulty to compute the minimum density, and hence mass loss, filling in this inner spine, which is required to recover the observed spectra.

## 7. Conclusion

A methodology to produce dynamic spectra of disk winds in post-AGB binaries using a self-consistent physical description
<figure><img src="image_17.png" /><figcaption>Fig. 14. Disk winds separated into sub- or super-Alfvénic. The synthetic dynamic spectra of the models only includes either sub- or super-Alfvénic flow. In the left panel, all the material that is above the Alfvénic point has been removed. In the right panel, all the material below the Alfvén point has been removed.</figcaption></figure>

we expect ejection to take place up to the companion surface, where Keplerian speed will be 440 km s⁻¹.

There are two different means of decreasing the rotation rate of an MHD wind. The first one is to simply decrease the rotation of the disk itself so that it becomes significantly sub-Keplerian. The deviation to the Keplerian rotation law at the disk midplane is

$$ \delta^2 = \frac{\Omega_o^2}{\Omega_{K_o}^2} = 1 - \left(\frac{5}{2} - \xi + \frac{m_s^2}{2}\right) \epsilon^2 - p\mu\epsilon, \quad (16) $$

where $p \sim B_r^+/B_z$ is a measure of the magnetic field curvature at the disk surface and is on the order of unity (Casse & Ferreira 2000a). As can be seen from this expression, rotation depends mostly on the disk aspect ratio $\epsilon$ and the disk magnetisation $\mu$. The thicker (hotter) the disk becomes, the slower the rotation gets. In order to reach values of $\delta \sim 0.5$ or less, the disk must become geometrically thick with $\epsilon \sim 0.3$ or so. This might, however, be problematic for purely centrifugally driven winds. Indeed, for such a driving mechanism, the magnetic field lines need to be bent by more than 30° at the disk surface (i.e. $p \ge 1$, Blandford & Payne 1982), which results in a decrease of the rotation speed of the plasma at the disk surface. But if that surface rotation becomes too small, the centrifugal push becomes inefficient to drive ejection. This is the reason why purely magnetically driven winds cannot be launched from thick accretion disks (Casse & Ferreira 2000a). However, the accretion disk in our case is irradiated by the post-AGB primary star. This certainly provides additional heating at the disk surface so that the wind is probably also thermally driven. The MHD solution used here already includes such a thermal effect, but it remained negligible (less than 1%, parameter $f$). Therefore, to capture the influence of the post-AGB star on the secondary accretion disk, one needs to seek thermomagnetic MHD winds with a large thermal input ($f$). We note that this condition is also consistent with the constraint of having a denser wind, as thermal effects are known to enhance the disk ejection efficiency $\xi$ (Casse & Ferreira 2000b).

The second way to decrease the wind rotation is by lowering the MHD efficiency of the wind acceleration. This is done by increasing the jet mass loading parameter $\kappa$ (Blandford & Payne 1982), which is related to the disk by Ferreira (1997), Jacquemin-Ide et al. (2019):

$$ \kappa \approx \frac{\xi m_s}{\mu} = \alpha_m p \xi \mu^{-1/2}. \quad (17) $$

Thus, the larger $\kappa$ is, the more mass that is loaded per magnetic flux unit, implying that the magnetic field becomes less able to shape the outflow. It can be seen right away that the more massive the outflow (i.e. $\xi$ large) and/or the smaller the disk magnetisation (i.e. $\mu$ small), the more the outflow becomes matter dominated. In the large $\kappa$ limit, the magnetic lever arm $\lambda$ decreases close to unity, and the field lines become initially more conical, leading to a drastic and fast decrease of the wind rotation velocity (Jacquemin-Ide et al. 2019).

To summarise all the above constraints, a new wind model setup would require ejection down to the stellar surface $R_2$ with a disk accretion rate two or three orders of magnitude smaller than that used here, requiring an MHD wind model that is deeply modified. These modifications include: a much thicker accretion disk with $\epsilon = 0.3$ or more; a much larger disk ejection efficiency, $\xi = 0.3$ or more, probably with a significant thermal content as well (larger $f$ value); and a lower disk magnetisation $\mu$. This exploration will be done in a subsequent work.

### 6.4. The inner axial wind

Our disk wind model has, by construction, a central empty space between the vertical axis and the innermost field line anchored at $r_{in}$. In order to assess the effect of only the disk wind, we arbitrarily set the density to zero in this inner region. This void can actually be seen in the synthetic dynamical spectra as a slightly redshifted zone of no absorption (see Fig. 8). However, observations do not show any evidence of such a feature. Moreover, our exploration of the effect of $r_{in}$ shows that even for a wind with $r_{in} = R_2$ launched near the companion star, our MHD model still displays this lack of absorption. Again, only an approach using a different MHD model will firmly assess whether this is systematic to disk winds or dependant on the specific model used. According to the exploration of many objects done by Bollen et al. (2022), this is systematic to disk winds. If this is so, then it means that an extra stellar wind component must also be included in the model. There would be no special difficulty to compute the minimum density, and hence mass loss, filling in this inner spine, which is required to recover the observed spectra.

## 7. Conclusion

A methodology to produce dynamic spectra of disk winds in post-AGB binaries using a self-consistent physical description
for the wind is now functional. This approach not only allows us
to study the behaviour of the disk wind, but it also puts strong
constraints on the disk itself. In this work, we used one particular
MHD model that has already been proven to explain most obser-
vational constraints in YSO jets. The resulting synthetic spec-
tra are able to cover velocities, phase coverages, and absorption
strengths very similar to those from observations of HD 52961.
This is very promising and indicates that magnetically driven jets
launched from accretion disks may be a good representation of
the actual winds detected in post-AGB systems.

However, observed velocity ranges require the jet emitting
zone to be quite far (≥20 R<sub>2</sub>) from the central companion star,
which is somewhat problematic. Moreover, observed absorp-
tion features require rather high wind densities that can only
be achieved at those distances by assuming a very large disk
accretion rate (as it is self-consistent with the wind mass-loss
rate). This need for a high mass accretion rate was already found
using parametric winds (Bollen 2021) and is also problematic,
as it questions the accretion disk lifetime. Moreover, such high
accretion rates should lead to the innermost disk regions being
much hotter than what is currently observed. Most remarkably,
the redshifted absorption produced by the rotation of the jet does
not appear in observations.

Nevertheless, we argue that the accretion disk in HD 52961
has most probably settled down to the surface of the companion
star. This is based on the conservative assumption that the stellar
magnetic field is solar-like and thus unable to truncate its circum-
stellar accretion disk and enforce material to accrete along the
magnetosphere. Although, accretion should lead to stellar spin-
up under such circumstances, we suspect that timescales are too
short to lead to a revitalised dynamo action and a stellar magnetic
field larger than solar-like.

As a consequence of our first use of a self-consistent
accretion-ejection model, we argue that disk winds require a
different MHD configuration: a larger disk ejection efficiency
to enhance wind density and lower velocities, a thicker (hot-
ter) accretion disk to lower disk and wind rotation, and prob-
ably a smaller disk magnetisation, along with some thermal
input at the disk surface (due to the post-AGB radiation field).
Such models already exist in the literature (Casse & Ferreira
2000b; Jacquemin-Ide et al. 2019), and the exploration of differ-
ent thermo-magnetic wind models will be performed in a forth-
coming paper.

Acknowledgements. O.V. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). H.V.W. acknowledges support of the FWO under grant G097619N. J.K. acknowledges support from FWO under the senior postdoctoral fellowship (1281121N). T.D.P. acknowledges support of the Research Foundation-Flanders (FWO) under grant 11P6I24N (Aspirant Fellowship). D.B. and H.V.W. acknowledge support from the Research Council of the KU Leuven under grant number C14/17/082. D.K. acknowledges the support

of the Australian Research Council (ARC) Discovery Early Career Research
Award (DECRA) grant (DE190100813). D.K. is also supported in part by the
Australian Research Council Centre of Excellence for All Sky Astrophysics
in 3 Dimensions (ASTRO 3D), through project number CE170100013. Based
on observations obtained with the HERMES spectrograph, which is supported
by the Research Foundation - Flanders (FWO), Belgium, the Research Coun-
cil of KU Leuven, Belgium, the Fonds National de la Recherche Scientifique
(F.R.S.-FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de
Genève, Switzerland and the Thüringer Landessternwarte Tautenburg, Germany.
Finally we would like to thank the anonymous referee to help streamline the
paper.

References

Balbus, S. A. 2003, ARA&A, 41, 555

Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155

Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883

Bollen, D. 2021, PhD Thesis, KU Leuven, Belgium, https://fys.kuleuven.be/ster/pub/phd-thesis-dylan-bollen/phd-thesis-dylan-bollen

Bollen, D., Van Winckel, H., & Kamath, D. 2017, A&A, 607, A60

Bollen, D., Kamath, D., Van Winckel, H., & De Marco, O. 2019, A&A, 631, A53

Bollen, D., Kamath, D., De Marco, O., Van Winckel, H., & Wardle, M. 2020, A&A, 641, A175

Bollen, D., Kamath, D., Van Winckel, H., et al. 2022, A&A, 666, A40

Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., & Van Winckel, H. 2015, A&A, 575, L7

Casse, F., & Ferreira, J. 2000a, A&A, 353, 1115

Casse, F., & Ferreira, J. 2000b, A&A, 361, 1178

Demircan, O., & Kahraman, G. 1991, Astrophys. Space Sci., 181, 313

Eggleton, P. P. 1983, ApJ, 268, 368

Ferreira, J. 1997, A&A, 319, 340

Ferreira, J., & Casse, F. 2013, MNRAS, 428, 307

Ferreira, J., & Pelletier, G. 1995, A&A, 295, 807

Ferreira, J., Petrucci, P. O., Henri, G., Saugé, L., & Pelletier, G. 2006, A&A, 447, 813

Gorlova, N., Van Winckel, H., Gielen, C., et al. 2012, A&A, 542, A27

Hillen, M., Menu, J., Van Winckel, H., et al. 2014, A&A, 568, A12

Hillen, M., Kluska, J., Bouquin, J.-B. L., et al. 2016, A&A, 588, L1

Ireland, L. G., Matt, S. P., & Zanni, C. 2022, ApJ, 929, 65

Jacquemin-Ide, J., Ferreira, J., & Lesur, G. 2019, MNRAS, 490, 3112

Jacquemin-Ide, J., Lesur, G., & Ferreira, J. 2021, A&A, 647, A192

Kamath, D., Wood, P. R., & Van Winckel, H. 2015, MNRAS, 454, 1468

Kluska, J., Hillen, M., Van Winckel, H., et al. 2018, A&A, 616, A153

Kluska, J., Van Winckel, H., Coppée, Q., et al. 2022, A&A, 658, A36

Manick, R., Van Winckel, H., Kamath, D., Hillen, M., & Escorza, A. 2017, A&A, 597, A129

Oomen, G.-M., Van Winckel, H., Pols, O., et al. 2018, A&A, 620, A85

Oomen, G.-M., Van Winckel, H., Pols, O., & Nelemans, G. 2019, A&A, 629, A49

Panoglou, D., Cabrit, S., Forêts, G. P. D., et al. 2012, A&A, 538, A2

Raskin, G. 2011, PhD Thesis, KU Leuven, Belgium, https://fys.kuleuven.be/ster/pub/phd-thesis-gert-raskin/phd-thesis-gert-raskin-1

Raskin, G., Van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69

Van Winckel, H. 2003, ARA&A, 41, 391

Van Winckel, H. 2007, Balt. Astron., 16, 112

Van Winckel, H. 2018, Camb. Astrophys. Ser., 54, 92

Zanni, C., & Ferreira, J. 2009, A&A, 508, 1117
for the wind is now functional. This approach not only allows us to study the behaviour of the disk wind, but it also puts strong constraints on the disk itself. In this work, we used one particular MHD model that has already been proven to explain most observational constraints in YSO jets. The resulting synthetic spectra are able to cover velocities, phase coverages, and absorption strengths very similar to those from observations of HD 52961. This is very promising and indicates that magnetically driven jets launched from accretion disks may be a good representation of the actual winds detected in post-AGB systems.

However, observed velocity ranges require the jet emitting zone to be quite far (≥20 $R_2$) from the central companion star, which is somewhat problematic. Moreover, observed absorption features require rather high wind densities that can only be achieved at those distances by assuming a very large disk accretion rate (as it is self-consistent with the wind mass-loss rate). This need for a high mass accretion rate was already found using parametric winds (Bollen 2021) and is also problematic, as it questions the accretion disk lifetime. Moreover, such high accretion rates should lead to the innermost disk regions being much hotter than what is currently observed. Most remarkably, the redshifted absorption produced by the rotation of the jet does not appear in observations.

Nevertheless, we argue that the accretion disk in HD 52961 has most probably settled down to the surface of the companion star. This is based on the conservative assumption that the stellar magnetic field is solar-like and thus unable to truncate its circumstellar accretion disk and enforce material to accrete along the magnetosphere. Although, accretion should lead to stellar spin-up under such circumstances, we suspect that timescales are too short to lead to a revitalised dynamo action and a stellar magnetic field larger than solar-like.

As a consequence of our first use of a self-consistent accretion-ejection model, we argue that disk winds require a different MHD configuration: a larger disk ejection efficiency to enhance wind density and lower velocities, a thicker (hotter) accretion disk to lower disk and wind rotation, and probably a smaller disk magnetisation, along with some thermal input at the disk surface (due to the post-AGB radiation field). Such models already exist in the literature (Casse & Ferreira 2000b; Jacquemin-Ide et al. 2019), and the exploration of different thermo-magnetic wind models will be performed in a forthcoming paper.

**Acknowledgements.** O.V. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). H.V.W. acknowledges support of the FWO under grant G097619N. J.K. acknowledges support from FWO under the senior postdoctoral fellowship (1281121N). T.D.P. acknowledges support of the Research Foundation-Flanders (FWO) under grant 11P6I24N (Aspirant Fellowship). D.B. and H.V.W. acknowledge support from the Research Council of the KU Leuven under grant number C14/17/082. D.K. acknowledges the support

of the Australian Research Council (ARC) Discovery Early Career Research Award (DECRA) grant (DE190100813). D.K. is also supported in part by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. Based on observations obtained with the HERMES spectrograph, which is supported by the Research Foundation - Flanders (FWO), Belgium, the Research Council of KU Leuven, Belgium, the Fonds National de la Recherche Scientifique (F.R.S.-FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de Genève, Switzerland and the Thüringer Landessternwarte Tautenburg, Germany. Finally we would like to thank the anonymous referee to help streamline the paper.

## References

* Balbus, S. A. 2003, ARA&A, 41, 555

* Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155

* Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883

* Bollen, D. 2021, PhD Thesis, KU Leuven, Belgium, https://fys.kuleuven.be/ster/pub/phd-thesis-dylan-bollen/phd-thesis-dylan-bollen

* Bollen, D., Van Winckel, H., & Kamath, D. 2017, A&A, 607, A60

* Bollen, D., Kamath, D., Van Winckel, H., & De Marco, O. 2019, A&A, 631, A53

* Bollen, D., Kamath, D., De Marco, O., Van Winckel, H., & Wardle, M. 2020, A&A, 641, A175

* Bollen, D., Kamath, D., Van Winckel, H., et al. 2022, A&A, 666, A40

* Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., & Van Winckel, H. 2015, A&A, 575, L7

* Casse, F., & Ferreira, J. 2000a, A&A, 353, 1115

* Casse, F., & Ferreira, J. 2000b, A&A, 361, 1178

* Demircan, O., & Kahraman, G. 1991, Astrophys. Space Sci., 181, 313

* Eggleton, P. P. 1983, ApJ, 268, 368

* Ferreira, J. 1997, A&A, 319, 340

* Ferreira, J., & Casse, F. 2013, MNRAS, 428, 307

* Ferreira, J., & Pelletier, G. 1995, A&A, 295, 807

* Ferreira, J., Petrucci, P. O., Henri, G., Saugé, L., & Pelletier, G. 2006, A&A, 447, 813

* Gorlova, N., Van Winckel, H., Gielen, C., et al. 2012, A&A, 542, A27

* Hillen, M., Menu, J., Van Winckel, H., et al. 2014, A&A, 568, A12

* Hillen, M., Kluska, J., Bouquin, J.-B. L., et al. 2016, A&A, 588, L1

* Ireland, L. G., Matt, S. P., & Zanni, C. 2022, ApJ, 929, 65

* Jacquemin-Ide, J., Ferreira, J., & Lesur, G. 2019, MNRAS, 490, 3112

* Jacquemin-Ide, J., Lesur, G., & Ferreira, J. 2021, A&A, 647, A192

* Kamath, D., Wood, P. R., & Van Winckel, H. 2015, MNRAS, 454, 1468

* Kluska, J., Hillen, M., Van Winckel, H., et al. 2018, A&A, 616, A153

* Kluska, J., Van Winckel, H., Coppée, Q., et al. 2022, A&A, 658, A36

* Manick, R., Van Winckel, H., Kamath, D., Hillen, M., & Escorza, A. 2017, A&A, 597, A129

* Oomen, G.-M., Van Winckel, H., Pols, O., et al. 2018, A&A, 620, A85

* Oomen, G.-M., Van Winckel, H., Pols, O., & Nelemans, G. 2019, A&A, 629, A49

* Panoglou, D., Cabrit, S., Forêts, G. P. D., et al. 2012, A&A, 538, A2

* Raskin, G. 2011, PhD Thesis, KU Leuven, Belgium, https://fys.kuleuven.be/ster/pub/phd-thesis-gert-raskin/phd-thesis-gert-raskin-1

* Raskin, G., Van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69

* Van Winckel, H. 2003, ARA&A, 41, 391

* Van Winckel, H. 2007, Balt. Astron., 16, 112

* Van Winckel, H. 2018, Camb. Astrophys. Ser., 54, 92

* Zanni, C., & Ferreira, J. 2009, A&A, 508, 1117