Showing 1–50 of 128 results for author: Roche, P

Polycyclic Aromatic Hydrocarbon Emission in Galaxies as seen with JWST

Authors: Dimitra Rigopoulou, Fergus R. Donnan, Ismael Garc��a-Bernete, Miguel Pereira-Santaella, Almudena Alonso-Herrero, Ric Davies, Leslie K. Hunt, Patrick F. Roche, Taro Shimizu

Abstract: We present a systematic study of mid-infrared spectra of galaxies including star-forming galaxies and Active Galactic Nuclei observed with JWST MIRI-MRS and NIRSpec-IFU. We focus on the relative variations of the 3.3, 6.2, 7.7, 11.3, 12.7 and 17 micron Polycyclic Aromatic Hydrocarbon (PAH) features within spatially resolved regions of galaxies including NGC 3256, NGC 7469, VV 114, II Zw96 and NGC… ▽ More We present a systematic study of mid-infrared spectra of galaxies including star-forming galaxies and Active Galactic Nuclei observed with JWST MIRI-MRS and NIRSpec-IFU. We focus on the relative variations of the 3.3, 6.2, 7.7, 11.3, 12.7 and 17 micron Polycyclic Aromatic Hydrocarbon (PAH) features within spatially resolved regions of galaxies including NGC 3256, NGC 7469, VV 114, II Zw96 and NGC 5728. Using theoretical PAH models and extending our earlier work, we introduce a new PAH diagnostic involving the 17 micron PAH feature. To determine the drivers of PAH band variations in galaxies, we compare observed PAH spectral bands to predictions from theoretical PAH models. We consider extinction, dehydrogenation and PAH size and charge as possible drivers of PAH band variations. We find a surprising uniformity in PAH size distribution among the spatially resolved regions of the galaxies studied here, with no evidence for preferential destruction of the smallest grains, contrary to earlier findings. Neither extinction nor dehydrogenation play a crucial role in setting the observed PAH bands. Instead, we find that PAH charge plays a significant role in PAH inter-band variations. We find a tight relation between PAH charge and the intensity of the radiation field as traced by the [NeIII]/[NeII] maps. In agreement with recent JWST results, we find a predominance of neutral PAH molecules in the nuclei of Active Galaxies and their outflows. Ionised PAHs are the dominant population in star-forming galaxies. We discuss the implications of our findings for the use of PAHs as ISM tracers in high redshift galaxies. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 16 pages, 13 figures, accepted for publication by MNRAS

Machine Learning for Exoplanet Detection in High-Contrast Spectroscopy: Revealing Exoplanets by Leveraging Hidden Molecular Signatures in Cross-Correlated Spectra with Convolutional Neural Networks

Authors: Emily O. Garvin, Markus J. Bonse, Jean Hayoz, Gabriele Cugno, Jonas Spiller, Polychronis A. Patapis, Dominique Petit Dit de la Roche, Rakesh Nath-Ranga, Olivier Absil, Nicolai F. Meinshausen, Sascha P. Quanz

Abstract: The new generation of observatories and instruments (VLT/ERIS, JWST, ELT) motivate the development of robust methods to detect and characterise faint and close-in exoplanets. Molecular mapping and cross-correlation for spectroscopy use molecular templates to isolate a planet's spectrum from its host star. However, reliance on signal-to-noise ratio (S/N) metrics can lead to missed discoveries, due… ▽ More The new generation of observatories and instruments (VLT/ERIS, JWST, ELT) motivate the development of robust methods to detect and characterise faint and close-in exoplanets. Molecular mapping and cross-correlation for spectroscopy use molecular templates to isolate a planet's spectrum from its host star. However, reliance on signal-to-noise ratio (S/N) metrics can lead to missed discoveries, due to strong assumptions of Gaussian independent and identically distributed noise. We introduce machine learning for cross-correlation spectroscopy (MLCCS); the method aims to leverage weak assumptions on exoplanet characterisation, such as the presence of specific molecules in atmospheres, to improve detection sensitivity for exoplanets. MLCCS methods, including a perceptron and unidimensional convolutional neural networks, operate in the cross-correlated spectral dimension, in which patterns from molecules can be identified. We test on mock datasets of synthetic planets inserted into real noise from SINFONI at K-band. The results from MLCCS show outstanding improvements. The outcome on a grid of faint synthetic gas giants shows that for a false discovery rate up to 5%, a perceptron can detect about 26 times the amount of planets compared to an S/N metric. This factor increases up to 77 times with convolutional neural networks, with a statistical sensitivity shift from 0.7% to 55.5%. In addition, MLCCS methods show a drastic improvement in detection confidence and conspicuity on imaging spectroscopy. Once trained, MLCCS methods offer sensitive and rapid detection of exoplanets and their molecular species in the spectral dimension. They handle systematic noise and challenging seeing conditions, can adapt to many spectroscopic instruments and modes, and are versatile regarding atmospheric characteristics, which can enable identification of various planets in archival and future data. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 27 pages, 24 figures. Submitted for publication in A&A January 2, 2024. After first iteration with the referee, resubmitted May 17, 2024

Peeling Back the Layers of Extinction of Dusty Galaxies in the Era of JWST: Modelling Joint NIRSpec + MIRI Spectra at rest-frame 1.5-28 $μ$m

Authors: F. R. Donnan, I. García-Bernete, D. Rigopoulou, M. Pereira-Santaella, P. F. Roche, A. Alonso-Herrero

Abstract: We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $μ$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical na… ▽ More We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $μ$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical nature of the dust in these environments. We show that differential extinction is necessary to reproduce the spectra of 4 highly obscured Luminous Infrared Galaxies observed with NIRSpec IFU and MIRI MRS, where simple screen or uniformly mixed dust distributions fail to fit the data. We additionally compare the extinction of HII regions in these galaxies via hydrogen recombination lines, the extinction of molecular gas via the H$_2$ lines, Polycyclic Aromatic Hydrocarbons via the 12.7/11.3 PAH ratio and the stellar continuum. We find that the molecular gas is deeply buried with the HII regions in star-forming regions, with a similar extinction to the hottest dust components. However we find the cooler dust to be less obscured, at a similar extinction to the stellar continuum and PAHs. The nuclei show a complex dust distribution with VV114 NE, NGC 3256 S, IIZw96 SW showing a deeply buried continuum source relative to the molecular gas/HII regions. Additionally, NGC 3256 S, NGC 7469 and VV114 SW show an isolated hot dust component, indicative of AGN heating, where NGC 3256 S and NGC 7469 are previously known AGN. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 19 pages, 13 figures. Accepted for publication in MNRAS

Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

Authors: Taylor J. Bell, Nicolas Crouzet, Patricio E. Cubillos, Laura Kreidberg, Anjali A. A. Piette, Michael T. Roman, Joanna K. Barstow, Jasmina Blecic, Ludmila Carone, Louis-Philippe Coulombe, Elsa Ducrot, Mark Hammond, João M. Mendonça, Julianne I. Moses, Vivien Parmentier, Kevin B. Stevenson, Lucas Teinturier, Michael Zhang, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Benjamin Charnay, Katy L. Chubb, Brice-Olivier Demory, Peter Gao , et al. (58 additional authors not shown)

Abstract: Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5… ▽ More Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions). △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 61 pages, 13 figures, 4 tables. This preprint has been submitted to and accepted in principle for publication in Nature Astronomy without significant changes

SAGE: A tool to constrain impacts of stellar activity on transmission spectroscopy

Authors: Hritam Chakraborty, Monika Lendl, Babatunde Akinsanmi, Dominique J. M. Petit dit de la Roche, Adrien Deline

Abstract: Transmission spectroscopy is a proven technique to study a transiting exoplanet's atmosphere. However, stellar surface inhomogeneities, spots and faculae, alter the observed transmission spectra: the stellar contamination effect. The variable nature of the stellar activity also makes it difficult to stitch together multi-epoch observations and evaluate any potential variability in the exoplanet's… ▽ More Transmission spectroscopy is a proven technique to study a transiting exoplanet's atmosphere. However, stellar surface inhomogeneities, spots and faculae, alter the observed transmission spectra: the stellar contamination effect. The variable nature of the stellar activity also makes it difficult to stitch together multi-epoch observations and evaluate any potential variability in the exoplanet's atmosphere. This paper introduces SAGE, a tool to correct for the time-dependent impact of stellar activity on transmission spectra. It uses a pixelation approach to model the stellar surface with spots and faculae, while fully accounting for limb-darkening and rotational line-broadening. The current version is designed for low to medium-resolution spectra. We used SAGE to evaluate stellar contamination for F to M-type hosts, testing various spot sizes and locations, and quantify the impact of limb-darkening. We find that limb-darkening enhances the importance of the spot location on the stellar disk, with spots close to the disk center impacting the transmission spectra more strongly than spots near the limb. Moreover, due to the chromaticity of limb darkening, the shape of the contamination spectrum is also altered. Additionally, SAGE can be used to retrieve the properties and distribution of active regions on the stellar surface from photometric monitoring. We demonstrate this for WASP-69 using TESS data, finding that two spots at mid-latitudes and a combined coverage fraction of $\sim$1% are favoured. SAGE allows us to connect the photometric variability to the stellar contamination of transmission spectra, enhancing our ability to jointly interpret transmission spectra obtained at different epochs. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 685, A173 (2024)

The Galaxy Activity, Torus, and Outflow Survey (GATOS). III: Revealing the inner icy structure in local AGN

Authors: I. García-Bernete, A. Alonso-Herrero, D. Rigopoulou, M. Pereira-Santaella, T. Shimizu, R. Davies, F. R. Donnan, P. F. Roche, O. González-Martín, C. Ramos Almeida, E. Bellocchi, P. Boorman, F. Combes, A. Efstathiou, D. Esparza-Arredondo, S. García-Burillo, E. González-Alfonso, E. K. S. Hicks, S. Hönig, A. Labiano, N. A. Levenson, E. López-Rodríguez, C. Ricci, C. Packham, D. Rouan , et al. (2 additional authors not shown)

Abstract: We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1… ▽ More We use JWST/MIRI MRS spectroscopy of a sample of six local obscured type 1.9/2 active galactic nuclei (AGN) to compare their nuclear mid-IR absorption bands with the level of nuclear obscuration traced by X-rays. This study is the first to use sub-arcsecond angular resolution data of local obscured AGN to investigate the nuclear mid-IR absorption bands with a wide wavelength coverage (4.9-28.1 $μ$m). All the nuclei show the 9.7 $μ$m silicate band in absorption. We compare the strength of the 9.7 and 18 $μ$m silicate features with torus model predictions. The observed silicate features are generally well explained by clumpy and smooth torus models. We report the detection of the 6 $μ$m dirty water ice band (i.e., a mix of water and other molecules such as CO and CO$_2$) at sub-arcsecond scales ($\sim$0.26 arcsec at 6 $μ$m; inner $\sim$50 pc) in a sample of local AGN with different levels of nuclear obscuration in the range log N$_{\rm H}^{\rm X-Ray}$(cm$^{-2}$)$\sim22-25$. We find a good correlation between the 6 $μ$m water ice optical depths and N$_{\rm H}^{\rm X-Ray}$. This result indicates that the water ice absorption might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The weak water ice absorption in less obscured AGN (log N$_H^{X-ray}$ (cm$^{-2}$)$\lesssim$23.0 cm$^{-2}$) might be related to the hotter dust temperature ($>$T$_{sub}^{H_2O}\sim$110 K) expected to be reached in the outer layers of the torus due to their more inhomogeneous medium. Our results suggest it might be necessary to include the molecular content, such as, H$_2$O, aliphatic hydrocarbons (CH-) and more complex PAH molecules in torus models to better constrain key parameters such as the torus covering factor (i.e. nuclear obscuration). △ Less

Submitted 4 December, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A. 13 pages, 12 Figures

