Statistical Shape Analysis of Galactic Hii Regions

Hii regions are diffuse nebulae of ionised hydrogen, excited by the extreme ultraviolet emission from massive stars. Due to the embedded nature of massive star formation, there are many observational difficulties involved when investigating such stars. Hii regions, however, are readily observed via their infrared and radio emission. As such, they highlight the location of their massive star sources. Furthermore, Hii region properties are directly resultant of their progenitors and environment. The overall aim of the work presented herein, is to determine whether statistical shape analysis of observational and numerically modelled Hii region data can be used to probe the associated astrophysical properties. Radio continuum and computer simulated synthetic images of Hii regions were analysed using the shape extraction and statistical comparison methods constructed in this work. For the radio data, six morphological groups were identified. Visual inspection and quantitative ordinance techniques confirmed that the shape analysis and grouping procedure were working as intended. It was found that in the first Galactic quadrant, location is mostly independent of group, with a small preference for regions of similar Galactic longitudes to share common morphologies. The shapes are homogeneously distributed across Galactocentric distance and latitude. One group contained regions that are all younger than 0.5 Myr and ionised by relatively low- to intermediate-mass sources. Those in another group are all driven by intermediate- to high-mass sources. One group was distinctly separated from the other five and contained regions at the surface brightness detection limit for the survey. The hierarchical procedure employed was most sensitive to the spatial sampling resolution used, which is determined for each region from its heliocentric distance. The numerical Hii region data was the result of photoionisation and feedback of a 34 M⊙ star, in a 1000 M⊙ cloud. Synthetic observations (SOs) were provided, comprising four evolutionary snapshots (0.1, 0.2, 0.4 and 0.6Myr), and multiple viewing projection angles. The shape analysis results provided conclusive evidence of the efficacy of the numerical simulations. When comparing the shapes of the synthetic regions to their observational counterparts, the SOs were grouped in amongst the Galactic Hii regions by the hierarchical procedure. There was also an association between the evolutionary distribution of regions of the respective samples. This suggested that this method could be further developed for classification of the observational regions by using the synthetic data, with its well defined parameters

Shape Analysis of HII Regions — II. Synthetic Observations

The statistical shape analysis method developed for probing the link between physical parameters and morphologies of Galactic HII regions is applied here to a set of synthetic observations (SOs) of a numerically modelled HII region. The systematic extraction of HII region shape, presented in the first paper of this series, allows for a quantifiable confirmation of the accuracy of the numerical simulation, with respect to the real observational counterparts of the resulting SOs. A further aim of this investigation is to determine whether such SOs can be used for direct interpretation of the observational data, in a future supervised classification scheme based upon HII region shape. The numerical HII region data was the result of photoionisation and radiation pressure feedback of a 34 M_\bigodot star, in a 1000 M_\bigodot cloud. The SOs analysed herein comprised four evolutionary snapshots (0.1, 0.2, 0.4 and 0.6 Myr), and multiple viewing projection angles. The shape analysis results provided conclusive evidence of the efficacy of the numerical simulations. When comparing the shapes of the synthetic regions to their observational counterparts, the SOs were grouped in amongst the Galactic HII regions by the hierarchical clustering procedure. There was also an association between the evolutionary distribution of regions and the respective groups. This suggested that the shape analysis method could be further developed for morphological classification of HII regions by using a synthetic data training set, with differing initial conditions of well-defined parameters

The STAR-MELT Python package for emission line analysis of YSOs

We introduce the STAR-MELT Python package that we developed to facilitate the analysis of time-resolved emission line spectroscopy of young stellar objects. STAR-MELT automatically extracts, identifies and fits emission lines. We summarise our analysis methods that utilises the time domain of high-resolution stellar spectra to investigate variability in the line profiles and corresponding emitting regions. This allows us to probe the innermost disc and accretion structures of YSOs. Local temperatures and densities can be determined using Boltzmann statistics, the Saha equation, and the Sobolev large velocity gradient approximation. STAR-MELT allows for new results to be obtained from archival data, as well as facilitating timely analysis of new data as it is obtained. We present the results of applying STAR-MELT to three YSOs, using spectra from UVES, XSHOOTER, FEROS, HARPS, and ESPaDOnS. We demonstrate what can be achieved for data with disparate time sampling, for stars with different inclinations and variability types. For EX Lupi, we confirm the presence of a localised and stable stellar-surface hotspot associated with the footprint of the accretion column. For GQ Lupi A, we find that the maximum infall rate from an accretion column is correlated with lines produced in the lowest temperatures. For CVSO109 we investigate the rapid temporal variability of a redshifted emission wing, indicative of rotating and infalling material in the inner disc. Our results show that STAR-MELT is a useful tool for such analysis, as well as other applications for emission lines.<br/

The Radial Distribution and Excitation of H2 around Young Stars in the HST-ULLYSES Survey

