Showing 1–50 of 114 results for author: Yates, J

Bayesian model reconstruction based on spectral line observations

Authors: Frederik De Ceuster, Thomas Ceulemans, Leen Decin, Taïssa Danilovich, Jeremy Yates

Abstract: Spectral line observations encode a wealth of information. A key challenge, therefore, lies in the interpretation of these observations in terms of models to derive the physical and chemical properties of the astronomical environments from which they arise. In this paper, we present pomme: an open-source Python package that allows users to retrieve 1D or 3D models of physical properties, such as c… ▽ More Spectral line observations encode a wealth of information. A key challenge, therefore, lies in the interpretation of these observations in terms of models to derive the physical and chemical properties of the astronomical environments from which they arise. In this paper, we present pomme: an open-source Python package that allows users to retrieve 1D or 3D models of physical properties, such as chemical abundance, velocity, and temperature distributions of (optically thin) astrophysical media, based on spectral line observations. We discuss how prior knowledge, for instance, in the form of a steady-state hydrodynamics model, can be used to guide the retrieval process, and demonstrate our methods both on synthetic and real observations of cool stellar winds. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ

Optimisation-based alignment of wideband integrated superconducting spectrometers for sub-mm astronomy

Authors: A. Moerman, K. Karatsu, S. J. C. Yates, R. Huiting, F. Steenvoorde, S. O. Dabironezare, T. Takekoshi, J. J. A. Baselmans, B. R. Brandl, A. Endo

Abstract: Integrated superconducting spectrometers (ISSs) for wideband sub-mm astronomy utilise quasi-optical systems for coupling radiation from the telescope to the instrument. Misalignment in these systems is detrimental to the system performance. The common method of using an optical laser to align the quasi-optical components requires accurate alignment of the laser to the sub-mm beam coming from the i… ▽ More Integrated superconducting spectrometers (ISSs) for wideband sub-mm astronomy utilise quasi-optical systems for coupling radiation from the telescope to the instrument. Misalignment in these systems is detrimental to the system performance. The common method of using an optical laser to align the quasi-optical components requires accurate alignment of the laser to the sub-mm beam coming from the instrument, which is not always guaranteed to a sufficient accuracy. We develop an alignment strategy for wideband ISSs directly utilising the sub-mm beam of the wideband ISS. The strategy should be applicable in both telescope and laboratory environments. Moreover, the strategy should deliver similar quality of the alignment across the spectral range of the wideband ISS. We measure misalignment in a quasi-optical system operating at sub-mm wavelengths using a novel phase and amplitude measurement scheme, capable of simultaneously measuring the complex beam patterns of a direct-detecting ISS across a harmonic range of frequencies. The direct detection nature of the MKID detectors in our device-under-test, DESHIMA 2.0, necessitates the use of this measurement scheme. Using geometrical optics, the measured misalignment, a mechanical hexapod, and an optimisation algorithm, we follow a numerical approach to optimise the positioning of corrective optics with respect to a given cost function. Laboratory measurements of the complex beam patterns are taken across a harmonic range between 205 and 391 GHz and simulated through a model of the ASTE telescope in order to assess the performance of the optimisation at the ASTE telescope. Laboratory measurements show that the optimised optical setup corrects for tilts and offsets of the sub-mm beam. Moreover, we find that the simulated telescope aperture efficiency is increased across the frequency range of the ISS after the optimisation. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 10 pages, 12 figures. Accepted for publication by A&A

Journal ref: A&A 684, A161 (2024)

ATOMIUM: Molecular inventory of 17 oxygen-rich evolved stars observed with ALMA

Authors: S. H. J. Wallstrom, T. Danilovich, H. S. P. Muller, C. A. Gottlieb, S. Maes, M. Van de Sande, L. Decin, A. M. S. Richards, A. Baudry, J. Bolte, T. Ceulemans, F. De Ceuster, A. de Koter, I. El Mellah, M. Esseldeurs, S. Etoka, D. Gobrecht, E. Gottlieb, M. Gray, F. Herpin, M. Jeste, D. Kee, P. Kervella, T. Khouri, E. Lagadec , et al. (13 additional authors not shown)

Abstract: The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but… ▽ More The dusty winds of cool evolved stars are a major contributor of the newly synthesised material enriching the Galaxy and future generations of stars. However, the details of the physics and chemistry behind dust formation and wind launching have yet to be pinpointed. Recent spatially resolved observations show the importance of gaining a more comprehensive view of the circumstellar chemistry, but a comparative study of the intricate interplay between chemistry and physics is still difficult because observational details such as frequencies and angular resolutions are rarely comparable. Aiming to overcome these deficiencies, ATOMIUM is an ALMA Large Programme to study the physics and chemistry of the circumstellar envelopes of a diverse set of oxygen-rich evolved stars under homogeneous observing conditions at three angular resolutions between ~0.02"-1.4". Here we summarize the molecular inventory of these sources, and the correlations between stellar parameters and molecular content. Seventeen oxygen-rich or S-type asymptotic giant branch (AGB) and red supergiant (RSG) stars have been observed in several tunings with ALMA Band 6, targeting a range of molecules to probe the circumstellar envelope and especially the chemistry of dust formation close to the star. We systematically assigned the molecular carriers of the spectral lines and measured their spectroscopic parameters and the angular extent of the emission of each line from integrated intensity maps. Across the ATOMIUM sample, we detect 291 transitions of 24 different molecules and their isotopologues. This includes several first detections in oxygen-rich AGB/RSG stars: PO v=1, SO2 v1=1 and v2=2, and several high energy H2O transitions. We also find several first detections in S-type AGB stars: vibrationally excited HCN v2=2,3 and SiS v=4,5,6, as well as first detections of the molecules SiC, AlCl, and AlF in W Aql... △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 19 pages plus appendices, forthcoming publication in Astronomy & Astrophysics

Applying Energy Absorption Interferometry to THz direct detectors using photomixers

Authors: Ian Veenendaal, Edgar Castillo-Dominguez, Stephen J. C. Yates, Bram Lap, Willem Jellema

Abstract: Detector requirements for far infrared astronomy generally result in devices which exhibit a few-moded response to incident radiation. The sensitivity and spatial form of the individual modes to which such a detector is sensitive can be determined with knowledge of the complex valued cross-spectral density of the system, which we label the detector response function (DRF). A matrix representing th… ▽ More Detector requirements for far infrared astronomy generally result in devices which exhibit a few-moded response to incident radiation. The sensitivity and spatial form of the individual modes to which such a detector is sensitive can be determined with knowledge of the complex valued cross-spectral density of the system, which we label the detector response function (DRF). A matrix representing the discretized cross-spectral density can be measured from the complex amplitudes of interference fringes generated by two identical sources as they are independently scanned through the field of view. We provide experimental verification of this technique using monochromatic THz beams generated by photomixers in which the relative phase is varied with fiber stretchers. We use this system to characterize the modal response of a single pixel from an array of microwave kinetic inductance detectors (MKIDs). △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: in IEEE Transactions on Terahertz Science and Technology. 2023

The VLT/SPHERE view of the ATOMIUM cool evolved star sample. I. Overview: Sample characterization through polarization analysis

Authors: M. Montargès, E. Cannon, A. de Koter, T. Khouri, E. Lagadec, P. Kervella, L. Decin, I. McDonald, W. Homan, L. B. F. M. Waters, R. Sahai, C. A. Gottlieb, J. Malfait, S. Maes, B. Pimpanuwat, M. Jeste, T. Danilovich, F. De Ceuster, M. Van de Sande, D. Gobrecht, S. H. J. Wallström, K. T. Wong, I. El Mellah, J. Bolte, F. Herpin , et al. (10 additional authors not shown)

Abstract: Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution dat… ▽ More Aims. Through the ATOMIUM project, based on an ALMA large program, we aim to present a consistent view of a sample of 17 nearby cool evolved stars (Aymptotic Giant Branch and red supergiant stars). Methods. Here we present VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function, and in or near the plane of the sky, with a total optical depth close to unity in the line of sight, representing only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03-0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of common dust species. The spatial structure of the DoLP shows a diverse set of shapes. Only for three sources do we note a correlation between the ALMA CO and SiO lines, which trace the gas density, and the DoLP, which traces the dust. Conclusion. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner circumstellar environment (CSE), dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the CSE. Except for $π^1$~Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 22 pages, 15 figures, 5 tables

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

Radiative Transfer as a Bayesian Linear Regression problem

Authors: Frederik De Ceuster, Thomas Ceulemans, Jon Cockayne, Leen Decin, Jeremy Yates

Abstract: Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptio… ▽ More Electromagnetic radiation plays a crucial role in various physical and chemical processes. Hence, almost all astrophysical simulations require some form of radiative transfer model. Despite many innovations in radiative transfer algorithms and their implementation, realistic radiative transfer models remain very computationally expensive, such that one often has to resort to approximate descriptions. The complexity of these models makes it difficult to assess the validity of any approximation and to quantify uncertainties on the model results. This impedes scientific rigour, in particular, when comparing models to observations, or when using their results as input for other models. We present a probabilistic numerical approach to address these issues by treating radiative transfer as a Bayesian linear regression problem. This allows us to model uncertainties on the input and output of the model with the variances of the associated probability distributions. Furthermore, this approach naturally allows us to create reduced-order radiative transfer models with a quantifiable accuracy. These are approximate solutions to exact radiative transfer models, in contrast to the exact solutions to approximate models that are often used. As a first demonstration, we derive a probabilistic version of the method of characteristics, a commonly-used technique to solve radiative transfer problems. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 16 pages, 1 figure, accepted to Monthly Notices of the Royal Astronomical Society (MNRAS)

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

