3D maps of the extinction density in the Galaxy can be built through the inversion of catalogues of distance-extinction pairs for individual target stars. The spatial resolution of the maps that can be achieved increases with the spatial density of the targets, and subsequently with the combination of catalogues. However, this requires their careful inter-calibration. Our aim is to develop methods of inter-calibration of two different catalogues. We used as reference a spectrophotometric catalogue. A principal component analysis was performed in G,GB,GR,J,H,K multi-colour space for the second catalog. The subspace constituted by the two first components was split into cells in which we estimated deviations from the reference. The deviations were computed using all targets from the reference located at a short spatial distance of each secondary target. Corrections and filtering were deduced for each cell in the multi-colour space. We applied the technique to two different datasets: on the one hand, the spectrophotometric catalogue, and, on the other hand, a catalogue of extinctions based on photometry of Gaia eDR3 and 2MASS. After calibration, the dispersion of the extinction among neighbouring points in the second catalogue is reduced, regardless of whether reference targets are present locally. Weak structures are then more apparent. The extinction of high Galactic latitude targets is more tightly correlated with the dust emission, a property acquired from the first catalogue. A hierarchical inversion technique was applied to the two merged inter-calibrated catalogues to produce 3D extinction density maps corresponding to different volumes and maximum spatial resolution. The maximum resolution is 10pc for a 3000x3000x800pc3 volume around the Sun, and the maximum size of the maps is 10x10x0.8 kpc3 for a resolution of 50pc. Maps can be downloaded or used by means of on-line tools.

State-of-the-art molecular databases and realistic global atmospheric models allow to predict accurate atmospheric transmittance spectra. The TAPAS online service provides atmospheric transmittance spectra of the most important species as well as Rayleigh extinction, adapted to any observing location, date, and direction. We describe recent updates, improvements, and additional tools. TAPAS is interpolating in location in the atmospheric profiles of temperature, pressure, H2O, O2 and O3 extracted from the ECMWF and produced by Arletty, a Data Terra/AERIS/ESPRI product, supplemented by auxiliary climatological models for additional species. The transmittance spectra are computed with the LBLRTM code. The resolution is chosen to ensure that the shapes of all absorption lines are reproduced for each species. Major improvements are: - the extension of the wavelength range in the near-UV down to 300 nm; in the near-IR up to 3500 nm; - the use of HITRAN2020; - the addition of NO2, to complement H2O, O2, O3, N2O, CO2 and CH4; - an increased accessibility; -the possibility to force the total H2O column to match the one measured at the observatory. We show O3 absorption in the near-UV and near-IR and NO2 absorption in the visible. We illustrate the quality of TAPAS by means of comparisons between models and ESO/VLT/CRIRES recorded spectra of a hot star with spectral resolution 130,000, in two intervals with strong H2O, NO2, CO2 and CH4 absorption. We describe the measurement of an instrumental Line Spread Function (LSF) based on TAPAS O2 lines and a method using the Singular Value Decomposition technique that can be made entirely automated. TAPAS can be used to improve the wavelength and to measure the LSF in regions where telluric features are present. The extended wavelength range will be particularly useful for future or recent spectrographs in the near-UV and in the near-infrared.

The planetary nebula NGC 6720, also known as the ``Ring Nebula", is one of the most iconic examples of nearby planetary nebulae whose morphologies present a challenge to our theoretical understanding of the processes that govern the deaths of most stars in the Universe that evolve on a Hubble time. We present new imaging with JWST of the central star of this planetary nebula (CSPN) and its close vicinity, in the near- to mid-IR wavelength range. We find the presence of a dust cloud around the CSPN, both from the spectral energy distribution at wavelengths >~5 micron, as well as radially-extended emission in the 7.7, 10 and 11.3 micron images. From modeling of these data, we infer that the CSPN has a luminosity of 310 Lsun and is surrounded by a dust cloud with a size of ~2600 au, consisting of relatively small amorphous silicate dust grains (radius ~0.01 micron) with a total mass of 1.9 x 10^(-6) M_earth. However, our best-fit model shows a significant lack of extended emission at 7.7 micron -- we show that such emission can arise from a smaller (7.3 x 10^(-7) M_earth) but uncertain mass of (stochastically-heated) ionized PAHs. However, the same energetic radiation also rapidly destroys PAH molecules, suggesting that these are most likely being continuously replenished, via the outgassing of cometary bodies and/or the collisional grinding of planetesimals. We also find significant photometric variability of the central source that could be due to the presence of a close dwarf companion of mass < ~0.1 Msun.

