Simulations of Protostar-Driven Photoionization in Herbig-Haro Jets
Pith reviewed 2026-05-23 18:53 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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
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
-
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
-
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
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
free parameters (1)
- X-ray luminosity =
typical for classical T-Tauri stars
axioms (2)
- standard math Axisymmetric ideal MHD equations govern jet launching and propagation
- domain assumption Photoionization and optically thin radiative losses can be treated as local source terms
Reference graph
Works this paper leans on
-
[1]
Bacciotti, F. and Eisl \"o ffel, J. (1999). Ionization and density along the beams of herbig-haro jets. A&A , 342:717--735
work page 1999
-
[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
work page 2002
-
[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
work page 2003
-
[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
work page 2007
-
[5]
Blandford, R. and Payne, D. (1982). Hydromagnetic flows from accretion discs and the production of radio jets. MNRAS , 199(4):883--903
work page 1982
-
[6]
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
work page 1990
-
[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
work page 2008
-
[8]
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
work page 2002
-
[9]
Cox, D. and Raymond, J. (1985). Preionization-dependent families of radiative shock waves. ApJ , 298:651--659
work page 1985
-
[10]
De Colle, F., Raga, A. C., and Esquivel, A. (2008). The dynamics of internal working surfaces in magnetohydrodynamic jets. ApJ , 689(1):302
work page 2008
-
[11]
Dopita, A. (1978). The herbig-haro objects in the gum nebula. A&A , 63:237--241
work page 1978
-
[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
work page 2000
-
[13]
Eisloeffel, J. and Mundt, R. (1992). Proper motion measurements in the hh34 jet and its associated bow shocks. A&A , 263:292--300
work page 1992
-
[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
work page 1994
-
[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
work page 2002
-
[16]
Fendt, C. (2006). Collimation of astrophysical jets: The role of the accretion disk magnetic field distribution. ApJ , 651(1):272
work page 2006
-
[17]
Fendt, C. (2009). Formation of protostellar jets as two-component outflows from star-disk magnetospheres. ApJ , 692(1):346
work page 2009
-
[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
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[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
work page 2019
-
[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
work page 1997
-
[21]
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
work page 1997
-
[22]
Hartigan, P., Edwards, S., and Ghandour, L. (1995). Disk accretion and mass loss from young stars. ApJ , 452:736
work page 1995
-
[23]
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
work page 1994
-
[24]
Lada, C. J. (1985). Cold outflows, energetic winds, and enigmatic jets around young stellar objects. Annu. Rev. Astron. Astrophys , 23(1):267--317
work page 1985
-
[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
work page 1997
-
[26]
Lavalley-Fouquet, C., Cabrit, S., and Dougados, C. (2000). Dg tau: A shocking jet. A&A , 356:L41--L44
work page 2000
-
[27]
Livio, M. (1999). Astrophysical jets: a phenomenological examination of acceleration and collimation. Physics Reports , 311(3-5):225--245
work page 1999
-
[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
work page 2005
-
[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
work page 1992
-
[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
work page 2009
-
[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
work page 2008
-
[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
work page 2012
-
[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
work page 2014
-
[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
work page 2006
-
[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
work page 2009
-
[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
work page 2000
-
[37]
Mignone, A. (2014). High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates. J. Comput. Phys. , 270:784--814
work page 2014
-
[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
work page 2007
-
[39]
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
work page 2010
-
[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
work page 2012
-
[41]
Moraghan, A., Smith, M. D., and Rosen, A. (2006). Velocity study of axisymmetric protostellar jets with molecular cooling. MNRAS , 371(3):1448--1458
work page 2006
-
[42]
Mundt, R. (1986). Jets from young stars: estimates of their physical parameters. Can.J.Phys. , 64(4):407--413
work page 1986
-
[43]
Ouyed, R. and Pudritz, R. E. (1997). Numerical simulations of astrophysical jets from keplerian disks. i. stationary models. ApJ , 482(2):712
work page 1997
-
[44]
Raga, A. and Cant \'o , J. (1995). The initial stages of an hh jet/cloud core collision. RevMexAA , 31:51--61
work page 1995
-
[45]
Ramsey, J. P. and Clarke, D. A. (2011). Simulating protostellar jets simultaneously at launching and observational scales. ApJ L , 728(1):L11
work page 2011
-
[46]
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
work page 1996
-
[47]
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
work page 1986
-
[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
work page 1997
-
[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
work page 2014
-
[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
work page 2002
-
[51]
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
work page 1994
-
[52]
Schwartz, R. (1975). T tauri nebulae and herbig-haro nebulae-evidence for excitation by a strong stellar wind. ApJ , 195:631--642
work page 1975
-
[53]
Shang, H., Glassgold, A. E., Shu, F. H., and Lizano, S. (2002). Heating and ionization of x-winds. ApJ , 564(2):853
work page 2002
-
[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
work page 1994
-
[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
work page 1985
-
[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
work page 2011
-
[57]
Stone, J. M., Xu, J., and Hardee, P. E. (1997). The stability of radiatively cooling jets. ii. nonlinear evolution. ApJ , 483(1):136
work page 1997
-
[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
work page 2009
-
[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
work page 2014
-
[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
work page 2008
-
[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
work page 2012
-
[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
work page 2009
-
[63]
Uchida, Y. and Shibata, K. (1985). Magnetodynamical acceleration of co and optical bipolar flows from the region of star formation. PASJ , 37:515--535
work page 1985
-
[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
work page 2015
-
[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
work page 2007
-
[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
work page 2019
-
[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
work page 2009
-
[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
work page 2008
-
[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
work page 2000
-
[70]
Mignone , A. (2014). High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates . J. Comput. Phys. , 270:784--814
work page 2014
-
[71]
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
work page 2010
-
[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
work page 2012
-
[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
work page 2002
-
[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
work page 2014
-
[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
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.