Authors: C. A. Gottlieb, L. Decin, A. M. S. Richards, F. De Ceuster, W. Homan, S. H. J. Wallstrom, T. Danilovich, T. J. Millar, M. Montarges, K. T. Wong, I. McDonald, A. Baudry, J. Bolte, E. Cannon, E. De Beck, A. de Koter, I. El Mellah, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, P. Kervella, T. Khouri, E. Lagadec , et al. (11 additional authors not shown)

Abstract: This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical e… ▽ More This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field. △ Less

Submitted 13 December, 2021; v1 submitted 8 December, 2021; originally announced December 2021.

Journal ref: A&A 660, A94 (2022)

DESHIMA 2.0: development of an integrated superconducting spectrometer for science-grade astronomical observations

Authors: Akio Taniguchi, Tom J. L. C. Bakx, Jochem J. A. Baselmans, Robert Huiting, Kenichi Karatsu, Nuria Llombart, Matus Rybak, Tatsuya Takekoshi, Yoichi Tamura, Hiroki Akamatsu, Stefanie Brackenhoff, Juan Bueno, Bruno T. Buijtendorp, Shahab Dabironezare, Anne-Kee Doing, Yasunori Fujii, Kazuyuki Fujita, Matthijs Gouwerok, Sebastian Hähnle, Tsuyoshi Ishida, Shun Ishii, Ryohei Kawabe, Tetsu Kitayama, Kotaro Kohno, Akira Kouchi , et al. (10 additional authors not shown)

Abstract: Integrated superconducting spectrometer (ISS) technology will enable ultra-wideband, integral-field spectroscopy for (sub)millimeter-wave astronomy, in particular, for uncovering the dust-obscured cosmic star formation and galaxy evolution over cosmic time. Here we present the development of DESHIMA 2.0, an ISS for ultra-wideband spectroscopy toward high-redshift galaxies. DESHIMA 2.0 is designed… ▽ More Integrated superconducting spectrometer (ISS) technology will enable ultra-wideband, integral-field spectroscopy for (sub)millimeter-wave astronomy, in particular, for uncovering the dust-obscured cosmic star formation and galaxy evolution over cosmic time. Here we present the development of DESHIMA 2.0, an ISS for ultra-wideband spectroscopy toward high-redshift galaxies. DESHIMA 2.0 is designed to observe the 220-440 GHz band in a single shot, corresponding to a redshift range of $z$=3.3-7.6 for the ionized carbon emission ([C II] 158 $μ$m). The first-light experiment of DESHIMA 1.0, using the 332-377 GHz band, has shown an excellent agreement among the on-sky measurements, the lab measurements, and the design. As a successor to DESHIMA 1.0, we plan the commissioning and the scientific observation campaign of DESHIMA 2.0 on the ASTE 10-m telescope in 2023. Ongoing upgrades for the full octave-bandwidth system include the wideband 347-channel chip design and the wideband quasi-optical system. For efficient measurements, we also develop the observation strategy using the mechanical fast sky-position chopper and the sky-noise removal technique based on a novel data-scientific approach. In the paper, we show the recent status of the upgrades and the plans for the scientific observation campaign. △ Less

Submitted 4 October, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: 9 pages, 4 figures, 1 table, accepted for publication in the Journal of Low Temperature Physics for the LTD19 special issue

ATOMIUM: Halide molecules around the S-type AGB star W Aquilae

Authors: T. Danilovich, M. Van de Sande, J. M. C. Plane, T. J. Millar, P. Royer, M. A. Amor, K. Hammami, L. Decock, C. A. Gottlieb, L. Decin, A. M. S. Richards, E. De Beck, A. Baudry, J. Bolte, E. Cannon, F. De Ceuster, A. de Koter, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, W. Homan, M. Jeste, P. Kervella, T. Khouri , et al. (14 additional authors not shown)

Abstract: S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed… ▽ More S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed for the first time towards an S-type AGB star, W Aql. In combination with models based on PACS observations, we aim to update our chemical kinetics network based on these results. We analyse ALMA observations towards W Aql of AlCl in the ground and first two vibrationally excited states and AlF in the ground vibrational state. Using radiative transfer models, we determine the abundances and spatial abundance distributions of Al$^{35}$Cl, Al$^{37}$Cl, and AlF. We also model HCl and HF emission and compare these models to PACS spectra to constrain the abundances of these species. AlCl is found in clumps very close to the star, with emission confined within 0.1$^{\prime\prime}$ of the star. AlF emission is more extended, with faint emission extending 0.2$^{\prime\prime}$ to 0.6$^{\prime\prime}$ from the continuum peak. We find peak abundances, relative to H$_2$, of $1.7\times 10^{-7}$ for Al$^{35}$Cl, $7\times 10^{-8}$ for Al$^{37}$Cl and $1\times 10^{-7}$ for AlF. From the PACS spectra, we find abundances of $9.7\times 10^{-8}$ and $\leq 10^{-8}$, relative to H$_2$, for HCl and HF, respectively. The AlF abundance exceeds the solar F abundance, indicating that fluorine synthesised in the AGB star has already been dredged up to the surface of the star and ejected into the circumstellar envelope. From our analysis of chemical reactions in the wind, we conclude that AlF may participate in the dust formation process, but we cannot fully explain the rapid depletion of AlCl seen in the wind. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 655, A80 (2021)

Atomium: The astounding complexity of the near circumstellar environment of the M-type AGB star R Hydrae. I. Morpho-kinematical interpretation of CO and SiO emission

Authors: Ward Homan, Bannawit Pimpanuwat, Fabrice Herpin, Taissa Danilovich, Iain McDonald, Sofia H. J. Wallström, Anita M. S. Richards, Alain Baudry, Raghvendra Sahai, Tom J. Millar, Alex de Koter, C. A. Gottlieb, Pierre Kervella, Miguel Montargès, Marie Van de Sande, Leen Decin, Albert Zijlstra, Sandra Etoka, Manali Jeste, Holger S. P. Müller, Silke Maes, Jolien Malfait, Karl Menten, John Plane, Kelvin Lee , et al. (14 additional authors not shown)

Abstract: Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we… ▽ More Evolved low- to intermediate-mass stars are known to shed their gaseous envelope into a large, dusty, molecule-rich circumstellar nebula which typically develops a high degree of structural complexity. Most of the large-scale, spatially correlated structures in the nebula are thought to originate from the interaction of the stellar wind with a companion. As part of the Atomium large programme, we observed the M-type asymptotic giant branch (AGB) star R Hydrae with ALMA. The morphology of the inner wind of R Hya, which has a known companion at ~3500 au, was determined from maps of CO and SiO obtained at high angular resolution. A map of the CO emission reveals a multi-layered structure consisting of a large elliptical feature at an angular scale of ~10'' that is oriented along the north-south axis. The wind morphology within the elliptical feature is dominated by two hollow bubbles. The bubbles are on opposite sides of the AGB star and lie along an axis with a position angle of ~115 deg. Both bubbles are offset from the central star, and their appearance in the SiO channel maps indicates that they might be shock waves travelling through the AGB wind. An estimate of the dynamical age of the bubbles yields an age of the order of 100 yr, which is in agreement with the previously proposed elapsed time since the star last underwent a thermal pulse. When the CO and SiO emission is examined on subarcsecond angular scales, there is evidence for an inclined, differentially rotating equatorial density enhancement, strongly suggesting the presence of a second nearby companion. The position angle of the major axis of this disc is ~70 deg in the plane of the sky. We tentatively estimate that a lower limit on the mass of the nearby companion is ~0.65 Msol on the basis of the highest measured speeds in the disc and the location of its inner rim at ~6 au from the AGB star. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 21 pages, 23 figures

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

Magritte, a modern software library for 3D radiative transfer: II. Adaptive ray-tracing, mesh construction and reduction

Authors: Frederik De Ceuster, Jan Bolte, Ward Homan, Silke Maes, Jolien Malfait, Leen Decin, Jeremy Yates, Peter Boyle, James Hetherington

Abstract: Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable… ▽ More Radiative transfer is a notoriously difficult and computationally demanding problem. Yet, it is an indispensable ingredient in nearly all astrophysical and cosmological simulations. Choosing an appropriate discretization scheme is a crucial part of the simulation, since it not only determines the direct memory cost of the model but also largely determines the computational cost and the achievable accuracy. In this paper, we show how an appropriate choice of directional discretization scheme as well as spatial model mesh can help alleviate the computational cost, while largely retaining the accuracy. First, we discuss the adaptive ray-tracing scheme implemented in our 3D radiative transfer library Magritte, that adapts the rays to the spatial mesh and uses a hierarchical directional discretization based on HEALPix. Second, we demonstrate how the free and open-source software library Gmsh can be used to generate high quality meshes that can be easily tailored for Magritte. In particular, we show how the local element size distribution of the mesh can be used to optimise the sampling of both analytically and numerically defined models. Furthermore, we show that when using the output of hydrodynamics simulations as input for a radiative transfer simulation, the number of elements in the input model can often be reduced by an order of magnitude, without significant loss of accuracy in the radiation field. We demonstrate this for two models based on a hierarchical octree mesh resulting from adaptive mesh refinement (AMR), as well as two models based on smoothed-particle hydrodynamics (SPH) data. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Comments: 12 pages, published in MNRAS