The spatial distribution and evolution of gas in the inner 10 au of protoplanetary disks form the basis for estimating the initial conditions of planet formation. Among the most important constraints derived from spectroscopic observations of the inner disk are the radial distributions of the major gas phase constituents, how the properties of the gas change with inner disk dust evolution, and how chemical abundances and excitation conditions are influenced by the high-energy radiation from the central star. We present a survey of the radial distribution, excitation, and evolution of inner disk molecular hydrogen (H$_{2}$) obtained as part of the $HST$/ULLYSES program. We analyze far-ultraviolet spectroscopy of 71 (63 accreting) pre-main sequence systems in the ULLYSES DR5 release to characterize the H$_{2}$ emission lines, H$_{2}$ dissociation continuum emission, and major photochemical/disk evolution driving UV emissions (Ly$\alpha$, UV continuum, and C IV). We use the widths of the H$_{2}$ emission lines to show that most fluorescent H$_{2}$ arises between 0.1 - 1.4 au from the parent star, and show positive correlations of the average emitting radius with the accretion luminosity and with the dust disk mass. We find a strong correlation between H$_{2}$ dissociation emission and both the accretion-dominated Ly$\alpha$ luminosity and the inner disk dust clearing, painting a picture where water molecules in the inner 3 au are exposed to and dissociated by strong Ly$\alpha$ emission as the opacity of the inner disk declines with time.Comment: 19 pages, 11 figures. Accepted to the Astronomical Journa

A Magnetically Driven Disk Wind in the Inner Disk of PDS 70

PDS 70 is so far the only young disk where multiple planets have been detected by direct imaging. The disk has a large cavity when seen at submillimeter and near-infrared wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disk shaped by two giant planets, and in particular to test whether disk winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O i] 6300 Å line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O i] line luminosity is high with respect to the accretion luminosity when compared to a large sample of disks with cavities in nearby star-forming regions. The FWHM and blueshifted peak of the [O i] line suggest an emission in a region very close to the star, favoring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He i 10830 Å line, which is unusual for low accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner-disk wind could also have had a significant contribution to clearing the inner disk.</p

A magnetically driven disc wind in the inner disc of PDS 70

PDS 70 is so far the only young disc where multiple planets have been detected by direct imaging. The disc has a large cavity when seen at sub-mm and NIR wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disc shaped by two giant planets, and in particular to test whether disc winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O I] 6300AA line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O I] line luminosity is high with respect to the accretion luminosity when compared to a large sample of discs with cavities in nearby star-forming regions. The FWHM and blue-shifted peak of the [O I] line suggest an emission in a region very close to the star, favouring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He I 10830AA line, which is unusual for low-accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner disc wind could also have had a significant contribution to clearing the inner-disc.<br/

A Magnetically Driven Disk Wind in the Inner Disk of PDS 70

PDS 70 is so far the only young disk where multiple planets have been detected by direct imaging. The disk has a large cavity when seen at submillimeter and near-infrared wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disk shaped by two giant planets, and in particular to test whether disk winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O i] 6300 Å line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O i] line luminosity is high with respect to the accretion luminosity when compared to a large sample of disks with cavities in nearby star-forming regions. The FWHM and blueshifted peak of the [O i] line suggest an emission in a region very close to the star, favoring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He i 10830 Å line, which is unusual for low accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner-disk wind could also have had a significant contribution to clearing the inner disk

Variability in IC5070: two young stars with deep recurring eclipses

Investigating the structure and properties of the innermost parts of protoplanetary accretion disks on sub-AU scales is currently only possible via indirect methods. One option to map the planet-forming zone is to search for occultations of the central young stellar object (YSO) by circumstellar material, e.g., warps or clumps in the inner disks. Such disk eclipses typically last hours to days (Cody et al. 2014) and have been identified in massive HAeBe stars such as UX Ori (Herbst & Shevchenko 1999) and lower mass objects such as AA Tau (Bouvier et al. 1999). Of particular interest are quasi-periodic dimming events. They allow distance determinations of the occulting material from the central star. In such cases the actual azimuthal physical extent of the material can be determined from the duration of the dimming event relative to the period. Observations over several periods enable investigations into temporal changes in the line of sight column density distribution, and multi-wavelength data allows us to probe the dust scattering properties. Our citizen science project HOYS-CAPS (Froebrich et al. 2018) aims to identify such periodic dimming events around YSOs. We used this data-set to search for periodic signatures in light-curves from YSOs in the Pelican nebula (IC 5070). For this field we have ~200 individual observations in the V, R, and I-band filters, distributed over ~800 days. Hence, the average cadence is 4 days, but the most frequent gap (30%) between subsequent observations is 2 days. Observations are usually taken as 8 × 2 minutes integrations in all filters to achieve a consistent S/N

Empirical determination of the lithium 6707.856 Å wavelength in young stars

Absorption features in stellar atmospheres are often used to calibrate photocentric velocities for the kinematic analysis of further spectral lines. The Li feature at ∼6708 Å is commonly used, especially in the case of young stellar objects, for which it is one of the strongest absorption lines. However, this complex line comprises two isotope fine-structure doublets. We empirically measured the wavelength of this Li feature in a sample of young stars from the PENELLOPE/VLT programme (using X-shooter, UVES, and ESPRESSO data) as well as HARPS data. For 51 targets, we fit 314 individual spectra using the STAR-MELT package, resulting in 241 accurately fitted Li features given the automated goodness-of-fit threshold. We find the mean air wavelength to be 6707.856 Å, with a standard error of 0.002 Å (0.09 km s−1), and a weighted standard deviation of 0.026 Å (1.16 km s−1). The observed spread in measured positions spans 0.145 Å, or 6.5 km s−1, which is higher by up to a factor of six than the typically reported velocity errors for high-resolution studies. We also find a correlation between the effective temperature of the star and the wavelength of the central absorption. We discuss that exclusively using this Li feature as a reference for photocentric velocity in young stars might introduce a systematic positive offset in wavelength to measurements of further spectral lines. If outflow tracing forbidden lines, such as [O