We report in this paper circularly and linearly polarized observations of the young active M dwarfs AU Mic and EV Lac with the near-infrared SPIRou spectropolarimeter at the Canada-France-Hawaii Telescope, collected from August to October 2023 over a few rotation cycles of both stars. Applying Least-Squares Deconvolution (LSD) to our spectra, we detected Zeeman signatures in circular (Stokes V) and linear (Stokes QU) polarization, and Zeeman broadening in unpolarized (Stokes I) LSD profiles, all exhibiting clear rotational modulation. Using the stellar surface tomographic technique of Zeeman-Doppler imaging on our sets of observations, along with a simple parametric description of how the small-scale and large-scale fields relate to each other, we recovered the magnetic topologies of AU Mic and EV Lac successively from LSD Stokes V, Stokes IV and Stokes IVQU profiles, to investigate how the reconstructed maps evolve as we provide more information, and ultimately infer reliable magnetic maps of both stars. We find that AU Mic hosts a fairly simple and mostly poloidal large-scale field aligned with the rotation axis within about 10deg, whereas that of EV Lac is more complex, stronger and less axisymmetric. Both stars feature intense small-scale fields, of about 4 kG for AU Mic and 6 kG for EV Lac when averaged over the whole stellar surface. Stokes QU Zeeman signatures allow one to reconstruct stellar magnetic fields more reliably, and are especially useful for stars with more complex fields and low vsini's like EV Lac.

More than 500 diffuse interstellar bands (DIBs) have been observed in astronomical spectra, and their signatures and correlations in different environments have been studied over the past decades to reveal clues about the nature of the carriers. We compare the equivalent widths of the DIBs, normalized to the amount of reddening, E_B-V, to search for anti-correlated DIB pairs using a data sample containing 54 DIBs measured in 25 sight lines. This data sample covers most of the strong and commonly detected DIBs in the optical region, and the sight lines probe a variety of ISM conditions. We find that 12.9% of the DIB pairs are anti-correlated, and the lowest Pearson correlation coefficient is r_norm ~ -0.7. We revisit correlation-based DIB families and are able to reproduce the assignments of such families for the well-studied DIBs by applying hierarchical agglomerative and k-means clustering algorithms. We visualize the dissimilarities between DIBs, represented by 1 - r_norm, using multi-dimensional scaling (MDS). With this representation, we find that the DIBs form a rather continuous sequence, which implies that some properties of the DIB carriers are changing gradually following this sequence. We also find at that least two factors are needed to properly explain the dissimilarities between DIBs. While the first factor may be interpreted as related to the ionization properties of the DIB carriers, a physical interpretation of the second factor is less clear and may be related to how DIB carriers interact with surrounding interstellar material.

Planetary nebulae are sites where ejected stellar material evolves into complex molecules, but the precise physical conditions and chemical routes that govern these processes are unclear. The presence of abundant carbon-rich molecules in O-rich environments poses particular challenges. Here we report the first detection of methyl cation (CH3+) in any planetary nebula, observed in the O-rich nebula NGC 6302 using JWST MIRI/MRS observations. CH3+ is a key driver of organic chemistry in UV-irradiated environments. Spatially resolved observations reveal that CH3+ is co-located with 12CO, H2, H II, HCO+, and Polycyclic aromatic hydrocarbons (PAHs). LTE modelling of the CH3+ emission yields excitation temperatures of 500-800K in the inner bubble and torus, rising to 1000-2000K in the outer bubble of NGC 6302, with column densities ranging from ~10^11 to 10^13 cm^-2. This detection demonstrates that hydrocarbon radical chemistry must be incorporated into planetary nebulae chemical models. Further near-IR observations are crucial to map different chemical networks operating in these environments.

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations.

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 photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding 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.