LeMMINGs. II. The e-MERLIN legacy survey of nearby galaxies. The deepest radio view of the Palomar sample on parsec scale

Authors: R. D. Baldi, D. R. A. Williams, I. M. McHardy, R. J. Beswick, E. Brinks, B. T. Dullo, J. H. Knapen, M. K. Argo, S. Aalto, A. Alberdi, W. A. Baan, G. J. Bendo, S. Corbel, D. M. Fenech, J. S. Gallagher, D. A. Green, R. C. Kennicutt, H. -R. Klöckner, E. Körding, T. J. Maccarone, T. W. B. Muxlow, C. G. Mundell, F. Panessa, A. B. Peck, M. A. Pérez-Torres , et al. (8 additional authors not shown)

Abstract: We present the second data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the LeMMINGs (Legacy e-MERLIN Multi-band Imaging of Nearby Galaxy Sample) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and… ▽ More We present the second data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the LeMMINGs (Legacy e-MERLIN Multi-band Imaging of Nearby Galaxy Sample) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and Seyfert) and inactive galaxies HII galaxies and Absorption Line Galaxies, ALG). This large program is the deepest radio survey of the local Universe, $\gtrsim$10$^{17.6}$ W Hz$^{-1}$, regardless of the host and nuclear type: we detect radio emission $\gtrsim$0.25 mJy beam$^{-1}$ for 125/280 galaxies (44.6 per cent) with sizes of typically $\lesssim$100 pc. Of those 125, 106 targets show a core which coincides within 1.2 arcsec with the optical nucleus. Although we observed mostly cores, around one third of the detected galaxies features jetted morphologies. The detected radio core luminosities of the sample range between $\sim$10$^{34}$ and 10$^{40}$ erg s$^{-1}$. LINERs and Seyferts are the most luminous sources, whereas HII galaxies are the least. LINERs show FRI-like core-brightened radio structures, while Seyferts reveal the highest fraction of symmetric morphologies. The majority of HII galaxies have single radio core or complex extended structures, which probably conceal a nuclear starburst and/or a weak active nucleus (seven of them show clear jets). ALGs, which are typically found in evolved ellipticals, although the least numerous, exhibit on average the most luminous radio structures, similar to LINERs. △ Less

Submitted 5 November, 2020; originally announced November 2020.

Comments: Accepted for publication on MNRAS (48 pages, 7 figures, 10 tables, Appendix)

Atomium: A high-resolution view on the highly asymmetric wind of the AGB star Pi1 Gruis. I. First detection of a new companion and its effect on the inner wind

Authors: Ward Homan, Miguel Montarges, Bannawit Pimpanuwat, Anita M. S. Richards, Sofia H. J. Wallstrom, Pierre Kervella, Leen Decin, Albert Zijlstra, Taissa Danilovich, Alex de Koter, Karl Menten, Raghvendra Sahai, John Plane, Kelvin Lee, Rens Waters, Alain Baudry, Ka Tat Wong, Tom J. Millar, Marie Van de Sande, Eric Lagadec, David Gobrecht, Jeremy Yates, Daniel Price, Emily Cannon, Jan Bolte , et al. (13 additional authors not shown)

Abstract: The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is tho… ▽ More The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is thought to harbour a second, closer-by (<10 au) companion, was observed with the Atacama Large Millimeter/submillimeter Array as part of the ATOMIUM Large programme. In this work, the brightest CO, SiO, and HCN molecular line transitions are analysed. The continuum map shows two maxima, separated by 0.04'' (6 au). The CO data unambiguously reveal that Pi1 Gru's circumstellar environment harbours an inclined, radially outflowing, equatorial density enhancement. It contains a spiral structure at an angle of 38+/-3 deg with the line-of-sight. The HCN emission in the inner wind reveals a clockwise spiral, with a dynamical crossing time of the spiral arms consistent with a companion at a distance of 0.04'' from the asymptotic giant branch star, which is in agreement with the position of the secondary continuum peak. The inner wind dynamics imply a large acceleration region, consistent with a beta-law power of ~6. The CO emission suggests that the spiral is approximately Archimedean within 5'', beyond which this trend breaks down as the succession of the spiral arms becomes less periodic. The SiO emission at scales smaller than 0.5'' exhibits signatures of gas in rotation, which is found to fit the expected behaviour of gas in the wind-companion interaction zone. An investigation of SiO maser emission reveals what could be a stream of gas accelerating from the surface of the AGB star to the companion. Using these dynamics, we have tentatively derived an upper limit on the companion mass to be ~1.1 Msol. △ Less

Submitted 12 October, 2020; originally announced October 2020.

Comments: 19 pages, 20 figures

Journal ref: A&A 644, A61 (2020)

(Sub)stellar companions shape the winds of evolved stars

Authors: L. Decin, M. Montargès, A. M. S. Richards, C. A. Gottlieb, W. Homan, I. McDonald, I. El Mellah, T. Danilovich, S. H. J. Wallström, A. Zijlstra, A. Baudry, J. Bolte, E. Cannon, E. De Beck, F. De Ceuster, A. de Koter, J. De Ridder, S. Etoka, D. Gobrecht, M. Gray, F. Herpin, M. Jeste, E. Lagadec, P. Kervella, T. Khouri , et al. (10 additional authors not shown)

Abstract: Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their… ▽ More Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe. △ Less

Submitted 28 September, 2020; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: 19 pages main journal, 97 pages Supplementary Information

Journal ref: Science, 2000, Vol 369, Issue 6510, p. 1497

Complex beam mapping and Fourier optics analysis of a wide field Microwave Kinetic Inductance Detector camera

Authors: S. J. C. Yates, K. K. Davis, W. Jellema, J. J. A. Baselmans, A. M. Baryshev

Abstract: For astronomical instruments, accurate knowledge of the optical pointing and coupling are essential to characterize the alignment and performance of (sub-)systems prior to integration and deployment. Ideally, this requires the phase response of the optical system, which for direct (phase insensitive) detectors was not previously accessible. Here we show development of the phase sensitive complex b… ▽ More For astronomical instruments, accurate knowledge of the optical pointing and coupling are essential to characterize the alignment and performance of (sub-)systems prior to integration and deployment. Ideally, this requires the phase response of the optical system, which for direct (phase insensitive) detectors was not previously accessible. Here we show development of the phase sensitive complex beam pattern technique using a dual optical source heterodyne technique for a large field of view Microwave Kinetic Inductance Detector camera at 350 GHz. We show here how you can analyze the measured data with Fourier optics, which allows integration into a telescope model to calculate the on sky beam pattern and telescope aperture efficiency prior to deployment at a telescope. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: This is a post-peer-review, pre-copy edit version of an article published in Journal of Low Temperature Physics. The final authenticated version is available online

Journal ref: Journal of Low Temperature Physics, 2020

COBRaS: The e-MERLIN 21 cm Legacy survey of Cygnus OB2

Authors: J. C. Morford, D. M. Fenech, R. K. Prinja, R. Blomme, J. A. Yates, J. J. Drake, S. P. S. Eyres, A. M. S. Richards, I. R. Stevens, N. J. Wright, J. S. Clark, S. Dougherty, J. M. Pittard, H. Smith, J. S. Vink

Abstract: The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar… ▽ More The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2. △ Less

Submitted 20 March, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 41 pages, 12 figures, accepted in A&A

Journal ref: A&A 637, A64 (2020)

DESHIMA on ASTE: On-sky Responsivity Calibration of the Integrated Superconducting Spectrometer

Authors: Tatsuya Takekoshi, Kenichi Karatsu, Junya Suzuki, Yoichi Tamura, Tai Oshima, Akio Taniguchi, Shin'ichiro Asayama, Tom J. L. C. Bakx, Jochem J. A. Baselmans, Sjoerd Bosma, Juan Bueno, Kah Wuy Chin, Yasunori Fujii, Kazuyuki Fujita, Robert Huiting, Soh Ikarashi, Tsuyoshi Ishida, Shun Ishii, Ryohei Kawabe, Teun M. Klapwijk, Kotaro Kohno, Akira Kouchi, Nuria Llombart, Jun Maekawa, Vignesh Murugesan , et al. (14 additional authors not shown)