Why every observatory needs a disco ball

Authors: Robert J. Cumming, Alexander G. M. Pietrow, Livia Pietrow, Maria Cavallius, Dominique Petit dit de la Roche, Casper Pietrow, Ilane Schroetter, Moa Skan

Abstract: Commercial disco balls provide a safe, effective and instructive way of observing the Sun. We explore the optics of solar projections with disco balls, and find that while sunspot observations are challenging, the solar disk and its changes during eclipses are easy and fun to observe. We explore the disco ball's potential for observing the moon and other bright astronomical phenomena. Commercial disco balls provide a safe, effective and instructive way of observing the Sun. We explore the optics of solar projections with disco balls, and find that while sunspot observations are challenging, the solar disk and its changes during eclipses are easy and fun to observe. We explore the disco ball's potential for observing the moon and other bright astronomical phenomena. △ Less

Submitted 27 September, 2023; v1 submitted 25 September, 2023; originally announced September 2023.

Comments: 6 pages, 7 figures. Submitted to Physics Education. Comments welcome

Detection of Carbon Monoxide in the Atmosphere of WASP-39b Applying Standard Cross-Correlation Techniques to JWST NIRSpec G395H Data

Authors: Emma Esparza-Borges, Mercedes López-Morales, Jéa I. Adams Redai, Enric Pallé, James Kirk, Núria Casasayas-Barris, Natasha E. Batalha, Benjamin V. Rackham, Jacob L. Bean, S. L. Casewell, Leen Decin, Leonardo A. Dos Santos, Antonio García Muñoz, Joseph Harrington, Kevin Heng, Renyu Hu, Luigi Mancini, Karan Molaverdikhani, Giuseppe Morello, Nikolay K. Nikolov, Matthew C. Nixon, Seth Redfield, Kevin B. Stevenson, Hannah R. Wakeford, Munazza K. Alam , et al. (8 additional authors not shown)

Abstract: Carbon monoxide was recently reported in the atmosphere of the hot Jupiter WASP-39b using the NIRSpec PRISM transit observation of this planet, collected as part of the JWST Transiting Exoplanet Community Early Release Science (JTEC ERS) Program. This detection, however, could not be confidently confirmed in the initial analysis of the higher resolution observations with NIRSpec G395H disperser. H… ▽ More Carbon monoxide was recently reported in the atmosphere of the hot Jupiter WASP-39b using the NIRSpec PRISM transit observation of this planet, collected as part of the JWST Transiting Exoplanet Community Early Release Science (JTEC ERS) Program. This detection, however, could not be confidently confirmed in the initial analysis of the higher resolution observations with NIRSpec G395H disperser. Here we confirm the detection of CO in the atmosphere of WASP-39b using the NIRSpec G395H data and cross-correlation techniques. We do this by searching for the CO signal in the unbinned transmission spectrum of the planet between 4.6 and 5.0 $μ$m, where the contribution of CO is expected to be higher than that of other anticipated molecules in the planet's atmosphere. Our search results in a detection of CO with a cross-correlation function (CCF) significance of $6.6 σ$ when using a template with only ${\rm ^{12}C^{16}O}$ lines. The CCF significance of the CO signal increases to $7.5 σ$ when including in the template lines from additional CO isotopologues, with the largest contribution being from ${\rm ^{13}C^{16}O}$. Our results highlight how cross-correlation techniques can be a powerful tool for unveiling the chemical composition of exoplanetary atmospheres from medium-resolution transmission spectra, including the detection of isotopologues. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Comments: Accepted for publication in The Astrophysical Journal Letters

AGN feedback in action in the molecular gas ring of the Seyfert galaxy NGC7172

Authors: A. Alonso-Herrero, S. Garcia-Burillo, M. Pereira-Santaella, T. Shimizu, F. Combes, E. K. S. Hicks, R. Davies, C. Ramos Almeida, I. Garcia-Bernete, S. F. Hoenig, N. A. Levenson, C. Packham, E. Bellocchi, L. K. Hunt, M. Imanishi, C. Ricci, P. Roche

Abstract: We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring w… ▽ More We present new ALMA observations of the CO(3-2) transition and 854micron continuum at 0.06-0.3" resolution, together with new VLT/SINFONI observations of NGC7172. This is a luminous (bolometric luminosity of ~10^44 erg/s) Seyfert galaxy that belongs to the Galaxy Activity, Torus, and Outflow Survey (GATOS). The CO(3-2) observations reveal the presence of a highly inclined cold molecular gas ring with an approximate radius of 3-4"~540-720 pc, which is likely associated with an inner Lindblad resonance of a putative stellar bar. There are noncircular motions in the VLT/SINFONI [SiVI]1.96micron and H2 at 2.12micron, and ALMA CO(3-2) velocity fields. After subtracting the stellar velocity field, we detected [SiVI] blueshifted velocities of a few hundred km/s to the south of the AGN. They trace outflowing ionized gas outside the plane of the galaxy and out to projected distances of ~200 pc. The CO(3-2) position-velocity diagram along the kinematic minor axis displays noncircular motions with observed velocities of up to ~150 km/s. Assuming that these are taking place in the disk of the galaxy, the observed velocity signs imply that the molecular gas ring is not only rotating but also outflowing. We derived an integrated cold molecular gas mass outflow rate of ~40 Msun/yr for the ring. Using the 854micron map, we resolved a 32 pc radius torus with a gas mass of 8x10^5 Msun. These torus properties are similar to other Seyfert galaxies in the GATOS sample. We measured a decreased cold molecular gas concentration in the nuclear-torus region relative to the circumnuclear region when compared to other less luminous Seyfert galaxies. We conclude that the effects of AGN feedback in NGC7172, which are likely caused by the AGN wind and/or the moderate luminosity radio jet, are seen as a large-scale outflowing molecular gas ring and accompanying redistribution of molecular gas in the nuclear regions. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted for publication to A&A

The role of grain size in AGN torus dust models

Authors: Omaira González-Martín, Cristina Ramos Almeida, Jacopo Fritz, Almudena Alonso-Herrero, Sebastian F. Hönig, Patrick F. Roche, Donaji Esparza-Arredondo, Ismael García-Bernete, Santiago García-Burillo, Natalia Osorio-Clavijo, Ulises Reyes-Amador, Marko Stalevski, César Victoria-Ceballos

Abstract: Fits the infrared spectra from the nuclear regions of AGN can place constraints on the dust properties, distribution, and geometry by comparison with models. However, none of the currently available models fully describe the observations of AGN currently available. Among the aspects least explored, here we focus on the role of dust grain size. We offer the community a new spectral energy distribut… ▽ More Fits the infrared spectra from the nuclear regions of AGN can place constraints on the dust properties, distribution, and geometry by comparison with models. However, none of the currently available models fully describe the observations of AGN currently available. Among the aspects least explored, here we focus on the role of dust grain size. We offer the community a new spectral energy distribution (SED) library, hereinafter [GoMar23] model, which is based on the two-phase torus model developed before with the inclusion of the grain size as a model parameter, parameterized by the maximum grain size Psize or equivalently the mass-weighted average grain size < P >. We created 691,200 SEDs using the SKIRT code, where the maximum grain size can vary within the range Psize = 0.01 - 10.0um ( < P >= 0.007 - 3.41um). We fit this new and several existing libraries to a sample of 68 nearby and luminous AGNs with Spitzer/IRS spectra dominated by AGN-heated dust. We find that the [GoMar23] model can adequately reproduce up to 85-88% of the spectra. The dust grain size parameter significantly improves the final fit in up to 90% of these spectra. Statistical tests indicate that the grain size is the third most important parameter in the fitting procedure (after the size and half opening angle of the torus). The requirement of a foreground extinction by our model is lower compared to purely clumpy models. We find that 41% of our sample requires that the maximum dust grain size is as large as Psize =10um (< P >= 3.41um). Nonetheless, we also remark that disk+wind and clumpy torus models are still required to reproduce the spectra of a non-negligible fraction of objects, suggesting the need for several dust geometries to explain the infrared continuum of AGN. This work provides tentative evidence for dust grain growth in the proximity of the AGN. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 26 pages, 14 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 676, A73 (2023)

Spatially resolving polycyclic aromatic hydrocarbons in Herbig Ae disks with VISIR-NEAR at the VLT

Authors: Gideon Yoffe, Roy van Boekel, Aigen Li, L. B. F. M Waters, Koen Maaskant, Ralf Siebenmorgen, Mario van den Ancker, D. J. M Petit dit de la Roche, Bruno Lopez, Alexis Matter, Jozsef Varga, M. R Hogerheijde, Gerd Weigelt, R. D Oudmaijer, Eric Pantin, M. R Meyer, Jean-Charles Augereau, Thomas Henning

Abstract: We use the long-slit spectroscopy mode of the VISIR-NEAR experiment to perform diffraction-limited observations of eight nearby Herbig Ae protoplanetary disks. We extract spectra for various locations along the slit with a spectral resolution of R = 300 and perform a compositional fit at each spatial location using spectral templates of silicates and the four PAH bands. This yields the intensity v… ▽ More We use the long-slit spectroscopy mode of the VISIR-NEAR experiment to perform diffraction-limited observations of eight nearby Herbig Ae protoplanetary disks. We extract spectra for various locations along the slit with a spectral resolution of R = 300 and perform a compositional fit at each spatial location using spectral templates of silicates and the four PAH bands. This yields the intensity vs. location profiles of each species. Results. We could obtain spatially-resolved intensity profiles of the PAH emission features in the N-band for five objects (AB Aurigae, HD 97048, HD 100546, HD 163296, and HD 169142). We observe two kinds of PAH emission geometry in our sample: centrally-peaked (HD 97048) and ring-like (AB Aurigae, HD 100546, HD 163296, and potentially HD 169142). Comparing the spatial PAH emission profiles with near-infrared scattered light images, we find a strong correlation in the disk sub-structure but a difference in radial intensity decay rate. The PAH emission shows a less steep decline with distance from the star. Finally, we find a correlation between the presence of (sub-) micron-sized silicate grains leading to the depletion of PAH emission within the inner regions of the disks. In this work, we find the following: (1) PAH emission traces the extent of Herbig Ae disks to a considerable radial distance. (2) The correlation between silicate emission within the inner regions of disks and the depletion of PAH emission can result from dust-mixing and PAH coagulation mechanisms and competition over UV photons. (3) For all objects in our sample, PAHs undergo stochastic heating across the entire spatial extent of the disk and are not saturated. (4) The difference in radial intensity decay rates between the PAHs and scattered-light profiles may be attributed to shadowing and dust-settling effects, which affect the scattering grains more than the PAHs. △ Less

Submitted 27 March, 2023; v1 submitted 12 March, 2023; originally announced March 2023.

Journal ref: A&A 674, A57 (2023)

Center-to-limb variation of spectral lines and continua observed with SST/CRISP and SST/CHROMIS

Authors: A. G. M. Pietrow, D. Kiselman, O. Andriienko, D. J. M. Petit dit de la Roche, C. J. Díaz Baso, F. Calvo

