{'fraction_non_alphanumeric': 0.036180432285684454, 'fraction_numerical': 0.028819045630155582, 'mean_word_length': 4.11810287241149, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The structure of the Orion Nebula in the direction of Theta 1 Ori C

We have used existing optical emission and absorption lines, [C II] emission lines, and H I absorption lines to create a new model for a Central Column of material near the Trapezium region of the Orion Nebula. This was necessary because recent high spectral resolution spectra of optical emission lines and imaging spectra in the [C II] 158 micron line have shown that there are new velocity systems associated with the foreground Veil and the material lying between Theta 1 Ori C and the Main Ionization Front of the nebula. When a family of models generated with the spectral synthesis code Cloudy were compared with the surface brightness of the emission lines and strengths of the Veil absorption lines seen in the Trapezium stars, distances from Theta 1 Ori C, were derived, with the closest, highest ionization layer being 1.3 pc. The line of sight distance of this layer is comparable with the size of the inner Huygens Region in the plane of the sky. These layers are all blueshifted with respect to the Orion Nebula Cluster of stars, probably because of the pressure of a hot central bubble created by Theta 1 Ori C's stellar wind. We find velocity components that are ascribed to both sides of this bubble. Our analysis shows that the foreground [C II] 158 micron emission is part of a previously identified layer that forms a portion of a recently discovered expanding shell of material covering most of the larger Extended Orion Nebula.


INTRODUCTION
As the closest star formation region that has a bright H II region, the Orion Nebula presents a unique testingground for theories of star formation and the physical processes governing astronomical plasmas. This opportunity has created a rich literature on these subjects (Ferland 2001;O'Dell 2001;Muench et al. 2008;O'Dell et al. 2008;Goicoechea et al. 2015;Kong et al. 2018). This article reports on a detailed study of material beyond and in the foreground of the dominant ionizing star (θ 1 Ori C) in the direction of the Trapezium stars. The brightest components (we shall call them layers as we study only a small area in the plane of the sky) lie beyond θ 1 Ori C and are largely associated with emission from the blister of ionized gas on the surface of the host molecular cloud. Less visible layers lie in the fore-ground of θ 1 Ori C and are collectively called the Veil. The Huygens Region falls near the northeast corner of a much lower surface brightness elliptical bubble known as the Extended Orion Nebula (EON, Güdel et al. (2008)).

Background of this study
Orion's Veil is a series of atomic and ionized layers, physically associated with the Orion star-forming environment. The Veil covers the brightest part of M42 (the Huygens Region) and is the primary source of extinction for the Nebula and for the stars ionizing M42, the Trapezium cluster. O'Dell & Yusef-Zadeh (2000) have mapped the extinction across the entire face of the Nebula, with about 1.6 mag of extinction towards the Trapezium. Additionally, van der Werf & Goss (1990), Troland, Heiles, & Goss (1989) and Troland et al. (2016) have mapped the magnetic field strength in the Veil through H I 21 cm absorption.
The simultaneous mapping of extinction and magnetic field for an environment with a well-characterized spectral energy distribution (SED) has made the Veil an ideal testbed to understand the physics of the interstellar medium (ISM), in particular the partitioning of magnetic, turbulent, and thermal energetics. In a series of works (Abel et al. (2004); Abel et al. (2006); Abel et al. (2016)) combined UV and optical absorption line data towards the Trapezium with the existing magnetic field and extinction maps to produce detailed calculations of the physical conditions in the two primary atomic layers seen in H I 21 cm absorption towards the Trapezium (typically referred to as Components A and B). Utilizing the spectral synthesis code Cloudy Ferland, et al. (2017), these works computed the density, temperature (gas and spin), and chemical composition for each component, then used this data to compare magnetic, thermal, and turbulent energies. Overall, the combination of observations and theoretical modeling showed Component B (2.0 pc) to be closer to the Trapezium than Component A (4.2 pc), with Component B also being somewhat warmer and denser than Component A. Finally, both layers appeared to be magnetically dominated, meaning the energy associated with the magnetic field of each component exceeds both the thermal and turbulent energies.
Although the series of Abel et al. papers focused primarily in the atomic Veil layers, an H + region must exist between the Trapezium and Component B. In the optical, 2S 3 He 0 absorption at 388.9 nm is observed towards the Trapezium and θ 2 Ori A. In the UV, S 2+ and P 2+ absorption is observed towards the Trapezium Abel et al. (2016). While the series of Abel et al. papers attempted to account for the ionized gas in their modeling efforts, the lack of emission features associated with this ionized component (due to the much brighter M42 on the other side of the Trapezium) have made it difficult to state anything definitive.
Two recent studies allow the possibility of understanding the ionized component between the Trapezium and Component B. O'Dell (2018) performed a detailed study of the velocity components from the inner portions of the Orion Nebula and found optical [N II] (658.4 nm) and [O III] (500.7 nm) emission features at similar velocities to the 2S 3 He 0 , S 2+ , and P 2+ absorption features seen in the optical and UV. These emission features are critical, as they allow the calculation of the physical properties of the ionized component, which can then be used to determine the effects of the ionized component model on the model results for the atomic layers. Another recent study, by O'Dell et al. (2017), calculated the relative effects of θ 1 Ori C and another nearby O star, θ 2 Ori A, on the ionization balance throughout the Nebula. This study found, in many regions of Orion, that θ 2 Ori A cannot be ignored as a source of ionizing photons. The effects of θ 2 Ori A therefore, should be explored in any calculation of the ionized component. Additionally, in the light of O'Dell et al. (2017), the degree to which θ 2 Ori A affects the neutral layers along the line of sight also requires investigation.

Outline of this paper
The format of this work is as follows. In Section 2, we highlight the key observations used in this study, emphasizing the recent work of O'Dell (2018) along with previous optical, UV, and radio studies of the region and the line of sight towards the Trapezium. Section 3 presents revised models of the Veil, taking into account the O'Dell (2018)  We end with a series of conclusions in Section 6.
Throughout this paper we designate radial velocities in the the Heliocentric system. These may be converted to LSR velocities by subtracting 18.1 km s −1 . We adopt the distance of 383±3 pc derived by Kounkel et al. (2017), which is in agreement with more recent results using Gaia DR2 (Großschedl et al. 2019). This distance gives a scale of 1.86×10 −3 pc/ .

OBSERVATIONS
We draw on previously published optical window observations, particularly those in the high spectral resolution Spectroscopic Atlas of Orion Spectra (García-Díaz et al. 2008) (henceforth 'the Atlas') prepared with spectra from Kitt Peak National Observatory, Cerro-Tololo Interamerican Observatory, and San Pedro Mártir Observatory. The Atlas was composed from north-south slit spectra at intervals of 2 in Right Ascension and have a velocity resolution of 10 km s −1 . The spatial resolution in Declination was seeing limited at about 2 . The [O III] observations of Doi et al. (2004) were used in creating the Atlas and a good illustration of a spectrum is given in Fig. 1 of that study. We also used the lower spectral resolution study made with the MUSE 2-D spectrograph at the European Southern Observatory's Paranal observatory (Weilbacher et al. 2015). In this case the velocity resolution was about 107 km s −1 and the spatial resolution was about 1 .
In the ultraviolet we employed the HST ultraviolet spectra of θ 1 Ori B at 2.5 km s −1 resolution described in detail in Abel et al. (2016). All of the lines of interest were seen in absorption against the stellar continuum of θ 1 Ori B.
We also used the results from Goicoechea et al. (2015). That study included all of the Huygens Region but did not extensively go into the EON. It reported Herschel space craft spectra of the [C II] 158 µm line at 0.2 km s −1 and 11. 4 resolution. The study also presented H41α observations with the IRAM-30m telescope. These had 0.65 km s −1 and 27 resolution. Their discussion also used CO 2-1 observations by Berné et al. (2014) at 0.4 km s −1 and 11 resolution.
We also employed calibrated emission line images made with the Hubble Space Telescope (O'Dell & Wong 1996;O'Dell 2009) and lower resolution ground-based camera images by Robert Gendler reproduced in O'Dell & Harris (2010). The latter cover the entire EON region and go beyond the mosaic of HST images first reproduced in Henney et al. (2007).

Surface Brightnesses and column densities of lines used in our study
Our analysis has required derivation of the extinction corrected surface brightness of the V low velocity component discussed in O'Dell (2018). A refined definition of the multiple velocity components is given in Table 1 and discussed in Section 2.2.
In  Wen (1992), and it also appears in the high resolution spectrum of Baldwin et al. (2000). Unfortunately none of these detections fall within the region used to obtain the surface brightnesses in [N II] and [O III]. However, they lie close to this region and the ratio of the signals of the V low and V mif components are quite similar, at 0.10±0.01. We then used the extinction-corrected images in the MUSE study (Weilbacher et al. 2015) to determine the surface brightness of the Main Ionization Front (MIF) component and then scaled this down by a factor of ten to find S(V low ,[O I]) = 1.6 ±0.53 ×10 −4 ergs s −1 cm −2 sr −1 .
The final line used in the analysis was the 388.9 nm line of He I. Most of the helium lies in the metastable 2 3 S state that is populated as ionized helium recombines. In the study of absorption lines in Trapezium stars ) a column density of 1.45×10 13 ions cm −2 was found in θ 1 Ori C.