Abstract: We are developing an ultra-wideband spectroscopic instrument, DESHIMA (DEep Spectroscopic HIgh-redshift MApper), based on the technologies of an on-chip filter-bank and Microwave Kinetic Inductance Detector (MKID) to investigate dusty star-burst galaxies in the distant universe at millimeter and submillimeter wavelength. An on-site experiment of DESHIMA was performed using the ASTE 10-m telescope.… ▽ More We are developing an ultra-wideband spectroscopic instrument, DESHIMA (DEep Spectroscopic HIgh-redshift MApper), based on the technologies of an on-chip filter-bank and Microwave Kinetic Inductance Detector (MKID) to investigate dusty star-burst galaxies in the distant universe at millimeter and submillimeter wavelength. An on-site experiment of DESHIMA was performed using the ASTE 10-m telescope. We established a responsivity model that converts frequency responses of the MKIDs to line-of-sight brightness temperature. We estimated two parameters of the responsivity model using a set of skydip data taken under various precipitable water vapor (PWV, 0.4-3.0 mm) conditions for each MKID. The line-of-sight brightness temperature of sky is estimated using an atmospheric transmission model and the PWVs. As a result, we obtain an average temperature calibration uncertainty of $1σ=4$%, which is smaller than other photometric biases. In addition, the average forward efficiency of 0.88 in our responsivity model is consistent with the value expected from the geometrical support structure of the telescope. We also estimate line-of-sight PWVs of each skydip observation using the frequency response of MKIDs, and confirm the consistency with PWVs reported by the Atacama Large Millimeter/submillimeter Array. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: Accepted for publication in the Journal of Low Temperature Physics for LTD-18 special issue

Ozone chemistry on tidally locked M dwarf planets

Authors: Jack S. Yates, Paul I. Palmer, James Manners, Ian Boutle, Krisztian Kohary, Nathan Mayne, Luke Abraham

Abstract: We use the Met Office Unified Model to explore the potential of a tidally locked M dwarf planet, nominally Proxima Centauri b irradiated by a quiescent version of its host star, to sustain an atmospheric ozone layer. We assume a slab ocean surface layer, and an Earth-like atmosphere of nitrogen and oxygen with trace amounts of ozone and water vapour. We describe ozone chemistry using the Chapman m… ▽ More We use the Met Office Unified Model to explore the potential of a tidally locked M dwarf planet, nominally Proxima Centauri b irradiated by a quiescent version of its host star, to sustain an atmospheric ozone layer. We assume a slab ocean surface layer, and an Earth-like atmosphere of nitrogen and oxygen with trace amounts of ozone and water vapour. We describe ozone chemistry using the Chapman mechanism and the hydrogen oxide (HO$_x$, describing the sum of OH and HO$_2$) catalytic cycle. We find that Proxima Centauri radiates with sufficient UV energy to initialize the Chapman mechanism. The result is a thin but stable ozone layer that peaks at 0.75 parts per million at 25 km. The quasi-stationary distribution of atmospheric ozone is determined by photolysis driven by incoming stellar radiation and by atmospheric transport. Ozone mole fractions are smallest in the lowest 15 km of the atmosphere at the sub-stellar point and largest in the nightside gyres. Above 15 km the ozone distribution is dominated by an equatorial jet stream that circumnavigates the planet. The nightside ozone distribution is dominated by two cyclonic Rossby gyres that result in localized ozone hotspots. On the dayside the atmospheric lifetime is determined by the HO$_x$ catalytic cycle and deposition to the surface, with nightside lifetimes due to chemistry much longer than timescales associated with atmospheric transport. Surface UV values peak at the substellar point with values of 0.01 W/m$^2$, shielded by the overlying atmospheric ozone layer but more importantly by water vapour clouds. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: Accepted for publication in MNRAS. 25 pages, 17 figures

Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling

Authors: Frederik De Ceuster, Ward Homan, Jeremy Yates, Leen Decin, Peter Boyle, James Hetherington

Abstract: Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of ray… ▽ More Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster. △ Less

Submitted 18 December, 2019; originally announced December 2019.

Comments: 16 pages, accepted in MNRAS

Developing a self-consistent AGB wind model: II. Non-classical, non-equilibrium polymer nucleation in a chemical mixture

Authors: Jels Boulangier, David Gobrecht, Leen Decin, Alex de Koter, Jeremy Yates

Abstract: Unravelling the composition and characteristics of gas and dust lost by asymptotic giant branch (AGB) stars is important as these stars play a vital role in the chemical life cycle of galaxies. The general hypothesis of their mass loss mechanism is a combination of stellar pulsations and radiative pressure on dust grains. However, current models simplify dust formation, which starts as a microscop… ▽ More Unravelling the composition and characteristics of gas and dust lost by asymptotic giant branch (AGB) stars is important as these stars play a vital role in the chemical life cycle of galaxies. The general hypothesis of their mass loss mechanism is a combination of stellar pulsations and radiative pressure on dust grains. However, current models simplify dust formation, which starts as a microscopic phase transition called nucleation. Various nucleation theories exist, yet all assume chemical equilibrium, growth restricted by monomers, and commonly use macroscopic properties for a microscopic process. Such simplifications for initial dust formation can have large repercussions on the type, amount, and formation time of dust. By abandoning equilibrium assumptions, discarding growth restrictions, and using quantum mechanical properties, we have constructed and investigated an improved nucleation theory in AGB wind conditions for four dust candidates, TiO$_2$, MgO, SiO and Al$_2$O$_3$. This paper reports the viability of these candidates as first dust precursors and reveals implications of simplified nucleation theories. Monomer restricted growth underpredicts large clusters at low temperatures and overpredicts formation times. Assuming the candidates are present, Al$_2$O$_3$ is the favoured precursor due to its rapid growth at the highest considered temperatures. However, when considering an initially atomic chemical mixture, only TiO$_2$-clusters form. Still, we believe Al$_2$O$_3$ to be the prime candidate due to substantial physical evidence in presolar grains, observations of dust around AGB stars at high temperatures, and its ability to form at high temperatures and expect the missing link to be insufficient quantitative data of Al-reactions. △ Less

Submitted 26 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS. 19 pages (68 incl. appendix)

Incorporating Astrochemistry into Molecular Line Modelling via Emulation

Authors: Damien de Mijolla, Serena Viti, Jonathan Holdship, Ioanna Manolopoulou, Jeremy Yates

Abstract: In studies of the interstellar medium in galaxies, radiative transfer models of molecular emission are useful for relating molecular line observations back to the physical conditions of the gas they trace. However, doing this requires solving a highly degenerate inverse problem. In order to alleviate these degeneracies, the abundances derived from astrochemical models can be converted into column… ▽ More In studies of the interstellar medium in galaxies, radiative transfer models of molecular emission are useful for relating molecular line observations back to the physical conditions of the gas they trace. However, doing this requires solving a highly degenerate inverse problem. In order to alleviate these degeneracies, the abundances derived from astrochemical models can be converted into column densities and fed into radiative transfer models. This enforces that the molecular gas composition used by the radiative transfer models be chemically realistic. However, because of the complexity and long running time of astrochemical models, it can be difficult to incorporate chemical models into the radiative transfer framework. In this paper, we introduce a statistical emulator of the UCLCHEM astrochemical model, built using neural networks. We then illustrate, through examples of parameter estimations, how such an emulator can be applied on real and synthetic observations. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: Accepted by A&A; Emulator soon to be publicly released

Journal ref: A&A 630, A117 (2019)

First light demonstration of the integrated superconducting spectrometer

Authors: Akira Endo, Kenichi Karatsu, Yoichi Tamura, Tai Oshima, Akio Taniguchi, Tatsuya Takekoshi, Shin'ichiro Asayama, Tom J. L. C. Bakx, Sjoerd Bosma, Juan Bueno, Kah Wuy Chin, Yasunori Fujii, Kazuyuki Fujita, Robert Huiting, Soh Ikarashi, Tsuyoshi Ishida, Shun Ishii, Ryohei Kawabe, Teun M. Klapwijk, Kotaro Kohno, Akira Kouchi, Nuria Llombart, Jun Maekawa, Vignesh Murugesan, Shunichi Nakatsubo , et al. (14 additional authors not shown)

Abstract: Ultra-wideband 3D imaging spectrometry in the millimeter-submillimeter (mm-submm) band is an essential tool for uncovering the dust-enshrouded portion of the cosmic history of star formation and galaxy evolution. However, it is challenging to scale up conventional coherent heterodyne receivers or free-space diffraction techniques to sufficient bandwidths ($\geq$1 octave) and numbers of spatial pix… ▽ More Ultra-wideband 3D imaging spectrometry in the millimeter-submillimeter (mm-submm) band is an essential tool for uncovering the dust-enshrouded portion of the cosmic history of star formation and galaxy evolution. However, it is challenging to scale up conventional coherent heterodyne receivers or free-space diffraction techniques to sufficient bandwidths ($\geq$1 octave) and numbers of spatial pixels (>$10^2$). Here we present the design and first astronomical spectra of an intrinsically scalable, integrated superconducting spectrometer, which covers 332-377 GHz with a spectral resolution of $F/ΔF \sim 380$. It combines the multiplexing advantage of microwave kinetic inductance detectors (MKIDs) with planar superconducting filters for dispersing the signal in a single, small superconducting integrated circuit. We demonstrate the two key applications for an instrument of this type: as an efficient redshift machine, and as a fast multi-line spectral mapper of extended areas. The line detection sensitivity is in excellent agreement with the instrument design and laboratory performance, reaching the atmospheric foreground photon noise limit on sky. The design can be scaled to bandwidths in excess of an octave, spectral resolution up to a few thousand and frequencies up to $\sim$1.1 THz. The miniature chip footprint of a few $\mathrm{cm^2}$ allows for compact multi-pixel spectral imagers, which would enable spectroscopic direct imaging and large volume spectroscopic surveys that are several orders of magnitude faster than what is currently possible. △ Less

Submitted 18 September, 2019; v1 submitted 24 June, 2019; originally announced June 2019.