Abstract: Observations of center-to-limb variations (CLV) of spectral lines and continua provide a good test for the accuracy of models of solar and stellar atmospheric structure and spectral-line formation. They are also widely used to constrain elemental abundances and are becoming increasingly important in atmospheric studies of exoplanets. However, only a few such data sets exist for chromospheric lines… ▽ More Observations of center-to-limb variations (CLV) of spectral lines and continua provide a good test for the accuracy of models of solar and stellar atmospheric structure and spectral-line formation. They are also widely used to constrain elemental abundances and are becoming increasingly important in atmospheric studies of exoplanets. However, only a few such data sets exist for chromospheric lines. We aim to create a set of standard profiles by means of mosaics made with the CRISP and CHROMIS instruments of the Swedish 1-m Solar Telescope (SST), as well as to explore the robustness of said profiles obtained using this method. For each spectral line we use a mosaic that ranges from the center to the limb. Each of these mosaics are averaged down to 50 individual spectral profiles, spaced by 0.02 in the $μ$ scale. These profiles are corrected for p-mode oscillations, and their line parameters (equivalent width, line shift, full-width at half-maximum, and line depth) are then compared against literature values where possible. We present a set of 50 average profiles that are spaced equidistantly along the cosine of the heliocentric angle ($μ$) by steps of 0.02 for five continuum points between 4001 and 7772 Å, as well as ten of the most commonly observed spectral lines at the SST (Ca II H & K, H$β$, Mg I 5173 Å, C I 5380 Å, Fe I 6173 Å, Fe I 6301 Å, H$α$, O I 7772 Å, and Ca II 8542 Å). Center-to-limb variation of line profiles and continua are shared in the CDS as machine-readable tables; providing a quantitative constraint on theoretical models that aim to model stellar atmospheres. △ Less

Submitted 19 March, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: Published in A&A. CLV data is available in the CDS

Journal ref: A&A 671, A130 (2023)

Early Release Science of the exoplanet WASP-39b with JWST NIRISS

Authors: Adina D. Feinstein, Michael Radica, Luis Welbanks, Catriona Anne Murray, Kazumasa Ohno, Louis-Philippe Coulombe, Néstor Espinoza, Jacob L. Bean, Johanna K. Teske, Björn Benneke, Michael R. Line, Zafar Rustamkulov, Arianna Saba, Angelos Tsiaras, Joanna K. Barstow, Jonathan J. Fortney, Peter Gao, Heather A. Knutson, Ryan J. MacDonald, Thomas Mikal-Evans, Benjamin V. Rackham, Jake Taylor, Vivien Parmentier, Natalie M. Batalha, Zachory K. Berta-Thompson , et al. (64 additional authors not shown)

Abstract: Transmission spectroscopy provides insight into the atmospheric properties and consequently the formation history, physics, and chemistry of transiting exoplanets. However, obtaining precise inferences of atmospheric properties from transmission spectra requires simultaneously measuring the strength and shape of multiple spectral absorption features from a wide range of chemical species. This has… ▽ More Transmission spectroscopy provides insight into the atmospheric properties and consequently the formation history, physics, and chemistry of transiting exoplanets. However, obtaining precise inferences of atmospheric properties from transmission spectra requires simultaneously measuring the strength and shape of multiple spectral absorption features from a wide range of chemical species. This has been challenging given the precision and wavelength coverage of previous observatories. Here, we present the transmission spectrum of the Saturn-mass exoplanet WASP-39b obtained using the SOSS mode of the NIRISS instrument on the JWST. This spectrum spans $0.6 - 2.8 μ$m in wavelength and reveals multiple water absorption bands, the potassium resonance doublet, as well as signatures of clouds. The precision and broad wavelength coverage of NIRISS-SOSS allows us to break model degeneracies between cloud properties and the atmospheric composition of WASP-39b, favoring a heavy element enhancement ("metallicity") of $\sim 10 - 30 \times$ the solar value, a sub-solar carbon-to-oxygen (C/O) ratio, and a solar-to-super-solar potassium-to-oxygen (K/O) ratio. The observations are best explained by wavelength-dependent, non-gray clouds with inhomogeneous coverage of the planet's terminator. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 48 pages, 12 figures, 2 tables. Under review at Nature

Early Release Science of the exoplanet WASP-39b with JWST NIRCam

Authors: Eva-Maria Ahrer, Kevin B. Stevenson, Megan Mansfield, Sarah E. Moran, Jonathan Brande, Giuseppe Morello, Catriona A. Murray, Nikolay K. Nikolov, Dominique J. M. Petit dit de la Roche, Everett Schlawin, Peter J. Wheatley, Sebastian Zieba, Natasha E. Batalha, Mario Damiano, Jayesh M Goyal, Monika Lendl, Joshua D. Lothringer, Sagnick Mukherjee, Kazumasa Ohno, Natalie M. Batalha, Matthew P. Battley, Jacob L. Bean, Thomas G. Beatty, Björn Benneke, Zachory K. Berta-Thompson , et al. (74 additional authors not shown)

Abstract: Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength covera… ▽ More Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $μ$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $μ$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 35 pages, 13 figures, 3 tables, Nature, accepted

Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H

Authors: Lili Alderson, Hannah R. Wakeford, Munazza K. Alam, Natasha E. Batalha, Joshua D. Lothringer, Jea Adams Redai, Saugata Barat, Jonathan Brande, Mario Damiano, Tansu Daylan, Néstor Espinoza, Laura Flagg, Jayesh M. Goyal, David Grant, Renyu Hu, Julie Inglis, Elspeth K. H. Lee, Thomas Mikal-Evans, Lakeisha Ramos-Rosado, Pierre-Alexis Roy, Nicole L. Wallack, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Zachory K. Berta-Thompson , et al. (67 additional authors not shown)

Abstract: Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the m… ▽ More Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the medium-resolution (R$\sim$600) transmission spectrum of an exoplanet atmosphere between 3-5 $μ$m covering multiple absorption features for the Saturn-mass exoplanet WASP-39b, obtained with JWST NIRSpec G395H. Our observations achieve 1.46x photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO$_2$ (28.5$σ$) and H$_2$O (21.5$σ$), and identify SO$_2$ as the source of absorption at 4.1 $μ$m (4.8$σ$). Best-fit atmospheric models range between 3 and 10x solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO$_2$, underscore the importance of characterising the chemistry in exoplanet atmospheres, and showcase NIRSpec G395H as an excellent mode for time series observations over this critical wavelength range. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 44 pages, 11 figures, 3 tables. Resubmitted after revision to Nature

Early Release Science of the exoplanet WASP-39b with JWST NIRSpec PRISM

Authors: Z. Rustamkulov, D. K. Sing, S. Mukherjee, E. M. May, J. Kirk, E. Schlawin, M. R. Line, C. Piaulet, A. L. Carter, N. E. Batalha, J. M. Goyal, M. López-Morales, J. D. Lothringer, R. J. MacDonald, S. E. Moran, K. B. Stevenson, H. R. Wakeford, N. Espinoza, J. L. Bean, N. M. Batalha, B. Benneke, Z. K. Berta-Thompson, I. J. M. Crossfield, P. Gao, L. Kreidberg , et al. (69 additional authors not shown)

Abstract: Transmission spectroscopy of exoplanets has revealed signatures of water vapor, aerosols, and alkali metals in a few dozen exoplanet atmospheres. However, these previous inferences with the Hubble and Spitzer Space Telescopes were hindered by the observations' relatively narrow wavelength range and spectral resolving power, which precluded the unambiguous identification of other chemical species… ▽ More Transmission spectroscopy of exoplanets has revealed signatures of water vapor, aerosols, and alkali metals in a few dozen exoplanet atmospheres. However, these previous inferences with the Hubble and Spitzer Space Telescopes were hindered by the observations' relatively narrow wavelength range and spectral resolving power, which precluded the unambiguous identification of other chemical species$-$in particular the primary carbon-bearing molecules. Here we report a broad-wavelength 0.5-5.5 $μ$m atmospheric transmission spectrum of WASP-39 b, a 1200 K, roughly Saturn-mass, Jupiter-radius exoplanet, measured with JWST NIRSpec's PRISM mode as part of the JWST Transiting Exoplanet Community Early Release Science Team program. We robustly detect multiple chemical species at high significance, including Na (19$σ$), H$_2$O (33$σ$), CO$_2$ (28$σ$), and CO (7$σ$). The non-detection of CH$_4$, combined with a strong CO$_2$ feature, favours atmospheric models with a super-solar atmospheric metallicity. An unanticipated absorption feature at 4$μ$m is best explained by SO$_2$ (2.7$σ$), which could be a tracer of atmospheric photochemistry. These observations demonstrate JWST's sensitivity to a rich diversity of exoplanet compositions and chemical processes. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 41 pages, 4 main figures, 10 extended data figures, 4 tables. Under review in Nature

A Detailed Look at the Most Obscured Galactic Nuclei in the Mid-Infrared

Authors: F. R. Donnan, D. Rigopoulou, I. García-Bernete, M. Pereira-Santaella, A. Alonso-Herrero, P. F. Roche, S. Aalto, A. Hernán-Caballero, H. W. W. Spoon

Abstract: Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromat… ▽ More Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromatic Hydrocarbon (PAH) features. Methods. We model Spitzer spectra by decomposing the continua into nuclear and star-forming components from which we then measure the nuclear optical depth, $τ_N$, of the $9.8 μ$m silicate absorption feature. We also use Spitzer spectral maps to investigate how PAH EW ratios vary with aperture size for objects hosting CONs. Results. We find that the nuclear optical depth, $τ_N$, strongly correlates with the HCN-vib emission line in the millimetre for CONs with a Pearson correlation coefficient of 0.91. We find the PAH EW ratios technique to be effective at selecting CONs and robust against highly inclined galaxies where strong dust lanes may mimic a CON like spectrum by producing a high $τ_N$. Our analysis of the Spitzer spectral maps showed that the efficacy of the PAH EW ratios to isolate CONs is reduced when there is a strong star-forming component from the host galaxy. In addition, we find that the use of the inferred nuclear optical depth is a reliable method to identify CONs in $36^{+8}_{-7}\%$ of ULIRGs and $17^{+3}_{-3}\%$ of LIRGs, consistent with previous work. Conclusions. We confirm mid-IR spectra to be a powerful diagnostic of CONs where the increased sensitivity of JWST will allow identification of CONs at cosmic noon revealing this extreme but hidden phase of galaxy evolution. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A. 23 pages, 15 figures, 3 tables

Journal ref: A&A 669, A87 (2023)

The Obscured Nucleus and Shocked Environment of VV 114E Revealed by JWST/MIRI Spectroscopy

Authors: F. R. Donnan, I. García-Bernete, D. Rigopoulou, M. Pereira-Santaella, A. Alonso-Herrero, P. F. Roche, A. Hernán-Caballero, H. W. W. Spoon

Abstract: Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28… ▽ More Compact Obscured Nuclei (CONs) potentially hide extreme supermassive black hole (SMBH) growth behind large column densities of gas/dust. We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the JWST Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28 $μ$m to extract Polycyclic Aromatic Hydrocarbon (PAH) emission and the underlying obscured continuum. We find that the NE nucleus (A) is highly obscured where the low PAH equivalent width (EW) ratio, EW(12.7)/EW(11.3), reveals a dust enshrouded continuum source. This is confirmed by decomposing the continuum into nuclear and star-forming where the nuclear component is found to be typical of CONs. The 11.3/6.2 PAH flux ratio is consistent with originating in star-forming regions rather than typical AGN. The second nucleus (B) is much less obscured, with PAH flux ratios also typical of star-forming regions. We do not detect any high ionisation lines such as [Ne V] or [Ne VI] which suggests that if an AGN is present it must be highly obscured. Additionally, we detect a shock front south of the secondary nucleus (B) in the [Fe II] (5.34 $μ$m) line and in warm molecular hydrogen. The 6.2 PAH emission does not spatially coincide with the low-J transitions of H$_2$ but rather appears strong at the shock front which may suggest destruction of the ionised PAHs in the post-shock gas behind the shock front. △ Less

Submitted 15 December, 2022; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS

Torus and polar dust dependence on AGN properties

Authors: I. García-Bernete, O. González-Martín, C. Ramos Almeida, A. Alonso-Herrero, M. Martínez-Paredes, M. J. Ward, P. F. Roche, J. A. Acosta-Pulido, E. López-Rodríguez, D. Rigopoulou, D. Esparza-Arredondo