Three-dimensional (3D) maps of Galactic interstellar dust are a tool for a wide range of uses. We aim to construct 3D maps of dust extinction in the Local Arm and surrounding regions. Gaia EDR3 photometric data were combined with 2MASS measurements to derive extinction towards stars with accurate photometry and relative uncertainties on parallaxes of less than 20%. We applied our hierarchical inversion algorithm adapted to inhomogeneous spatial distributions of target stars to this catalogue of extinctions. We present the updated 3D dust extinction distribution and provide an estimate of the error on integrated extinctions from the Sun to each area in the 3D map. The computational area is similar to the one of the previous DR2 map, a 6 kpc x 6 kpc x 0.8 kpcAstrophysics volume around the Sun. Due to the addition of fainter target stars, the volume in which the clouds can be reconstructed has increased. Due to the improved accuracy of the parallaxes and photometric data in EDR3, extinctions among neighbouring targets are more consistent, allowing one to reach an increased contrast in the dense areas, while cavity contours are more regular. We show several comparisons with recent results on dust and star distributions. The wavy pattern around the Plane of the dust concentrations is better seen and exists over large regions. Its mean vertical peak-to-peak amplitude is of the order of 300 pc; interestingly, it is similar to the vertical period of the spectacular snail-shaped stellar kinematical pattern discovered in Gaia data. The Gaia EDR3 catalogue allows for a significant improvement of the extinction maps to be made and the hierarchical technique confirms its efficiency for massive datasets. Future comparisons between 3D maps of interstellar matter and stellar distributions may help to understand which mergers or internal perturbations have shaped the Galaxy within the first 3 kpc.

Stellar atmosphere modelling predicts the luminosity and temperature of a star, together with parameters such as the effective gravity and the metallicity, by reproducing the observed spectral energy distribution. Most observational data comes from photometric surveys, using a variety of passbands. We herein present the Python Stellar Spectral Energy Distribution (PySSED) routine, designed to combine photometry from disparate catalogues, fit the luminosity and temperature of stars, and determine departures from stellar atmosphere models such as infrared or ultraviolet excess. We detail the routine's operation, and present use cases on both individual stars, stellar populations, and wider regions of the sky. PySSED benefits from fully automated processing, allowing fitting of arbitrarily large datasets at the rate of a few seconds per star.

M dwarfs are known to host magnetic fields, impacting exoplanet studies and playing a key role in stellar and planetary formation and evolution. Observations revealed the long-term evolution of the large-scale magnetic field reconstructed with Zeeman-Doppler imaging, and a diversity of their topologies. These large-scale magnetic fields only account for a small amount of the unsigned magnetic flux that can be probed by directly modeling the Zeeman broadening of spectral lines in unpolarized spectra. We aim at investigating the long-term behavior of the average small-scale magnetic field of M dwarfs with time, and assess our ability to detect rotational modulation from time series of field measurements derived from unpolarized spectra. We perform fits of synthetic spectra computed with ZeeTurbo to near-infrared high-resolution spectra recorded with SPIRou between 2019 and 2024 in the context of the SLS and SPICE large programs. The analysis is performed on the spectra of 2 partially convective (AD Leo, DS Leo) and 3 fully convective (PM J18482+0741, CN Leo, Barnard star) M dwarfs, along with EV Lac whose mass is close to the fully-convective limit. Our analysis provides measurements of the average small-scale magnetic field, which are compared to longitudinal magnetic field and temperature variation measurements (d$Temp$) obtained from the same data. We were able to detect the rotation period in the small-scale magnetic field series for 4 of the 6 stars in our sample. We find that the average magnetic field can vary by up to 0.3 kG throughout the year (e.g., CN Leo), or of up to 1 kG across rotation phases. The rotation periods retrieved from longitudinal and small-scale magnetic fields are found in agreement within error bars. d$Temp$ measurements are found to anti-correlate with small-scale magnetic field measurements for three stars (EV Lac, DS Leo and Barnard's star).

Context: There have been many attempts to identify families of diffuse interstellar bands (DIBs) with perfectly correlating band strengths. Although major efforts have been made to classify broadly based DIB families and important insights have been gained, no family has been identified with sufficient accuracy or statistical significance to prove that a series of selected DIBs originates from the same carrier. This can be attributed in part to the exclusive use of equivalent widths to establish DIB families. Aims: In a change of strategy, we search for DIBs that are highly correlated in both band strength and profile shape. This approach increases the chance of correlating DIBs being members of one family and originating from the same carrier molecule. We also search for correlations between DIB profile families and atomic interstellar lines, with the goal of further chemically constraining possible DIB carriers. Methods: We adapted the well-known method of time-series alignment to perform a spectral alignment; that is, DIB alignment. In a second step, we analysed the alignment results using a clustering analysis. This method required a statistically significant data set of DIB sight lines. The ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES) data were perfectly suited for this application. Results: We report eight DIB families with correlating strengths and profiles, as well as four previously unreported DIBs in the visual range, found using DIB alignment. All profile family members show Pearson correlation coefficients in band strength higher than 0.9. In particular, we report the 6614 - 6521 AA DIB pair, in which both DIBs show the same triple-peak substructure and an unprecedented band strength Pearson correlation coefficient of 0.9935. The presented approach opens up new perspectives that can guide the laboratory search for DIB carriers.

Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The highly excited OH infrared rotational lines are the tell-tale signs of H2O destruction by FUV. The OH infrared ro-vibrational lines are attributed to chemical excitation via the key reaction O+H=OH+H which seeds the formation of water in the gas-phase. We infer that the equivalent of the Earth ocean's worth of water is destroyed per month and replenished. These results show that under warm and irradiated conditions water is destroyed and efficiently reformed via gas-phase reactions. This process, assisted by diffusive transport, could reduce the HDO/H2O ratio in the warm regions of planet-forming disks.