Comments: Published in Nature Astronomy. SharedIt Link to the full published paper: https://rdcu.be/bM2FN

Wideband on-chip terahertz spectrometer based on a superconducting filterbank

Authors: Akira Endo, Kenichi Karatsu, Alejandro Pascual Laguna, Behnam Mirzaei, Robert Huiting, David J. Thoen, Vignesh Murugesan, Stephen J. C. Yates, Juan Bueno, Nuri van Marrewijk, Sjoerd Bosma, Ozan Yurduseven, Nuria Llombart, Junya Suzuki, Masato Naruse, Pieter J. de Visser, Paul P. van der Werf, Teun M. Klapwijk, Jochem J. A. Baselmans

Abstract: Terahertz spectrometers with a wide instantaneous frequency coverage for passive remote sensing are enormously attractive for many terahertz applications, such as astronomy, atmospheric science and security. Here we demonstrate a wide-band terahertz spectrometer based on a single superconducting chip. The chip consists of an antenna coupled to a transmission line filterbank, with a microwave kinet… ▽ More Terahertz spectrometers with a wide instantaneous frequency coverage for passive remote sensing are enormously attractive for many terahertz applications, such as astronomy, atmospheric science and security. Here we demonstrate a wide-band terahertz spectrometer based on a single superconducting chip. The chip consists of an antenna coupled to a transmission line filterbank, with a microwave kinetic inductance detector behind each filter. Using frequency division multiplexing, all detectors are read-out simultaneously creating a wide-band spectrometer with an instantaneous bandwidth of 45 GHz centered around 350 GHz. The spectrometer has a spectral resolution of $F/ΔF=380$ and reaches photon-noise limited sensitivity. We discuss the chip design and fabrication, as well as the system integration and testing. We confirm full system operation by the detection of an emission line spectrum of methanol gas. The proposed concept allows for spectroscopic radiation detection over large bandwidths and resolutions up to $F/ΔF\sim1000$, all using a chip area of a few $\mathrm{cm^2}$. This will allow the construction of medium resolution imaging spectrometers with unprecedented speed and sensitivity. △ Less

Submitted 21 June, 2019; v1 submitted 21 January, 2019; originally announced January 2019.

Comments: 21 pages, 7 figures

Journal ref: Journal of Astronomical Telescopes, Instruments, and Systems, 5(3), 035004 (2019). (Open Access)

Analysis techniques for complex-field radiation pattern measurements

Authors: Kristina K. Davis, Willem Jellema, Stephen J. C. Yates, Christopher E. Groppi, Jochem J. A. Baselmans, Andrey M. Baryshev

Abstract: Complex field measurements are increasingly becoming the standard for state-of-the-art astronomical instrumentation. Complex field measurements have been used to characterize a suite of ground, airborne, and space-based heterodyne receiver missions [1], [2], [3], [4], [5], [6], and a description of how to acquire coherent field maps for direct detector arrays was demonstrated in Davis et. al. 2017… ▽ More Complex field measurements are increasingly becoming the standard for state-of-the-art astronomical instrumentation. Complex field measurements have been used to characterize a suite of ground, airborne, and space-based heterodyne receiver missions [1], [2], [3], [4], [5], [6], and a description of how to acquire coherent field maps for direct detector arrays was demonstrated in Davis et. al. 2017. This technique has the ability to determine both amplitude and phase radiation patterns from individual pixels on an array. Phase information helps to better characterize the optical performance of the array (as compared to total power radiation patterns) by constraining the fit in an additional plane [4]. Here we discuss the mathematical framework used in an analysis pipeline developed to process complex field radiation pattern measurements. This routine determines and compensates misalignments of the instrument and scanning system. We begin with an overview of Gaussian beam formalism and how it relates to complex field pattern measurements. Next we discuss a scan strategy using an offset in z along the optical axis that allows first-order optical standing waves between the scanned source and optical system to be removed in post-processing. Also discussed is a method by which the co- and cross-polarization fields can be extracted individually for each pixel by rotating the two orthogonal measurement planes until the signal is the co-polarization map is maximized (and the signal in the cross-polarization field is minimized). We detail a minimization function that can fit measurement data to an arbitrary beam shape model. We conclude by discussing the angular plane wave spectral (APWS) method for beam propagation, including the near-field to far-field transformation. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Journal ref: Proceedings Volume 10708, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy IX; 107082R (2018)

Bayesian Inference of Reaction Rates in Icy Mantles

Authors: Jonathan Holdship, Niall Jeffrey, Antonios Makrymallis, Serena Viti, Jeremy Yates

Abstract: Grain surface chemistry and its treatment in gas-grain chemical models is an area of large uncertainty. Whilst laboratory experiments are making progress, there is still much that is unknown about grain surface chemistry. Further, the results and parameters produced by experiment are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. It is poss… ▽ More Grain surface chemistry and its treatment in gas-grain chemical models is an area of large uncertainty. Whilst laboratory experiments are making progress, there is still much that is unknown about grain surface chemistry. Further, the results and parameters produced by experiment are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. It is possible that statistical methods can reduce the uncertainty in grain surface chemical networks. In this work, a simple model of grain surface chemistry in a molecular cloud is developed and a Bayesian inference of the reactions rates is performed through MCMC sampling. Using observational data of the solid state abundances of major chemical species in molecular clouds, the posterior distributions for the rates of seven reactions producing CO, CO$_2$, CH$_3$OH and H$_2$O are calculated, in a form that is suitable for rate equation models. This represents a vital first step in the development of a method to infer reaction rates from observations of chemical abundances in astrophysical environments. △ Less

Submitted 25 September, 2018; v1 submitted 17 September, 2018; originally announced September 2018.

Comments: Accepted for publication in ApJ

Eliminating stray radiation inside large area imaging arrays

Authors: Stephen J. C. Yates, Simon Doyle, Peter Barry, Andrey M. Baryshev, Juan Bueno, Lorenza Ferrari, Nuria Llombart, Vignesh Murugesan, David J. Thoen, Ozan Yurduseven, Jochem J. A. Baselmans

Abstract: With increasing array size, it is increasingly important to control stray radiation inside the detector chips themselves. We demonstrate this effect with focal plane arrays of absorber coupled Lumped Element microwave Kinetic Inductance Detectors (LEKIDs) and lens-antenna coupled distributed quarter wavelength Microwave Kinetic Inductance Detectors (MKIDs). In these arrays the response from a poin… ▽ More With increasing array size, it is increasingly important to control stray radiation inside the detector chips themselves. We demonstrate this effect with focal plane arrays of absorber coupled Lumped Element microwave Kinetic Inductance Detectors (LEKIDs) and lens-antenna coupled distributed quarter wavelength Microwave Kinetic Inductance Detectors (MKIDs). In these arrays the response from a point source at the pixel position is at a similar level to the stray response integrated over the entire chip area. For the antenna coupled arrays, we show that this effect can be suppressed by incorporating an on-chip stray light absorber. A similar method should be possible with the LEKID array, especially when they are lens coupled. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: arXiv admin note: substantial text overlap with arXiv:1707.02142

Journal ref: Proc. SPIE 10708, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy IX, 107081X (9 July 2018)

Solar System Ice Giants: Exoplanets in our Backyard

Authors: Abigail Rymer, Kathleen Mandt, Dana Hurley, Carey Lisse, Noam Izenberg, H. Todd Smith, Joseph Westlake, Emma Bunce, Christopher Arridge, Adam Masters, Mark Hofstadter, Amy Simon, Pontus Brandt, George Clark, Ian Cohen, Robert Allen, Sarah Vine, Kenneth Hansen, George Hospodarsky, William Kurth, Paul Romani, Laurent Lamy, Philippe Zarka, Hao Cao, Carol Paty , et al. (88 additional authors not shown)

Abstract: Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of s… ▽ More Future remote sensing of exoplanets will be enhanced by a thorough investigation of our solar system Ice Giants (Neptune-size planets). What can the configuration of the magnetic field tell us (remotely) about the interior, and what implications does that field have for the structure of the magnetosphere; energy input into the atmosphere, and surface geophysics (for example surface weathering of satellites that might harbour sub-surface oceans). How can monitoring of auroral emission help inform future remote observations of emission from exoplanets? Our Solar System provides the only laboratory in which we can perform in-situ experiments to understand exoplanet formation, dynamos, systems and magnetospheres. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: Exoplanet Science Strategy White Paper, submitted to the National Academies of Sciences, Engineering and Medicine, Space Studies Board, 9 March 2018

Journal ref: 2020n3i176

The JCMT Gould Belt Survey: A First Look at the Auriga-California Molecular Cloud with SCUBA-2

Authors: H. Broekhoven-Fiene, B. C. Matthews, P. Harvey, H. Kirk, M. Chen, M. J. Currie, K. Pattle, J. Lane, J. Buckle, J. Di Francesco, E. Drabek-Maunder, D. Johnstone, D. S. Berry, M. Fich, J. Hatchell, T. Jenness, J. C. Mottram, D. Nutter, J. E. Pineda, C. Quinn, C. Salji, S. Tisi, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien , et al. (35 additional authors not shown)

Abstract: We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on… ▽ More We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency. △ Less

Submitted 24 January, 2018; originally announced January 2018.

Comments: 39 pages (54 including the Appendix), 9 figures, 3 tables, in press in ApJ