Abstract: We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and c… ▽ More We present a statistical analysis of the properties of the obscuring material around active galactic nuclei (AGN). This study represents the first of its kind for an ultra-hard X-ray (14-195keV; Swift/BAT) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample) using high angular resolution infrared data and various torus models: smooth, clumpy and two-phase torus models and clumpy disc+wind models. We find that the smooth, clumpy and two-phase torus models (i.e. without including the polar dusty wind component) and disc+wind models provide best fits for a comparable number of galaxies, 8/24 (33.3%) and 9/24 (37.5%), respectively. We find that the best-fitted models depend on the hydrogen column density (NH), which is related to the X-ray (unobscured/obscured) and/or optical (Sy1/Sy2) classification. In particular, smooth, clumpy and two-phase torus models best reproduce the infrared (IR) emission of AGN with relatively high hydrogen column density (median value of log (NH)=23.5+-0.8; i.e. Sy2). However, clumpy disc+wind models provide the best fits to the nuclear IR spectral energy distributions (SEDs) of Sy1/1.8/1.9 (median value of log (NH)=21.0+-1.0), specifically in the near-IR (NIR) range. The success of the disc+wind models in fitting the NIR emission of Sy1 galaxies is due to the combination of adding large graphite grains to the dust composition and self-obscuration effects caused by the wind at intermediate inclinations. In general, we find that the Seyfert galaxies having unfavourable (favourable) conditions, i.e. nuclear hydrogen column density and Eddington ratio, for launching IR dusty polar outflows are best-fitted with smooth, clumpy and two-phase torus (disk+wind) models confirming the predictions from simulations. Therefore, our results indicate that the nature of the inner dusty structure in AGN depend on the intrinsic AGN properties. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: Accepted for publication in A&A, 43 pages, 43 figures, 16 tables

Journal ref: A&A 667, A140 (2022)

Probing computational methodologies in predicting mid-infrared spectra for large polycyclic aromatic hydrocarbons

Authors: Boutheïna Kerkeni, Ismael García-Bernete, Dimitra Rigopoulou, David P. Tew, Patrick F. Roche, David C. Clary

Abstract: We extend the prediction of vibrational spectra to large sized polycyclic aromatic hydrocarbon (PAH) molecules comprising up to \sim 1500 carbon atoms by evaluating the efficiency of several computational chemistry methodologies. We employ classical mechanics methods (Amber and Gaff) with improved atomic point charges, semi-empirical (PM3, and density functional tight binding), and density functio… ▽ More We extend the prediction of vibrational spectra to large sized polycyclic aromatic hydrocarbon (PAH) molecules comprising up to \sim 1500 carbon atoms by evaluating the efficiency of several computational chemistry methodologies. We employ classical mechanics methods (Amber and Gaff) with improved atomic point charges, semi-empirical (PM3, and density functional tight binding), and density functional theory (B3LYP) and conduct global optimizations and frequency calculations in order to investigate the impact of PAH size on the vibrational band positions. We primarily focus on the following mid-infrared emission bands 3.3, 6.2, 7.7, 8.6, 11.3, 12.7, and 17.0 microns. We developed a general Frequency Scaling Function (FSF) to shift the bands and to provide a systematic comparison versus the three methods for each PAH. We first validate this procedure on IR scaled spectra from the NASA Ames PAH Database, and extend it to new large PAHs. We show that when the FSF is applied to the Amber and Gaff IR spectra, an agreement between the normal mode peak positions with those inferred from the B3LYP/4-31G model chemistry is achieved. As calculations become time intensive for large sized molecules Nc > 450, this proposed methodology has advantages. The FSF has enabled extending the investigations to large PAHs where we clearly see the emergence of the 17.0 microns feature, and the weakening of the 3.3 microns one. We finally investigate the trends in the 3.3 microns 17.0 microns PAH band ratio as a function of PAH size and its response following the exposure to fields of varying radiation intensities. △ Less

Submitted 26 September, 2022; originally announced October 2022.

Comments: 19 pages, 12 figures

Journal ref: 513, 2022, 3663

Identification of carbon dioxide in an exoplanet atmosphere

Authors: The JWST Transiting Exoplanet Community Early Release Science Team, Eva-Maria Ahrer, Lili Alderson, Natalie M. Batalha, Natasha E. Batalha, Jacob L. Bean, Thomas G. Beatty, Taylor J. Bell, Björn Benneke, Zachory K. Berta-Thompson, Aarynn L. Carter, Ian J. M. Crossfield, Néstor Espinoza, Adina D. Feinstein, Jonathan J. Fortney, Neale P. Gibson, Jayesh M. Goyal, Eliza M. -R. Kempton, James Kirk, Laura Kreidberg, Mercedes López-Morales, Michael R. Line, Joshua D. Lothringer, Sarah E. Moran, Sagnick Mukherjee , et al. (107 additional authors not shown)