NGC 6302 is a spectacular bipolar planetary nebula (PN) whose spectrum exhibits fast outflows and highly ionized emission lines, indicating the presence of a very hot central star (~220,000 K). Its infrared spectrum reveals a mixed oxygen and carbon dust chemistry, displaying both silicate and polycyclic aromatic hydrocarbon (PAH) features. Using the JWST Mid-Infrared Instrument (MIRI) and Medium Resolution Spectrometer, a mosaic map was obtained over the core of NGC 6302, covering the wavelength range of 5--28 micron and spanning an area of ~18.5 arcsec x15 arcsec. The spatially resolved spectrum reveals ~200 molecular and ionized lines from species requiring ionisation potentials of up to 205 eV. The spatial distributions highlight a complex structure at the nebula's centre. Highly ionized species such as Mg VII and Si VII show compact structures, while lower-ionization species such as H^+ extend much farther outwards, forming filament-defined rims that delineate a bubble. Within the bubble, the H^+ and H_2 emission coincide, while the PAH emission appears farther out, indicating an ionization structure distinct from typical photodissociation regions, such as the Orion Bar. This may be the first identification of a PAH formation site in a PN. This PN appears to be shaped not by a steady, continuous outflow, but by a series of dynamic, impulsive bubble ejections, creating local conditions conducive to PAH formation. A dusty torus surrounds the core, primarily composed of large (micron-sized) silicate grains with crystalline components. The long-lived torus contains a substantial mass of material, which could support an equilibrium chemistry and a slow dust-formation process.

Although there are some real world applications where the use of variable length representation (VLR) in Evolutionary Algorithm is natural and suitable, an academic framework is lacking for such representations. In this work we propose a family of tunable fitness landscapes based on VLR of genotypes. The fitness landscapes we propose possess a tunable degree of both neutrality and epistasis; they are inspired, on the one hand by the Royal Road fitness landscapes, and the other hand by the NK fitness landscapes. So these landscapes offer a scale of continuity from Royal Road functions, with neutrality and no epistasis, to landscapes with a large amount of epistasis and no redundancy. To gain insight into these fitness landscapes, we first use standard tools such as adaptive walks and correlation length. Second, we evaluate the performances of evolutionary algorithms on these landscapes for various values of the neutral and the epistatic parameters; the results allow us to correlate the performances with the expected degrees of neutrality and epistasis.

The carbonaceous macromolecules imprinting in astronomical spectra the numerous absorptions called Diffuse Interstellar Bands (DIBs) are omnipresent in the Galaxy and beyond and represent a considerable reservoir of organic matter. However, their chemical formulae, formation and destruction sites remain open questions. Their spatial distribution and the local relation to other interstellar species is paramount to unravel their role in the lifecycle of organic matter. Volume density maps bring local instead of line-of-sight distributed information, and allow for new diagnostics. We present the first large-scale volume (3D) density map of a DIB carrier and compare it with an equivalent map of interstellar dust. The DIB carrier map is obtained through hierarchical inversion of about 202,000 measurements of the 8621 nm DIB obtained with the Gaia-RVS instrument. It covers about 4000 pc around the Sun in the Galactic plane. A dedicated interstellar dust map is built based on extinction towards the same target stars. At the 50 pc resolution of the maps, the 3D DIB distribution is found remarkably similar in shape to the 3D distribution of dust. On the other hand, the DIB-to-dust local density ratio increases in low-dust areas. It is also increasing away from the disk, however, the minimum ratio is found to be shifted above the Plane to Z = +50pc. Finally, the average ratio is also surprisingly found to increase away from the Galactic Center. We suggest that the three latter trends may be indications of a dominant contribution of material from the carbon-rich category of dying giant stars to the formation of the carriers. Our suggestion is based on recent catalogues of AGB stars and estimates of their mass fluxes of C-rich and O-rich ejecta.