Journal ref: The Astrophysical Journal, Volume 852, Issue 2, article id. 73, 20 pp. (2018)

Surface wave control for large arrays of microwave kinetic inductance detectors

Authors: Stephen J. C. Yates, Andrey M. Baryshev, Ozan Yurduseven, Juan Bueno, Kristina K. Davis, Lorenza Ferrari, Willem Jellema, Nuria Llombart, Vignesh Murugesan, David J. Thoen, Jochem J. A. Baselmans

Abstract: Large ultra-sensitive detector arrays are needed for present and future observatories for far infra-red, submillimeter wave (THz), and millimeter wave astronomy. With increasing array size, it is increasingly important to control stray radiation inside the detector chips themselves, the surface wave. We demonstrate this effect with focal plane arrays of 880 lens-antenna coupled Microwave Kinetic I… ▽ More Large ultra-sensitive detector arrays are needed for present and future observatories for far infra-red, submillimeter wave (THz), and millimeter wave astronomy. With increasing array size, it is increasingly important to control stray radiation inside the detector chips themselves, the surface wave. We demonstrate this effect with focal plane arrays of 880 lens-antenna coupled Microwave Kinetic Inductance Detectors (MKIDs). Presented here are near field measurements of the MKID optical response versus the position on the array of a reimaged optical source. We demonstrate that the optical response of a detector in these arrays saturates off-pixel at the $\sim-30$ dB level compared to the peak pixel response. The result is that the power detected from a point source at the pixel position is almost identical to the stray response integrated over the chip area. With such a contribution, it would be impossible to measure extended sources, while the point source sensitivity is degraded due to an increase of the stray loading. However, we show that by incorporating an on-chip stray light absorber, the surface wave contribution is reduced by a factor $>$10. With the on-chip stray light absorber the point source response is close to simulations down to the $\sim-35$ dB level, the simulation based on an ideal Gaussian illumination of the optics. In addition, as a crosscheck we show that the extended source response of a single pixel in the array with the absorbing grid is in agreement with the integral of the point source measurements. △ Less

Submitted 7 July, 2017; originally announced July 2017.

Comments: accepted for publication in IEEE Transactions on Terahertz Science and Technology

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

The JCMT Gould Belt Survey: A First Look at IC 5146

Authors: D. Johnstone, S. Ciccone, H. Kirk, S. Mairs, J. Buckle, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, J. Hatchell, T. Jenness, J. C. Mottram, K. Pattle, S. Tisi J. Di Francesco, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, D. Bresnahan, H. Butner, M. Chen, A. Chrysostomou, S. Coude, C. J. Davis, E. Drabek-Maunder, A. Duarte-Cabral, M. Fich , et al. (31 additional authors not shown)

Abstract: We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual… ▽ More We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: 30 pages, 11 figures, 8 tables, accepted by ApJ

Atmospheric Habitable Zones in Cool Y Dwarf Atmospheres

Authors: Jack S. Yates, Paul I. Palmer, Beth Biller, Charles S. Cockell

Abstract: We use a simple organism lifecycle model to explore the viability of an atmospheric habitable zone (AHZ), with temperatures that could support Earth-centric life, which sits above an environment that does not support life. To illustrate our model we use a cool Y dwarf atmosphere, such as $\mathrm{WISE~J}085510.83-0714442.5$ whose $4.5-5.2$ micron spectrum shows absorption features consistent with… ▽ More We use a simple organism lifecycle model to explore the viability of an atmospheric habitable zone (AHZ), with temperatures that could support Earth-centric life, which sits above an environment that does not support life. To illustrate our model we use a cool Y dwarf atmosphere, such as $\mathrm{WISE~J}085510.83-0714442.5$ whose $4.5-5.2$ micron spectrum shows absorption features consistent with water vapour and clouds. We allow organisms to adapt to their atmospheric environment (described by temperature, convection, and gravity) by adopting different growth strategies that maximize their chance of survival and proliferation. We assume a constant upward vertical velocity through the AHZ. We found that the organism growth strategy is most sensitive to the magnitude of the atmospheric convection. Stronger convection supports the evolution of more massive organisms. For a purely radiative environment we find that evolved organisms have a mass that is an order of magnitude smaller than terrestrial microbes, thereby defining a dynamical constraint on the dimensions of life that an AHZ can support. Based on a previously defined statistical approach we infer that there are of order $10^9$ cool Y brown dwarfs in the Milky Way, and likely a few tens of these objects are within ten parsecs from Earth. Our work also has implications for exploring life in the atmospheres of temperate gas giants. Consideration of the habitable volumes in planetary atmospheres significantly increases the volume of habitable space in the galaxy. △ Less

Submitted 20 December, 2016; v1 submitted 28 November, 2016; originally announced November 2016.

Comments: Revised version after comments received. Re-accepted to ApJ on 20/12/2016

The JCMT Gould Belt Survey: First results from SCUBA-2 observations of the Cepheus Flare Region

Authors: Kate Pattle, Derek Ward-Thompson, Jason M. Kirk, James Di Francesco, Helen Kirk, Joseph C. Mottram, Jared Keown, Jane Buckle, Sylvie F. Beaulieu, David S. Berry, Hannah Broekhoven-Fiene, Malcolm J. Currie, Michel Fich, Jenny Hatchell, Tim Jenness, Doug Johnstone, David Nutter, Jaime E. Pineda, Ciera Quinn, Carl Salji, Sam Tisi, Samantha Walker-Smith, Michiel R. Hogerheijde, Pierre Bastien, David Bresnahan , et al. (35 additional authors not shown)

Abstract: We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Sur… ▽ More We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields. △ Less

Submitted 12 October, 2016; originally announced October 2016.

Comments: Accepted for publication by MNRAS, 29 pages, 16 figures, 8 tables

A kilo-pixel imaging system for future space based far-infrared observatories using microwave kinetic inductance detectors

Authors: J. J. A. Baselmans, J. Bueno, S. J. C. Yates, O. Yurduseven, N. Llombart, K. Karatsu, A. M. Baryshev, L. Ferrari, A. Endo, D. J. Thoen, P. J. de Visser, R. M. J. Janssen, V. Murugesan, E. F. C. Driessen, G. Coiffard, J. Martin-Pintado, P. Hargrave, M. Griffin

Abstract: Future astrophysics and cosmic microwave background space missions operating in the far-infrared to millimetre part of the spectrum will require very large arrays of ultra-sensitive detectors in combination with high multiplexing factors and efficient low-noise and low-power readout systems. We have developed a demonstrator system suitable for such applications. The system combines a 961 pixel ima… ▽ More Future astrophysics and cosmic microwave background space missions operating in the far-infrared to millimetre part of the spectrum will require very large arrays of ultra-sensitive detectors in combination with high multiplexing factors and efficient low-noise and low-power readout systems. We have developed a demonstrator system suitable for such applications. The system combines a 961 pixel imaging array based upon Microwave Kinetic Inductance Detectors (MKIDs) with a readout system capable of reading out all pixels simultaneously with only one readout cable pair and a single cryogenic amplifier. We evaluate, in a representative environment, the system performance in terms of sensitivity, dynamic range, optical efficiency, cosmic ray rejection, pixel-pixel crosstalk and overall yield at at an observation centre frequency of 850 GHz and 20% fractional bandwidth. The overall system has an excellent sensitivity, with an average detector sensitivity NEPdet=3x10^-19 W/rt(Hz) measured using a thermal calibration source. At a loading power per pixel of 50fW we demonstrate white, photon noise limited detector noise down to 300 mHz. The dynamic range would allow the detection of 1 Jy bright sources within the field of view without tuning the readout of the detectors. The expected dead time due to cosmic ray interactions, when operated in an L2 or a similar far-Earth orbit, is found to be <4%. Additionally, the achieved pixel yield is 83% and the crosstalk between the pixels is <-30dB. This demonstrates that MKID technology can provide multiplexing ratios on the order of a 1000 with state-of-the-art single pixel performance, and that the technology is now mature enough to be considered for future space based observatories and experiments. △ Less

Submitted 1 March, 2017; v1 submitted 7 September, 2016; originally announced September 2016.

Comments: Resubmitted manuscript after review A&A (2016) 16 pages, 10 figures

Journal ref: A&A 601, A89 (2017)

Measurements and analysis of optical crosstalk in a microwave kinetic inductance detector array

Authors: L. Bisigello, S. J. C. Yates, L. Ferrari, J. J. A. Baselmans, A. M. Baryshev

Abstract: The main advantage of Microwave Kinetic Inductance Detector arrays (MKID) is their multiplexing capability, which allows for building cameras with a large number of pixels and good sensitivity, particularly suitable to perform large blank galaxy surveys. However, to have as many pixels as possible it is necessary to arrange detectors close in readout frequency. Consequently KIDs overlap in frequen… ▽ More The main advantage of Microwave Kinetic Inductance Detector arrays (MKID) is their multiplexing capability, which allows for building cameras with a large number of pixels and good sensitivity, particularly suitable to perform large blank galaxy surveys. However, to have as many pixels as possible it is necessary to arrange detectors close in readout frequency. Consequently KIDs overlap in frequency and are coupled to each other producing crosstalk. Because crosstalk can be only minimised by improving the array design, in this work we aim to correct for this effect a posteriori. We analysed a MKID array consisting of 880 KIDs with readout frequencies at 4-8 GHz. We measured the beam patterns for every detector in the array and described the response of each detector by using a two-dimensional Gaussian fit. Then, we identified detectors affected by crosstalk above -30 dB level from the maximum and removed the signal of the crosstalking detectors. Moreover, we modelled the crosstalk level for each KID as a function of the readout frequency separation starting from the assumption that the transmission of a KID is a Lorenztian function in power. We were able to describe the general crosstalk level of the array and the crosstalk of each KID within 5 dB, so enabling the design of future arrays with the crosstalk as a design criterion. In this work, we demonstrate that it is possible to process MKID images a posteriori to decrease the crosstalk effect, subtracting the response of each coupled KID from the original map. △ Less

