pith. sign in

arxiv: 2410.11692 · v2 · submitted 2024-10-15 · 🌌 astro-ph.SR · astro-ph.GA· astro-ph.HE

Simulations of Protostar-Driven Photoionization in Herbig-Haro Jets

Pith reviewed 2026-05-23 18:53 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GAastro-ph.HE
keywords Herbig-Haro jetsphotoionizationX-ray emissionMHD simulationsprotostarsionization fractionT-Tauri starsjet launching
0
0 comments X

The pith

X-ray photoionization from the central star ionizes 10 to 20 percent of the jet close to the protostar.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Young stellar object jets require a substantial pre-shock ionization fraction to produce the line emissions seen in shocks. This work tests whether X-rays emitted near the central star can supply that ionization through photoionization that is then carried outward. Axisymmetric MHD simulations add photoionization and optically thin cooling to the equations and track the resulting ionization levels in the launched jet. The results show that typical T-Tauri X-ray luminosities achieve 10 to 20 percent ionization near the star.

Core claim

The simulations demonstrate that for typical X-ray luminosities in classical T-Tauri stars, photoionization is responsible for ionizing 10 to 20 percent of the jet close to the star. This ionized material can be advected at large distances along the jet.

What carries the argument

Axisymmetric MHD jet launching simulations that incorporate X-ray photoionization and optically thin cooling terms.

If this is right

  • The pre-ionized gas supplies the ionization fraction needed to explain observed shock line emissions in Herbig-Haro objects.
  • Ionization occurs near the launch region and is carried outward by the flow.
  • The effect appears for standard T-Tauri X-ray output without requiring extreme stellar activity.
  • Optically thin cooling is included alongside photoionization in the MHD equations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Stellar X-ray variability could produce time-dependent changes in jet ionization and emission properties.
  • This pre-ionization channel may lessen reliance on other mechanisms such as cosmic-ray ionization in jet models.
  • Three-dimensional extensions could reveal whether ionization is uniform across the jet cross-section.

Load-bearing premise

X-ray emission originates in the immediate vicinity of the central source and the added photoionization and cooling terms accurately capture the microphysics without dominant numerical or missing-physics errors.

What would settle it

Spatially resolved measurement of the ionization fraction within a few hundred AU of a classical T-Tauri star whose X-ray luminosity is known; values consistently outside the 10-20 percent range for typical luminosities would challenge the result.

Figures

Figures reproduced from arXiv: 2410.11692 by A. Bouldjderi, A. Mignone, C. Zanni, S. Massaglia, Z. Ahmane.