Abstract: Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres… ▽ More Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 μm in wavelength and show a prominent CO2 absorption feature at 4.3 μm (26σ significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 μm that is not reproduced by these models. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

A high angular resolution view of the PAH emission in Seyfert galaxies using JWST/MRS data

Authors: I. García-Bernete, D. Rigopoulou, A. Alonso-Herrero, F. R. Donnan, P. F. Roche, M. Pereira-Santaella, A. Labiano, L. Peralta de Arriba, T. Izumi, C. Ramos Almeida, T. Shimizu, S. Hönig, S. García-Burillo, D. J. Rosario, M. J. Ward, E. Bellocchi, E. K. S. Hicks, L. Fuller, C. Packham

Abstract: Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work, we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming regions. This study represents the first of its kind using sub-arcsecond angular resolution data of local luminous Seyferts (Lb… ▽ More Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules that are ubiquitous in a variety of astrophysical objects and environments. In this work, we use JWST/MIRI MRS spectroscopy of three Seyferts to compare their nuclear PAH emission with that of star-forming regions. This study represents the first of its kind using sub-arcsecond angular resolution data of local luminous Seyferts (Lbol>10^44.46 erg/s) on a wide wavelength coverage (4.9-28.1 micron). We present an analysis of their nuclear PAH properties by comparing the observed ratios with PAH diagnostic model grids, derived from theoretical spectra. Our results show that a suite of PAH features is present in the innermost parts (~0.45 arcsec at 12 micron; in the inner ~142-245 pc) of luminous Seyfert galaxies. We find that the nuclear regions of AGN lie at different positions of the PAH diagnostic diagrams, whereas the SF regions are concentrated around the average values of SF galaxies. In particular, we find that the nuclear PAH emission mainly originates in neutral PAHs. In contrast, PAH emission originating in the SF regions favours ionised PAH grains. The observed PAH ratios in the nuclear region of AGN-dominated galaxy NGC 6552 indicate the presence of larger-sized PAH molecules compared with those of the SF regions. Therefore, our results provide evidence that the AGN have a significant impact on the ionization state (and probably the size) of the PAH grains on scales of ~142-245 pc. △ Less

Submitted 28 September, 2022; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: 12 pages, 10 figures. Small changes to match the accepted for publication version in A&A (official date of acceptance: 20/09/2022)

Journal ref: A&A 666, L5 (2022)

A technique to select the most obscured galaxy nuclei

Authors: I. García-Bernete, D. Rigopoulou, S. Aalto, H. W. W. Spoon, A. Hernán-Caballero, A. Efstathiou, P. F. Roche, S. König

Abstract: Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS… ▽ More Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS spectra of a representative sample of local ULIRGs (z<0.27), we examine their PAH and underlying continuum emission ratios. For deeply embedded sources, we find that the 9.7 micron silicate absorption band has a particularly pronounced effect on the 11.3 micron PAH feature. The low flux level in the nuclear silicate absorption band enhances the 11.3 micron PAH feature contrast (high PAH equivalent width) compared to that of the other PAH features. The technique has been extended to include the use of the continuum ratios. However, the latter are affected both by the extinction coming from the host galaxy as well as the nuclear region, whereas the foreground extinction is cancelled out when using the PAH equivalent width ratios. We apply our method to the HERUS and GOALS samples and classify as CON candidates 14 ULIRGs and 10 LIRGs, corresponding to 30% of ULIRGs and 7% of LIRGs from these samples. We find that the observed continuum ratios of CON-dominated sources can be explained by assuming torus models with a tapered disk geometry and a smooth dust distribution. This suggests that the nuclear dusty structure of CONs has an extremely high dust coverage. We also demonstrate that the use of mid-IR color-color diagrams is an effective way to select CON-dominated sources at different redshifts. In particular, the combination of filters of the JWST/MIRI will enable the selection of CONs out to z~1.5. This will allow extending the selection of CONs to high redshifts where U/LIRGs are more numerous. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 13 pages, 11 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 663, A46 (2022)

Physical Characterization of Main-Belt Comet (248370) 2005 QN173

Authors: Henry H. Hsieh, Colin O. Chandler, Larry Denneau, Alan Fitzsimmons, Nicolas Erasmus, Michael S. P. Kelley, Matthew M. Knight, Tim A. Lister, Jana Pittichova, Scott S. Sheppard, Audrey Thirouin, Chadwick A. Trujillo, Helen Usher, Edward Gomez, Joey Chatelain, Sarah Greenstreet, Tony Angel, Richard Miles, Paul Roche, Ben Wooding

Abstract: We report results from new and archival observations of the newly discovered active asteroid (248370) 2005 QN_137, which has been determined to be a likely main-belt comet based on a subsequent discovery that it is recurrently active near perihelion. From archival data analysis, we estimate g'-, r'-, i'-, and z'-band absolute magnitudes for the nucleus of H_g=16.62+/-0.13, H_r=16.12+/-0.10, H_i=16… ▽ More We report results from new and archival observations of the newly discovered active asteroid (248370) 2005 QN_137, which has been determined to be a likely main-belt comet based on a subsequent discovery that it is recurrently active near perihelion. From archival data analysis, we estimate g'-, r'-, i'-, and z'-band absolute magnitudes for the nucleus of H_g=16.62+/-0.13, H_r=16.12+/-0.10, H_i=16.05+/-0.11, and H_z=15.93+/-0.08, corresponding to nucleus colors of g'-r'=0.50+/-0.16, r'-i'=0.07+/-0.15, and i'-z'=0.12+/-0.14, an equivalent V-band absolute magnitude of H_V=16.32+/-0.08, and a nucleus radius of r_n=1.6+/-0.2 km (using a V-band albedo of p_V=0.054+/-0.012). Meanwhile, we find mean near-nucleus coma colors when 248370 was active of g'-r'=0.47+/-0.03, r'-i'=0.10+/-0.04, and i'-z'=0.05+/-0.05, and similar mean dust tail colors, suggesting that no significant gas coma is present. We find approximate ratios between the scattering cross-sections of near-nucleus dust (within 5000 km of the nucleus) and the nucleus of A_d/A_n=0.7+/-0.3 on 2016 July 22, and 1.8<A_d/A_n<2.9 in 2021 July and August. During the 2021 observation period, the coma declined in intrinsic brightness by ~0.35 mag (or ~25%) in 37 days, while the surface brightness of the dust tail remained effectively constant over the same period. Constraints derived from the sunward extent of the coma suggest that terminal velocities of ejected dust grains are extremely slow (~1 m/s for 1 micron particles), indicating that the observed dust emission may have been aided by rapid rotation of the nucleus lowering the effective escape velocity. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 12 pages, 2 figures, accepted for publication in ApJ Letters (submitted version posted here; please contact first author for final accepted version)

The Galaxy Activity, Torus and Outflow Survey (GATOS): II. Torus and polar dust emission in nearby Seyfert galaxies

Authors: A. Alonso-Herrero, S. García-Burillo, S. F. Hoenig, I. García-Bernete, C. Ramos Almeida, O. González-Martín, E. López-Rodríguez, P. G. Boorman, A. J. Bunker, L. Burtscher, F. Combes, R. Davies, T. Díaz-Santos, P. Gandhi, B. García-Lorenzo, E. K. S. Hicks, L. K. Hunt, K. Ichikawa, M. Imanishi, T. Izumi, A. Labiano, N. A. Levenson, C. Packham, M. Pereira-Santaella, C. Ricci , et al. (8 additional authors not shown)

Abstract: We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densitie… ▽ More We compare mid-IR and ALMA far-IR images of 12 nearby Seyferts selected from GATOS. The mid-IR unresolved emission contributes more than 60% of the nuclear emission in most galaxies. By contrast, the ALMA 870micron continuum emission is mostly resolved and typically along the torus equatorial direction (Paper I, Garcia-Burillo et al. 2021). The Eddington ratios and nuclear hydrogen column densities NH of half the sample are favorable to launching polar and/or equatorial dusty winds, according to simulations. Six show mid-IR extended emission in the polar direction as traced by the NLR and perpendicular to the ALMA emission. In a few, the nuclear NH might be too high to uplift large quantities of dusty material along the polar direction. Five galaxies have low NH and/or Eddington ratios and thus polar dusty winds are not likely. We generate new CAT3D-WIND disk-wind model images. At low wind-to-disk cloud ratios the far-IR model images have disk- and ring-like morphologies. The X-shape associated with dusty winds is seen better in the far-IR at intermediate-high inclinations for the extended-wind configurations. In most models, the mid-IR emission comes from the inner part of the disk/cone. Extended bi-conical and one-sided polar mid-IR emission is seen in extended-wind configurations and high wind-to-disk cloud ratios. When convolved to our resolution, the model images reproduce qualitative aspects of the observed morphologies. Low-intermediate wind-to-disk ratios are required to account for the large fractions of unresolved mid-IR emission. This work and Paper I provide observational support for the torus+wind scenario. The wind component is more relevant at high Eddington ratios and/or AGN luminosities, and polar dust emission is predicted at NH of up to $10^{24}$cm$^{-2}$. The torus/disk component, on the other hand, prevails at low luminosities and/or Eddington ratios. (Abridged) △ Less

Submitted 2 July, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 27 pages, 27 figures. Accepted for publication in A&A

Constraints on the nearby exoplanet $ε$ Ind Ab from deep near/mid-infrared imaging limits

Authors: Gayathri Viswanath, Markus Janson, Carl-Henrik Dahlqvist, Dominique Petit dit de la Roche, Matthias Samland, Julien Girard, Prashant Pathak, Markus Kasper, Fabo Feng, Michael Meyer, Anna Boehle, Sascha P. Quanz, Hugh R. A. Jones, Olivier Absil, Wolfgang Brandner, Anne-Lise Maire, Ralf Siebenmorgen, Michael Sterzik, Eric Pantin

Abstract: The past decade has seen increasing efforts in detecting and characterising exoplanets by high contrast imaging in the near/mid-infrared, which is the optimal wavelength domain for studying old, cold planets. In this work, we present deep AO imaging observations of the nearby Sun-like star $ε$ Ind A with NaCo ($L^{\prime}$) and NEAR (10-12.5 microns) instruments at VLT, in an attempt to directly d… ▽ More The past decade has seen increasing efforts in detecting and characterising exoplanets by high contrast imaging in the near/mid-infrared, which is the optimal wavelength domain for studying old, cold planets. In this work, we present deep AO imaging observations of the nearby Sun-like star $ε$ Ind A with NaCo ($L^{\prime}$) and NEAR (10-12.5 microns) instruments at VLT, in an attempt to directly detect its planetary companion whose presence has been indicated from radial velocity (RV) and astrometric trends. We derive brightness limits from the non-detection of the companion with both instruments, and interpret the corresponding sensitivity in mass based on both cloudy and cloud-free atmospheric and evolutionary models. For an assumed age of 5 Gyr for the system, we get detectable mass limits as low as 4.4 $M_{\rm J}$ in NaCo $L^{\prime}$ and 8.2 $M_{\rm J}$ in NEAR bands at 1.5$\arcsec$ from the central star. If the age assumed is 1 Gyr, we reach even lower mass limits of 1.7 $M_{\rm J}$ in NaCo $L^{\prime}$ and 3.5 $M_{\rm J}$ in NEAR bands, at the same separation. However, based on the dynamical mass estimate (3.25 $M_{\rm J}$) and ephemerides from astrometry and RV, we find that the non-detection of the planet in these observations puts a constraint of 2 Gyr on the lower age limit of the system. NaCo offers the highest sensitivity to the planetary companion in these observations, but the combination with the NEAR wavelength range adds a considerable degree of robustness against uncertainties in the atmospheric models. This underlines the benefits of including a broad set of wavelengths for detection and characterisation of exoplanets in direct imaging studies. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 13 pages, 9 figures, 2 tables Accepted for publication in section 10. Planets and planetary systems of Astronomy and Astrophysics

Journal ref: A&A 651, A89 (2021)

First detection of a disk free of volatile elements around a young A-type star: A sign of collisions between rocky planets?

Authors: M. E. van den Ancker, N. P. Gentile Fusillo, T. J. Haworth, C. F. Manara, P. A. Miles-Páez, R. D. Oudmaijer, O. Panic, D. J. M. Petit dit de la Roche, M. G. Petr-Gotzens, M. Vioque

Abstract: Aims. We present the first detailed analysis of the astrophysical parameters of the poorly studied Sco-Cen member HD 152384 and its circumstellar environment. Methods. We analyze newly obtained optical-near-IR XSHOOTER spectra, as well as archival TESS data, of HD 152384. In addition, we use literature photometric data to construct a detailed spectral energy distribution (SED) of the star. Results… ▽ More Aims. We present the first detailed analysis of the astrophysical parameters of the poorly studied Sco-Cen member HD 152384 and its circumstellar environment. Methods. We analyze newly obtained optical-near-IR XSHOOTER spectra, as well as archival TESS data, of HD 152384. In addition, we use literature photometric data to construct a detailed spectral energy distribution (SED) of the star. Results. The photospheric absorption lines in the spectrum of HD 152384 are characteristic of a A0 V star, for which we derive a stellar mass of 2.1 +/- 0.1 M_sun and a stellar age > 4.5 Myr. Superimposed on the photospheric absorption, the optical spectrum also displays double-peaked emission lines of Ca II, Fe I, Mg I and Si I, typical of circumstellar disks. Notably, all Hydrogen and Helium lines appear strictly in absorption. A toy model shows that the observed emission line profiles can be reproduced by emission from a compact (radius < 0.3 au) disk seen at an inclination of ~24 degrees. Further evidence for the presence of circumstellar material comes from the detection of a moderate infrared excess in the SED, similar to those found in extreme debris disk systems. Conclusions. We conclude that HD 152384 is surrounded by a tenuous circumstellar disk which, although rich in refractory elements, is highly depleted of volatile elements. To the best of our knowledge such a disk is unique within the group of young stars. However, it is reminiscent of the disks seen in some white dwarfs, which have been attributed to the disruption of rocky planets. We suggest that the disk around HD 152384 may have a similar origin and may be due to collisions in a newly formed planetary system. △ Less

Submitted 25 May, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 5 pages, 5 figures, accepted for publication in A&A Letters

Journal ref: A&A 651, L11 (2021)

High contrast imaging at 10 microns, a search for exoplanets around: Eps Indi A, Eps Eri, Tau Ceti, Sirius A and Sirius B

Authors: P. Pathak, D. J. M. Petit dit de la Roche, M. Kasper, M. Sterzik, O. Absil, A. Boehle, F. Feng, V. D. Ivanov, M. Janson, H. R. A. Jones, A. Kaufer, H. -U. Käufl, A. -L. Maire, M. Meyer, E. Pantin, R. Siebenmorgen, M. E. van den Ancker, G. Viswanath

Abstract: The direct imaging of rocky exoplanets is one of the major science goals for upcoming large telescopes. The contrast requirement for imaging such planets is challenging. However, the mid-IR (InfraRed) regime provides the optimum contrast to directly detect the thermal signatures of exoplanets in our solar neighbourhood. We aim to exploit novel fast chopping techniques newly developed for astronomy… ▽ More The direct imaging of rocky exoplanets is one of the major science goals for upcoming large telescopes. The contrast requirement for imaging such planets is challenging. However, the mid-IR (InfraRed) regime provides the optimum contrast to directly detect the thermal signatures of exoplanets in our solar neighbourhood. We aim to exploit novel fast chopping techniques newly developed for astronomy with the aid of adaptive optics to look for thermal signatures of exoplanets around bright stars in the solar neighbourhood. We use the upgraded VISIR (Very Large Telescope Imager and Spectrometer for the mid-InfraRed) instrument with high contrast imaging (HCI) capability optimized for observations at 10~$μ$m to look for exoplanets around five nearby ($d$ < 4 pc) stars. The instrument provides an improved signal-to-noise (S/N) by a factor of $\sim$4 in the N-band compared to standard VISIR for a given S/N and time. In this work we achieve a detection sensitivity of sub-mJy, which is sufficient to detect few Jupiter mass planets in nearby systems. Although no detections are made we achieve most sensitive limits within $<2''$ for all the observed targets compared to previous campaigns. For $ε$ Indi A and $ε$ Eri we achieve detection limits very close to the giant planets discovered by RV, with the limits on $ε$ Indi A being the most sensitive to date. Our non-detection therefore supports an older age for $ε$ Indi A. The results presented here show the promise for high contrast imaging and exoplanet detections in the mid-IR regime. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 652, A121 (2021)

New Mid-Infrared Imaging Constraints on Companions and Protoplanetary Disks around six Young Stars

Authors: D. J. M. Petit dit de la Roche, N. Oberg, M. E. van den Ancker, I. Kamp, R. van Boekel, D. Fedele, V. D. Ivanov, M. Kasper, H. U. Käufl, M. Kissler-Patig, P. A. Miles-Páez, E. Pantin, S. P. Quanz, Ch. Rab, R. Siebenmorgen, L. B. F. M. Waters

Abstract: Mid-infrared imaging traces the sub-micron and micron sized dust grains in protoplanetary disks and it offers constraints on the geometrical properties of the disks and potential companions, particularly if those companions have circumplanetary disks. We use the VISIR instrument and its upgrade NEAR on the VLT to take new mid-infrared images of five (pre-)transition disks and one circumstellar dis… ▽ More Mid-infrared imaging traces the sub-micron and micron sized dust grains in protoplanetary disks and it offers constraints on the geometrical properties of the disks and potential companions, particularly if those companions have circumplanetary disks. We use the VISIR instrument and its upgrade NEAR on the VLT to take new mid-infrared images of five (pre-)transition disks and one circumstellar disk with proposed planets and obtain the deepest resolved mid-infrared observations to date in order to put new constraints on the sizes of the emitting regions of the disks and the presence of possible companions. We derotate and stack the data to find the disk properties. Where available we compare the data to ProDiMo (Protoplanetary Disk Model) radiation thermo-chemical models to achieve a deeper understanding of the underlying physical processes within the disks. We apply the circularised PSF subtraction method to find upper limits on the fluxes of possible companions and model companions with circumplanetary disks. We resolve three of the six disks and calculate position angles, inclinations and (upper limits to) sizes of emission regions in the disks, improving upper limits on two of the unresolved disks. In all cases the majority of the mid-IR emission comes from small inner disks or the hot inner rims of outer disks. We refine the existing ProDiMo HD 100546 model SED fit in the mid-IR by increasing the PAH abundance relative to the ISM, adopting coronene as the representative PAH, and increase the outer cavity radius to 22.3 AU. We produce flux estimates for putative planetary-mass companions and circumplanetary disks, ruling out the presence of planetary-mass companions with $L > 0.0028 L_{\odot}$ for $a > 180$ AU in the HD 100546 system. Upper limits of 0.5 mJy-30 mJy are obtained at 8 $μ$m-12 $μ$m for potential companions in the different disks. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: Accepted in A&A

Journal ref: A&A 648, A92 (2021)

Mid-infrared photometry of the T Tauri triple system with kernel phase interferometry

Authors: J. Kammerer, M. Kasper, M. J. Ireland, R. Köhler, R. Laugier, F. Martinache, R. Siebenmorgen, M. E. van den Ancker, R. van Boekel, T. M. Herbst, E. Pantin, H. -U. Käufl, D. J. M. Petit dit de la Roche, V. D. Ivanov

Abstract: T Tauri has long been the prototypical young pre-main-sequence star. However, it has now been decomposed into a triple system with a complex disk and outflow geometry. We aim to measure the brightness of all three components of the T Tauri system (T Tau N, T Tau Sa, T Tau Sb) in the mid-infrared in order to obtain photometry around the $\sim 9.7~μm$ silicate feature. This allows us to study their… ▽ More T Tauri has long been the prototypical young pre-main-sequence star. However, it has now been decomposed into a triple system with a complex disk and outflow geometry. We aim to measure the brightness of all three components of the T Tauri system (T Tau N, T Tau Sa, T Tau Sb) in the mid-infrared in order to obtain photometry around the $\sim 9.7~μm$ silicate feature. This allows us to study their variability and to investigate the distribution of dust and the geometry of circumstellar and circumbinary disks in this complex system. We observe T Tauri with the VLT/VISIR-NEAR instrument. With kernel phase interferometry post-processing of the data, and using the astrometric positions of all three components from VLT/SPHERE, we measure the three components' individual brightnesses (including the southern binary at an angular separation down to $\sim 0.2~λ/D$) and obtain their photometry. In order to validate our methods, we simulate and recover mock data of the T Tauri system using the observed reference point-spread function of HD 27639. We find that T Tau N is rather stable and shows weak silicate emission, while T Tau Sa is highly variable and shows prominent silicate absorption. T Tau Sb became significantly fainter compared to data from 2004 and 2006, suggesting increased extinction by dust. The precision of our photometry is limited by systematic errors, which is consistent with previous studies using kernel phase interferometry. Our results confirm the complex scenario of misaligned disks in the T Tauri system that had been observed previously, and they are in agreement with the recently observed dimming of T Tau Sb in the near-infrared. Our mid-infrared photometry supports the interpretation that T Tau Sb has moved behind the dense region of the Sa-Sb circumbinary disk on its tight orbit around Sa, therefore suffering increased extinction. △ Less

Submitted 4 February, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

Comments: 14 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

Polycyclic Aromatic Hydrocarbons in Seyfert and star-forming galaxies

Authors: I. García-Bernete, D. Rigopoulou, A. Alonso-Herrero, M. Pereira-Santaella, P. F. Roche, B. Kerkeni

Abstract: Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules resulting from the union of aromatic rings and related species, which are likely responsible for strong infrared emission features (3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 microns). In this work, using a sample of 50 Seyfert galaxies (DL<100 Mpc) we compare the circumnuclear (inner kpc) PAH emission of AGN to that of a control sample of… ▽ More Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules resulting from the union of aromatic rings and related species, which are likely responsible for strong infrared emission features (3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 microns). In this work, using a sample of 50 Seyfert galaxies (DL<100 Mpc) we compare the circumnuclear (inner kpc) PAH emission of AGN to that of a control sample of star-forming galaxies (22 luminous infrared galaxies and 30 HII galaxies), and investigate the differences between central and extended PAH emission. Using Spitzer/InfraRed Spectrograph spectral data of Seyfert and star-forming galaxies and newly developed PAH diagnostic model grids, derived from theoretical spectra, we compare the predicted and observed PAH ratios. We find that star-forming galaxies and AGN-dominated systems are located in different regions of the PAH diagnostic diagrams. This suggests that not only are the size and charge of the PAH molecules different, but also the nature and hardness of the radiation field that excite them. We find tentative evidence that PAH ratios in AGN-dominated systems are consistent with emission from larger PAH molecules (Nc>300-400) as well as neutral species. By subtracting the spectrum of the central source from the total, we compare the PAH emission in the central vs extended region of a small sample of AGN. In contrast to the findings for the central regions of AGN-dominated systems, the PAH ratios measured in the extended regions of both type 1 and type 2 Seyfert galaxies can be explained assuming similar PAH molecular size distribution and ionized fractions of molecules to those seen in central regions of star-forming galaxies (100<Nc<300). △ Less

Submitted 26 October, 2021; v1 submitted 21 November, 2020; originally announced November 2020.

Comments: 29 pages, 21 figures, 6 tables, MNRAS, accepted 2021 October 22

The properties of Polycyclic Aromatic Hydrocarbons in galaxies: constraints on PAH sizes, charge and radiation fields

Authors: Dimitra Rigopoulou, Marie Barale, David Clary, Xiao Shan, Almudena Alonso-Herrero, Ismael García-Bernete, Leslie Hunt, Boutheina Kerkeni, Miguel Pereira-Santaella, Patrick Roche

Abstract: Based on theoretical spectra computed using Density Functional Theory we study the properties of Polycyclic Aromatic Hydrocarbons (PAH). In particular using bin-average spectra of PAH molecules with varying number of carbons we investigate how the intensity of the mid-infrared emission bands, 3.3, 6.2, 7.7 and 11.3 microns, respond to changes in the number of carbons, charge of the molecule, and t… ▽ More Based on theoretical spectra computed using Density Functional Theory we study the properties of Polycyclic Aromatic Hydrocarbons (PAH). In particular using bin-average spectra of PAH molecules with varying number of carbons we investigate how the intensity of the mid-infrared emission bands, 3.3, 6.2, 7.7 and 11.3 microns, respond to changes in the number of carbons, charge of the molecule, and the hardness of the radiation field that impinges the molecule. We confirm that the 6.2/7.7 band ratio is a good predictor for the size of the PAH molecule (based on the number of carbons present). We also investigate the efficacy of the 11.3/3.3 ratio to trace the size of PAH molecules and note the dependence of this ratio on the hardness of the radiation field. While the ratio can potentially also be used to trace PAH molecular size, a better understanding of the impact of the underlying radiation field on the 3.3 microns feature and the effect of the extinction on the ratio should be evaluated. The newly developed diagnostics are compared to band ratios measured in a variety of galaxies observed with the Infrared Spectrograph on board the Spitzer Space Telescope. We demonstrate that the band ratios can be used to probe the conditions of the interstellar medium in galaxies and differentiate between environments encountered in normal star forming galaxies and Active Galactic Nuclei. Our work highlights the immense potential that PAH observations with the James Webb Space Telescope will have on our understanding of the PAH emission itself and of the physical conditions in galaxies near and far. △ Less

Submitted 3 May, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 14 pages; 12 figures. Final version accepted for publication in MNRAS, accepted 2021 March 31

Seeing the Bigger Picture: The Rosetta Mission Amateur Observing Campaign and Lessons for the Future

Authors: Helen Usher, Colin Snodgrass, Simon F. Green, Andrew Norton, Paul Roche

Abstract: Amateur astronomers can make useful contributions to the study of comets. They add temporal coverage and multi-scale observations which can aid the study of fast-changing, and large-scale comet features. We document and review the amateur observing campaign set up to complement the Rosetta space mission, including the data submitted to date, and consider the campaign's effectiveness in the light o… ▽ More Amateur astronomers can make useful contributions to the study of comets. They add temporal coverage and multi-scale observations which can aid the study of fast-changing, and large-scale comet features. We document and review the amateur observing campaign set up to complement the Rosetta space mission, including the data submitted to date, and consider the campaign's effectiveness in the light of experience from previous comet amateur campaigns. We report the results of surveys of campaign participants, the amateur astronomy community, and schools who participated in a comet 46P observing campaign. We draw lessons for future campaigns which include the need for: clarity of objectives; recognising the wider impact campaigns can have on increasing science capital; clear, consistent, timely and tailored guidance; easy upload procedures with in-built quality control; and, regular communication, feedback and recognition. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 32 pages, 8 figures To be published in The AAS Planetary Science Journal A'Hearn Symposium Focus Issue

Superorbital Modulation in the High-Mass X-ray Binary 4U 1538-52, and Possible Modulation in IGR J16393-4643

Authors: Robin H. D. Corbet, Joel B. Coley, Hans A. Krimm, Katja Pottschmidt, Paul Roche

Abstract: Hard X-ray observations with the Neil Gehrels Swift Observatory Burst Alert Telescope (BAT) reveal superorbital modulation in the wind-accreting supergiant high-mass X-ray binary (HMXB) 4U 1538-52 at a period of 14.9130 +/- 0.0026 days that is consistent with four times the 3.73 day orbital period. These periods agree with a previously suggested correlation between superorbital and orbital periods… ▽ More Hard X-ray observations with the Neil Gehrels Swift Observatory Burst Alert Telescope (BAT) reveal superorbital modulation in the wind-accreting supergiant high-mass X-ray binary (HMXB) 4U 1538-52 at a period of 14.9130 +/- 0.0026 days that is consistent with four times the 3.73 day orbital period. These periods agree with a previously suggested correlation between superorbital and orbital periods in similar HMXBs. During the ~14 years of observations the superorbital modulation changes amplitude, and since ~MJD 57,650 it was no longer detected in the power spectrum, although a peak near the second harmonic of this was present for some time. Measurements of the spin period of the neutron star in the system with the Fermi Gamma-ray Burst Monitor show a long-term spin-down trend which halted towards the end of the light curve, suggesting a connection between dP(spin)/dt and superorbital modulation, as proposed for 2S 0114+650. However, an earlier torque reversal from INTEGRAL observations was not associated with superorbital modulation changes. B and V band photometry from the Las Cumbres Observatory reveals orbital ellipsoidal photometric variability, but no superorbital optical modulation. However the photometry was obtained when the 14.9130 day period was no longer detected in the BAT power spectrum. We revisit possible superorbital modulation in BAT observations of IGR J16393-4643 but cannot conclusively determine whether this is present, although is not persistent. We consider superorbital modulation mechanisms, and suggest that the Corotating Interaction Region model, with small deviations from orbital synchronization, appears promising. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: Accepted for publication in the Astrophysical Journal. 15 pages, 16 figures

Cold molecular gas and PAH emission in Seyfert galaxies

Authors: A. Alonso-Herrero, M. Pereira-Santaella, D. Rigopoulou, I. Garcia-Bernete, S. Garcia-Burillo, A. J. Dominguez-Fernandez, F. Combes, R. I. Davies, T. Diaz-Santos, D. Esparza-Arredondo, O. Gonzalez-Martin, A. Hernan-Caballero, E. K. S. Hicks, S. F. Hoenig, N. A. Levenson, C. Ramos Almeida, P. F. Roche, D. Rosario

Abstract: We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO… ▽ More We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,μ$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,μ$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A

New constraints on the HR 8799 planetary system from mid-infrared direct imaging

Authors: D. J. M. Petit dit de la Roche, M. E. van den Ancker, M. Kissler-Patig, V. D. Ivanov, D. Fedele

Abstract: Direct imaging is a tried and tested method of detecting exoplanets in the near infrared, but has so far not been extended to longer wavelengths. New data at mid-IR wavelengths (8-20μm) canprovide additional constraints on planetary atmospheric models. We use the VISIR instrumenton the VLT to detect or set stringent limits on the 8.7μm flux of the four planets surrounding HR8799, and to search for… ▽ More Direct imaging is a tried and tested method of detecting exoplanets in the near infrared, but has so far not been extended to longer wavelengths. New data at mid-IR wavelengths (8-20μm) canprovide additional constraints on planetary atmospheric models. We use the VISIR instrumenton the VLT to detect or set stringent limits on the 8.7μm flux of the four planets surrounding HR8799, and to search for additional companions. We use a novel circularised PSF subtractiontechnique to reduce the stellar signal and obtain instrument limited background levels andobtain optimal flux limits. The BT SETTL isochrones are then used to determine the resultingmass limits. We find flux limits between 0.7 and 3.3 mJy for the J8.9 flux of the differentplanets at better than5σlevel and derive a new mass limit of 30 MJupfor any objects beyond40 AU. While this work has not detected planets in the HR 8799 system at 8.7μm, it has foundthat an instrument with the sensitivity of VISIR is sufficient to detect at least 4 known hotplanets around close stars, including\b{eta}Pictoris b (1700 K, 19 pc), with more than5σcertaintyin 10 hours of observing time in the mid-IR. △ Less

Submitted 13 November, 2019; v1 submitted 12 November, 2019; originally announced November 2019.

Comments: Published in MNRAS, 5 pages, 5 figures

High angular resolution ALMA images of dust and molecules in the SN 1987A ejecta

Authors: Phil Cigan, Mikako Matsuura, Haley L. Gomez, Remy Indebetouw, Fran Abellán, Michael Gabler, Anita Richards, Dennis Alp, Tim Davis, Hans-Thomas Janka, Jason Spyromilio, M. J. Barlow, David Burrows, Eli Dwek, Claes Fransson, Bryan Gaensler, Josefin Larsson, P. Bouchet, Peter Lundqvist, J. M. Marcaide, C. -Y. Ng, Sangwook Park, Pat Roche, Jacco Th. van Loon, J. C. Wheeler , et al. (1 additional authors not shown)

Abstract: We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall t… ▽ More We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$α$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 32 pages, containing 19 figures and three appendices

An X-ray and optical study of the outbursting behaviour of the SMC Be X-ray binary SXP 91.1

Authors: Itumeleng Monageng, Malcolm Coe, Jamie Kennea, Lee Townsend, David Buckley, Vanessa McBride, Andrzej Udalski, Phil Evans, Paul Roche

Abstract: In this paper we report on the optical and X-ray behaviour of the Be X-ray binary, SXP 91.1, during a recent type I outburst. We monitored the outburst using the Neil Gehrels Swift Observatory. These data were supported by optical data from the Southern African Large Telescope (SALT) and the Optical Gravitational Lensing Experiment (OGLE) to show the circumstellar disc activity. Matter from this d… ▽ More In this paper we report on the optical and X-ray behaviour of the Be X-ray binary, SXP 91.1, during a recent type I outburst. We monitored the outburst using the Neil Gehrels Swift Observatory. These data were supported by optical data from the Southern African Large Telescope (SALT) and the Optical Gravitational Lensing Experiment (OGLE) to show the circumstellar disc activity. Matter from this disc accretes onto the neutron star, giving rise to the X-ray outburst as seen in the synchronous evolution of the optical and X-ray lightcurves. Using data taken with OGLE we show that the circumstellar disc has exhibited stable behaviour over two decades. A positive correlation is seen between the colour and magnitude from the OGLE and MACHO observations, which indicates that the disc is orientated at relatively low inclination angles. From the OGLE and Swift data, we demonstrate that the system has shown relative phase offsets that have persisted for many years. The spin period derivative is seen to be at maximum spin-up at phases when the mass accretion rate is at maximum. We show that the neutron star in SXP 91.1 is an unusual member of its class in that it has had a consistent spin period derivative over many years, with the average spin-up rate being one of the highest for known SMC pulsars. The most recent measurements of the spin-up rate reveal higher values than the global trend, which is attributed to the recent mass accretion event leading to the current outburst. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS

A nuclear molecular outflow in the Seyfert galaxy NGC3227

Authors: A. Alonso-Herrero, S. Garcia-Burillo, M. Pereira-Santaella, R. I. Davies, F. Combes, M. Vestergaard, S. I. Raimundo, A. Bunker, T. Diaz-Santos, P. Gandhi, I. Garcia-Bernete, E. K. S. Hicks, S. F. Hönig, L. K. Hunt, M. Imanishi, T. Izumi, N. A. Levenson, W. Maciejewski, C. Packham, C. Ramos Almeida, C. Ricci, D. Rigopoulou, P. F. Roche, D. Rosario, M. Schartmann , et al. (2 additional authors not shown)

Abstract: We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emi… ▽ More We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$μ$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged) △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: Accepted for publication in Astronomy and Astrophysics, 17 pages, 12 figures

Torus model properties of an ultra-hard X-ray selected sample of Seyfert galaxies

Authors: I. García-Bernete, C. Ramos Almeida, A. Alonso-Herrero, M. J. Ward, J. A. Acosta-Pulido, M. Pereira-Santaella, A. Hernán-Caballero, A. Asensio Ramos, O. González-Martín, N. A. Levenson, S. Mateos, F. J. Carrera, C. Ricci, P. Roche, I. Marquez, C. Packham, J. Masegosa, L. Fuller

Abstract: We characterize for the first time the torus properties of an ultra-hard X-ray (14-195 keV) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample). The sample was selected from the Swift/BAT nine month catalog. We use high angular resolution nuclear infrared (IR) photometry and N-band spectroscopy, the CLUMPY torus models and a Bayesian tool to characterize the properties of… ▽ More We characterize for the first time the torus properties of an ultra-hard X-ray (14-195 keV) volume-limited (DL<40 Mpc) sample of 24 Seyfert (Sy) galaxies (BCS40 sample). The sample was selected from the Swift/BAT nine month catalog. We use high angular resolution nuclear infrared (IR) photometry and N-band spectroscopy, the CLUMPY torus models and a Bayesian tool to characterize the properties of the nuclear dust. In the case of the Sy1s we estimate the accretion disk contribution to the subarcsecond resolution nuclear IR SEDs (~0.4'') which is, on average, 46+-28, 23+-13 and 11+-5% in the J-, H- and K-bands, respectively. This indicates that the accretion disk templates that assume a steep fall for longer wavelengths than 1 micron might underestimate its contribution to the near-IR emission. Using both optical (broad vs narrow lines) and X-ray (unabsorbed vs absorbed) classifications, we compare the global posterior distribution of the torus model parameters. We confirm that Sy2s have larger values of the torus covering factor (CT~0.95) than Sy1s (CT~0.65) in our volume-limited Seyfert sample. These findings are independent of whether we use an optical or X-ray classification. We find that the torus covering factor remains essentially constant within the errors in our luminosity range and there is no clear dependence with the Eddington ratio. Finally, we find tentative evidence that even an ultra hard X-ray selection is missing a significant fraction of highly absorbed type 2 sources with very high covering factor tori. △ Less