The domain adaptation of satellite images has recently gained an increasing attention to overcome the limited generalization abilities of machine learning models when segmenting large-scale satellite images. Most of the existing approaches seek for adapting the model from one domain to another. However, such single-source and single-target setting prevents the methods from being scalable solutions, since nowadays multiple source and target domains having different data distributions are usually available. Besides, the continuous proliferation of satellite images necessitates the classifiers to adapt to continuously increasing data. We propose a novel approach, coined DAugNet, for unsupervised, multi-source, multi-target, and life-long domain adaptation of satellite images. It consists of a classifier and a data augmentor. The data augmentor, which is a shallow network, is able to perform style transfer between multiple satellite images in an unsupervised manner, even when new data are added over the time. In each training iteration, it provides the classifier with diversified data, which makes the classifier robust to large data distribution difference between the domains. Our extensive experiments prove that DAugNet significantly better generalizes to new geographic locations than the existing approaches.

Three-dimensional (3D) kinetic maps of the Milky Way interstellar medium are an essential tool in studies of its structure and of star formation. We aim to assign radial velocities to Galactic interstellar clouds now spatially localized based on starlight extinction and star distances from Gaia and stellar surveys. We developed an automated search for coherent projections on the sky of clouds isolated in 3D extinction density maps on the one hand, and regions responsible for CO radio emissions at specific Doppler shifts on the other hand. The discrete dust structures were obtained by application of the Fellwalker algorithm to a recent 3D extinction density map. For each extinction cloud, a technique using a narrow sliding spectral window moved along the contour-bounded CO spectrum and geometrical criteria was used to select the most likely velocity interval. We applied the new contour-based technique to the 3D extinction density distribution within the volume encompassing the Taurus, Auriga, Perseus and California molecular complexes. From the 45 clouds issued from the decomposition, 42 were assigned a velocity. The remaining structures correspond to very weak CO emission or extinction. We used the non-automated assignments of radial velocities to clouds of the same region presented in paper I and based on KI absorption spectra as a validation test. The new fully automated determinations were found in good agreement with these previous measurements. Our results show that an automated search based on cloud contour morphology can be efficient and that this novel technique may be extended to wider regions of the Milky Way and at larger distance. We discuss its limitations and potential improvements after combination with other techniques.

Forty years ago it was proposed that gas phase organic chemistry in the interstellar medium was initiated by the methyl cation CH3+, but hitherto it has not been observed outside the Solar System. Alternative routes involving processes on grain surfaces have been invoked. Here we report JWST observations of CH3+ in a protoplanetary disk in the Orion star forming region. We find that gas-phase organic chemistry is activated by UV irradiation.

We used 3D maps of 862nm DIB equivalent width (EW) and extinction, DIB catalogues, and measured parameters of dust extinction law and dust emission to study relationships between DIB and extinction level, total-to-selective extinction ratio Rv, dust emission spectral index beta. We revisited the link between several DIBs and the 220nm absorption bump. The ratio, DIBn862, between the 862nm DIB carrier density and the extinction density is increasing in low density clouds, confirming with local values the line-of-sight data. A fitted power law ranks this DIB in the high increase range among the 20 bands measured toward SDSS targets. Using map-integrated 862nm DIB EWs and extinctions along the paths to APOGEE targets with proxies R'v for Rv, we found that DIBn862 increases with R'v for low to moderate extinctions (Av<2-3 mag). Based on stars outside the thin disk, DIBn862 is found to be globally anti-correlated with the Planck opacity spectral index beta. In the light of a recent result on the variability of the carbon/silicate ratio in dust grains as a source of the Rv-beta anti-correlation, it suggests that DIBn862 increases with this ratio, in agreement with the carbonaceous nature of carriers and recent evidences for a spatial correlation between DIBn862 and carbon-rich ejecta of AGBs. At higher Av, both trends disappear. We found that two factors explain the absence of clear results on the link between the UV absorption bump height and DIBs: the correlation disappears when we move from sigma- to zeta-type DIBs and/or from single-cloud lines of sight to paths crossing multiple clouds distant from each other. We show examples of simple models of the bump height based on DIBs. We found an anti-correlation between DIBn and the bump width, similarly disappearing from sigma- to zeta-type DIBs. This suggests that a fraction of the bump is generated outside the dense molecular clouds.