PDRs4All VI: Probing the Photochemical Evolution of PAHs in the Orion Bar Using Machine Learning Techniques

[Abridged] JWST observations of the Orion Bar have shown the incredible richness of PAH bands and their variation on small scales. We aim to probe the photochemical evolution of PAHs across the key zones of the photodissociation region (PDR) that is the Orion Bar using unsupervised machine learning. We use NIRSpec and MIRI IFU data from the JWST ERS Program PDRs4All. We lever bisecting k-means clustering to generate detailed spatial maps of the spectral variability in several wavelength regions. We discuss the variations in the cluster profiles and connect them to the local physical conditions. We interpret these variations with respect to the key zones: the HII region, the atomic PDR zone, and the three dissociation fronts. The PAH emission exhibits spectral variation that depends strongly on spatial position in the PDR. We find the 8.6um band to behave differently than all other bands which vary systematically with one another. We find uniform variation in the 3.4-3.6um bands and 3.4/3.3 intensity ratio. We attribute the carrier of the 3.4-3.6um bands to a single side group attached to very similarly sized PAHs. Cluster profiles reveal a transition between characteristic profiles classes of the 11.2um feature from the atomic to the molecular PDR zone. We find the carriers of each of the profile classes to be independent, and reason the latter to be PAH clusters existing solely deep in the molecular PDR. Clustering also reveals a connection between the 11.2 and 6.2um bands; and that clusters generated from variation in the 10.9-11.63um region can be used to recover those in the 5.95-6.6um region. Clustering is a powerful tool for characterizing PAH variability on both spatial and spectral scales. For individual bands as well as global spectral behaviours, we find UV-processing to be the most important driver of the evolution of PAHs and their spectral signatures in the Orion Bar.Comment: 28 pages, 20 figures, submitted to A&

A far-ultraviolet-driven photoevaporation flow observed in a protoplanetary disk.

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photodissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, which affects planet formation within the disks. We report James Webb Space Telescope and Atacama Large Millimeter Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modeling their kinematics and excitation allowed us to constrain the physical conditions within the gas. We quantified the mass-loss rate induced by the FUV irradiation and found that it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk

Observational calibration of the projection factor of Cepheids: IV. Period-projection factor relation of Galactic and Magellanic Cloud Cepheids

International audienceContext. The Baade-Wesselink (BW) method, which combines linear and angular diameter variations, is the most common method to determine the distances to pulsating stars. However, the projection factor, p-factor, used to convert radial velocities into pulsation velocities, is still poorly calibrated. This parameter is critical on the use of this technique, and often leads to 5–10% uncertainties on the derived distances.Aims. We focus on empirically measuring the p-factor of a homogeneous sample of 29 LMC and 10 SMC Cepheids for which an accurate average distances were estimated from eclipsing binary systems.Methods. We used the SPIPS algorithm, which is an implementation of the BW technique. Unlike other conventional methods, SPIPS combines all observables, i.e. radial velocities, multi-band photometry and interferometry into a consistent physical modelling to estimate the parameters of the stars. The large number and their redundancy insure its robustness and improves the statistical precision.Results. We successfully estimated the p-factor of several Magellanic Cloud Cepheids. Combined with our previous Galactic results, we find the following P−p relation: −0.08± 0.04(log P−1.18) + 1.24± 0.02. We find no evidence of a metallicity dependent p-factor. We also derive a new calibration of the period-radius relation, log R = 0.684± 0.007(log P−0.517) + 1.489± 0.002, with an intrinsic dispersion of 0.020. We detect an infrared excess for all stars at 3.6 μm and 4.5 μm, which might be the signature of circumstellar dust. We measure a mean offset of Δm3.6 = 0.057 ± 0.006 mag and Δm4.5 = 0.065 ± 0.008 mag.Conclusions. We provide a new P−p relation based on a multi-wavelength fit that can be used for the distance scale calibration from the BW method. The dispersion is due to the LMC and SMC width we took into account because individual Cepheids distances are unknown. The new P−R relation has a small intrinsic dispersion: 4.5% in radius. This precision will allow us to accurately apply the BW method to nearby galaxies. Finally, the infrared excesses we detect again raise the issue of using mid-IR wavelengths to derive period-luminosity relation and to calibrate the Hubble constant. These IR excesses might be the signature of circumstellar dust, and are never taken into account when applying the BW method at those wavelengths. Our measured offsets may give an average bias of ~ 2.8% on the distances derived through mid-IR P−L relations

Consistent radial velocities of classical Cepheids from the cross-correlation technique

International audienc

VizieR Online Data Catalog: Classical Cepheids consistent radial velocities (Borgniet+, 2019)

VizieR On-line Data Catalog: J/A+A/631/A37. Originally published in: 2019A&A...631A..37BOur sample is made up of 64 Classical Galactic Cepheids with pulsation periods in the range 2 to 68 days. We provide the full details of our sample in Table A.1. We present both the six correlation templates built to derive tailored cross-correlation functions (CCFs), the CCFs themselves, and the derived radial velocity and line profile observable time series. All data are provided with the corresponding Cepheid name, the spectrograph used to observe the target, the Modified Julian Day (MJD) of the observation and the observation program identifier. The data are provided over three different wavelength ranges: "blue" from 390 to 498nm; "green" from 450 to 680nm; and "red" from 570 to 880nm, corresponding to the letters "b", "g", and "r" (first identifier within the file names). Depending on this range, the data are computed based on one (for the "blue" and "red" ranges) or four (for the "green" range) different correlation templates, corresponding to different sets of spectral lines selected based on their relative depth: weak, medium, deep and all lines. The corresponding identifiers are the letters "w", "m", "d", and "a" (second identifier within the file names). Finally, the radial velocities and other line profile observables are computed in three different ways for each star, each spectrograph, each wavelength range, and each template: first based on the CCF first moment or centroid, second based on a Gaussian fit of the CCF, and third based on a BiGaussian fit of the CCF. The corresponding identifiers are the letters "c", "g" and "b" (third identifier within the rv file names). (28 data files)