Submitted 7 April, 2019; originally announced April 2019.

Comments: 28 pages, 14 figures and 14 tables. Monthly Notices of the Royal Astronomical Society, Accepted 2019 April 4. Received 2019 March 14; in original form 2019 January 18, in press

Optical precursors to X-ray binary outbursts

Authors: David M. Russell, D. M. Bramich, F. Lewis, A. AlMannaei, T. Al Qaissieh, A. Al Qasim, A. Al Yazeedi, M. C. Baglio, F. Bernardini, N. Elgalad, A. Gabuya, J. P. Lasota, A. Palado, P. Roche, H. Shivkumar, S. Udrescu, G. Zhang

Abstract: Disc instability models predict that for X-ray binaries in quiescence, there should be a brightening of the optical flux prior to an X-ray outburst. Tracking the X-ray variations of X-ray binaries in quiescence is generally not possible, so optical monitoring provides the best means to measure the mass accretion rate variability between outbursts, and to identify the beginning stages of new outbur… ▽ More Disc instability models predict that for X-ray binaries in quiescence, there should be a brightening of the optical flux prior to an X-ray outburst. Tracking the X-ray variations of X-ray binaries in quiescence is generally not possible, so optical monitoring provides the best means to measure the mass accretion rate variability between outbursts, and to identify the beginning stages of new outbursts. With our regular Faulkes Telescope/Las Cumbres Observatory (LCO) monitoring we are routinely detecting the optical rise of new X-ray binary outbursts before they are detected by X-ray all-sky monitors. We present examples of detections of an optical rise in X-ray binaries prior to X-ray detection. We also present initial optical monitoring of the new black hole transient MAXI J1820+070 (ASASSN-18ey) with the Faulkes, LCO telescopes and Al Sadeem Observatory in Abu Dhabi, UAE. Finally, we introduce our new real-time data analysis pipeline, the 'X-ray Binary New Early Warning System (XB-NEWS)' which aims to detect and announce new X-ray binary outbursts within a day of first optical detection. This will allow us to trigger X-ray and multi-wavelength campaigns during the very early stages of outbursts, to constrain the outburst triggering mechanism. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Accepted to Astronomische Nachrichten/Astronomical Notes (proceedings of XMM-Newton workshop 'Time-Domain Astronomy: A High Energy View', ESAC, Madrid, Spain, 13-15 June 2018)