Characteristic Velocity Systems
In the study of the brightest part of the the Huygens Region O'Dell (2018) averaged the Atlas spectra over areas of 10 ×10 , thus increasing the signal to noise ratio. This grid is shown in Figure 1. In order to identify systematic patterns in the [N II] and [O III] velocity components we defined an area within this grid of Atlas spectra, as shown in Figure 1, designated as the NE Region. It was selected to include the dominant ionizing star θ 1 Ori C and to be free of major shocks and jets. We have examined each of the 10 ×10 averaged spectra within the NE-Region, assigning velocity components according to the criteria given in Table 1. A good illustration of one of these [N II] spectra and its deconvolution is shown in Fig. A1 of O'Dell (2018). Sometimes there were two components to the red or blue of the V mif or V wide components, making the assignments easy. In the case of there being only one, then the signal of that component as compared with the corresponding V mif or V wide was used for assignment. These criteria produced an almost unambiguous assignment of each component. The results are shown in Figure 2, with the assigned errors being based on the scatter of the values. The accuracy of very weak component lying on the shoulder of the strong V mif or V wide component is discussed in the Appendix of O'Dell (2018). In this section the velocity groups are arranged in order of descending radial velocities, whereas in Section 5 and Table 2 they are arranged according to their probable distance from the observer. All of these optical window lines arise from the MIF which is composed of material photo-evaporating towards the observer from the Photon Dominated Region (PDR). The PDR lies immediately beyond the MIF and is a shock compressed region lying along the surface of the Orion Molecular Cloud (OMC).
The V scat components in the Huygens Region are usually attributed to backscattering by dust in the PDR. and absorption line velocities gives further confimation of the ionized Component I as real.
The V blue components are the most enigmatic, with no current explanation. They are weaker than the V low components, even after factoring in the classification criteria. The velocities are more negative than the V low components and have a greater range. They are most common (23 of 32 samples) in the NE-Region [N II] spectra.