Figure 1
Figure 1. Figure 1: Left panel: Jet integrals along the magnetic field line rooted in the innermost disk. The jet specific angular momentum Ω, the mass load k, the field angular velocity L, the jet specific energy E, and the specific entropy Q. Central panel: Jet specific energy contributions along the vertical direction z. Different energy components are indicated by colors: kinetic (green), magnetic (red), gravitational (ma… view at source ↗
Figure 2
Figure 2. Figure 2: From top to bottom, we show four snapshots of the temporal evolution corre￾sponding to t = 25, 100, 175 and 250 yrs. In the first panel (A) shows the logarithmic maps of the temperature ( in units of 104 K) for the adiabatic case. The second three panels (B, C, D) compares the temperature maps when cooling is included (upper half) and the adiabatic case (lower half). The jet propagates from left to right. … view at source ↗
Figure 3
Figure 3. Figure 3: Top panel: velocity along the z direction in km/s profile - adiabatic and cooling cases, - at different heights. Bottom panel: temperature profile (in K) - adiabatic and cooling cases, - at different heights, corresponding to t = 175 yrs. The effect of energy losses on the jet dynamics is shown in [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Top panel: From top to bottom, we show two snapshots of the logarithmic ratio of the absolute value of the toroidal magnetic field over the poloidal, the temporal evolution corresponding to t = 125 and 250 yrs. The solid lines in (red, black and green) indicate the critical surfaces i.e. Alfv´enic, fast magnetosonic and slow magnetosonic respectively. Bottom panel: Zoomed image of the jet base, with logari… view at source ↗
Figure 5
Figure 5. Figure 5: SNEq (Simplified Non Equilibrium Cooling): Ionization fraction of the jet material when only cooling is included (no heating). Snapshots are shown at t = 100, 150 and t = 250 yrs. In the first case, only ∼ 10% of the hydrogen becomes ionized close to the central object and recombination takes place as the flow propagates outward, eventually lowering this value to less than 4% at large distances. This situa… view at source ↗
Figure 6
Figure 6. Figure 6: (SNEq+Heating) Bottom panel: Maps of the total ionization fraction of the jet material for different luminosity LX = 1028 erg/s, LX = 1032 erg/s at different times t = 100, 150 and t = 250 yrs. Top panel: The corresponding zoomed image of total hydrogen ionization fraction at time t = 100 yrs. 0 20 40 60 80 100 120 R(AU) 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 HII Ionisation fraction (t = 250yrs) LX =… view at source ↗
Figure 7
Figure 7. Figure 7: Ionisation fraction profiles for different luminosity LX = 1028 erg/s, LX = 1032 erg/s at different heights, corresponding to t = 250 yrs. ionization degree, lowering to ∼ 10% at large distances (see the zoomed image below). The ionization area is now also much larger extending to R = 120 AU with 10% of ionisation fraction, reaching the outer radial boundary at z = 2000 AU see ( [PITH_FULL_IMAGE:figures/f… view at source ↗
Figure 8
Figure 8. Figure 8: Left panel: Zoomed logarithmic 2D map of the cooling Λ [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
read the original abstract

Recent studies showed that observations of line emission from shocks in YSO jets require a substantial amount of ionization of the pre-shock matter. Photoionization from X-ray emitted close to the central source may be responsible for the initial ionization fraction. The aim of our work is to study the effect of X-ray photoionization, coming from the vicinity of the central star, on the ionization fraction inside the jet that can be advected at large distances. For this purpose we have performed axisymmetric MHD jet launching simulations including photoionization and optically thin losses using PLUTO. For typical X-ray luminosities in classical T-Tauri stars, we see that the photoionization is responsible for ionizing to 10 % -20 % the jet close to the star.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper presents axisymmetric MHD simulations of jet launching from young stellar objects using the PLUTO code, with added source terms for X-ray photoionization from the central protostar and optically thin radiative cooling. The central claim is that, for typical X-ray luminosities of classical T-Tauri stars, photoionization produces an ionization fraction of 10-20% in the jet material close to the star that can be advected outward.

Significance. If the added microphysical terms are shown to be accurate, the result would supply a concrete mechanism linking stellar X-ray emission to the pre-shock ionization fractions required by observations of line emission in Herbig-Haro jets. The self-consistent inclusion of MHD jet launching with photoionization is a positive step toward connecting stellar and jet physics.

major comments (2)
  1. [Abstract] Abstract: the headline 10-20% ionization result is stated without any reported convergence tests, resolution studies, or error analysis on the ionization fraction in the launching region; because the result is generated entirely by the added photoionization/cooling source terms, this validation is load-bearing for the central claim.
  2. [Numerical methods] Numerical methods (implementation of source terms): the manuscript provides no benchmark or test problem demonstrating that the X-ray propagation (attenuation vs. unattenuated 1/r^2), spectrum, and microphysical rates (recombination, cooling) reproduce known analytic or tabulated solutions; without this, it is unclear whether numerical advection or missing processes dominate the reported ionization fraction.
minor comments (1)
  1. [Abstract] Abstract: the phrasing 'ionizing to 10 % -20 % the jet' is grammatically awkward and should be revised to 'ionizing the jet to 10-20%'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We agree that additional validation is required to support the central claims and will revise the paper accordingly.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the headline 10-20% ionization result is stated without any reported convergence tests, resolution studies, or error analysis on the ionization fraction in the launching region; because the result is generated entirely by the added photoionization/cooling source terms, this validation is load-bearing for the central claim.

    Authors: We concur that the absence of convergence tests and error analysis weakens the presentation of the 10-20% ionization result. In the revised manuscript, we will conduct and report resolution studies focused on the launching region, including the dependence of the ionization fraction on grid resolution. Error bars or ranges will be provided based on these tests, and the abstract will be updated to reflect this validation. revision: yes

  2. Referee: [Numerical methods] Numerical methods (implementation of source terms): the manuscript provides no benchmark or test problem demonstrating that the X-ray propagation (attenuation vs. unattenuated 1/r^2), spectrum, and microphysical rates (recombination, cooling) reproduce known analytic or tabulated solutions; without this, it is unclear whether numerical advection or missing processes dominate the reported ionization fraction.

    Authors: The referee correctly identifies the lack of benchmark tests for the photoionization and cooling source terms. We will add a new section or subsection detailing benchmark problems. These will verify X-ray attenuation against analytic expectations, ionization equilibrium solutions, and cooling functions against standard tables. This will confirm the implementation before presenting the jet simulation results. revision: yes

Circularity Check

0 steps flagged

No circularity; simulation outputs are independent numerical results

full rationale

The paper performs axisymmetric MHD simulations in PLUTO with added photoionization and optically thin cooling source terms. The central claim (10-20% ionization fraction near the star for typical T-Tauri X-ray luminosities) is reported as a direct output of the numerical integration under the stated assumptions. No equations or steps reduce this fraction to a fitted parameter, self-defined quantity, or self-citation chain by construction. The work is a forward simulation whose result is not equivalent to its inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

Central claim rests on choice of a representative X-ray luminosity drawn from T-Tauri observations and on the standard MHD framework with added source terms for photoionization and cooling. No invented entities. Limited detail available from abstract.

free parameters (1)
  • X-ray luminosity = typical for classical T-Tauri stars
    Ionization fraction depends on adopting a typical observed value for classical T-Tauri stars; exact numerical value not stated in abstract.
axioms (2)
  • standard math Axisymmetric ideal MHD equations govern jet launching and propagation
    Invoked as the base dynamical model for the simulations.
  • domain assumption Photoionization and optically thin radiative losses can be treated as local source terms
    Standard modeling choice in astrophysical plasma codes.

pith-pipeline@v0.9.0 · 5676 in / 1290 out tokens · 29121 ms · 2026-05-23T18:53:15.133860+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

75 extracted references · 75 canonical work pages · 1 internal anchor

  1. [1]

    and Eisl \"o ffel, J

    Bacciotti, F. and Eisl \"o ffel, J. (1999). Ionization and density along the beams of herbig-haro jets. A&A , 342:717--735

  2. [2]

    P., Mundt, R., Eisl \"o ffel, J., and Solf, J

    Bacciotti, F., Ray, T. P., Mundt, R., Eisl \"o ffel, J., and Solf, J. (2002). Hubble space telescope/stis spectroscopy of the optical outflow from dg tauri: Indications for rotation in the initial jet channel. ApJ , 576(1):222

  3. [3]

    Bally, J., Feigelson, E., and Reipurth, B. (2003). X-rays from the vicinity of the protostar l1551 irs 5: Reflection or fast shocks? ApJ , 584(2):843

  4. [4]

    Bally, J., Reipurth, B., and Davis, C. J. (2007). Observations of jets and outflows from young stars. Protostars and Planets V , pages 215--230

  5. [5]

    and Payne, D

    Blandford, R. and Payne, D. (1982). Hydromagnetic flows from accretion discs and the production of radio jets. MNRAS , 199(4):883--903

  6. [6]

    E., and Strom, K

    Cabrit, S., Edwards, S., Strom, S. E., and Strom, K. M. (1990). Forbidden-line emission and infrared excesses in t tauri stars-evidence for accretion-driven mass loss? ApJ , 354:687--700

  7. [7]

    J., Shang, H., Lin, H.-H., and Shu, F

    Cai, M. J., Shang, H., Lin, H.-H., and Shu, F. H. (2008). X-winds in action. ApJ , 672(1):489

  8. [8]

    and Keppens, R

    Casse, F. and Keppens, R. (2002). Magnetized accretion-ejection structures: 2.5-dimensional magnetohydrodynamic simulations of continuous ideal jet launching from resistive accretion disks. ApJ , 581(2):988

  9. [9]

    and Raymond, J

    Cox, D. and Raymond, J. (1985). Preionization-dependent families of radiative shock waves. ApJ , 298:651--659

  10. [10]

    C., and Esquivel, A

    De Colle, F., Raga, A. C., and Esquivel, A. (2008). The dynamics of internal working surfaces in magnetohydrodynamic jets. ApJ , 689(1):302

  11. [11]

    Dopita, A. (1978). The herbig-haro objects in the gum nebula. A&A , 63:237--241

  12. [12]

    Dougados, C., Cabrit, S., Lavalley, C., and M \'e nard, F. (2000). T tauri stars microjets resolved by adaptive optics. A&A , 357:L61--L64

  13. [13]

    and Mundt, R

    Eisloeffel, J. and Mundt, R. (1992). Proper motion measurements in the hh34 jet and its associated bow shocks. A&A , 263:292--300

  14. [14]

    Eisloeffel, J., Mundt, R., and Bohm, K.-H. (1994). Structure and proper motions in herbig-haro objects 1 and 2. ApJ , 108:1042--1055

  15. [15]

    Favata, F., Fridlund, C., Micela, G., Sciortino, S., and Kaas, A. (2002). Discovery of x-ray emission from the protostellar jet l1551 irs5 (hh 154). A&A , 386(1):204--210

  16. [16]

    Fendt, C. (2006). Collimation of astrophysical jets: The role of the accretion disk magnetic field distribution. ApJ , 651(1):272

  17. [17]

    Fendt, C. (2009). Formation of protostellar jets as two-component outflows from star-disk magnetospheres. ApJ , 692(1):346

  18. [18]

    Numerical simulation of surface brightness of astrophysical jets

    Gardner, C. L., Jones, J. R., and Vargas, P. B. (2017). Numerical simulation of surface brightness of astrophysical jets. arXiv preprint arXiv:1706.03819

  19. [19]

    M., Bacciotti, F., Fedele, D., Frasca, A., Harutyunyan, A., Munari, U., Rigliaco, E., and Vitali, F

    Giannini, T., Nisini, B., Antoniucci, S., Biazzo, K., Alcal \'a , J. M., Bacciotti, F., Fedele, D., Frasca, A., Harutyunyan, A., Munari, U., Rigliaco, E., and Vitali, F. (2019). GIARPS High-Resolution Observations of T Tauri Stars (GHOsT). I. Jet Line Emission . Astronomy & Astrophysics , 631:A44

  20. [20]

    x-ray ionization of protoplanetary disks

    Glassgold, A., Najita, J., and Igea, J. (1997). Erratum:" x-ray ionization of protoplanetary disks"(apj, 480, 344 [1997]). ApJ , 485(2):920

  21. [21]

    E., Stone , J

    Hardee , P. E., Stone , J. M., and Xu , J. (1997). Theory and Simulation of Asymmetrically Perturbed Radiatively Cooled Jets . In Wickramasinghe , D. T., Bicknell , G. V., and Ferrario , L., editors, IAU Colloq. 163: Accretion Phenomena and Related Outflows , volume 121 of ASP Conference Series , page 575

  22. [22]

    Hartigan, P., Edwards, S., and Ghandour, L. (1995). Disk accretion and mass loss from young stars. ApJ , 452:736

  23. [23]

    A., and Raymond, J

    Hartigan, P., Morse, J. A., and Raymond, J. (1994). Mass-loss rates, ionization fractions, shock velocities, and magnetic fields of stellar jets. ApJ , 436:125--143

  24. [24]

    Lada, C. J. (1985). Cold outflows, energetic winds, and enigmatic jets around young stellar objects. Annu. Rev. Astron. Astrophys , 23(1):267--317

  25. [25]

    Lavalley, C., Cabrit, S., Dougados, C., Ferruit, P., and Bacon, R. (1997). Sub-arcsecond morphology and kinematics of the dg tauri jet in the [oi] lambda6300 line. A&A , 327:671--680

  26. [26]

    Lavalley-Fouquet, C., Cabrit, S., and Dougados, C. (2000). Dg tau: A shocking jet. A&A , 356:L41--L44

  27. [27]

    Livio, M. (1999). Astrophysical jets: a phenomenological examination of acceleration and collimation. Physics Reports , 311(3-5):225--245

  28. [28]

    Massaglia, S., Mignone, A., and Bodo, G. (2005). Time-dependent mhd shocks and line emission: the case of the dg tau jet. A&A , 442(2):549--554

  29. [29]

    Massaglia , S., Trussoni , E., Bodo , G., Rossi , P., and Ferrari , A. (1992). Radiative unstable modes in the jets of young stellar objects. A&A , 260:243--249

  30. [30]

    Matsakos, T., Massaglia, S., Trussoni, E., Tsinganos, K., Vlahakis, N., Sauty, C., and Mignone, A. (2009). Two-component jet simulations. ii. combining analytical disk and stellar mhd outflow solutions. A&A , 502(1):217--229

  31. [31]

    Matsakos, T., Tsinganos, K., Vlahakis, N., Massaglia, S., Mignone, A., and Trussoni, E. (2008). Two-component jet simulations. i. topological stability of analytical mhd outflow solutions. A&A , 477(2):521--533

  32. [32]

    P., Massaglia , S., Trussoni , E., and Mignone , A

    Matsakos , T., Vlahakis , N., Tsinganos , K., Karampelas , K., Sauty , C., Cayatte , V., Matt , S. P., Massaglia , S., Trussoni , E., and Mignone , A. (2012). Velocity asymmetries in young stellar object jets. Intrinsic and extrinsic mechanisms . A&A , 545:A53

  33. [33]

    Maurri, L., Bacciotti, F., Podio, L., Eisl \"o ffel, J., Ray, T., Mundt, R., Locatelli, U., and Coffey, D. (2014). Physical properties of the jet from dg tauri on sub-arcsecond scales with hst/stis. A&A , 565:A110

  34. [34]

    Meliani, Z., Casse, F., and Sauty, C. (2006). Two-component magnetohydrodynamical outflows around young stellar objects-interplay between stellar magnetospheric winds and disc-driven jets. A&A , 460(1):1--14

  35. [35]

    Y., Eisl \"o ffel, J., Bacciotti, F., Woitas, J., and Ray, T

    Melnikov, S. Y., Eisl \"o ffel, J., Bacciotti, F., Woitas, J., and Ray, T. P. (2009). Hst/stis observations of the rw aurigae bipolar jet: mapping the physical parameters close to the source. A&A , 506(2):763--777

  36. [36]

    Micono, M., Bodo, G., Massaglia, S., Rossi, P., Ferrari, A., and Rosner, R. (2000). Kelvin-helmholtz instability in three dimensional radiative jets. A&A , 360:795--808

  37. [37]

    Mignone, A. (2014). High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates. J. Comput. Phys. , 270:784--814

  38. [38]

    Mignone, A., Bodo, G., Massaglia, S., Matsakos, T., Tesileanu, O., Zanni, C., and Ferrari, A. (2007). Pluto: A numerical code for computational astrophysics. Astrophys. J. Suppl. Ser. , 170:228--242

  39. [39]

    and Tzeferacos, P

    Mignone, A. and Tzeferacos, P. (2010). A second-order unsplit godunov scheme for cell-centered mhd: The ctu-glm scheme. J. Comput. Phys. , 229(6):2117--2138

  40. [40]

    Mignone, A., Zanni, C., Tzeferacos, P., van Straalen, B., Colella, P., and Bodo, G. (2012). The pluto code for adaptive mesh computations in astrophysical fluid dynamics. ApJS , 198(1):7

  41. [41]

    D., and Rosen, A

    Moraghan, A., Smith, M. D., and Rosen, A. (2006). Velocity study of axisymmetric protostellar jets with molecular cooling. MNRAS , 371(3):1448--1458

  42. [42]

    Mundt, R. (1986). Jets from young stars: estimates of their physical parameters. Can.J.Phys. , 64(4):407--413

  43. [43]

    and Pudritz, R

    Ouyed, R. and Pudritz, R. E. (1997). Numerical simulations of astrophysical jets from keplerian disks. i. stationary models. ApJ , 482(2):712

  44. [44]

    and Cant \'o , J

    Raga, A. and Cant \'o , J. (1995). The initial stages of an hh jet/cloud core collision. RevMexAA , 31:51--61

  45. [45]

    Ramsey, J. P. and Clarke, D. A. (2011). Simulating protostellar jets simultaneously at launching and observational scales. ApJ L , 728(1):L11

  46. [46]

    P., Mundt, R., Dyson, J

    Ray, T. P., Mundt, R., Dyson, J. E., Falle, S. A., and Raga, A. C. (1996). Hubble space telescope observations of jets from young stars. Astrophys. J. Lett. , 468(2):L103

  47. [47]

    P., and Zealey, W

    Reipurth, B., Bally, J., Graham, J., Lane, A. P., and Zealey, W. (1986). The jet and energy source of hh 34. A&A , 164:51--66

  48. [48]

    Rossi, P., Bodo, G., Massaglia, S., and Ferrari, A. (1997). Evolution of kelvin-helmholtz instabilities in radiative jets. ii. shock structure and entrainment properties. A&A , 321:672--684

  49. [49]

    Rubini, F., Maurri, L., Inghirami, G., Bacciotti, F., and Del Zanna, L. (2014). Numerical simulations of stellar jets and comparison between synthetic and observed maps: clues to the launch mechanism. A&A , 567:A13

  50. [50]

    Sauty, C., Trussoni, E., and Tsinganos, K. (2002). Nonradial and nonpolytropic astrophysical outflows. v. acceleration and collimation of self-similar winds. A&A , 389:1068--1085

  51. [51]

    and Tsinganos , K

    Sauty , C. and Tsinganos , K. (1994). Nonradial and nonpolytropic astrophysical outflows III. A criterion for the transition from jets to winds . A&A , 287:893--926

  52. [52]

    Schwartz, R. (1975). T tauri nebulae and herbig-haro nebulae-evidence for excitation by a strong stellar wind. ApJ , 195:631--642

  53. [53]

    E., Shu, F

    Shang, H., Glassgold, A. E., Shu, F. H., and Lizano, S. (2002). Heating and ionization of x-winds. ApJ , 564(2):853

  54. [54]

    Shu, F., Najita, J., Ostriker, E., Wilkin, F., Ruden, S., and Lizano, S. (1994). Magnetocentrifugally driven flows from young stars and disks. 1: A generalized model. ApJ , 429:781--796

  55. [55]

    Shull, J. M. and van Steenberg, M. E. (1985). X-ray secondary heating and ionization in quasar emission-line clouds. ApJ , 298:268--274

  56. [56]

    L., Audard, M., and G \"u del, M

    Skinner, S. L., Audard, M., and G \"u del, M. (2011). Chandra evidence for extended x-ray structure in ry tau. ApJ , 737(1):19

  57. [57]

    M., Xu, J., and Hardee, P

    Stone, J. M., Xu, J., and Hardee, P. E. (1997). The stability of radiatively cooling jets. ii. nonlinear evolution. ApJ , 483(1):136

  58. [58]

    Te s ileanu, O., Massaglia, S., Mignone, A., Bodo, G., and Bacciotti, F. (2009). Time-dependent mhd shocks and line intensity ratios in the hh 30 jet: a focus on cooling function and numerical resolution. A&A , 507(2):581--588

  59. [59]

    Te s ileanu, O., Matsakos, T., Massaglia, S., Trussoni, E., Mignone, A., Vlahakis, N., Tsinganos, K., Stute, M., Cayatte, V., and Sauty, C. (2014). Young stellar object jet models: From theory to synthetic observations. A&A , 562

  60. [60]

    Te s ileanu, O., Mignone, A., and Massaglia, S. (2008). Simulating radiative astrophysical flows with the pluto code: a non-equilibrium, multi-species cooling function. A&A , 488(1):429--440

  61. [61]

    Te s ileanu, O., Mignone, A., Massaglia, S., and Bacciotti, F. (2012). Numerical simulations of radiative magnetized herbig--haro jets: The influence of pre-ionization from x-rays on emission lines. ApJ , 746(1):96

  62. [62]

    Tzeferacos, P., Ferrari, A., Mignone, A., Zanni, C., Bodo, G., and Massaglia, S. (2009). On the magnetization of jet-launching discs. MNRAS , 400(2):820--834

  63. [63]

    and Shibata, K

    Uchida, Y. and Shibata, K. (1985). Magnetodynamical acceleration of co and optical bipolar flows from the region of star formation. PASJ , 37:515--535

  64. [64]

    Vaidya, B., Mignone, A., Bodo, G., and Massaglia, S. (2015). Astrophysical fluid simulations of thermally ideal gases with non-constant adiabatic index: numerical implementation. A&A , 580

  65. [65]

    Zanni, C., Ferrari, A., Rosner, R., Bodo, G., and Massaglia, S. (2007). Mhd simulations of jet acceleration from keplerian accretion disks-the effects of disk resistivity. A&A , 469(3):811--828

  66. [66]

    Giannini , T., Nisini , B., Antoniucci , S., Biazzo , K., Alcal \'a , J., Bacciotti , F., Fedele , D., Frasca , A., Harutyunyan , A., Munari , U., Rigliaco , E., and Vitali , F. (2019). GIARPS High-resolution Observations of T Tauri stars (GHOsT). I. Jet line emission . , 631:A44

  67. [67]

    Matsakos , T., Massaglia , S., Trussoni , E., Tsinganos , K., Vlahakis , N., Sauty , C., and Mignone , A. (2009). Two-component jet simulations. II. Combining analytical disk and stellar MHD outflow solutions . A&A , 502(1):217--229

  68. [68]

    Matsakos , T., Tsinganos , K., Vlahakis , N., Massaglia , S., Mignone , A., and Trussoni , E. (2008). Two-component jet simulations. I. Topological stability of analytical MHD outflow solutions . A&A , 477(2):521--533

  69. [69]

    Micono , M., Bodo , G., Massaglia , S., Rossi , P., Ferrari , A., and Rosner , R. (2000). Kelvin-Helmholtz instability in three dimensional radiative jets . A&A , 360:795--808

  70. [70]

    Mignone , A. (2014). High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates . J. Comput. Phys. , 270:784--814

  71. [71]

    and Tzeferacos , P

    Mignone , A. and Tzeferacos , P. (2010). A second-order unsplit Godunov scheme for cell-centered MHD: The CTU-GLM scheme . J. Comput. Phys. , 229(6):2117--2138

  72. [72]

    Mignone , A., Zanni , C., Tzeferacos , P., van Straalen , B., Colella , P., and Bodo , G. (2012). The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics . , 198(1):7

  73. [73]

    Sauty , C., Trussoni , E., and Tsinganos , K. (2002). Nonradial and nonpolytropic astrophysical outflows. V. Acceleration and collimation of self-similar winds . A&A , 389:1068--1085

  74. [74]

    Te s ileanu , O., Matsakos , T., Massaglia , S., Trussoni , E., Mignone , A., Vlahakis , N., Tsinganos , K., Stute , M., Cayatte , V., and Sauty , C. (2014). Young stellar object jet models: From theory to synthetic observations . A&A , 562:A117

  75. [75]

    Vaidya , B., Mignone , A., Bodo , G., and Massaglia , S. (2015). Astrophysical fluid simulations of thermally ideal gases with non-constant adiabatic index: numerical implementation . A&A , 580:A110