The SMC X-ray binary SXP4.78 : a new Type II outburst and the identification and study of the optical counterpart

Authors: Itumeleng Monageng, Malcolm Coe, Lee Townsend, David Buckley, Vanessa McBride, Paul Roche, Jamie Kennea, Andrzej Udalski, Phil Evans

Abstract: SXP4.78 was originally discovered in 2000 as a pulsar in the Small Magellanic Cloud (SMC) by the Rossi X-ray Timing Explorer (RXTE) but it was not spatially located at that time. A new detection in 2018 with the Neil Gehrels Swift Observatory during a Type II outburst permitted its position to be accurately located and its optical counterpart identified. We report X-ray and optical monitoring cove… ▽ More SXP4.78 was originally discovered in 2000 as a pulsar in the Small Magellanic Cloud (SMC) by the Rossi X-ray Timing Explorer (RXTE) but it was not spatially located at that time. A new detection in 2018 with the Neil Gehrels Swift Observatory during a Type II outburst permitted its position to be accurately located and its optical counterpart identified. We report X-ray and optical monitoring covering epochs before and during the outburst. Using photometric data we show the long-term variability of the Be disc where we present flux and colour changes associated with the disc growth and decay over a period of ~6000 days. We show evidence of disc growth during the recent outburst through an increase in the H-alpha equivalent width and photometric flux. Period analysis was performed using both optical photometric and spectroscopic data, but with no significant detection of an orbital period. A modest periodic signature of 2.65 days was detected from the OGLE I band data, however, but we attribute that to the non-radial pulsations (NRPs) of the Be star. We also obtained a blue spectrum from the Southern African Large Telescope (SALT) which permits us to classify the spectral type as B0.5 IV-V. △ Less

Submitted 11 March, 2019; v1 submitted 2 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS after minor revision

Molecule mapping of HR8799b using OSIRIS on Keck: Strong detection of water and carbon monoxide, but no methane

Authors: D. J. M. Petit dit de la Roche, H. J. Hoeijmakers, I. A. G. Snellen

Abstract: Context. In 2015, Barman et al. (ApJ, 804, 61) presented detections of absorption from water, carbon monoxide, and methane in the atmosphere of the directly imaged exoplanet HR8799b using integral field spectroscopy (IFS) with OSIRIS on the Keck II telescope. We recently devised a new method to analyse IFU data, called molecule mapping, searching for high-frequency signatures of particular molecul… ▽ More Context. In 2015, Barman et al. (ApJ, 804, 61) presented detections of absorption from water, carbon monoxide, and methane in the atmosphere of the directly imaged exoplanet HR8799b using integral field spectroscopy (IFS) with OSIRIS on the Keck II telescope. We recently devised a new method to analyse IFU data, called molecule mapping, searching for high-frequency signatures of particular molecules in an IFU data cube. Aims. The aim of this paper is to use the molecule mapping technique to search for the previously detected spectral signatures in HR8799b using the same data, allowing a comparison of molecule mapping with previous methods. Methods. The medium-resolution H- and K-band pipeline-reduced archival data were retrieved from the Keck archive facility. Telluric and stellar lines were removed from each spectrum in the data cube, after which the residuals were cross-correlated with model spectra of carbon monoxide, water, and methane. Results. Both carbon monoxide and water are clearly detected at high signal-to-noise, however, methane is not retrieved. Conclusions. Molecule mapping works very well on the OSIRIS data of exoplanet HR8799b. However, it is not evident why methane is detected in the original analysis, but not with the molecule mapping technique. Possible causes could be the presence of telluric residuals, different spectral filtering techniques, or the use of different methane models. We do note that in the original analysis methane was only detected in the K-band, while the H-band methane signal could be expected to be comparably strong. More sensitive observations with the JWST will be capable of confirming or disproving the presence of methane in this planet at high confidence. △ Less

