The gasphase abundance of carbon xSUBCSUB  CHSUBgasSUB  xSUBCSUPSUPSUB  xSUBCSUP0SUPSUB  xSUBCOSUB    and its depletion factors are essential parameters for understanding the gas and solid compositions that are ultimately incorporated into exoplanets. The majority of protoplanetary disks are born in clusters and as a result are exposed to external farultraviolet FUV radiation. These FUV photons potentially affect the disks evolution chemical composition and line excitation. We present the first detection of the C I 609 m finestructure SUP3SUPPSUB1SUBSUP3SUPPSUB0SUB line of neutral carbon CSUP0SUP achieved with ALMA toward one of these disks d203506 in the Orion Nebula Cluster. We also report the detection of C I forbidden and C I permitted lines from electronically excited states up to 10 eV observed with JWST in the nearinfrared NIR. These lines trace the irradiated outer disk and photoevaporative wind. Contrary to the common belief that these NIR lines are CSUPSUP recombination lines we find that they are dominated by FUVpumping of CSUP0SUP followed by fluorescence cascades. They trace the transition from atomic to molecular gas and their intensities scale with GSUB0SUB. The lack of outstanding NIR O I fluorescent emission however implies a sharper attenuation of external FUV radiation with E  12 eV   Lyman. This is related to a lower effective FUV dust absorption cross section compared to that of interstellar grains implying a more prominent role for FUV shielding by the CSUP0SUP photoionization continuum. The C I 609 m line intensity is proportional to NCSUP0SUP and can be used to infer xSUBCSUB. We derive xSUBCSUB  1.4  10SUP4SUP. This implies that there is no major depletion of volatile carbon compared to xSUBCSUB measured in the natal cloud hinting at a young disk. We also show that external FUV radiation impacts the outer disk and wind by vertically shifting the water freezeout depth which likely results in less efficient grain growth and settling. This shift leads to nearly solar gasphase CO abundance ratios in these irradiated layers.
['protoplanetary disks', 'ISM: abundances', 'photon-dominated region (PDR)', 'Astrophysics - Astrophysics of Galaxies', 'Astrophysics - Earth and Planetary Astrophysics', 'Astrophysics - Solar and Stellar Astrophysics']
PDRs4All X. ALMA and JWST detection of neutral carbon in the externally irradiated disk d203506 Undepleted gasphase carbon
{'No Header': 'L etter to the E ditor', 'PDRs4All. X. ALMA and JWST detection of neutral carbon in the externally irradiated disk d203-506: Undepleted gas-phase carbon': 'Javier R. Goicoechea 1 , J. Le Bourlot 2 , J. H. Black 3 , F. Alarcón 4 , E. A. Bergin 4 , O. Berné 5 , E. Bron 6 , A. Canin 5 , E. Chapillon 7 , R. Chown 8 , 9 , 10 , E. Dartois 11 , M. Gerin 6 , E. Habart 12 , T. J. Haworth 13 , C. Joblin 5 , O. Kannavou 12 , F. Le Petit 6 , T. Onaka 14 , E. Peeters 8 , 9 , 15 , J. Pety 6 , 7 , E. Roue ff 6 , A. Sidhu 8 , 9 , I. Schroetter 5 , B. Tabone 12 , A. G. G. M. Tielens 16 , 17 , B. Trahin 12 , D. Van De Putte 18 , 8 , 9 , S. Vicente 19 , and M. Zannese 12  \n(A ffi liations can be found after the references)  \nReceived 4 June 2024 / Accepted 7 August 2024', 'ABSTRACT': "The gas-phase abundance of carbon, x C = [C / H]gas = x C + + x C 0 + x CO + ... , and its depletion factors are essential parameters for understanding the gas and solid compositions that are ultimately incorporated into (exo)planets. The majority of protoplanetary disks are born in clusters and, as a result, are exposed to external far-ultraviolet (FUV) radiation. These FUV photons potentially a ff ect the disk's evolution, chemical composition, and line excitation. We present the first detection of the [C i ] 609 µ mfine-structure ( 3 P13 P0) line of neutral carbon (C 0 ), achieved with ALMA, toward one of these disks, d203-506, in the Orion Nebula Cluster. We also report the detection of [C i ] forbidden and C i permitted lines (from electronically excited energy levels up to ∼ 10 eV) observed with JWST in the near-infrared (NIR). These lines trace the irradiated outer disk and photo-evaporative wind. Contrary to the common belief that these NIR lines are C + recombination lines, we find that they are dominated by FUV-pumping of C 0 followed by fluorescence cascades. They trace the transition from atomic to molecular gas, and their intensities scale with G 0. The lack of outstanding NIR O i fluorescent emission, however, implies a sharper attenuation of external FUV radiation with E ≳ 12 eV ( λ ≲ Lymanβ ). This is related to a lower e ff ective FUV dust absorption cross section compared to that of interstellar grains, implying a more prominent role for FUV shielding by the C 0 photoionization continuum. The [C i ] 609 µ mline intensity is proportional to N (C 0 ) and can be used to infer x C. We derive x C ≃ 1.4 × 10 -4 . This implies that there is no major depletion of volatile carbon compared to x C measured in the natal cloud, hinting at a young disk. We also show that external FUV radiation impacts the outer disk and wind by vertically shifting the water freeze-out depth, which likely results in less e ffi cient grain growth and settling. This shift leads to nearly solar gas-phase C / O abundance ratios in these irradiated layers.  \nKey words. Protoplanetary disks - ISM: abundances - photon-dominated region (PDR)", '1. Introduction': "Tracing the carbon reservoir in protoplanetary disks is essential for understanding the variations in the gaseous elemental carbon-to-oxygen abundance 1 ratio, which connects the atmospheric composition of a planet with its formation site (e.g., Öberg et al. 2011; Madhusudhan 2012; Öberg et al. 2023; Miotello et al. 2019; Bosman et al. 2021; Mah et al. 2023). In addition, carbon plays a key role in the formation of the complex organics that serve as the foundation for prebiotic chemistry and can be incorporated into planetary systems as they form.  \nObservations show that the gas-phase CO abundance in 1-10 M yr old disks can be depleted by a factor of 10-100 compared to molecular clouds, even after accounting for freeze-out and photodissociation. This suggests that something happens to CO by ∼ 1Myr (Zhang et al. 2020; Bergner et al. 2020). Grain growth, facilitated by the formation of water-rich ice mantles, followed by settling to the midplane and radial drift may explain these depletion factors (Bergin et al. 2016; Krijt et al. 2020). However, the majority of observational studies target 'isolated' disks in nearby low-mass star-forming regions such as Taurus, Lupus, Ophiuchus, and Chameleon (at d ∼ 150 pc; e.g., ALMA Partnership et al. 2015; Öberg et al. 2021). Yet, most  \nlow-mass stars (and thus planets) are born in clusters that contain one or more high-mass stars (e.g., Lada & Lada 2003). These OB stars emit intense ultraviolet (UV) radiation that illuminates the disk populations that gradually emerge from the natal cloud (e.g., O'Dell et al. 1993). External UV illumination leads to disk mass-loss driven by photo-evaporation (Johnstone et al. 1998; Störzer & Hollenbach 1999; Winter & Haworth 2022) and potentially impacts the properties of planets forming within those disks (Winter et al. 2022; Qiao et al. 2023). However, due to the greater distance to clusters, spectroscopic studies are scarcer (Boyden & Eisner 2020, 2023; Maucó et al. 2023), and it is still a matter of debate whether environmental UV radiation modifies the chemistry of 'externally irradiated' disks (Walsh et al. 2013; Ramírez-Tannus et al. 2023; Díaz-Berríos et al. 2024).  \nThe d203-506 system, in the Orion Nebula Cluster (ONC; at d ∼ 400pc; see Habart et al. 2024), is a remarkable example of a disk around a low-mass star that is externally irradiated by far-ultraviolet (FUV) radiation ( ∼ 6 < E < 13.6 eV) with no evidence for the comet-shaped ionization front seen in proplyds (e.g., O'Dell et al. 1993; Ricci et al. 2008). This implies a shielding from ionizing extreme-ultraviolet (EUV) radiation ( E > 13.6 eV). The estimated FUV flux from the O-type stars θ 1 Ori C and θ 2 Ori A is ∼ 10 4 times higher than around the profusely studied isolated disks in Taurus-like clouds. Thus, d203-506 may be more representative of the initial conditions of the proto-Solar System disk (Fatuzzo & Adams 2008; Adams 2010; Bergin et al. 2023). This flared and nearly edge-on disk, located toward the  \nb)  \nFigure 2 shows the velocity-resolved [C i ] 609 µ m and C 18 O 3-2 line profiles observed by ALMA. The [C i ] 609 µ m line does not show the characteristic double-peak profile of an  \n<!-- image -->  \n<!-- image -->  \nFig. 1. Images of d203-506 (same field of view). a): C 18 O J = 3-2 and [C i ] 609 µ mintegrated line emission (spatially unresolved). Contours show 10 σ and 20 σ intensity levels for C 18 O, and 5 σ and 10 σ for [C i ] 609 µ m. Synthesized beam sizes are shown at the bottom corner. The apparent larger extent of the C 18 O emission is due to the larger beam and faint background emission. b): [C i ] 0.9853 µ memission from the photo-evaporative wind observed with NIRSpec (colored image). Orange contours show the 0.8 mm dust continuum observed by ALMA (from Berné et al. 2024).  \n<!-- image -->  \nOrion Bar (Champion et al. 2017, but see our Appendix A), was first detected by the Hubble Space Telescope in silhouette against the optical nebular background (Bally et al. 2000). More recent Keck, James Webb Space Telescope (JWST), and Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations unveiled the extended nature of the FUV-pumped vibrationally excited H2 (hereafter H ∗ 2 ) and [C i ] 0.8729 µ memissions (Habart et al. 2023; Haworth et al. 2023; Berné et al. 2024).  \nThe JWST PDRs4All team constrained the mass of the host star ( ∼ 0.3 M ⊙ ), the disk's mass ( ∼ 10 MJup ≃ 10 -2 M ⊙ ) and radius ( ∼ 100 au ≃ 0.25 '' ), and a high mass-loss rate due to external photo-evaporation of (0.1-4.6) × 10 -6 M ⊙ yr -1 (Berné et al. 2024). The infrared (IR) spectrum shows the presence of highly excited rotational lines of OH probing the ongoing photodissociation of a hidden reservoir of water vapor (Zannese et al. 2024). It also shows bright CH + and CH + 3 rovibrational emission (which was first detected in space toward this disk by Berné et al. 2023 and Changala et al. 2023). These ions form due to the high temperatures and the enhanced reactivity of H ∗ 2 with C + ions. Their abundances scale with the FUV flux (Goicoechea et al. 2019). These photochemical signatures di ff er from what is commonly observed in isolated disks (e.g., van Dishoeck et al. 2023; Kamp et al. 2023; Perotti et al. 2023; Grant et al. 2023), suggesting that the outer disk gas reservoir in d203-506 is significantly reprocessed by the external FUV field.  \nA powerful tracer of this photo-processing is the groundstate 3 P1 3 P0 fine-structure line of neutral carbon (C 0 ), the [C i ] 609 µ m line, which can be observed by the Atacama Large Millimeter / submillimeter Array (ALMA) at subarcsecond and sub-km s -1 resolutions. However, low-angular-resolution searches toward proplyds in the NGC 1977 cluster (illuminated by the B1V-type star 43 Ori, north of the ONC) yielded nondetections, which were interpreted as being due to either very low-mass disks or a depletion of carbon in the outer disk (Haworth et al. 2022). Here we present its first detection toward an externally irradiated disk, complemented with the JWST detection of electronically excited C 0 lines in the near-infrared (NIR).", '2. ALMA and JWST observations of d203-506': "We used ALMA to observe the protoplanetary disk d203-506, located at α (2000) = 5 h 35 m 20.32 s , δ (2000) = -5 · 25 ' 05.55 '' . These observations, a ∼ 40 '' × 40 '' mosaic using 47 ALMA 12 m antennas, are part of an imaging program of the Orion Bar (2021.1.01369.S, P.I.: J. R. Goicoechea). Here we present the detection of the [C i ] 609 µ m(492GHz, in band 8) and C 18 O J = 3-2  \n(329 GHz, in band 7) lines. We used correlators that provide ∼ 282kHz and ∼ 564kHz resolution, respectively. We binned all spectra to a common velocity resolution of 0.4 km s -1 . The total observation times with the ALMA 12 m array were ∼ 4.6 h (492 GHz) and 2.7 h (329 GHz). The final synthesized beams are 0.52 '' × 0.38 '' at position angle PA = 110 · (492 GHz) and 0.77 '' × 0.60 '' at PA = 64 · (329.3 GHz). The achieved rms noises are 25 mJy and 10 mJy per velocity channel, respectively. The complete data cubes and calibration strategy will be described in an accompanying paper that focuses on the Bar.  \nWe also observed d203-506 as part of the PDRs4All JWST Early Release Science program (ID # 1288; Berné et al. 2022). Our study focuses on NIRSpec's spectral cube from 0.97 to ∼ 5 µ m, observed with grating dispersers at R = λ/ ∆ λ ∼ 2700 resolution. The angular resolution is ∼ 0.1 '' ( ∼ 40au at the distance to Orion). Peeters et al. (2024) describe the data reduction in detail. We extracted a NIR spectrum in two apertures: one toward d203-506 (ON measurement) and the other near the disk, providing the OFF reference. The intrinsic spectrum of d203-506 is the ON f · OFF measurement (see Appendices B and C).", '3. Results': 'Figure 1 shows the detection of [C i ] 609 µ m emission in d203-506 as well as extended [C i ] 0.9853 µ mforbidden line emission, which connects the first electronic excited state 1 D2 with the ground state. To our knowledge, this is the first detection of these lines toward an externally photoevaporating disk 2 . Only the NIRSpec [C i ] 0.9853 µ mobservation spatially resolves the outer disk and photo-evaporative wind. The [C i ] 0.9853 µ memission is similar to that of H ∗ 2 (see Fig. B.1), and significantly more extended than the submillimeter dust continuum emission from the inner, dense Keplerian disk (orange contours in Fig. 1b).  \n1  \n[CI]  \nD  \n-  \n2  \n3  \nP  \n0.9853  \n𝝁  \nm  \n2  \nFig. 2. Neutral carbon toward d203-506. Upper: [C i ] 609 µ m and C 18 O J = 3-2 line profiles. The gray-shaded spectrum shows HCO + J = 4-3 (Berné et al. 2024) scaled and convolved to the angular resolution of the [C i ] 609 µ mobservations. Lower: NIR detection of [C i ] forbidden lines at 0.9827 and 0.9853 µ m (top) and C i permitted line multiplets at ∼ 1.069, ∼ 1.176, and ∼ 1.355 µ m (bottom). Magenta lines represent the position and relative local thermodynamic equilibrium (LTE) intensity of each component.  \n<!-- image -->  \nedge-on rotating disk. Instead, it is similar to that of HCO + 4-3 (convolved to the same angular resolution), which is dominated by emission from the outer disk and the quasi-spherical photoevaporative wind (Berné et al. 2024). Since our observation does not spatially resolve d203-506, the resulting line profile is more or less Gaussian. The remaining blueshifted [C i ] 609 µ m emission may originate from a slow wind not traced by HCO + . A two-Gaussian fit to the two [C i ] 609 µ m components shows a velocity centroid di ff erence of only 2.3 ± 0.4 km s -1 , which is inconsistent with emission from the fast jet that emerges from the object (e.g., Haworth et al. 2023). In contrast, C 18 O 3-2 shows a clear double-peak profile, which traces deeper layers (we infer N H ≳ 10 23 cm -2 ) of the Keplerian disk (Miotello et al. 2016).  \nFigure 2 also shows NIR excited lines from neutral carbon: [C i ] (2 p 2 ) 1 D2 - (2 p 2 ) 3 P forbidden lines at 0.9827 and 0.9853 µ m ( E u ≃ 1.3 eV), and the highly excited (3 p ) 3 D-(3 s ) 3 P o , (3 d ) 3 F o - (3 p ) 3 D, and (4 s ) 3 P o - (3 p ) 3 S permitted line multiplets at ∼ 1.069, ∼ 1.176, and ∼ 1.355 µ m, respectively (with E u ≃ 10 eV). In addition, Haworth et al. (2023) imaged the [C i ] (2 p ) 1 S 0 - (2 p ) 1 D 2 forbidden line emission at 0.8729 µ m( E u ≃ 2.7 eV) with VLT / MUSE. Table I.1 summarizes the observed line intensities. As we show here, all these lines arise  \nFig. 3. Models with G 0 = 2 × 10 4 and di ff erent values of n H and x C . Red and blue markers show the predicted C 0 column density and [C i ] 609 µ m intensity, respectively, from AV = 0 to 10 mag into the wind and disk system. The gas density derived from H2 observations is ≈ 10 7 cm -3 (Berné et al. 2024). The horizontal line marks the observed intensity ( ± 1 σ ). Circles and triangles represent models with small depletion factors of two and three, respectively.  \n<!-- image -->  \nfrom neutral photon-dominated region (PDR) gas (not from H + ionized gas), which is consistent with the low ionization potential (IP) of C 0 (11.26 eV) and with the lack of cometary ionization fronts in d203-506.', '4. Analysis: Thermochemistry and line excitation': 'Neutral carbon lines are expected to trace warm molecular gas at the C + / C 0 / CO transition zone of the outer disk and inner photoevaporative wind (e.g., Haworth & Owen 2020). At the high gas densities relative to the FUV flux in a disk ( n H / G 0 > 10 2 cm -3 ), H2 and CO line self-shielding shifts the atomic to molecular gas transitions close to the irradiated surface, where the FUV flux is strong and the H2-emitting gas is heated to T k ≳ 1000K (Johnstone et al. 1998; Störzer & Hollenbach 1999; Champion et al. 2017). To derive the beam-averaged column density of C 0 in this zone, N (C 0 ), and the (total) gas-phase abundance of carbon, x C = x C + + x C 0 + x CO + ... , we used the Meudon code (Le Petit et al. 2006), a fully benchmarked PDR model (see Röllig et al. 2007). We simulated the [C i ] 609 µ m- and NIR-emitting zones as a 1D stationary slab of constant-density gas. This model involves solving the λ -dependent attenuation of external FUV photons (considering dust extinction and line self-shielding; Goicoechea & Le Bourlot 2007), C 0 photoionization, C + recombination, and CO photodissociation. Our detailed treatment of the penetration of FUV radiation, H2 excitation, thermal balance, and chemistry allows a precise determination of the H / H2 and C + / C / CO transition layers in the slab, but neglects the wind dynamics (e.g., Haworth & Owen 2020). We implicitly assumed that the contribution to the [C i ] 609 µ mintensity from the inner disk close to the central star (Gressel et al. 2020) is negligible and beam-diluted. We treated the λ -dependent absorption and scattering of FUV photons by dust, providing an e ff ective absorption cross section of 7 × 10 -22 cm 2 H -1 at 1000 Å. These grains are bigger than interstellar medium (ISM) grains but consistent with the modest grain growth expected in the upper layers of irradiated disks (Störzer & Hollenbach 1999; Birnstiel et al. 2018).  \nThe [C i ] 609 and 370 µ mlines involve forbidden transitions with very low Einstein coe ffi cients for spontaneous emission. In dense and warm gas, and because of their low excitation requirements - low critical density ( n cr ≲ 10 3 cm -3 ) and low level energy separation compared to T k ( ∆ E / k B = 23.6 K ≪ T k) - the [C i ] 609 µ memission is optically thin 3 , collisionally excited, and thermalized ( T ex = T k). In this regime, I ([C i ] 609 µ m) is proportional to N (C 0 ) irrespective of the physical conditions. In addition, N (C 0 ) is not very sensitive to variations in G 0 (Appendix G). Figure 3 shows the N (C 0 ) and I ([C i ] 609 µ m), integrating up to an equivalent extinction depth of AV = 10mag, predicted by models of varying gas densities, with n H = n (H) + 2 n (H2), and G 0 = 2 × 10 4 . Our models encompass the [C i ] 609 µ m-emitting layers and thus can be used, for a given n H, to constrain N (C 0 ) and x C. Figure 3 shows that I ([C i ] 609 µ m), a surface tracer, scales with N (C 0 ) and x C. Models marked with squares have x C = 1.4 × 10 -4 , which is the gas-phase abundance of interstellar carbon measured 4 in the line of sight toward the Trapezium stars (Sofia et al. 2004). Our main result is that the observed I ([C i ] 609 µ m) implies high x C abundances, roughly the same as those in the natal cloud. The lack of significant volatile carbon depletion, which would evolve with time (e.g., Kama et al. 2016; Krijt et al. 2020), suggests that this protoplanetary disk is young.  \nWe adopted n H ≈ 10 7 cm -3 as the gas density in the [C i ] 609 µ m-emitting zone, as determined from H2 observations with JWST (Berné et al. 2024). At these densities, the spatial scale of the [C i ] 609 µ memission is a few tens of au, a fraction of an arcsecond, in agreement with the observed size of d203-506. Models with n H ≃ 10 7 cm -3 fit the observed I ([C i ] 609 µ m) with N (C 0 ) ≃ 2 × 10 17 cm -2 and x C ≃ 1.4 × 10 -4 , with an uncertainty of less than a factor of 2. Figure 4 (upper panel) shows the predicted structure of the [C i ] 609 µ m-emitting layer, roughly representing the vertical structure of the outer disk and inner wind. The middle panel shows the normalized [C i ] 609 µ m, H2 1-0 S (1), and CO 3-2 line emissivities. [C i ] 609 µ mpeaks just beyond the H2 dissociation front. The lower panel shows the C + , C 0 , CO, water ice, and CO ice abundance profiles. The predicted gas temperature in the H2-emitting zone is T k ≳ 1000-1200 K, consistent with the observed H2 rotational temperatures (Berné et al. 2024), and drops to T k ≃ 700 K at the [C i ] 609 µ mintensity peak, where C 0 becomes the most abundant carbon species.  \nTurning back to the NIR carbon lines, we note that the (3 p ) 3 D-(3 s ) 3 P o (1.069 µ m) and (3 d ) 3 F-(3 p ) 3 D (1.176 µ m) multiplets are the brightest permitted lines predicted by recombination theory (C + + e -→ C ∗ + line cascade; Escalante & Victor 1990). However, the observed I (1.069 µ m) / I (1.176 µ m) and I (0.984 µ m) / I (1.069 µ m) intensity ratios do not match the ratios predicted by this theory (see Appendix D). In analogy with the NIR O i fluorescent lines detected very close behind the ionization fronts of interstellar PDRs (Walmsley et al. 2000; Lucy 2002; Peeters et al. 2024), the NIR C 0 lines could form via de-excitation cascades following FUV-pumping from the ground-state 3 P to the high-energy triplets (3 s ) 3 P o , (3 d ) 3 D o , (3 d ) 3 F o , (4 s ) 3 P o , and so on (i.e., C + h ν 0 → C ∗ + line cascade). The FUV-pumping lines  \nFig. 4. Model of the outer disk and inner wind ( n H = 10 7 cm -3 and G 0 = 2 × 10 4 ). Upper panel : H2 and H densities and gas and dust temperatures as a function of depth into the disk. Middle : Normalized line emissivity profiles, including NIR C 0 lines. Lower : Abundance profiles.  \n<!-- image -->  \nof these transitions lie at 1656 Å, 1277 Å, 1279 Å, and 1280 Å, respectively, that is, FUV photons with energies h ν 0 = 7.5-9.6 eV (see Fig. E.1 for a Gotrian diagram). Here we included FUVpumping followed by radiative cascades in the Meudon PDR code (Appendix E). This mechanism greatly increases the intensities of all the observed NIR C 0 lines (both forbidden and permitted) compared to models that consider collisional excitation and C + recombination alone, thereby fitting the observations (Table I.1). NIR C 0 line intensities scale with G 0 (see Appendix G) and are a powerful diagnostic of external FUV irradiation (we derive G 0 ≃ (1-2) × 10 4 in d203-506). Figure 4 shows their emissivity profiles, which peak very close to the H2 1-0 S (1) emission, as observed. Thus, NIR C 0 lines trace the transition from neutral atomic gas to molecular gas. The lack of bright NIR O i fluorescent lines (other than a faint O i 1.317 µ memission; Appendix C) implies that d203-506 is pervaded by a reduced flux of FUV photons at the pumping transitions of these lines: photons with E ≳ 12 eV ( λ ≲ Lymanβ , Walmsley et al. 2000). Furthermore, grain sizes and gas-to-dust mass ratios in the outer disk and wind are likely larger than those in the ISM (e.g., Winter & Haworth 2022). This results in lower e ff ective FUV dust absorption cross sections, and implies that the attenuation of these FUV photons is dominated by H i and H2 absorption lines, and by the C 0 photoionization continuum at λ< 1101 Å (Appendix F). In contrast, the attenuation of FUV photons at the pumping transitions of C 0 ( λ ≃ 1300 Å) is less pronounced and is dominated by dust grains. The lack of outstanding oxygen fluorescent lines but detectable carbon lines in d203-506 agrees with this scenario.', '5. Discussion: Photo-processing in d203-506': "The [C i ] 609 µ m line fluxes observed in nearby isolated disks (i.e., disks that are not externally irradiated) such as TW Hya suggest carbon depletion by 1-2 orders of magnitude in these disks (e.g., Kama et al. 2016). In d203-506, however, we determine a high gas-phase x C abundance. Furthermore, the outer layers of externally FUV-irradiated disks are heated to high gas and dust temperatures, T k up to 1000 K and T d up to ≃ 100 K (red curves in the upper panel of Fig. 4). These warm dust temperatures are on the order of the freezing temperature of H2O and significantly higher than that of CO ( ∼ 25K; e.g., Öberg et al. 2011), thus preventing their freeze-out in the outer disk and wind.  \nAstrong external FUV flux leads to high ice photo-desorption rates and high T d, and shifts the freeze-out depth (e.g., in disk height) where most oxygen becomes trapped in ice mantles (see Fig. 4 and Appendix H). Thus, we expect the vertical distribution of the x C / x O abundance ratio in irradiated disks to di ff er from that of isolated disks. In particular, grain growth, settling, and subsequent radial drift (Krijt et al. 2020; Zhang et al. 2020) will be less e ffi cient due to the lack of thick ice mantles in the strongly irradiated outer disk. This will result in lower depletions of gaseous CO compared to more evolved, ≳ 1Myr, disks (e.g., Bergner et al. 2020). Interestingly, the NIRSpec spectrum of d203-506 shows bright CO ro-vibrational emission at ∼ 4.7 µ m (see Berné et al. 2024), which we attribute to the externally irradiated layers. All in all, the outer gas layers of d203-506, and presumably of other externally irradiated disks (e.g., Boyden & Eisner 2023), must be 'oxygen rich,' with x C / x O very close to the solar ratio (assuming gaseous oxygen abundances as in the Orion cloud 4 ). In d203-506, this is further demonstrated by the presence of abundant gaseous OH, with x OH + x H2O ≳ 10 -5 (Zannese et al. 2024), and by the detection of bright and extended [O i ] 6302 Å emission (Haworth et al. 2023), where the excited O( 1 D) atoms in PDR gas (close to the H / H2 dissociation front) are produced by OH photodissociation (van Dishoeck & Dalgarno 1983; Störzer & Hollenbach 1998). In these irradiated layers, water undergoes photo-processing and cycles between gas-phase formation and photodissociation. Other inherited ice mantles will undergo photo-desorption and photo-processing as well. Does this mean that strongly irradiated disks develop a less rich carbon chemistry than isolated disks? On the contrary, we expect that hydrocarbons will be abundant in these O-rich disk layers, formed by gas-phase reactions triggered by the external FUV radiation, high temperatures, and enhanced abundances of C + , C ∗ , C 0 , and H ∗ 2 (see also Agúndez et al. 2008; Berné et al. 2023). These FUV-induced processes, facilitated by vertical transport, likely lead to di ff erent chemical compositions - similar to interstellar PDRs (e.g., rich in hydrocarbons, Cuadrado et al. 2015) - in the outer disk regions where gas-giant planets form. Additionally, the absence of carbon depletion in other proplyds would imply low disk masses for the non-detections in NGC 1977 by Haworth et al. (2022). Still, even higher angular resolution will be needed to map radial and vertical chemical abundances in Orion's disks.  \nAcknowledgements. We made used of ADS / JAO.ALMA#2021.1.01369.S data. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSTC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI / NRAO and NAOJ. The JWST data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program #1288. Support for program #1288 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127. We thank  \nour referee for a constructive report. JRG thanks the Spanish MCINN for funding support under grants PID2019-106110GB-I00 and PID2023-146667NB-I00. TJH acknowledges funding from a Royal Society Dorothy Hodgkin Fellowship and UKRI guaranteed funding for a Horizon Europe ERC consolidator grant (EP / Y024710 / 1). We thank the Programme National 'Physique et Chimie du Milieu Interstellaire' (PCMI) of CNRS / INSU with INC / INP, co-funded by CEA and CNES. EP acknowledges support from the University of Western Ontario, the Institute for Earth and Space Exploration, the Canadian Space Agency (CSA, 22JWGO1-16), and the Natural Sciences and Engineering Research Council of Canada. TO is supported by JSPS Bilateral Program, Grant Number 120219939.", 'References': "Abel, N. P., Ferland, G. J., & O'Dell, C. R. 2019, ApJ, 881, 130  \nAbrahamsson, E., Krems, R. V., & Dalgarno, A. 2007, ApJ, 654, 1171 Adams, F. C. 2010, ARA&A, 48, 47 Agúndez, M., Cernicharo, J., & Goicoechea, J. R. 2008, A&A, 483, 831 Alarcón, F., Bergin, E. A., & Teague, R. 2022, ApJ, 941, L24 ALMA Partnership, Brogan, C. L., Pérez, L. M., et al. 2015, ApJ, 808, L3 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 Badnell, N. R. 2006, ApJS, 167, 334 Bally, J., O'Dell, C. R., & McCaughrean, M. J. 2000, AJ, 119, 2919 Bergin, E. A., Alexander, C., Drozdovskaya, M., Gounelle, M., & Pfalzner, S.  \n- 2023, arXiv e-prints, arXiv:2301.05212\n- Bergin, E. A., Du, F., Cleeves, L. I., et al. 2016, ApJ, 831, 101\n- Bergner, J. B., Öberg, K. I., Bergin, E. A., et al. 2020, ApJ, 898, 97\n- Berné, O., Habart, É., Peeters, E., et al. 2022, PASP, 134, 054301\n- Berné, O., Martin-Drumel, M.-A., Schroetter, I., et al. 2023, Nature, 621, 56\n- Berné, O., Habart, E., Peeters, E., et al. 2024, Science, 383, 988\n- Birnstiel, T., Dullemond, C. P., Zhu, Z., et al. 2018, ApJ, 869, L45  \nBosman, A. D., Alarcón, F., Bergin, E. A., et al. 2021, ApJS, 257, 7  \nBowen, I. S. 1928, ApJ, 67, 1  \n- Boyden, R. D. & Eisner, J. A. 2020, ApJ, 894, 74\n- Boyden, R. D. & Eisner, J. A. 2023, ApJ, 947, 7\n- Cartledge, S. I. B., Meyer, D. M., Lauroesch, J. T., & Sofia, U. J. 2001, ApJ, 562, 394\n- Champion, J., Berné, O., Vicente, S., et al. 2017, A&A, 604, A69\n- Changala, P. B., Chen, N. L., Le, H. L., et al. 2023, A&A, 680, A19\n- Cuadrado, S., Goicoechea, J. R., Pilleri, P., et al. 2015, A&A, 575, A82\n- Díaz-Berríos, J. K., Guzmán, V. V., Walsh, C., et al. 2024, ApJ, 969, 165\n- Escalante, V., Sternberg, A., & Dalgarno, A. 1991, ApJ, 375, 630\n- Escalante, V. & Victor, G. A. 1990, ApJS, 73, 513\n- Fatuzzo, M. & Adams, F. C. 2008, ApJ, 675, 1361\n- Goicoechea, J. R., Aguado, A., Cuadrado, S., et al. 2021, A&A, 647, A10\n- Goicoechea, J. R. & Le Bourlot, J. 2007, A&A, 467, 1\n- Goicoechea, J. R., Pabst, C. H. M., Kabanovic, S., et al. 2020, A&A, 639, A1\n- Goicoechea, J. R., Santa-Maria, M. G., Bron, E., et al. 2019, A&A, 622, A91\n- Gordon, K. D., Clayton, G. C., Decleir, M., et al. 2023, ApJ, 950, 86\n- Grant, S. L., van Dishoeck, E. F., Tabone, B., et al. 2023, ApJ, 947, L6\n- Gressel, O., Ramsey, J. P., Brinch, C., et al. 2020, ApJ, 896, 126\n- Guan, L., Jiang, P., Zhang, G., et al. 2021, A&A, 647, A127\n- Habart, E., Le Gal, R., Alvarez, C., et al. 2023, A&A, 673, A149\n- Habart, E., Peeters, E., Berné, O., et al. 2024, A&A, 685, A73\n- Haris, K. & Kramida, A. 2017, ApJS, 233, 16\n- Haworth, T. J., Kim, J. S., Qiao, L., et al. 2022, MNRAS, 512, 2594\n- Haworth, T. J. & Owen, J. E. 2020, MNRAS, 492, 5030\n- Haworth, T. J., Reiter, M., O'Dell, C. R., et al. 2023, MNRAS, 525, 4129\n- Henney, W. J. 2021, MNRAS, 502, 4597\n- Hollenbach, D., Kaufman, M. J., Bergin, E. A., & Melnick, G. J. 2009, ApJ, 690, 1497\n- Jenkins, E. B. 2009, ApJ, 700, 1299\n- Johnson, C. T., Burke, P. G., & Kingston, A. E. 1987, Journal of Physics B Atomic Molecular Physics, 20, 2553\n- Johnstone, D., Hollenbach, D., & Bally, J. 1998, ApJ, 499, 758\n- Kama, M., Bruderer, S., Carney, M., et al. 2016, A&A, 588, A108\n- Kamp, I., Henning, T., Arabhavi, A. M., et al. 2023, Faraday Discuss., 245, 112\n- Krijt, S., Bosman, A. D., Zhang, K., et al. 2020, ApJ, 899, 134\n- Lada, C. J. & Lada, E. A. 2003, ARA&A, 41, 57\n- Le Petit, F., Nehmé, C., Le Bourlot, J., & Roue ff , E. 2006, ApJS, 164, 506 Lucy, L. B. 2002, MNRAS, 329, 787\n- Madhusudhan, N. 2012, ApJ, 758, 36\n- Mah, J., Bitsch, B., Pascucci, I., & Henning, T. 2023, A&A, 677, L7\n- Maucó, K., Manara, C. F., Ansdell, M., et al. 2023, A&A, 679, A82\n- McClure, M. K. 2019, A&A, 632, A32\n- McClure, M. K., Dominik, C., & Kama, M. 2020, A&A, 642, L15\n- Mendoza, C. 1983, in Planetary Nebulae, ed. L. H. Aller, Vol. 103, 143-172  \nMiotello, A., Facchini, S., van Dishoeck, E. F., et al. 2019, A&A, 631, A69  \nMiotello, A., van Dishoeck, E. F., Kama, M., & Bruderer, S. 2016, A&A, 594, A85  \n- Nieva, M. F. & Przybilla, N. 2012, A&A, 539, A143\n- Öberg, K. I., Facchini, S., & Anderson, D. E. 2023, ARA&A, 61, 287\n- Öberg, K. I., Guzmán, V. V., Walsh, C., et al. 2021, ApJS, 257, 1\n- Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011, ApJ, 743, L16\n- O'Dell, C. R., Abel, N. P., & Ferland, G. J. 2020, ApJ, 891, 46\n- O'Dell, C. R., Wen, Z., & Hu, X. 1993, ApJ, 410, 696\n- Pabst, C. H. M., Goicoechea, J. R., Teyssier, D., et al. 2020, A&A, 639, A2\n- Pascucci, I., Skinner, B. N., Deng, D., et al. 2023, ApJ, 953, 183\n- Peeters, E., Habart, E., Berné, O., et al. 2024, A&A, 685, A74\n- Pequignot, D. & Aldrovandi, S. M. V. 1976, A&A, 50, 141\n- Perotti, G., Christiaens, V., Henning, T., et al. 2023, Nature, 620, 516\n- Plomp, V., Wang, X.-D., Kłos, J., et al. 2023, arXiv e-prints, arXiv:2312.00203\n- Putaud, T., Michaut, X., Le Petit, F., Roue ff , E., & Lis, D. C. 2019, A&A, 632, A8  \nQiao, L., Coleman, G. A. L., & Haworth, T. J. 2023, MNRAS, 522, 1939  \nRamírez-Tannus, M. C., Bik, A., Cuijpers, L., et al. 2023, ApJ, 958, L30  \n- Ricci, L., Robberto, M., & Soderblom, D. R. 2008, AJ, 136, 2136\n- Rogers, C., de Marchi, G., & Brandl, B. 2024, A&A, 684, L8\n- Röllig, M., Abel, N. P., Bell, T., et al. 2007, A&A, 467, 187\n- Savage, B. D. & Sembach, K. R. 1996, ARA&A, 34, 279  \nSimón-Díaz, S. & Stasi'nska, G. 2011, A&A, 526, A48  \nSofia, U. J., Lauroesch, J. T., Meyer, D. M., & Cartledge, S. 2004, ApJ, 605, 272  \n- Staemmler, V. & Flower, D. R. 1991, JPhB., 24, 2343\n- Störzer, H. & Hollenbach, D. 1998, ApJ, 502, L71\n- Störzer, H. & Hollenbach, D. 1999, ApJ, 515, 669\n- Sturm, J. A., McClure, M. K., Harsono, D., et al. 2022, A&A, 660, A126\n- Temmink, M., Booth, A. S., van der Marel, N., & van Dishoeck, E. F. 2023, A&A, 675, A131\n- Tsukagoshi, T., Momose, M., Saito, M., et al. 2015, ApJ, 802, L7\n- van der Werf, P. P., Goss, W. M., & O'Dell, C. R. 2013, ApJ, 762, 101\n- van Dishoeck, E. F. & Black, J. H. 1988, ApJ, 334, 771\n- van Dishoeck, E. F. & Dalgarno, A. 1983, J. Chem. Phys., 79, 873\n- van Dishoeck, E. F., Grant, S., Tabone, B., et al. 2023, Faraday Discuss., 245, 52 van Dishoeck, E. F., Kristensen, L. E., Mottram, J. C., et al. 2021, A&A, 648, A24\n- Walmsley, C. M., Natta, A., Oliva, E., & Testi, L. 2000, A&A, 364, 301 Walsh, C., Millar, T. J., & Nomura, H. 2013, ApJ, 766, L23\n- Winter, A. J., Clarke, C. J., Rosotti, G. P., Hacar, A., & Alexander, R. 2019, MNRAS, 490, 5478\n- Winter, A. J. & Haworth, T. J. 2022, European Physical Journal Plus, 137, 1132 Winter, A. J., Haworth, T. J., Coleman, G. A. L., & Nayakshin, S. 2022, MNRAS, 515, 4287\n- Zannese, M., Tabone, B., Habart, E., et al. 2024, Nature Astronomy  \nZhang, K., Schwarz, K. R., & Bergin, E. A. 2020, ApJ, 891, L17  \n- 1 Instituto de Física Fundamental (CSIC). Calle Serrano 121-123, 28006, Madrid, Spain.e-mail: javier.r.goicoechea@csic.es\n- 2 Université Paris Cité, France.\n- 3 Chalmers University of Technology, Onsala Space Observatory, Onsala, Sweden.\n- 4 Dept. of Astronomy, University of Michigan, Ann Arbor, MI, USA.\n- 5 Institut de Recherche en Astrophysique et Planétologie, Université Toulouse III - Paul Sabatier, CNRS, CNES, Toulouse, France.\n- 6 LERMA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, Paris-Meudon, France.\n- 7 IRAM, 300 rue de la Piscine, 38406 Saint Martin d'Hères, France.\n- 8 Department of Physics and Astronomy, University of Western Ontario, London, Ontario, Canada.\n- 9 Institute for Earth and Space Exploration, University of Western Ontario, London, Ontario, Canada.\n- 10 Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA.\n- 11 Institut des Sciences Moléculaires d'Orsay, CNRS, Université ParisSaclay, Orsay, France.\n- 12 Université Paris-Saclay, CNRS, Institut d'Astrophysique Spatiale, Orsay, France.\n- 13 School of Physics and Astronomy, Queen Mary University of London, London E1 4NS, UK.\n- 14 Department of Astronomy, Graduate School of Science, University of Tokyo, Tokyo, Japan.\n- 15 Carl Sagan Center, SETI Institute, Mountain View, CA, USA.\n- 16 Leiden Observatory, Leiden University, The Netherlands.\n- 17 Astronomy Department, Univ. of Maryland, College Park, MD, USA.\n- 18 Space Telescope Science Institute, Baltimore, MD, USA.\n- 19 Instituto de Astrofísica e Ciências do Espaço, Lisbon, Portugal.", 'Appendix A: On the location of d203-506': "The NIR [C i ] lines observed by JWST in d203-506 are 3-5 times brighter than toward the Orion Bar PDR (Peeters et al. 2024). The submm [C i ] 609 µ m( 3 P13 P0) line integrated flux is 605 mJy km s -1 , slightly above the CO 3-2 line fluxes measured by Boyden & Eisner (2020, 2023) in other disks of the ONC. The velocity centroid of the main [C i ] 609 µ m line emission, v LSR (C i ) = 9.5 ± 0.1 km s -1 (see Fig. A.1), is in the range of LSR velocities of these disks. This velocity would be consistent with d203-506 being embedded at the rim of the Orion Bar where only FUV photons with E < 13 . 6eV penetrate inside. The amount of extinction internal to the Bar PDR can be estimated using the observed line intensity ratio of H2 lines that arise from the same upper v and J state. The inferred extinction toward d203-506 is nearly zero (Fig. 11 of Peeters et al. 2024), suggesting that the disk is actually not embedded in the Bar. Indeed, d203-506 shows a silhouette appearance at the wavelengths of several visible and NIR lines (H α , [S ii ] 6732 Å, Paschenα , etc.) arising from ionized gas in the background H ii region. This lead Haworth et al. (2023) to conclude that the disk lies in front of a structure of ionized gas known as the near ionization layer (NIL, Abel et al. 2019; O'Dell et al. 2020). Thus, located closer to the observer than the Bar, the Huygens H ii region, and the Trapezium stars. In this picture, d203-506 is in the foreground, well separated from the Bar, at least by ∼ 0.4 pc. The NIL would absorb any EUV radiation emitted by the massive stars θ 1 Ori C (in the Trapezium) and θ 2 Ori A, thus leaving the disk illuminated only by FUV radiation. Still, the LSR velocity of the ionized gas emission in the NIL, -12 ± 2kms -1 (O'Dell et al. 2020) is very di ff erent to that of d203-506, so this layer is likely not related with the birth place of the disk. Yet, the strong external FUV field and high gas-phase carbon abundance derived in d203-506 is very similar to that in the Orion cloud, suggesting a young disk, but separated from the main star-forming molecular core (OMC-1). Plausible scenarios are: i) the disk has recently emerged from a molecular globule, such as those detected by Goicoechea et al. (2020) in the Veil bubble of Orion, or, ii) d203-506 belongs to a recent star formation episode followed by radial migration (e.g., Störzer & Hollenbach 1999; Winter et al. 2019). Irrespective of its exact location, the important point is that d203-506 is externally illuminated by FUV radiation, but no EUV.  \nFig. A.1. [C i ] 609 µ m line profile toward d203-506 and the twocomponent Gaussian fit. The velocity centroid di ff erence of the two components is 2.3 ± 0.4 km s -1 and suggests that the blueshifted emission arises from a slow wind not seen in HCO + 4-3 emission (see Fig. 2).  \n<!-- image -->", 'H 2 * 2.12 𝝁 m (NIRCam)': 'Fig. B.1. [C i ] 0.9853 µ m observed with NIRSpec (orange contours at 3 × 10 -4 and 6 × 10 -4 erg s -1 cm -2 sr -1 ) over a NIRCam image of the H ∗ 2 1-0 S (1) emission. Cyan contours represent the H ∗ 2 1-0 S (1) emission observed by NIRSpec (7 × 10 -4 and 2 × 10 -3 erg s -1 cm -2 sr -1 ) at a similar angular resolution and spatial sampling as the [C i ] 0.9853 µ memission.  \n<!-- image -->', 'Appendix B: Spatial distribution of the [C i ] 0.9853 µ memission in d203-506': "Figure B.1 compares the [C i ] 0.9853 µ m(orange contours) and H ∗ 2 1-0 S (1) (cyan contours) emission observed with NIRSpec at similar angular resolution and spatial sampling. These observations have worse resolution and sampling than the H ∗ 2 1-0 S (1) image obtained by NIRCam (greenish colors). They show that both line emissions have very similar spatial distributions (at the ∼ 0.1 '' ≃ 40 au resolution of these images). Since in a PDR the intensity of FUV-pumped H ∗ 2 lines peaks slightly ahead of the H / H2 transition, the observed NIR C 0 lines also trace this transition from neutral atomic to molecular gas (as confirmed by our model predictions in Fig. 4), that is, the FUV-irradiated gas layers where the H2 abundance sharply increases. We note that spatially resolving the H ∗ 2 and NIR [C i ] line emissivity gradients will require an order of magnitude higher angular resolution.", 'Appendix C: Extraction of NIR spectra, background line emission, and FUV-pumping photons': "We first extracted the NIRSpec spectrum toward d203-506 (ON source measurement). As in our previous studies, we extracted the ON spectrum from an elliptical aperture centered on α (2000) = 5 h 35 m 20.357 s , δ (2000) = -5 · 25 ' 05.81 '' with dimension length l = 0.52 '' , height h = 0.38 '' , and a PA = 33 o East of North (Berné et al. 2024). This aperture includes both the inner, dense disk detectable in millimeter dust continuum emission (Berné et al. 2024) and also the extended, optically thinner, outer disk and photo-evaporative wind. Among other atomic ion lines, the ON spectrum contains a plethora of bright He i lines. These permitted NIR lines arise from very high-energy levels, up to ∼ 24 eV, thus they cannot be FUV-pumped lines of neutral helium. Instead, they are produced by EUV photoionization of He atoms, IP(He) = 24.6 eV, followed by radiative recombination cascades  \nFig. B.2. Selected JWST / NIRSpec spectra toward d203-506. The red spectrum corresponds to the ON source measurement, which includes emission from the background H ii region and the Bar PDR. The intrinsic spectrum of d203-506 is the ON f · OFF measurement (with f = 0.95 in blue and f = 1 in cyan). The scaling factor, f , is determined from the nebular He i recombination lines (see the main text). These spectra show NIR carbon lines (permitted and forbidden), [Fe ii ], and vibrationally excited H2 and OH lines (Zannese et al. 2024), among other species (see Berné et al. 2024). However, NIR O i fluorescent lines are not detected in the ON-OFF spectrum. These FUV-pumped O i fluorescent lines appear in the ON spectrum and arise from predominantly atomic gas close behind the ionization front of the Bar (the ∆ feature; Haworth et al. 2023; Peeters et al. 2024).  \n<!-- image -->  \nof He + . These He i recombination lines are nebular lines (Rogers et al. 2024) from the background H ii region (Peeters et al. 2024). This implies that toward the ON aperture, d203-506 is not opaque to the background NIR line emission. In order to obtain the intrinsic NIR spectrum of d203-506, one needs to correct for any background H ii region and Bar PDR line emission. As in our previous studies, we extracted a reference spectrum very close to the disk and providing the OFF-target correction. The OFF spectrum was obtained in a circular aperture of radius r = 0.365 '' centered α (2000) = 5 h 35 m 20.37 s , δ (2000) = -5 · 25 ' 04.97 '' , that is, at (0.20 '' , 0.84 '' ) from the disk (see Berné et al. 2024). The adopted intrinsic spectrum of d203-506 (Fig. 2) is the ON f · OFF spectrum, where f is the scaling factor needed to have no He i recombination line emission from d203-506 itself. We find f ≃ 0 . 95 -1, similar to Rogers et al. (2024) for other disks. This range provides a ∼ 5 % uncertainty in the line intensities of the NIR carbon lines, which we derived after baseline subtraction and Gaussian line fitting of the ON f · OFF spectrum.  \nThe ON source spectrum also shows several O i permitted lines: (3 d ) 3 D o - (3 p ) 3 P at 1.129 µ m, (4 s ) 3 S o - (3 p ) 3 P at 1.317 µ m, and (4 p ) 3 P-(3 s ) 3 S o at 2.894 µ m(see Fig. B.2), that are not seen in the spectrum of d203-506 (ON-OFF). The above lines arise from the background, from the predominantly neutral gas ( x H 0 > x H + ) that borders close behind the ionization front of the Bar PDR (Walmsley et al. 2000; Peeters et al. 2024). These NIR O i lines are not O + recombination lines -consistent with the  \nslightly higher IP of oxygen (13.618 eV) compared to hydrogen (13.598 eV)-, as there are no EUV photons capable of ionizing neutral oxygen atoms inside PDRs. The excitation of the NIR O i lines in the Bar PDR is attributed to FUV-pumping followed by fluorescence (Walmsley et al. 2000; Lucy 2002; Henney 2021; Peeters et al. 2024). As for neutral carbon, the ground-state of neutral oxygen is a triplet 3 P state, so that FUV excitations occur for high triplet states. However, the main pumping lines of the NIR O i lines at 1.129, 1.317, and 2.894 µ m lie at short FUV wavelengths: respectively at 1027 Å (12.1 eV; around Lymanβ ), 1040Å (11.9 eV), and multiple lines between 918 and 979 Å (13.5-12.7 eV). In a neutral PDR, these high-energy FUV photons (significantly more energetic than the 7.5-9.6 eV photons needed to pump C 0 ) are only available close to the ionization front, where the gas is mostly atomic.  \nThe lack of remarkable NIR O i emission in d203-506 (other than a faint O i 1.317 µ m emission in the ON- 0.95 · OFF spectrum; Fig. B.2) is not due to a reduced gas-phase elemental abundance of oxygen. Instead, it is caused by a reduced flux of FUV photons with energies capable of pumping the O i levels that are permeating the disk. This deficiency leads to very faint NIR O i fluorescent emission, likely confined to a very thin layer on the outermost surface of the wind. Hence, we conclude that d203-506 is not pervaded by a high flux of FUV photons with energies above ∼ 12eV: both FUV continuum photons and  \nLymanβ photons at 1026 Å (accidental resonance with the main FUV-pumping line of the NIR O i 1.317 µ m emission line in a Bowen-fluorescence mechanism). These photons are absorbed close to the wind surface layers. Indeed, at high n H / G 0 ratios ( ≳ 500cm -3 ) and relatively low e ff ective dust FUV absorption cross sections, these photons are removed by gas absorption: C 0 photoionization continuum (for λ< 1101 Å; E > 11.26 eV), Lyman H i , and Lyman-Werner H2 lines (see Fig. E.2), which develop broad line-absorption wings as gas column densities sharply increase (e.g., Goicoechea & Le Bourlot 2007). In Appendix E, we explicitly model these processes for the physical conditions in d203-506. In addition, a small fraction of ≳ 12 eV photons might be absorbed by any dust and neutral gas component situated between d203-506 and the massive stars θ 1 Ori C and θ 2 Ori A (e.g., van der Werf et al. 2013; O'Dell et al. 2020; Pabst et al. 2020).  \nThe above reasoning is also supported by the lack of fluorescent O i (3 p ) 3 P-(3 s ) 3 S o line emission at 0.845 µ m(VLT / MUSE images, see Haworth et al. 2023), where the upper energy level of the O i 0.845 µ mtransition is the lower-energy level of the NIR O i 1.129 µ mand 1.317 µ mlines targeted by NIRSpec (i.e., the O i 0.845 µ mline is pumped by the same ≳ 12eV FUV-photons than the O i lines targeted by NIRSpec, Bowen 1928; Walmsley et al. 2000; Henney 2021). Instead, the disk is seen in silhouette against the background O i 0.845 µ mline emission from the rim of the Bar PDR (the ∆ -feature, Haworth et al. 2023). Hence, most of the NIR O i fluorescent emission in the ON spectrum arises from this background feature (Fig. 8 of Peeters et al. 2024).  \nA reduced number of ≳ 12eV photons also implies that CO (X 1 Σ + ) photodissociation will mainly produce neutral carbon in the ground 3 P state (dissociation energy of 11.1 eV, e.g., van Dishoeck & Black 1988). Indeed, experiments show that FUV photons with energies ≥ 12.4 eV are needed to produce neutral carbon atoms directly in the first electronically excited 1 D state (Guan et al. 2021). Thus, this mechanism will be less dominant in d203-506.", 'Appendix D: Recombination theory applied to the observed NIR carbon lines': 'Figure D.1 is analogous to Fig. 14 of Walmsley et al. (2000). It shows the I (0.984 µ m) / I (1.069 µ m) (bottom) and I (1.069 µ m) / I (1.176 µ m) (top) line intensity ratios predicted by the recombination theory (C + + e -→ C ∗ + line cascade; Escalante & Victor 1990; Escalante et al. 1991) in two limiting cases: Case A (optically thin gas in all transitions) and the more realistic Case B (optically thick gas to photons produced by permitted transitions, and direct radiative recombinations to the ground state). Here, the I (0.984 µ m) refers to the intensities of the observed 0.9827 + 0.9853 µ mforbidden lines, and I (1.069 µ m) and I (1.176 µ m) are the sum of all permitted line intensities in these multiplets (see Table I.1).  \nIn the recombination theory, the I (1.069 µ m) / I (1.176 µ m) ratio depends mainly on the opacity of the FUV-resonant transitions. However, the observed I (1.069 µ m) / I (1.176 µ m) ratio in d203-506 (9.8 ± 6.0; blue-shaded areas) is significantly above the ratio predicted by this theory. Furthermore, the I (0.984 µ m) / I (1.069 µ m) line intensity ratio depends on the electron temperature. Without any extinction correction, the observed ratio only fits the predictions of Case A, which is less realistic. The upper observational value of the I (0.984 µ m) / I (1.069 µ m) ratio in Fig. D.1 includes a foreground extinction correction of 1.5 mag (Peeters et al. 2024) adopting the RV = 5.5 extinction law of Gordon et al. (2023). In any case, the absolute intensities of  \nFig. D.1. NIR carbon line intensity ratios predicted by the recombination theory of Escalante & Victor (1990). Curves represent the expected ratios under Case B (continuous) and Case A (dashed) conditions, taken from Walmsley et al. (2000). The blue-shaded area shows the line ratios in d203-506 that are compatible with the observations. See the main text.  \n<!-- image -->  \nthe NIR C 0 lines observed in d203-506 are stronger (see the next section and Table I.1) than those predicted by C + recombination theory (e.g., Escalante et al. 1991). This comparison implies that the involved C 0 levels should be predominantly populated by means other than recombination.', 'Appendix E: Inclusion of FUV-pumping for C 0 levels and C + +e -recombination followed by radiative cascades in the Meudon PDR code': 'Prompted by the impossibility of explaining the observed NIR carbon lines within the recombination theory, and drawing an analogy with the FUV-pumping mechanism of NIR O i lines at 1.129, 1.317, and 2.894 µ mseen toward the ionization fronts of interstellar PDRs (e.g., Bowen 1928; Walmsley et al. 2000; Lucy 2002; Henney 2021), we incorporated the FUV radiative excitation of C 0 in the Meudon PDR code (Le Petit et al. 2006). We specifically included 401 levels of C 0 (electronic states and their fine-structure levels), with electronic configurations up to 29 d (at ∼ 11.25 eV) and involving a total of 1572 radiative transitions. Figure E.1 shows a reduced Gotrian diagram of C 0 , showing up to the electronic levels associated with the NIR lines detected with NIRSpec. We took the main spectroscopic parameters from the NIST data base 5 . The external FUV field ( λ> 912Å) penetrating d203-506 is self-consistently computed as in Goicoechea & Le Bourlot (2007). The attenuation of FUV continuum and line photons impinging in d203-506 includes absorption and anisotropic scattering by dust grains as well as gas FUV-line absorption and photoionization continuum. The disk-position and λ -dependent FUV-field couples to C 0 by the resonant absorption of FUV photons from the 3 P ground-state to high-energy electronic triplet states through several pumping transitions: 1656 Å, 1277Å, 1279 Å, 1280 Å, etc. (see Fig. E.2). Subsequent visible and NIR fluorescent line emission cascades populate the lowerenergy C 0 levels. The brightest NIR carbon lines predicted by the model are precisely the ones detected by JWST / NIRSpec and VLT / MUSE: the forbidden lines at ∼ 0.984 and 0.8729 µ m, as well as the higher-energy permitted lines at ∼ 1.069, ∼ 1.176 , and ∼ 1.355 µ m, in decreasing order of intensity. The model also predicts that other NIR C i permitted line multiplets should be detectable (e.g., at ∼ 0.910, ∼ 0.966, and ∼ 1.455 µ m). However, the NIRSpec observation of d203-506 does not cover these specific wavelength ranges.  \nIn addition, we included a reduced description of C + + e -radiative recombination, including cascades toward the lowest 254 levels of C 0 (with core electrons in one of the ground-state 2 P3 / 2 or 2 P1 / 2 fine-structure levels of C + ). We adopted the recombination rates from Badnell (2006). However, for the physical conditions of the NIR-emitting gas layers of d203-506 ( G 0 of few 10 4 and n H ≃ 10 7 cm -3 ), we find that C + + e -recombinations only increase the intensity of the observed NIR carbon lines by ≲ 5 %. Therefore, in d203-506, this is not the dominant excitation process of these NIR lines.  \nOur excitation model also incorporates inelastic collisions of C 0 with e -, from the ground-state 3 P to the first and second excited C 0 electronic levels, 1 D2 and 1 S0 (rate coe ffi cients from Pequignot & Aldrovandi 1976; Mendoza 1983). Nevertheless, we find that collisional excitation to these high-energy levels plays a negligible role compared to FUV-pumping. Within the groundstate 3 P J fine-structure levels, the model also includes inelastic collisions with e -(Johnson et al. 1987), with H (Abrahamsson et al. 2007), with He (Staemmler & Flower 1991), and with o -/ p -H2 (Plomp et al. 2023, and extrapolated to high temperatures as in Goicoechea et al. in prep.). However, at the high densities in disks ( n H ≫ n cr , [CI]609 ≃ 10 3 cm -3 ), the [C i ]609 µ mfine-structure line is thermalized ( T ex = T k), and its excitation is independent on the collisional rate coe ffi cients. [C i ] 370 and 609 µ m line intensities are not a ff ected by FUV-pumping either.  \nFig. E.1. Reduced Gotrian diagram of C 0 showing up to the electronic levels associated with the NIR lines discussed in this study. The figure also includes the main FUV-pumping lines (in blue) of the detected NIR C 0 lines, as well as the submm fine-structure lines within the ground state detectable with ALMA (in green; energy level splitting is exaggerated).  \n<!-- image -->  \nFigure 1:  \nC  \n3d  \nradiative transitions.  \n2  \nOur PDR models including FUV-pumping and fluorescent deexcitation match the intensity level of the observed NIR carbon lines. Without this radiative excitation mechanism, the predicted NIR line intensities would be up to two orders of magnitude lower than observed (Table I.1). FUV-pumping also produces higher values of the I (1.069 µ m) / I (1.176 µ m) ratio than those predicted by recombination theory. The reference model predicts that the observed NIR and [C i ] 609 µ mlines are optically thin, with τ ≃ 10 -6 and τ ≃ 10 -1 , respectively. In contrast, the main C i FUV-pumping lines at 1656 Å, 1277 Å, 1279 Å, and 1280 Å have very high opacities, ranging from τ ≃ 200 to ≃ 5000. The FUV continuum is optically thick, τ ≃ 30-40 (see Fig. E.2).  \nIn the Orion Bar PDR, Peeters et al. (2024) found high I (0.984 µ m) / I (1.069 µ m) line ratios and I (1.069 µ m) / I (1.176 µ m) line ratios that are inconsistent with recombination theory too (their Fig. G.1). This suggests that recombination alone may not drive the excitation of NIR carbon lines in strongly irradiated interstellar PDRs either. Adding the role of FUV-pumping will likely relax the hot temperatures (several thousand K) and densities ( ∼ 10 8 cm -3 ) determined by Peeters et al. (2024). That is, there will be no need to invoke the presence of very high-pressure small clumps, which are otherwise not seen in the JWST images.', 'Appendix F: FUV attenuation by gas lines and dust': 'In this section we provide evidence that in irradiated protoplanetary disks, FUV photons at the pumping transitions of the O i fluorescent lines (photons with E ≳ 12eV) can be more quickly removed than photons pumping neutral carbon ( E ≲ 9.6 eV). Compared to interstellar PDRs, where small dust grains dominate the absorption of FUV photons, gas absorption becomes very important in the upper layers of irradiated  \nFig. E.2. External FUV radiation field (in energy density) propagating inside the photoevaporating wind and outer disk (reference model of d203-506 with G 0 = 2 × 10 4 and n H = 10 7 cm -3 ). The upper panels (blue spectra) are for a slab position at AV = 0.1 mag (close to the PDR surface), whereas the lower panels (red spectra) are for a position deeper inside the PDR, at AV = 0.5 mag. The disk-depth-dependent FUV field results from a FUV radiative transfer calculation that includes gas absorption, H2 fluorescent emission, and grain absorption and scattering. The resulting FUV continuum and line emission are modulated by strong H i absorption lines (Lyα and Lyβ lines are indicated in this figure) and by a forest of H2 absorption and fluorescent emission lines. The black vertical lines show the wavelength position of the main C 0 pumping transitions at 1277, 1279, and 1280 Å ( ∼ 9.6 eV). The vertical dashed gray lines show the wavelength position of the main O 0 pumping transitions at 1027 and 1040 Å ( ∼ 12 eV). These photons are quickly absorbed by small column densities of gas, mainly via C 0 photoionization below 1101 Å for AV < 1 mag (Fig. F.1). Lower-energy FUV photons propagate deeper inside the disk and are attenuated by dust grains.  \n<!-- image -->  \nprotoplanetary disks. This is due to the reduced e ff ective FUV dust absorption cross section, caused by larger grain sizes and possibly higher gas-to-dust mass ratios than in the ISM. We determine that, close to the PDR surface of d203-506 (from AV ≃ 0.05 to ≃ 0.8 mag), C 0 photoionization dominates the absorption of FUV photons with λ< 1101Å ( > 11.26 eV). For longer photon wavelengths and for PDR depths higher than AV ≃ 0.8 mag, dust absorption dominates. To be more quantitative, Fig. F.1 shows the evolution of the ratio of the FUV gas absorption coe ffi cient over the dust extinction coe ffi cient as a function of AV for two representative continuum wavelengths: 1100.73 Å (below the C 0 ionization threshold; in black) and 1102.50 Å (above this threshold, i.e., not producing C 0 ionization; in red). The black curve shows that this ratio reaches nearly a factor of 3 (that is, C 0 photoionization dominates over FUV dust absorption), thus producing enhanced absorption of < 1101Å photons compared to photons at longer wavelengths. This is exemplified by the sharp step, delineated by the green vertical line, in the FUV spectrum of Fig. E.2. This plot shows that FUV photons at the pumping transitions of neutral oxygen atoms ( ∼ 1000Å) are scarcer than those at the pumping transitions of neutral carbon ( ∼ 1300 Å).', 'Appendix G: Variations in I ([C i ] 609 µ m) and NIR carbon lines with FUV flux and gas density': 'In this study, we assume that at the physical scales probed by our observations, the density distribution is smooth across the outer disk and inner wind. Therefore, as a first approximation, we adopt constant density models. In PDR environments, the C 0 column density, and thus the intensity of the [C i ] 609 µ m line, are not strong functions of the external FUV flux or gas density. In general, N (C 0 ) somewhat decreases with increasing n H and slightly increases with increasing G 0. For the irradiation conditions in the outer layers of d203-506, G 0 of a few 10 4 , and gas densities varying between n H = 10 6 cm -3 and 10 8 cm -3 , the predicted line intensity I ([C i ] 609 µ m) only changes by a factor of 3. Adopting n H = 10 7 cm -3 in the outer disk and inner wind, I ([C i ] 609 µ m) only increases by a factor of ≲ 2 when the external FUV flux changes from G 0 = 10 3 and 10 5 (see Fig. G.1). Therefore, with a reasonable knowledge of G 0 and n H, the I ([C i ] 609 µ m) line intensity can be used to accurately determine the column density of C 0 and, through modeling, the (total) gas-phase carbon abundance, x C = x C + + x C 0 + x CO + ... , in the [C i ] 609 µ m-emitting layers.  \nFig. F.1. Evolution of the ratio between the wavelength-dependent FUV gas absorption coe ffi cient and the dust extinction coe ffi cient (absorption plus scattering) as a function of AV for two representative continuum wavelengths: one below the C 0 photoionization threshold (black curve) and one above (red curve). For positions where the black curve is greater than 1, C 0 photoionization dominates the absorption of FUV photons with λ< 1101 Å. This process reduces the amount of FUV photons that pump the NIR oxygen lines close to the PDR surface.  \n<!-- image -->  \nIn our models, however, the NIR carbon lines follow a different emission trend as they are mostly sensitive to the flux of FUV-pumping photons (i.e., to G 0; see Fig. G.2). This is a consequence of the negligible collisional (de)excitation of the NIR C 0 lines ( n H ≪ n cr) and the attenuation of the FUV-pumping photons at high gas densities. At high G 0, this makes the column density of C 0 that is sensitive to FUV-pumping, N (C 0 pump ), roughly proportional to G 0 / n H, and the intensity of the NIR C 0 fluorescent lines approximately proportional to n H N (C 0 pump ) ≃ G 0. This conclusion applies to both the high-energy C i permitted lines and also the [C i ] forbidden lines at 0.8729 µ m and ∼ 0.984 µ m (we recall the later are not observed in the innermost regions, < 1 au, of isolated protoplanetary disks; McClure 2019, but are bright in externally irradiated disks). Therefore, the NIR C 0 line intensities depend less on the gas density; they only change by ∼ 30% if the gas density changes by an order of magnitude compared to the reference model. For the same reason, a reduction of an order of magnitude in the total gas-phase carbon abundance, x C, only changes the intensity of the NIR carbon lines by a factor of ≲ 2. Thus, we conclude that the extended emission from these NIR lines is a powerful probe of external FUV radiation fields impinging on protoplanetary disks in cluster environments.  \nFig. G.1. Grid of constant-density PDR models for varying external FUV radiation fields but fixed n H = 10 7 cm -3 and x C = 1.4 × 10 -4 values. The red markers show the predicted column density of C 0 . The blue markers show the predicted [C i ] 609 µ mline intensity, integrating from AV = 0 to 10 mag into the wind and disk system. The horizontal lines mark the observed line intensity ( ± σ ) in d203-506.  \n<!-- image -->  \nFig. G.2. Grid of constant-density PDR models for varying external FUV radiation fields but fixed n H = 10 7 cm -3 and x C = 1.4 × 10 -4 . Colored curves show the predicted line intensities of the neutral carbon lines discussed in the text relative to the line intensity reference model ( G 0 = 2 × 10 4 ). While the NIR carbon line intensities increase with G 0, the [C i ] 609 µ mline intensity is much less dependent on G 0.  \n<!-- image -->', 'Appendix H: Little H 2 O freeze-out in the outer disk and wind layers traced by [C i ] 609 µ m': 'Water ice is the first ice mantle to condense when dust temperatures cool below T d ≲ 130K (depending on the assumed gas density, e.g., van Dishoeck et al. 2021). However, water ice photo-desorption and thermal desorption replenishes the gasphase oxygen abundance and reduces the gas-phase x C / x O ratio. Figure H.1 shows a simple analytical model of the predicted shift of the water freeze out depth ( AV , f ) at which most undepleted (i.e., not in refractory material) gas-phase oxygen is incorporated into water ice mantles. This plot shows AV , f , in mag of visual extinction from the edge of the wind / disk system. We plot AV , f as a function of the external FUV photon flux, in G 0 units, for a typical gas density in the outer disk and inner wind, n H = 10 7 cm -3 , and T k = T d. This analytical determination of the H2O freeze out depth implies equating the rate at which undepleted gas-phase O atoms stick on grains, R gr , O (first step in the grain surface formation of water ice) and the photo-desorption rate of water molecules from icy grain surfaces, R photodes , H2O. These rates go as R gr , O ∝ n H T 1 / 2 k x O and R photodes , H2O ∝ Y G 0 exp( -bAV ), respectively, where x O is the undepleted gas-phase abundance of O atoms, Y is the number of desorbed molecules per incident FUV photon, and b is a dust-related FUV field absorption factor (for the basic formalism and standard parameters see e.g., Hollenbach et al. 2009). The water freeze out takes place very close to the irradiated disk surface for high values of the gas density over the FUV flux ratio, n H / G 0 ≳ 10 5 cm -3 . However, for strong external FUV radiation fields ( G 0 > 10 2 in our representative example), the n H / G 0 ratio decreases, and water ice only becomes an abundant oxygen-reservoir deeper inside the disk. In this toy model of the outer disk layers of d203-506, the water freeze out depth shifts by ≃ 4.5 mag (tens of au) compared to isolated disks, that is, beyond the [C i ] 609 µ m emitting layers (Fig. 4). The blue curve in Fig. H.1 also shows the analytical estimation of the CO freeze-out depth using the rate at which gas-phase CO molecules stick on grains. Nearly all volatile oxygen and carbon in the disk zone traced by [C i ] 609 µ m and NIR C 0 line observations is in the gas-phase, and the gaseous elemental x C / x O abundance ratio is low, close to the solar value in the case of Orion.  \nOur more detailed thermo-photochemical PDR models in Fig. 4 include gas-phase and simple gas-grain exchanges for O, OH, H2O, O2, and CO. These species adsorb on dust grains as temperatures drop, are photo-desorbed by FUV photons, desorb via cosmic-ray impacts, and thermally sublimate (see Putaud et al. 2019, for the inclusion of these processes in the Meudon PDR code). Only for water ice formation we included the grain surface reactions s-H + s-O → s-OH and s-H + s-OH → s-H2O, where srefers to the species in the solid (see Goicoechea et al. 2021). Thus, while the water ice abundance and H2O freeze-out depths are self-consistently calculated, more detailed surface reactions and models will be needed to accurately describe the growth and (photo)processing of carbon-bearing ices under strong external FUV irradiation.', 'Appendix I: Spectroscopic parameters of neutral atomic carbon lines': 'Figure H.2 shows a reduced energy-level diagram of neutral carbon including fine-structure splittings. Table I.1 summarizes the spectroscopic parameters and line intensities of the carbon lines discussed in the text (including PDR model predictions).  \nFig. H.1. Shift in the H2O and CO freeze-out depth (simulating the outer disk) with increasing external FUV radiation field. This analytical model assumes n H = 10 7 cm -3 , as in d203-506. The vertical lines show the G 0 values compatible with the observed NIR carbon line intensities.  \n<!-- image -->  \nFig. H.2. Lowest-energy states (in K units) of neutral carbon, including the fine-structure levels and brightest NIR and submillimeter lines discussed in the text. Fine-structure splittings are exaggerated for clarity.  \n<!-- image -->  \nTable I.1. Spectroscopic parameters (from NIST 5 and Haris & Kramida 2017) of the neutral carbon lines discussed in the text, intensities of the lines detected by JWST and ALMA in d203-506, and model predictions.  \nNotes. a Wavelengths in a vacuum. † Observed with MUSE (Haworth et al. 2023) and quoted also as the [C i ] 8729 Å line. b In units of erg cm -2 s -1 sr -1 . c Adopting AV = 1.5 mag of foreground extinction (see text). d Line intensity predicted by the reference PDR model (including C 0 FUV-pumping and C + recombination cascades) with G 0 ≃ (1-2) × 10 4 , n H = 10 7 cm -3 , and x C = 1.4 × 10 -4 . e Same but neglecting FUV-pumping of C 0 lines.'}