A thin shell of ionized gas as the explanation for infrared excess among classical Cepheids

International audienceContext. The infrared (IR) excess of classical Cepheids is seldom studied and poorly understood despite observational evidence and the potential for its contribution to induce systematics on the period-luminosity (PL) relation used in the calibration of the extragalactic distance scale. Aims: This study aims to understand the physical origin of the IR excess found in the spectral energy distribution (SED) of 5 Cepheids: RS Pup (P = 41.46d), ζ Gem (P = 10.15d), η Aql (P = 7.18d), V Cen (P = 5.49d) and SU Cyg (P = 3.85d). Methods: A time series of atmospheric models along the pulsation cycle were fitted to a compilation of data, including optical and near-IR photometry, Spitzer spectra (secured at a specific phase), interferometric angular diameters, effective temperature estimates, and radial velocity measurements. Herschel images in two bands were also analyzed qualitatively. In this fitting process, based on the SPIPS algorithm, a residual was found in the SED, whatever the pulsation phase, and for wavelengths larger than about 1.2 μm, which corresponds to the so-determined infrared excess of Cepheids. This IR excess was then corrected from interstellar medium absorption in order to infer the presence (or absence) of dust shells and was, ultimately, used in order to fit a model for a shell of ionized gas. Results: For all Cepheids, we find a continuum IR excess increasing up to approximately -0.1 magnitudes at 30 μm, which cannot be explained by a hot or cold dust model of CircumStellar Environment (CSE). However, a weak but significant dust emission at 9.7 μm is found for ζ Gem, η Aql and RS Pup, while clear interstellar clouds are seen in the Herschel images for V Cen and RS Pup. We show, for the first time, that the IR excess of Cepheids can be explained by free-free emission from a thin shell of ionized gas, with a thickness of ≃15% of the star radius, a mass of 10-9-10-7M⊙ and a temperature ranging between 3500 and 4500 K. Conclusions: The presence of a thin shell of ionized gas around Cepheids must be tested with interferometers operating in the visible or mid-IR, or using radio telescopes. The impact of such CSEs of ionized gas on the PL relation of Cepheids also calls for further investigation

CRIRES high-resolution infrared spectroscopy of the long-period Cepheid l Carinae

International audienceContext. The dynamical structure of the atmosphere of Cepheids has been well studied in the optical. Several authors have found very interesting spectral features in the J band, but little data have been secured beyond 1.6 μm. However, such observations can probe different radial velocities and line asymmetry regimes, and are able to provide crucial insights into stellar physics.Aims. Our goal was to investigate the infrared line-forming region in the K-band domain, and its impact on the projection factor and the k-term of Cepheids.Methods. We secured CRIRES observations for the long-period Cepheid l Car, with a focus on the unblended spectral line NaI 2208.969 nm. We measured the corresponding radial velocities (by using the first moment method) and the line asymmetries (by using the bi-Gaussian method). These quantities are compared to the HARPS visible spectra we previously obtained on l Car.Results. The optical and infrared radial velocity curves show the same amplitude (only about 3% of difference), with a slight radial velocity shift of about 0.5 ± 0.3 km s−1 between the two curves. Around the minimum radius (phase ≃ 0.9) the visible radial velocity curve is found in advance compared to the infrared one (phase lag), which is consistent with an infrared line forming higher in the atmosphere (compared to the visible line) and with a compression wave moving from the bottom to the top of the atmosphere during maximum outward velocity. The asymmetry of the K-band line is also found to be significantly different from that of the optical line

PDRs4All VI: Probing the Photochemical Evolution of PAHs in the Orion Bar Using Machine Learning Techniques

International audience[Abridged] JWST observations of the Orion Bar have shown the incredible richness of PAH bands and their variation on small scales. We aim to probe the photochemical evolution of PAHs across the key zones of the photodissociation region (PDR) that is the Orion Bar using unsupervised machine learning. We use NIRSpec and MIRI IFU data from the JWST ERS Program PDRs4All. We lever bisecting k-means clustering to generate detailed spatial maps of the spectral variability in several wavelength regions. We discuss the variations in the cluster profiles and connect them to the local physical conditions. We interpret these variations with respect to the key zones: the HII region, the atomic PDR zone, and the three dissociation fronts. The PAH emission exhibits spectral variation that depends strongly on spatial position in the PDR. We find the 8.6um band to behave differently than all other bands which vary systematically with one another. We find uniform variation in the 3.4-3.6um bands and 3.4/3.3 intensity ratio. We attribute the carrier of the 3.4-3.6um bands to a single side group attached to very similarly sized PAHs. Cluster profiles reveal a transition between characteristic profiles classes of the 11.2um feature from the atomic to the molecular PDR zone. We find the carriers of each of the profile classes to be independent, and reason the latter to be PAH clusters existing solely deep in the molecular PDR. Clustering also reveals a connection between the 11.2 and 6.2um bands; and that clusters generated from variation in the 10.9-11.63um region can be used to recover those in the 5.95-6.6um region. Clustering is a powerful tool for characterizing PAH variability on both spatial and spectral scales. For individual bands as well as global spectral behaviours, we find UV-processing to be the most important driver of the evolution of PAHs and their spectral signatures in the Orion Bar