Submitted 31 August, 2018; originally announced August 2018.

Comments: 5 pages, 5 figures and 2 tables, accepted by A&A

Journal ref: D. J. M. Petit dit de la Roche, H. J. Hoeijmakers and I. A. G. Snellen, 2018, A&A, 616, A146

High Resolution Imaging of the Magnetic Field in the central parsec of the Galaxy

Authors: P. F. Roche, E. Lopez-Rodriguez, C. M. Telesco, R. Schodel, C. Packham

Abstract: We discuss a high resolution (FWHM~ 0.45 arcsec) image of the emissive polarization from warm dust in the minispiral in the Galactic Centre and discuss the implications for the magnetic field in the dusty filaments. The image was obtained at a wavelength of 12.5 microns with the CanariCam multimode mid-infrared imager on the Gran Telescopio Canarias. It confirms the results obtained from previous… ▽ More We discuss a high resolution (FWHM~ 0.45 arcsec) image of the emissive polarization from warm dust in the minispiral in the Galactic Centre and discuss the implications for the magnetic field in the dusty filaments. The image was obtained at a wavelength of 12.5 microns with the CanariCam multimode mid-infrared imager on the Gran Telescopio Canarias. It confirms the results obtained from previous observations but also reveals new details of the polarization structures. In particular, we identify regions of coherent magnetic field emission at position angles of ~45 deg to the predominantly north--south run of field lines in the Northern Arm which may be related to orbital motions inclined to the general flow of the Northern Arm. The luminous stars that have been identified as bow-shock sources in the Northern Arm do not disrupt or dilute the field but are linked by a coherent field structure, implying that the winds from these objects may push and compress the field but do not overwhelm it. The magnetic field in the the low surface brightness regions in the East-West Bar to the south of SgrA* lies along the Bar, but the brighter regions generally have different polarization position angles, suggesting that they are distinct structures. In the region of the Northern Arm sampled here, there is only a weak correlation between the intensity of the emission and the degree of polarization. This is consistent with saturated grain alignment where the degree of polarization depends on geometric effects, including the angle of inclination of the field to the line of sight and superposition of filaments with different field directions, rather than the alignment efficiency. △ Less

Submitted 31 July, 2018; originally announced July 2018.

Comments: 9 pages, 3 figures, Proceedings of Cosmic Dust X, held in Mitaka, Japan in August 2017

New Cepheid variables in the young open clusters Berkeley 51 and Berkeley 55

Authors: M. E. Lohr, I. Negueruela, H. M. Tabernero, J. S. Clark, F. Lewis, P. Roche

Abstract: As part of a wider investigation of evolved massive stars in Galactic open clusters, we have spectroscopically identified three candidate classical Cepheids in the little-studied clusters Berkeley 51, Berkeley 55 and NGC 6603. Using new multi-epoch photometry, we confirm that Be 51 #162 and Be 55 #107 are bona fide Cepheids, with pulsation periods of 9.83+/-0.01 d and 5.850+/-0.005 d respectively,… ▽ More As part of a wider investigation of evolved massive stars in Galactic open clusters, we have spectroscopically identified three candidate classical Cepheids in the little-studied clusters Berkeley 51, Berkeley 55 and NGC 6603. Using new multi-epoch photometry, we confirm that Be 51 #162 and Be 55 #107 are bona fide Cepheids, with pulsation periods of 9.83+/-0.01 d and 5.850+/-0.005 d respectively, while NGC 6603 star W2249 does not show significant photometric variability. Using the period-luminosity relationship for Cepheid variables, we determine a distance to Be 51 of 5.3(+1.0,-0.8) kpc and an age of 44(+9,-8) Myr, placing it in a sparsely-attested region of the Perseus arm. For Be 55, we find a distance of 2.2+/-0.3 kpc and age of 63(+12,-11) Myr, locating the cluster in the Local arm. Taken together with our recent discovery of a long-period Cepheid in the starburst cluster VdBH222, these represent an important increase in the number of young, massive Cepheids known in Galactic open clusters. We also consider new Gaia (data release 2) parallaxes and proper motions for members of Be 51 and Be 55; the uncertainties on the parallaxes do not allow us to refine our distance estimates to these clusters, but the well-constrained proper motion measurements furnish further confirmation of cluster membership. However, future final Gaia parallaxes for such objects should provide valuable independent distance measurements, improving the calibration of the period-luminosity relationship, with implications for the distance ladder out to cosmological scales. △ Less

Submitted 16 May, 2018; v1 submitted 15 May, 2018; originally announced May 2018.

Comments: 8 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society. Author details corrected

The Magnetic Field in the central parsec of the Galaxy

Authors: PF Roche, E Lopez-Rodriguez, CM Telesco, R Schodel, C Packham

Abstract: We present a polarisation map of the warm dust emission from the minispiral in the central parsec of the Galactic centre. The observations were made at a wavelength of 12.5 microns with CanariCam mounted on the 10.4-m Gran Telescopio Canarias. The magnetic field traced by the polarised emission from aligned dust grains is consistent with previous observations, but the increased resolution of the p… ▽ More We present a polarisation map of the warm dust emission from the minispiral in the central parsec of the Galactic centre. The observations were made at a wavelength of 12.5 microns with CanariCam mounted on the 10.4-m Gran Telescopio Canarias. The magnetic field traced by the polarised emission from aligned dust grains is consistent with previous observations, but the increased resolution of the present data reveals considerably more information on the detailed structure of the B field and its correspondence with the filamentary emission seen in both mid-infrared continuum emission and free-free emission at cm wavelengths. The magnetic field appears to be compressed and pushed by the outflows from luminous stars in the Northern Arm, but it is not disordered by them. We identify some magnetically coherent filaments that cross the Northern Arm at a Position Angle of ~45 degrees, and which may trace orbits inclined to the primary orientation of the Northern Arm and circumnuclear disk. In the East-West bar, the magnetic fields implied by the polarization in the lower intensity regions lie predominantly along the bar at a Position Angle of 130 - 140 degrees. In contrast to the Northern Arm, the brighter regions of the bar tend to have lower degrees of polarization with a greater divergence in position angle compared to the local diffuse emission. It appears that the diffuse emission in the East-West bar traces the underlying field and that the bright compact sources are unrelated objects presumably projected onto the bar and with different field orientations △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: 12 Pages, 5 figures, 1 Table. To be published in MNRAS

A mid-infrared statistical investigation of clumpy torus model predictions

Authors: J. García-González, A. Alonso-Herrero, S. F. Hönig, A. Hernán-Caballero, C. Ramos Almeida, N. A. Levenson, P. F. Roche, O. González-Martín, C. Packham, M. Kishimoto

Abstract: We present new calculations of the CAT3D clumpy torus models, which now include a more physical dust sublimation model as well as AGN anisotropic emission. These new models allow graphite grains to persist at temperatures higher than the silicate dust sublimation temperature. This produces stronger near-infrared emission and bluer mid-infrared (MIR) spectral slopes. We make a statistical compariso… ▽ More We present new calculations of the CAT3D clumpy torus models, which now include a more physical dust sublimation model as well as AGN anisotropic emission. These new models allow graphite grains to persist at temperatures higher than the silicate dust sublimation temperature. This produces stronger near-infrared emission and bluer mid-infrared (MIR) spectral slopes. We make a statistical comparison of the CAT3D model MIR predictions with a compilation of sub-arcsecond resolution ground-based MIR spectroscopy of 52 nearby Seyfert galaxies (median distance of 36 Mpc) and 10 quasars. We focus on the AGN MIR spectral index $α_{MIR}$ and the strength of the 9.7 $μ$m silicate feature $S_{Sil}$. As with other clumpy torus models, the new CAT3D models do not reproduce the Seyfert galaxies with deep silicate absorption ($S_{Sil}<-1$). Excluding those, we conclude that the new CAT3D models are in better agreement with the observed $α_{MIR}$ and $S_{Sil}$ of Seyfert galaxies and quasars. We find that Seyfert 2 are reproduced with models with low photon escape probabilities, while the quasars and the Seyfert 1-1.5 require generally models with higher photon escape probabilities. Quasars and Seyfert 1-1.5 tend to show steeper radial cloud distributions and fewer clouds along an equatorial line-of-sight than Seyfert 2. Introducing AGN anisotropic emission besides the more physical dust sublimation models alleviates the problem of requiring inverted radial cloud distributions (i.e., more clouds towards the outer parts of the torus) to explain the MIR spectral indices of type 2 Seyferts. △ Less

Submitted 23 June, 2017; v1 submitted 22 June, 2017; originally announced June 2017.

Comments: 23 pages, 17 figures, 5 tables, Accepted for publication in MNRAS

Very Deep Inside the SN 1987A Core Ejecta: Molecular Structures Seen in 3D

Authors: F. J. Abellán, R. Indebetouw, J. M. Marcaide, M. Gabler, C. Fransson, J. Spyromilio, D. N. Burrows, R. Chevalier, P. Cigan, B. M. Gaensler, H. L. Gomez, H. -Th. Janka, R. Kirshner, J. Larsson, P. Lundqvist, M. Matsuura, R. McCray, C. -Y. Ng, S. Park, P. Roche, L. Staveley-Smith, J. Th. Van Loon, J. C. Wheeler, S. E. Woosley

Abstract: Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct p… ▽ More Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: ApJL accepted, 9 pages, 5 figures

ALMA spectral survey of Supernova 1987A --- molecular inventory, chemistry, dynamics and explosive nucleosynthesis

Authors: M. Matsuura, R. Indebetouw, S. Woosley, V. Bujarrabal, F. J. Abellan, R. McCray, J. Kamenetzky, C. Fransson, M. J. Barlow, H. L. Gomez, P. Cigan, I De Looze, J. Spyromilio, L. Staveley-Smith, G. Zanardo, P. Roche, J. Larsson, S. Viti, J. Th. van Loon, J. C. Wheeler, M. Baes, R. Chevalier, P. Lundqvist, J. M. Marcaide, E. Dwek , et al. (4 additional authors not shown)

Abstract: We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta mo… ▽ More We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae. △ Less

Submitted 7 April, 2017; originally announced April 2017.

Comments: Accepted by MNRAS, 16 pages

Infrared polarimetry of Mrk 231: Scattering off hot dust grains in the central core

Authors: E. Lopez-Rodriguez, C. Packham, T. J. Jones, R. Siebenmorgen, P. F. Roche, N. A. Levenson, A. Alonso-Herrero, E. Perlman, K. Ichikawa, C. Ramos Almeida, O. Gonzalez-Martin, R. Nikutta, M. Martinez-Paredes, D. Shenoy, M. S. Gordon, C. M. Telesco

Abstract: We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $μ$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $μ$m, 10.3 $μ$m, and 11.6 $μ$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $μ$m total and polarized nuclear spectral energy distributio… ▽ More We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $μ$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $μ$m, 10.3 $μ$m, and 11.6 $μ$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $μ$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $μ$m to 0.31 $\pm$ 0.15 per cent at 11.6 $μ$m and follows a power-law function, $λ^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $μ$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions. △ Less

Submitted 28 September, 2016; originally announced September 2016.

Comments: 9 pages, 2 figures. Accepted for publication in MNRAS