Submitted 1 September, 2016; originally announced September 2016.

Comments: 8 pages, 8 figures, Proc. SPIE 9914, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII, 99143L (July 20, 2016)

The mass-loss rates of red supergiants at low metallicity: Detection of rotational CO emission from two red supergiants in the Large Magellanic Cloud

Authors: Mikako Matsuura, B. Sargent, Bruce Swinyard, Jeremy Yates, P. Royer, M. J. Barlow, Martha Boyer, L. Decin, Theo Khouri, Margaret Meixner, Jacco Th. van Loon, Paul M. Woods

Abstract: Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars i… ▽ More Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars in the LMC. Their CO line intensities are as strong as those of the Galactic RSG, VY CMa. Modelling the CO lines and the spectral energy distribution results in an estimated mass-loss rate for IRAS 05280-6910 of 3x10^-4 Msun per yr. The model assumes a gas-to-dust ratio and a CO-to-H2 abundance ratio is estimated from the Galactic values scaled by the LMC metallicity ([Fe/H]~-0.3), i.e., that the CO-to-dust ratio is constant for Galactic and LMC metallicities within the uncertainties of the model. The key factor determining the CO line intensities and the mass-loss rate found to be the stellar luminosity. △ Less

Submitted 4 August, 2016; originally announced August 2016.

Comments: 11 pages; accepted by MNRAS

e-MERLIN 21cm constraints on the mass-loss rates of OB stars in Cyg OB2

Authors: Jack Morford, Danielle Fenech, Raman Prinja, Ronnie Blomme, Jeremy Yates