Computational Details
We use c17 of the spectral synthesis code Cloudy ) to find a self-consistent model for each component observed along the line of sight to the Trapezium. The assumed gas and dust abundances,infrared, optical, and UV stellar continuum (SED), and grain properties are all identical to Abel et al. (2016), which in turn follows the general logic of the previous modeling efforts of the Veil presented in (Abel et al. 2004(Abel et al. , 2006 and Lykins et al. (2010). The primary difference between the calculations presented here and the one presented in previous works is the modeling of two ionized components along the line of sight (which must lie between the components observed in 21 cm absorption and the Trapezium) and the way we propagate the radiation field to model subsequent components along the line of sight. These components are seen in Figure 3 and are described below: I. The ionized gas component closest to the Trapezium, with observed values for the surface brightness of [N II] 658.3 nm and [O III] 500.7 nm. In addition, the 2S 3 He 0 (henceforth He I*) S 2+ and P 2+ column densities are known from previous absorption studies in the optical and UV Abel et al. 2006).
II III. The H I 21 cm component designated as Component B in previous Veil studies, whose velocity is associated with excited state H 2 observed in UV absorption. Column densities for a variety of elements are derived in Abel et al. (2006), and is the source of about two-thirds of the total H I column density along the line of sight. For clarity this component will frequently be designated as III (B).
IV. The H I 21 cm component, previously designated as Component A, with similar column densities observed in this component to that of Component III, with the exception being the lack of molecular hydrogen. The remaining one-third of the total H I column density lies in this component. For clarity this component will frequently be designated as IV (A).
The velocities of the individual components make it likely that none of the four components listed are physically adjacent to each other. However, since all components are affected by the same SED, each component is coupled through radiative transfer. The radiation which makes it through the first component will be the incident continuum impingent upon the 2 nd component, which continues until the radiation field which makes it through H I 21 cm component A, which ends up heading towards the observer. Our calculations account for this effect. Once we find the best-fitting model which accounts for the observed [N II], [O III], and 2S 3 He 0 column density, we save the transmitted SED for the last zone of the calculation, and use this transmitted SED as the incident SED for the next component. For each component, we vary the distance to the Trapezium and the total hydrogen density of the layer (assuming constant density, where "total" is defined as the sum of the hydrogen atoms per cm −3 in ionized, atomic, and molecular form). For component I, we stop all calculations at the observed N (S 2+ ) column density of 10 15 cm -2 . For component II, each calculation went until the H + /H 0 ionization fraction reached 10 -2 . Finally, for components III (B) and IV (A), we stopped the calculation at the deduced H 0 column density of 10 21.5 cm -2 (component B) and 10 21.2 cm -2 (component A).  Figure 6 is a similar figure for 2S 3 He 0 column density, where we assume an error bar of ± 0.3 dex since O' Dell et al. (1993) only report the average column density of 10 13.3 cm -2 . The best fit to the observations (shown by the dot in Figure  6) yield S([N II]) = 3.97x10 -3 ergs s −1 cm −2 sr −1 ,S([O III]) = 1.6x10 -2 ergs s −1 cm −2 sr −1 , and a HeI* column density of 10 12.98 cm -2 , corresponding to a density of around 10 2.65 cm -3 , and a distance of 10 18.6 cm (1.3 pc).

Component II
Given the density/distance combination deduced for Component I, we saved the transmitted continuum for the last zone of this calculation and used this as the input SED to model Component II, which represents the [O I] emitting region. Calculations for this component spanned a distance of log R (cm) = 18.7 -19.3, in increments of 0.1 dex, and a density range of log n (cm −3 ) = 2 -5 . The lower range of distance was chosen so that it is farther from the Trapezium than component I. Figure 7 shows the overall results, that is, both the calculated values and the range that agree with the observations. A broad range of distances exists where the [O I] observations are reproduced. However, we can exclude some  (2006) showed that if the neutral Veil layers are much farther than 10 19 cm from the Trapezium, then significantly more molecular hydrogen will form, contrary to the low H 2 abundance observed in the Veil. Therefore, we have set as an upper limit for the [O I] distance as 10 18.9 cm, so that in the modeling of component III (B), the closest initial distance of to the Trapezium would not exceed 10 19 cm. This leads to three distance/density combinations for the [O I] emitting layer (at a resolution of 0.25 dex in density, 0.1 dex in distance): I. log n (cm −3 ) = 2; log R (cm) = 18.7 II. log n (cm −3 ) = 2.25; log R (cm) = 18.7 -18.9 III. log n (cm −3 ) = 2.5; log R (cm) = 18.8 -18.9 For each of these combinations, we saved the transmitted continuum for use in modeling Component III (B) in Section 3.4.

Component III (B)
Using the transmitted SED through the [O I] emitting layer, we computed the predicted SED for component III (B). As with the first two components, we varied the distance and density away from the Trapezium, with the closest distance being 0.1 dex farther than the distances listed for the six combinations in Section 3.3. The upper limit to the distance for all calculations is 10 19.3 cm, which assures all calculations go beyond the upper limit to the total H 2 column density of 10 17.35 cm -2 . The density was varied from log n (cm −3 ) = 2 -4, in increments of 0.2 dex, which covers about 1 dex on either side of the density derived in Abel et al. (2006)  layer cannot be more than 10 18.9 cm away, or else Component III (B) will produce more H 2 than observed.
The combination of distance-density that best reproduces Component III (B) is around a density of log n (cm −3 ) = 3 -3.2, and log R (cm) = 18.8, which is nearly identical to previous modeling efforts of Component III (B) in Abel et al. (2006Abel et al. ( , 2016. Figure 8 shows a comparison of the predicted H 2 spectrum (red) for a density of log n (cm −3 ) = 3, log R (cm) = 18.8 to the derived H 2 column densities from STIS observations (blue). Other calculations, such as CI, CI * , and CI * * , are unchanged from the published contour plots from Abel et al. (2016).

Component IV (A)
Given the results of calculations for Veil Component III (B), and the similarities between the Abel et al.
(2016) paper and the current work, calculations for Component IV (A) are unnecessary, as given the same density, distance, and SED, Component IV (A) will produce the same deduced physical conditions as in the 2016 paper. In Abel et al. (2016), we deduced a density of log n (cm −3 ) = 2.4, and a distance of 4.2 parsecs. As the density/distance controls other deduced physical conditions, such as the balance between magnetic, thermal, and turbulent energies, these deduced properties are also unchanged from Abel et al. (2016).

Observational Constraints
In our models of the several components of the Veil we have only considered contributions of the Trapezium stars, dominated in the EUV by θ 1 Ori C, with FUV contributions by the later spectral type cooler members of the Trapezium. The models presented do not assume any contribution by the O9.5V star θ 2 Ori A, but it is worth examining the effects of that star, which has much less EUV luminosity than θ 1 Ori C, but relatively more FUV luminosity. In order to be important in photoionizing Component I, θ 2 Ori A would have to lie on the side of the Veil away from the observer, although this location does not guarantee that it is important.
It is hard to determine where θ 2 Ori A lies along the line of sight relative to θ 1 Ori C. The first constraint on this comes from its placement relative to the Veil. The visual extinction of the Trapezium stars is about 1.55 and that of θ 2 Ori A 0.7 (Goudis 1982), i.e. θ 2 Ori A has much less extinction. Since the extinction largely arises from within the Veil ), this appears to argue that θ 2 Ori A lies nearer to the observer than the Veil. However, the H I 21 cm absorption map, which shows the Veil to be thinning to the SW of the Dark Bay (O'Dell & Yusef-Zadeh 2000), which could account for the lower extinction of θ 2 Ori A even if it lies on the far side of the Veil. The O'Dell et al. (1993) study of very high resolution spectra of HeI gives total column densities of the two (-2.9 km s −1 and 5.2 km s −1 ) velocity components of θ 2 Ori A of 2.94×10 13 cm −2 , (individual values 1.70×10 13 cm −2 and 1.24×10 13 cm −2 respectively) and for θ 1 Ori C 1.45×10 13 cm −2 . The average of the four Trapezium stars is 1.60×10 13 ±0.34×10 13 cm −2 .
These numbers indicate that the Trapezium HeI components are very similar to those of the redder component in θ 2 Ori A. This suggest, but does not argue strongly for θ 2 Ori A being on the far side. However, the column density in θ 1 Ori C's HeI is high for the amount of extinction, when compared with the Trapezium stars, which may be due to different conditions of ionization.
In a recent study O'Dell et al. (2017) of the ionization across the Huygens Region, they concluded that most of the MIF ionization northwest of the Bright Bar was caused by θ 1 Ori C, and that to the southeast by θ 2 Ori A. However, they found that proplyds near the Bright Bar were ionized by θ 2 Ori A and some as close as half the distance in the plane of the sky were ionized by both θ 2 Ori A and θ 1 Ori C (one sees ionization fronts in the same proplyd facing both stars ). Since the emission line proplyds clearly lie within the cavity between the MIF and the Veil, these objects become a strong argument that θ 2 Ori A lies on the far side of the Veil and its radiation can affect our Component I. It is prudent to investigate the sensitivity of our modeling efforts to θ 2 Ori A.

Models considering both θ 1 Ori C and θ 2 Ori A
We used the best-fitting model to Component I (R (cm)= 10 18.6 , n (cm −3 ) = 10 2.65 ) and reran the model using an SED consistent with θ 2 Ori A. We assume θ 2 Ori A is a typical O9.5V star with a stellar temperature of 34,600K and Q(H) = 3.63×10 48 photons s −1 (O'Dell et al. 2017). Since θ 2 Ori A is about 0.41 parsecs away from the Trapezium in the plane of the sky, the distance of θ 2 Ori A from the line of sight to the Trapezium is highly uncertain. We assume θ 2 Ori A is between the ionized layer and the MIF (just like the Trapezium). Given this assumption, the closest distance θ 2 Ori A could be from the line of sight connecting the Earth to the Trapezium is 0.41 parsecs (if θ 2 Ori A were touching Component I). If θ 2 Ori A were the same distance from Component I as the Trapezium (1.0 parsecs, for the best model), then θ 2 Ori A, the Trapezium, and the line of sight would form a right angle, and the θ 2 Ori A distance to the line of sight would be 1.08 parsecs.
In addition to the distance, the direction the ionizing flux impinges Component I must also be considered. Since radiation from θ 2 Ori A strikes the ionized layer at an angle (relative to the straight line connecting Trapezium to Earth), the flux will be reduced by the cosine of the angle between the Trapezium-Earth line of sight and a line connecting θ 2 Ori A to the illuminated face of Component I directly in front of the Trapezium. We therefore calculate the ionizing photon flux, φ(H) = Q(H)/4πR 2 , where R 2 = 0.41 2 + y 2 is the distance of θ 2 Ori A from the illumination point and y is the perpendicular distance of θ 2 Ori A to Component I, allowed to vary between 0 and 1 parsec. We then multiplied by the ionizing flux by cos(θ), where cos(θ) = y/R. Given these geometric assumptions, the minimum and maximum values for the contribution of θ 2 Ori A to φ(H) comes out to 10 10.33 and 10 10.83 photons cm −2 s −1 .
Given the range of φ(H) considered, we found a small, but ultimately inconsequential change to the model results for Component I. Including θ 2 Ori A reduced the [N II] surface brightness from S([N II]) = 3.97x10 -3 ergs s −1 cm −2 sr −1 (without θ 2 Ori A) to (1.8-2.7)x10 -3 ergs s −1 cm −2 sr −1 . These slightly smaller values are still well within the observed range. The [O III] surface brightness is identical to before, to within one percent, and the HeI* column density changed by less than 0.01 dex. Overall, the best model for Component I is not significantly affected by including ionizing flux from θ 2 Ori A. Since the model for Component II (the [O I] emitting layer) is based on the results of Component I, both ionized layers are relatively insensitive to θ 2 Ori A.
One way in which the presence of θ 2 Ori A could affect the Veil is the assumed distance of the neutral H I layers (Component III (B) and IV (A)). Including θ 2 Ori A increases the FUV flux, G 0 , by a factor of 2 -4. Increasing G 0 means we would need to increase the distance of Component III and IV by a factor of 1.4 to 2 in order to reproduce the same physical conditions deduced from UV and radio observations (Abel et al. 2016). Overall, these results are expected, given the spectral type of θ 2 Ori A and we would expect the star to be a relatively minor contributor to the overall hydrogen ionizing flux, but perhaps a more important contributor to the FUV.

PROPERTIES OF A CENTRAL COLUMN PASSING FROM THE OBSERVER TO THE ORION MOLECULAR CLOUD
In Table 2 we present the velocities for multiple velocity components as derived in Section 2.2 and other studies cited in column 3. The assignment to components is by similar velocities and expected states of ionization or illumination by θ 1 Ori C. The Backscattered MIF component is not a physical layer, rather, it is the observed velocity of MIF emission scattered in the PDR. We refer to this collection as representative of a Central Column because multiple stars are used for the absorption lines and the sampled emission line regions vary slightly between different studies. Goicoechea et al. (2015) The study of Goicoechea et al. (2015) produced velocity images in CO 2-1, H41α, and [C II] 158 µm. The data relevant to the spatial column discussed in this section are nicely presented in the Trapezium panel of velocities in their Figure 3. There one sees that the H41α line is wide, as caused by thermal broadening, and is centered on 16 km s −1 (2.5 km s −1 LSR). They assign this line to the MIF. The much narrower CO 2-1 line peaks at 28 km s −1 (10 km s −1 LSR) and is assigned to the background PDR and OMC. The strongest [C II] 158 µm line peaks at 28 km s −1 (10 km s −1 LSR) and is assigned to the MIF.

Velocity components from
There is a weaker [C II] 158 µm line component at 19 km s −1 (0.5 km s −1 LSR). This velocity component shows up on their spectrum for the entire Huygens Region as a distinct shoulder on the blue side of the 28 km s −1 line. In Section 3.3.1 of Goicoechea et al. (2015) those authors assign this weaker blue [C II] 158 µm component to a foreground layer that is illuminated by θ 1 Ori C's strong FUV photodissociating continuum, based on of the absence of a similar CO 2-1 component, linking it to the H I foreground layers. We assign it to Component III(B) because of the good agreement in velocity with the other members of this component.
The detection of [C II] 158 µm emission is consistent with the observed UV absorption spectra and our theoretical modeling of the neutral Veil layers. Sofia et al. (2004) detected the C + 232.5 nm line in absorption towards the Trapezium, deducing a total column density along the line of sight of about 10 17.8 cm −2 . The velocity of C + corresponds with the velocity of Veil Component III (B) and IV (A), indicating the observed C + is associated with the neutral layers along the line of sight. Given the total H I column density towards the Trapezium (10 21.68 cm −2 , the ratio of ionized carbon to hydrogen is about 1.4x10 −4 , which is consistent with the observed ratio of C/H in the ISM. Therefore, almost all carbon along this sightline is in the form of C + , and is associated with the neutral Veil layers. The models computed in previous studies of the Veil (Abel et al. (2016) and references therein) assumes this ratio of C/H and stopped the calculation at the observed HI column density.

Possible velocity components associated with a central shocked bubble created by θ 1 Ori C's stellar wind
The V blue and V new components are assigned to the boundaries of a 'Central Bubble' of hot shocked gas created by the strong stellar wind from θ 1 Ori C. The boundary closer to the observer is the V blue component and the further boundary is the V new component. This V blue component is moving away from θ 1 Ori C at about 34±6 km s −1 , while the V new component on the far side is moving away from θ 1 Ori C by 10±2 km s −1 for [N II] and is almost stationary with respect to θ 1 Ori C in [O III]. The lower velocities in the far side of the Central Bubble are probably due to the shocked gas moving into higher density gas. The boundaries in the plane of the sky are probably outlined by the High Ionization Arc that encircles the Trapezium region except in the direction of the Orion-S Cloud O'Dell 2018).   Table 2 shows the Heliocentric Velocities for the various ionization layers within the MIF). Distances along the line-of-sight are not to scale, but they are indicated in Table 3. Figure 9. The relative velocity of the MIF components (shown as -11 km s −1 ) is fully consistent with evaporation of material away the PDR (whose -1 km s −1 value is probably uncertain to 1 km s −1 ). The relative velocities of the boundaries of the Wind Blown Cavity are consistent with expectation, with the expansion velocity on further side having been slowed by interaction with the dense MIF material and the nearer side showing a much higher relative velocity. This large velocity approaching velocity is contrary to the expectation that this inner cavity is essentially stationary.
It is not surprising that the well defined foreground Veil components all have negative velocities with respect to the OMC because the ionization of the immediate environment of θ 1 Ori C and the stellar wind from that star will provide an outward push. However, Components III(B) and IV(A) may have velocities inherently different because of their distance from the OMC. As noted before (Abel et al. 2016), the negative velocity differences mean that Components I and II are going to collide with component III(B) at a time determined by the separation of the layers. A detailed study of this collision timescale may tell us about the lifetime of the stellar wind and FUV flux of θ 1 Ori C. How a line of site correlation can exist presents a puzzle since Component I lies about 1.3 pc in the foreground of θ 1 Ori C and the MIF about one-fourth pc beyond θ 1 Ori C.
Could the X-ray emitting hot gas discovered by Güdel et al. (2008) affect Component I along with the other Veil layers? In planetary nebulae, where more complete observations exist, the hot and nebular phases are in rough pressure equilibrium. There are no observations of the hot gas for the regions close to the Trapezium that we study here due to extinction of soft X-rays by the Veil. Arthur (2012) made hydrodynamical simulations of the hot gas with the aim of matching the Güdel observations. Her Figure 2 suggests a gas pressure of nT ∼ 10 6 K cm −3 close to the stars. This is very uncertain since the calculations were not intended to reproduce the properties of the hot gas close to the star cluster, never-the-less they do serve as a rough guide to what might be present. This pressure is similar to the gas pressure in the two Veil components (Component I 10 6.5 K cm −3 and Component 2 10 6.0 K cm −3 ). It is possible that the Veil layers we observe are entrained with the hot gas flow. By comparison, the MIF has a much higher gas pressure, 10 7.5 K cm −3 , so would not be affected by the pressure from the X-ray gas. The MIF pressure would push on the hot bubble, which could in turn push the foreground Veil layers. This would account for the large scale correlation between the two velocities. This process would not explain a V mif and V low correlation on a small scale in the plane of the sky.
6. HOW DO THESE RESULTS RELATE TO THE STRUCTURE OF THE EON? In Table 3 we summarize the distances of the various major components within our Central Column. Most of these are derived from predictions from ionization modeling with observations. The exception to this is the Nearer Central Bubble value of 0.1pc. In this case we have adopted one-half of the diameter of the High Ionization Arc (106 or 0.20 pc diameter).
For reference, we note that the north-south axis of the EON optical window image is about 34. 7 (3.9 pc) and the east-west axis is about 30. 1 (3.3 pc). If the EON is one-half of a spherical bubble, its near side would be about 1.8 pc in front of θ 1 Ori C, a value very near the distance of Veil Component III(B). However, there is little to support the assumption of that geometry and θ 1 Ori C is not near the center of the EON, lying only 0.4 pc from the north-east corner of the EON.
A very relevant recent paper Pabst et al. (2019) examines the EON structure using [C II] 158 µm high velocity resolution images obtained with the airborne observatory SOFIA. That paper does not use the term EON, rather, they call it the Veil Bubble because of the interpretation they apply to their important new observations. Of particular interest is their conclusion that the EON actually extends as far as 2500 , which corresponds to a diameter of 4.6 pc at the distance we adopt (383±3 pc, (Kounkel et al. 2017;Großschedl et al. 2019) or 5.0 pc at the distance (414±7, (Menten et al. 2007)) they adopt. Position-velocity [C II] 158 µm emission plots across the EON (both north-south and east-west) indicate an irregular and weak [C II] 158 µm component with a velocity of expansion of 13 km s −1 away from the MIF component. In the vicinity of the Huygens region this component is displaced about 8 km s −1 , in agreement with our Central Column [C II] 158 µm displacement of 9 km s −1 derived from the Goicoechea et al. (2015) study. Their interpretation is that the displaced [C II] 158 µm component is the boundary of a 2 pc radius shell encompassing a high temperature region that is the source of the X-ray emission. A smaller region (which they do not detect but include in their model) contains a sphere of gas evacuated by the fast stellar wind of θ 1 Ori C. They do not use the 21 cm H I data because of the complexity of the emission and absorption results in this velocity region.
In our study of the Central Column near θ 1 Ori C we show strong evidence for linking the expanding [C II] 158 µm material to Component III(B), which is seen in H I, H 2 , Ca II, and Na I. In addition, we establish that Component IV(A) must lie outside of Component III(B) and that the well defined Component I must lie between both Component III(B) and the ionized Component II. We also argue that the High Ionization Arc is the edge of the sphere of low density material swept up by the stellar wind from θ 1 Ori C. This not to argue that the Pabst et al. (2019) model is wrong, rather, that the real situation is more complex than pictured from only the [C II] 158 µm line and complementary infrared data. The evidence for Component III(B) being part of a curved shell of about the same spatial size as the EON is strong, while we cannot comment on the form of Component I, Component II, and Component IV(A).

CONCLUSIONS
The foreground Veil associated with the Huygens Region of the Orion Nebula has recently been shown to have multiple new components.
Comparison of the observed and calculated surface brightnesses for [N II], [O III], and [O I] indicate the conditions within the two Veil layers closest to the Trapezium stars. The strength of the He I absorption lines in the Trapezium stars allow further refinement of the conditions within the Veil.
We establish that there are two faint but real ionized layers on either side of Trapezium and these probably represent the two sides of the hot central cavity caused by the strong stellar wind from θ 1 Ori C. The high pressure of this bubble is the driving force that causes the Veil components to be moving with a blueshift relative to the host Orion Molecular Cloud. This also provides the explanation for the radial velocities of the ionized Veil layer closest to θ 1 Ori C being correlated to the velocities of the MIF.
We calculated the effects of θ 2 Ori A on the deduced physical conditions, and find a negligible contribution to the ionized layers, but the FUV contribution to the overall SED could lead to a factor of two increase in the distance of the neutral layers from the Trapezium. We confirm that one layer (Component III(B)) in our Central Column is probably a part of a recently recognized expanding shell of material covering the entire EON.