Abstract: We present e-MERLIN 21 cm (L-band) observations of single luminous OB stars in the Cygnus OB2 association, from the COBRaS Legacy programme. The radio observations potentially offer the most straightforward, least model-dependent, determinations of mass-loss rates, and can be used to help resolve current discrepancies in mass-loss rates via clumped and structured hot star winds. We report here tha… ▽ More We present e-MERLIN 21 cm (L-band) observations of single luminous OB stars in the Cygnus OB2 association, from the COBRaS Legacy programme. The radio observations potentially offer the most straightforward, least model-dependent, determinations of mass-loss rates, and can be used to help resolve current discrepancies in mass-loss rates via clumped and structured hot star winds. We report here that the 21 cm flux densities of O3 to O6 supergiant and giant stars are less than ~ 70 microJy. These fluxes may be translated to `smooth' wind mass-loss upper limits of ~ 4.4 - 4.8 x 10^(-6) M_sol/yr for O3 supergiants and < 2.9 x 10^(-6) M_sol/yr for B0 to B1 supergiants. The first ever resolved 21 cm detections of the hypergiant (and LBV candidate) Cyg OB2 #12 are discussed; for multiple observations separated by 14 days, we detect a ~ 69% increase in its flux density. Our constraints on the upper limits for the mass-loss rates of evolved OB stars in Cyg OB2 support the model that the inner wind region close to the stellar surface (where Halpha forms) is more clumped than the very extended geometric region sampled by our radio observations. △ Less

Submitted 29 July, 2016; originally announced July 2016.

The JCMT Gould Belt Survey: A First Look at Southern Orion A with SCUBA-2

Authors: Steve Mairs, D. Johnstone, H. Kirk, J. Buckle, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, S. Graves, J. Hatchell, T. Jenness, J. C. Mottram, D. Nutter, K. Pattle, J. E. Pineda, C. Salji, J. Di Francesco, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, D. Bresnahan, H. Butner, M. Chen, A. Chrysostomou, S. Coudé , et al. (30 additional authors not shown)

Abstract: We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($δ\leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $μ$m continuum maps obtained… ▽ More We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($δ\leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $μ$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $μ$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: 31 Pages, 19 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The JCMT Gould Belt Survey: Evidence for Dust Grain Evolution in Perseus Star-forming Clumps

Authors: Michael Chun-Yuan Chen, J. Di Francesco, D. Johnstone, S. Sadavoy, J. Hatchell, J. C. Mottram, H. Kirk, J. Buckle, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, T. Jenness, D. Nutter, K. Pattle, J. E. Pineda, C. Quinn, C. Salji, S. Tisi, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, D. Bresnahan, H. Butner, A. Chrysostomou , et al. (34 additional authors not shown)

Abstract: The dust emissivity spectral index, $β$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $β$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $β$, dust temperature, and optical depth maps of… ▽ More The dust emissivity spectral index, $β$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $β$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $β$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $μ$m, 250 $μ$m, 350 $μ$m, 500 $μ$m, and 850 $μ$m bands. Most of the derived $β$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $β$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $β$ variations within individual clumps and find low $β$ regions correlate with local temperature peaks, hinting at the possible origins of low $β$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape. △ Less

Submitted 19 May, 2016; originally announced May 2016.

Comments: Accepted to ApJ, 15 figures, 3 tables

The JCMT Gould Belt Survey: Evidence for radiative heating and contamination in the W40 complex

Authors: D. Rumble, J. Hatchell, K. Pattle, H. Kirk, T. Wilson, J. Buckle, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, T. Jenness, D. Johnstone, J. C. Mottram, D. Nutter, J. E. Pineda, C. Quinn, C. Salji, S. Tisi, S. Walker-Smith, J. Di Francesco, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, D. Bresnahan, H. Butner , et al. (33 additional authors not shown)

Abstract: We present SCUBA-2 450μm and 850μm observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, β = 1.8, and a beam convolution kernel to achieve a co… ▽ More We present SCUBA-2 450μm and 850μm observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, β = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 μm images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850μm band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 μm or 12 per cent at 850 μm in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west. △ Less

Submitted 16 May, 2016; originally announced May 2016.

Comments: 27 pages, 25 figures, 7 tables, 3 online catalogues

Calibration scheme for large Kinetic Inductance Detector Arrays based on Readout Frequency Response

Authors: L. Bisigello, S. J. C. Yates, V. Murugesan, J. J. A. Baselmans, A. M. Baryshev

Abstract: Microwave kinetic inductance detector (MKID) provides a way to build large ground based sub-mm instruments such as NIKA and A-MKID. For such instruments, therefore, it is important to understand and characterize the response to ensure good linearity and calibration over wide dynamic range. We propose to use the MKID readout frequency response to determine the MKID responsivity to an input optical… ▽ More Microwave kinetic inductance detector (MKID) provides a way to build large ground based sub-mm instruments such as NIKA and A-MKID. For such instruments, therefore, it is important to understand and characterize the response to ensure good linearity and calibration over wide dynamic range. We propose to use the MKID readout frequency response to determine the MKID responsivity to an input optical source power. A signal can be measured in a KID as a change in the phase of the readout signal with respect to the KID resonant circle. Fundamentally, this phase change is due to a shift in the KID resonance frequency, in turn due to a radiation induced change in the quasiparticle number in the superconducting resonator. We show that shift in resonant frequency can be determined from the phase shift by using KID phase versus frequency dependence using a previously measured resonant frequency. Working in this calculated resonant frequency, we gain near linearity and constant calibration to a constant optical signal applied in a wide range of operating points on the resonance and readout powers. This calibration method has three particular advantages: first, it is fast enough to be used to calibrate large arrays, with pixel counts in the thousand of pixels; second, it is based on data that are already necessary to determine KID positions; third, it can be done without applying any optical source in front of the array. △ Less

Submitted 1 February, 2016; originally announced February 2016.

Comments: Accepted to Journal of Low Temperature Physics LTD16 Special Issue, Low Temperature Detector 16 Conference Proceedings,manuscript number: #JLTP-D-15-00356R1, 6 pages, 5 figures

Journal ref: Journal of Low Temperature Physics, 2006, 1-6

The JCMT Gould Belt Survey: A First Look at Dense Cores in Orion B

Authors: H. Kirk, J. Di Francesco, D. Johnstone, A. Duarte-Cabral, S. Sadavoy, J. Hatchell, J. C. Mottram, J. Buckle, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, T. Jenness, D. Nutter, K. Pattle, J. E. Pineda, C. Quinn, C. Salji, S. Tisi, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, D. Bresnahan, H. Butner, M. Chen , et al. (32 additional authors not shown)

Abstract: We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes,… ▽ More We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars. △ Less

Submitted 2 December, 2015; originally announced December 2015.

Comments: Accepted for publication in ApJ; 28 pages, 19 figures. Data associated with the paper, including the full table 4, can be found at https://doi.org/10.11570/16.0003

The Physics of Water Masers observable with ALMA and SOFIA: Model Predictions for Evolved Stars

Authors: M. D. Gray, A. Baudry, A. M. S. Richards, E. M. L. Humphreys, A. M. Sobolev, J. A. Yates

Abstract: We present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space i… ▽ More We present the results of models that were designed to study all possible water maser transitions in the frequency range 0-1.91THz, with particular emphasis on maser transitions that may be generated in evolved-star envelopes and observed with the ALMA and SOFIA telescopes. We used tens of thousands of radiative transfer models of both spin species of H2O, spanning a considerable parameter space in number density, kinetic temperature and dust temperature. Results, in the form of maser optical depths, have been summarized in a master table, Table 6. Maser transitions identified in these models were grouped according to loci of inverted regions in the density/kinetic temperature plane, a property clearly related to the dominant mode of pumping. A more detailed study of the effect of dust temperature on maser optical depth enabled us to divide the maser transitions into three groups: those with both collisional and radiative pumping schemes (22,96,209,321,325,395,941 and 1486\,GHz), a much larger set that are predominantly radiatively pumped, and another large group with a predominantly collisional pump. The effect of accelerative and decelerative velocity shifts of up to 5km/s was found to be generally modest, with the primary effect of reducing computed maser optical depths. More subtle asymmetric effects, dependent on line overlap, include maximum gains offset from zero shift by >1km/s, but these effects were predominantly found under conditions of weak amplification. These models will allow astronomers to use multi-transition water maser observations to constrain physical conditions down to the size of individual masing clouds (size of a few astronomical units). △ Less

Submitted 21 October, 2015; originally announced October 2015.

TORUS-3DPDR: A self-consistent code treating three-dimensional photoionization and photodissociation regions

Authors: T. G. Bisbas, T. J. Haworth, M. J. Barlow, S. Viti, T. J. Harries, T. Bell, J. A. Yates

Abstract: The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium eith… ▽ More The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations. △ Less

Submitted 16 September, 2015; originally announced September 2015.

Comments: 17 pages, 10 Figures, accepted for publication in MNRAS. Comments welcome

ALMA reveals sunburn: CO dissociation around AGB stars in the globular cluster 47 Tucanae

Authors: Iain McDonald, Albert A. Zijlstra, Eric Lagadec, Gregory C. Sloan, Martha L. Boyer, Mikako Matsuura, Rowan J. Smith, Christina L. Smith, Jeremy A. Yates, Jacco Th. van Loon, Olivia C. Jones, Sofia Ramstedt, Adam Avison, Kay Justtanont, Hans Olofsson, Joris A. D. L. Blommaert, Steven R. Goldman, Martin A. T. Groenewegen

Abstract: ALMA observations show a non-detection of carbon monoxide around the four most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tucanae. Stellar evolution models and star counts show that the mass-loss rates from these stars should be ~1.2-3.5 x 10^-7 solar masses per year. We would naively expect such stars to be detectable at this distance (4.5 kpc). By modelling the ultra… ▽ More ALMA observations show a non-detection of carbon monoxide around the four most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tucanae. Stellar evolution models and star counts show that the mass-loss rates from these stars should be ~1.2-3.5 x 10^-7 solar masses per year. We would naively expect such stars to be detectable at this distance (4.5 kpc). By modelling the ultraviolet radiation field from post-AGB stars and white dwarfs in 47 Tuc, we conclude CO should be dissociated abnormally close to the stars. We estimate that the CO envelopes will be truncated at a few hundred stellar radii from their host stars and that the line intensities are about two orders of magnitude below our current detection limits. The truncation of CO envelopes should be important for AGB stars in dense clusters. Observing the CO (3-2) and higher transitions and targeting stars far from the centres of clusters should result in the detections needed to measure the outflow velocities from these stars. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: 14 pages, 4 figures, accepted MNRAS

The non-equilibrium response of a superconductor to pair-breaking radiation measured over a broad frequency band

Authors: P. J. de Visser, S. J. C. Yates, T. Guruswamy, D. J. Goldie, S. Withington, A. Neto, N. Llombart, A. M. Baryshev, T. M. Klapwijk, J. J. A. Baselmans

Abstract: We have measured the absorption of terahertz radiation in a BCS superconductor over a broad range of frequencies from 200 GHz to 1.1 THz, using a broadband antenna-lens system and a tantalum microwave resonator. From low frequencies, the response of the resonator rises rapidly to a maximum at the gap edge of the superconductor. From there on the response drops to half the maximum response at twice… ▽ More We have measured the absorption of terahertz radiation in a BCS superconductor over a broad range of frequencies from 200 GHz to 1.1 THz, using a broadband antenna-lens system and a tantalum microwave resonator. From low frequencies, the response of the resonator rises rapidly to a maximum at the gap edge of the superconductor. From there on the response drops to half the maximum response at twice the pair-breaking energy. At higher frequencies, the response rises again due to trapping of pair-breaking phonons in the superconductor. In practice this is the first measurement of the frequency dependence of the quasiparticle creation efficiency due to pair-breaking in a superconductor. The efficiency, calculated from the different non-equilibrium quasiparticle distribution functions at each frequency, is in agreement with the measurements. △ Less

Submitted 22 May, 2015; originally announced May 2015.

Journal ref: Applied Physics Letters 106, 252602 (2015)

Herschel observations of extreme OH/IR stars - the isotopic ratios of oxygen as a sign-post for the stellar mass

Authors: K. Justtanont, M. J. Barlow, J. Blommaert, L. Decin, F. Kerschbaum, M. Matsuura, H. Olofsson, P. Owen, P. Royer, B. Swinyard, D. Teyssier, L. B. F. M. Waters, J. Yates

Abstract: Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the st… ▽ More Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy. △ Less

Submitted 2 June, 2015; v1 submitted 21 May, 2015; originally announced May 2015.

Comments: Accepted for publication in A&A. 50 pages including online material (28 figures)

The JCMT Gould Belt Survey: SCUBA-2 observations of circumstellar disks in L 1495

Authors: J. V. Buckle, E. Drabek-Maunder, J. Greaves, J. S. Richer, B. C. Matthews, D. Johnstone, H. Kirk, S. F. Beaulieu, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, J. Hatchell, T. Jenness, J. C. Mottram, D. Nutter, K. Pattle, J. E. Pineda, C. Salji, S. Tisi, J. Di Francesco, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, H. Butner , et al. (29 additional authors not shown)

Abstract: We present 850$μ$m and 450$μ$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk… ▽ More We present 850$μ$m and 450$μ$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M. △ Less

Submitted 27 February, 2015; originally announced February 2015.

Comments: 11 pages, 10 figures plus Appendix. MNRAS accepted

The JCMT Gould Belt Survey: First results from the SCUBA-2 observations of the Ophiuchus molecular cloud and a virial analysis of its prestellar core population

Authors: K. Pattle, D. Ward-Thompson, J. M. Kirk, G. J. White, E. Drabek-Maunder, J. Buckle, S. F. Beaulieu, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, J. Hatchell, H. Kirk, T. Jenness, D. Johnstone, J. C. Mottram, D. Nutter, J. E. Pineda, C. Quinn, C. Salji, S. Tisi, S. Walker-Smith, J. Di Francesco, M. R. Hogerheijde, Ph. André , et al. (37 additional authors not shown)

Abstract: In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We… ▽ More In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region. △ Less

Submitted 20 February, 2015; originally announced February 2015.

Comments: 34 pages, 20 figures, 5 tables. Accepted for publication in MNRAS

Mapping CS in Starburst Galaxies: Disentangling and Characterising Dense Gas

Authors: G. Kelly, S. Viti, E. Bayet, R. Aladro, J. Yates

Abstract: Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense gas varies across the circumnuclear regions in starburst galaxies. Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2-1) and CS(3-2) lines in the circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected the formaldeh… ▽ More Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense gas varies across the circumnuclear regions in starburst galaxies. Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2-1) and CS(3-2) lines in the circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected the formaldehyde (H2CO) and methanol (CH3OH) in both galaxies. We marginally detect the isotopologue C34S. Results. We calculate column densities under LTE conditions for CS and CH3OH. Using the detections accumulated here to guide our inputs, we link a time and depth dependent chemical model with a molecular line radiative transfer model; we reproduce the observations, showing how conditions where CS is present are likely to vary away from the galactic centres. Conclusions. Using the rotational diagram method for CH3OH, we obtain a lower limit temperature of 14 K. In addition to this, by comparing the chemical and radiative transfer models to observations, we determine the properties of the dense gas as traced by CS (and CH3OH). We also estimate the quantity of the dense gas. We find that, provided that there are a between 10^5 and 10^6 dense cores in our beam, for both target galaxies, emission of CS from warm (T = 100 - 400 K), dense (n(H2) = 10^5-6 cm-3) cores, possibly with a high cosmic ray ionisation rate (zeta = 100 zeta0) best describes conditions for our central pointing. In NGC 6946, conditions are generally cooler and/or less dense further from the centre, whereas in NGC 3079, conditions are more uniform. The inclusion of shocks allows for more efficient CS formation, leading to an order of magnitude less dense gas being required to replicate observations in some cases. △ Less

Submitted 12 March, 2015; v1 submitted 16 February, 2015; originally announced February 2015.

Comments: 14 pages, 10 figures, accepted to A&A

Journal ref: A&A 578, A70 (2015)