Showing 1–50 of 231 results for author: Müller, H

Millimeter-millimeter-wave double-modulation double-resonance spectroscopy

Authors: Oliver Zingsheim, Luis Bonah, Frank Lewen, Sven Thorwirth, Holger S. P. Müller, Stephan Schlemmer

Abstract: A new millimeter- to millimeter-wave double-modulation double-resonance (MMW-MMW DM-DR) scheme has been applied to record spectra of two astronomically relevant complex organic molecules (COMs), propanal (C2H5CHO) and ethyl cyanide (C2H5CN), to demonstrate advantages of the DM-DR experimental technique. The DR technique helps to identify target transitions in a forest of lines and the implementati… ▽ More A new millimeter- to millimeter-wave double-modulation double-resonance (MMW-MMW DM-DR) scheme has been applied to record spectra of two astronomically relevant complex organic molecules (COMs), propanal (C2H5CHO) and ethyl cyanide (C2H5CN), to demonstrate advantages of the DM-DR experimental technique. The DR technique helps to identify target transitions in a forest of lines and the implementation of a DM procedure (modulation of the pump and probe source) allows for confusion- and baseline-free spectra containing only the line(s) of interest. In particular the unambiguous assignment of weak and blended transitions in very dense MMW spectra is highlighted. Details of the observed Autler-Townes line splitting and possible future applications, such as automated analyses and adaptions of DM-DR methods to other experimental setups, are discussed. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Investigation of the rotational spectrum of CH$_3$$^{17}$OH and its tentative detection toward Sagittarius B2(N)

Authors: Holger S. P. Müller, Vadim V. Ilyushin, Arnaud Belloche, Frank Lewen, Stephan Schlemmer

Abstract: Methanol is an abundant molecule in space. The column density of CH$_3^{18}$OH is in some star-forming regions so high that the search for CH$_3^{17}$OH is promising. But only very few transition frequencies of CH$_3^{17}$OH with a microwave accuracy have been published thus far. We recorded the rotational spectrum of CH$_3^{17}$OH between 38 and 1095 GHz employing a methanol sample enriched in… ▽ More Methanol is an abundant molecule in space. The column density of CH$_3^{18}$OH is in some star-forming regions so high that the search for CH$_3^{17}$OH is promising. But only very few transition frequencies of CH$_3^{17}$OH with a microwave accuracy have been published thus far. We recorded the rotational spectrum of CH$_3^{17}$OH between 38 and 1095 GHz employing a methanol sample enriched in $^{17}$O to 20\%. A torsion-rotation Hamiltonian model based on the rho-axis method was employed to fit the data, as in our previous studies. We searched for rotational transitions of CH$_3^{17}$OH in the imaging spectral line survey ReMoCA obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the high-mass star-forming region Sgr B2(N). The observed spectra were modeled under the assumption of local thermodynamic equilibrium (LTE). The assignments cover $0 \le J \le 45$, $K_a \le 16$, and mainly the $v_ t = 0$ and 1 torsional states. The Hamiltonian model describes our data well. The model was applied to derive a line list for radio-astronomical observations. We report a tentative detection of CH$_3^{17}$OH along with secure detections of the more abundant isotopologs of methanol toward Sgr B2(N2b). The derived column densities yield isotopic ratios $^{12}$C/$^{13}$C = 25, $^{16}$O/$^{18}$O = 240, and $^{18}$O/$^{17}$O = 3.3, which are consistent with values found earlier for other molecules in Sgr B2. The agreement between the $^{18}$O/$^{17}$O isotopic ratio that we obtained for methanol and the $^{18}$O/$^{17}$O ratios reported in the past for other molecules in Sgr B2(N) strongly supports our tentative interstellar identification of CH$_3^{17}$OH. The accuracy of the derived line list is sufficient for further radio astronomical searches for this methanol isotopolog toward other star-forming regions. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 16 pages, including tables, figures, references and appendices. Abstract slightly shortened. Astronomy and Astrophysics, accepted for publication

Interstellar Neutral Hydrogen in the Heliosphere: New Horizons Observations in the Context of Models

Authors: P. Swaczyna, M. Bzowski, K. Dialynas, L. Dyke, F. Fraternale, A. Galli, J. Heerikhuisen, M. Z. Kornbleuth, D. Koutroumpa, I. Kowalska-Leszczyńska, M. A. Kubiak, A. T. Michael, H. -R. Müller, M. Opher, F. Rahmanifard

Abstract: Interstellar neutral (ISN) hydrogen is the most abundant species in the outer heliosheath and the very local interstellar medium (VLISM). Charge exchange collisions in the outer heliosheath result in filtration, reducing the ISN hydrogen density inside the heliosphere. Additionally, these atoms are intensively ionized close to the Sun, resulting in a substantial reduction of their density within a… ▽ More Interstellar neutral (ISN) hydrogen is the most abundant species in the outer heliosheath and the very local interstellar medium (VLISM). Charge exchange collisions in the outer heliosheath result in filtration, reducing the ISN hydrogen density inside the heliosphere. Additionally, these atoms are intensively ionized close to the Sun, resulting in a substantial reduction of their density within a few au from the Sun. The products of this ionization - pickup ions (PUIs) - are detected by charged particle detectors. The Solar Wind Around Pluto (SWAP) instrument on New Horizons provides, for the first time, PUI observations from the distant heliosphere. We analyze the observations collected between 22 and 52 au from the Sun to find the ISN hydrogen density profile and compare the results with predictions from global heliosphere models. We conclude that the density profile derived from the observations is inconsistent with steady-state model predictions. This discrepancy is not explained by time variations close to the Sun and thus may be related to the temporal evolution of the outer boundaries or VLISM conditions. Furthermore, we show that the cold and hot models of ISN hydrogen distribution are not a good approximation closer to the termination shock. Therefore, we recommend a new fiduciary point based on the available New Horizons observations at 40 au from the Sun, at ecliptic direction (285.62°, 1.94°), where the ISN hydrogen density is 0.11 cm$^{-3}$. The continued operation of New Horizons should give better insight into the source of the discussed discrepancy. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Author Accepted Manuscript. Accepted for publication in The Astrophysical Journal Letters. 11 pages, 4 figures, 1 table

Rotational spectroscopy of CH$_3$OD with a reanalysis of CH$_3$OD toward IRAS 16293$-$2422

Authors: V. V. Ilyushin, H. S. P. Müller, M. N. Drozdovskaya, J. K. Jørgensen, S. Bauerecker, C. Maul, R. Porohovoi, E. A. Alekseev, O. Dorovskaya, O. Zakharenko, F. Lewen, S. Schlemmer, R. M. Lees

Abstract: We have started a measurement campaign of numerous methanol isotopologs in low-lying torsional states in order to provide extensive line lists for radio astronomical observations from an adequate spectroscopic model and to investigate how the intricate vibration-torsion-rotation interactions manifest themselves in the spectra of different isotopic species. After CD$_3$OH and CD$_3$OD, we turn our… ▽ More We have started a measurement campaign of numerous methanol isotopologs in low-lying torsional states in order to provide extensive line lists for radio astronomical observations from an adequate spectroscopic model and to investigate how the intricate vibration-torsion-rotation interactions manifest themselves in the spectra of different isotopic species. After CD$_3$OH and CD$_3$OD, we turn our focus to CH$_3$OD, which is an important species for studying deuteration in prestellar cores and envelopes that enshroud protostars. Notably, deuteration is frequently viewed as a diagnostic tool for star formation. The measurements used in this study were obtained in two spectroscopic laboratories and cover large fractions of the 34 GHz--1.35 THz range. As done in previous studies, we employed a torsion-rotation Hamiltonian model for our analysis that is based on the rho-axis method. The resulting model describes the ground and first excited torsional states of CH$_3$OD well up to quantum numbers $J \leqslant 51$ and $K_a \leqslant 18$. We derived a line list for radio astronomical observations from this model that is accurate up to at least 1.35~THz and should be sufficient for all types of radio astronomical searches for this methanol isotopolog in these two lowest torsional states. This line list was applied to a reinvestigation of CH$_3$OD in data from the Protostellar Interferometric Line Survey of IRAS 16293--2422 obtained with the Atacama Large Millimeter/submillimeter Array. The new accurately determined value for the column density of CH$_3$OD implies that the deuteration in methanol differs in its two functional groups by a factor of $\sim$7.5. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 12 pages including figures, tables, references and appendix. Astron. Astrophys., in press

Swarm intelligence for full Stokes dynamic imaging reconstruction of interferometric data

Authors: Alejandro Mus, Hendrik Müller, Andrei Lobanov

Abstract: In very long baseline interferometry (VLBI) the combination of multiple antennas permits the synthesis of a virtual telescope with a larger diameter and consequently higher resolution than the individual antennae. Yet, due to the sparse nature of the array, recovering an image from the observed data is a challenging ill-posed inverse problem. The VLBI community is interested in not only recovering… ▽ More In very long baseline interferometry (VLBI) the combination of multiple antennas permits the synthesis of a virtual telescope with a larger diameter and consequently higher resolution than the individual antennae. Yet, due to the sparse nature of the array, recovering an image from the observed data is a challenging ill-posed inverse problem. The VLBI community is interested in not only recovering an image in total intensity from interferometric data, but also to obtain results in the polarimetric and the temporal domain. Only a few algorithms are able to work in all these domains simultaneously. In particular, the algorithms based on optimization that consider various penalty terms specific to static total intensity imaging, time-variability and polarimetry are restricted to grids the domain of the objective function. In this work we present a novel algorithm, multiobjective particle swarm optimization, that is able to recover the optimal weights without any space-gridding, and to obtain the marginal contribution of each the playing terms. To this end, we utilize multiobjective optimization together with particle swarm metaheuristics. We let the swarm of weights to converge together to the best position. We evaluate our algorithm with representative synthetic data sets focused on the instrumental configuration of the Event Horizon Telescope Collaboration and its planned successors. We successfully recover the polarimetric, static and time-dynamic signature of the ground truth movie, even with relative sparsity, and a set of realistic data corruptions. This is a novel, fast, weighting space gridding-free algorithm that successfully recovers static and dynamic polarimetric reconstructions. Compared to Regularized Maximum Likelihood methods, it avoids the need for parameter surveys, and it is not limited to the number of pixels such as recently proposed multiobjective imaging algorithms. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: Accepted in A&A

The rotation-tunneling spectrum of 3-hydroxypropenal and confirmation of its detection toward IRAS 16293$-$2422 B

Authors: H. S. P. Müller, A. Coutens, J. K. Jørgensen, L. Margulès, R. A. Motiyenko, J. -C. Guillemin

Abstract: 3-Hydroxypropenal (HOCHCHCHO) is the lower energy tautomer of malonaldehyde which displays a complex rotation-tunneling spectrum. It was detected tentatively toward the solar-type protostar IRAS 16293$-$2422 B with ALMA in the framework of the Protostellar Interferometric Line Survey (PILS). Several transitions, however, had large residuals, preventing not only their detection, but also the excita… ▽ More 3-Hydroxypropenal (HOCHCHCHO) is the lower energy tautomer of malonaldehyde which displays a complex rotation-tunneling spectrum. It was detected tentatively toward the solar-type protostar IRAS 16293$-$2422 B with ALMA in the framework of the Protostellar Interferometric Line Survey (PILS). Several transitions, however, had large residuals, preventing not only their detection, but also the excitation temperature of the species from being determined unambiguously. We want to extend the existing rotational line list of 3-hydroxypropenal to shed more light on the recent observational results and to facilitate additional radio astronomical searches for this molecule. We analyzed the rotation-tunneling spectrum of 3-hydroxypropenal in the frequency regions between 150 and 330 GHz and between 400 and 660 GHz. Transitions were searched for in the PILS observations of IRAS 16293$-$2422. Local thermodynamic equilibrium (LTE) models were carried out and compared to the observations to constrain the excitation temperature. Additional transitions were searched for in other ALMA archival data of the same source to confirm the presence of 3-hydroxypropenal. More than 11500 transitions were assigned in the course of our investigation with quantum numbers $2 \le J \le 100$, $K_a \le 59$, and $K_c \le 97$, resulting in a greatly improved set of spectroscopic parameters. The comparison between the LTE models and the observations yields an excitation temperature of 125 K with a column density $N = 1.0 \times 10^{15}$ cm$^{-2}$ for this species. We identified seven additional lines of 3-hydroxypropenal that show a good agreement with the model in the ALMA archive data. The calculated rotation-tunneling spectrum of 3-hydroxypropenal has sufficient accuracy for radio astronomical searches. The detection of 3-hydroxypropenal toward IRAS 16293$-$2422 B is now secure. △ Less

Submitted 12 June, 2024; v1 submitted 2 May, 2024; originally announced May 2024.

Comments: 11 pages with figures, tables, and references. Abstract here slightly abbreviated. Accepted at Astron. Astrophys

Journal ref: A&A 687, A8 (2024)

Effects of an Explicit Time-Dependent Radiation Pressure Force on Trajectories of Primary Neutral Hydrogen in the Heliosphere

Authors: Lucas Dyke, Hans-Reinhard Müller

Abstract: Radiation pressure exerted by solar photon output is salient to the motion of primary neutral hydrogen atoms streaming into the inner heliosphere directly from the Local Interstellar Medium (LISM). The action of a time-dependent radiation pressure force, when coupled with the usual gravitational force, changes the characteristic velocities, and therefore energies, of the atoms when they reach regi… ▽ More Radiation pressure exerted by solar photon output is salient to the motion of primary neutral hydrogen atoms streaming into the inner heliosphere directly from the Local Interstellar Medium (LISM). The action of a time-dependent radiation pressure force, when coupled with the usual gravitational force, changes the characteristic velocities, and therefore energies, of the atoms when they reach regions in which explorer probes are present. A study is presented that uses a 2D code to backtrace neutral hydrogen trajectories from representative target points located 1 au from the Sun. It makes use of both a radiation pressure function and a function for the photoionization rate at 1 au that oscillate with time based on measurements over a typical solar cycle, as well as a time-independent charge exchange ionization rate at 1 au. Assuming a Maxwellian distribution in the distant upwind direction, phase space data is calculated at the target points, at different moments in time. The dependence of the force on the radial particle velocity has been omitted in the analysis, such that the emphasis is on the effects of the global solar UV intensity variations through the solar cycle. This process allows for analysis of direct and indirect Maxwellian components through time and space in the time-dependent force environment. Additionally, pseudo-bound orbits caused by energy losses associated with this force environment are observed and their properties evaluated with the aim of determining their effects on potential measurements by explorer probes. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 15 pages, 20 figures, to be published in The Astrophysical Journal (ApJ)

Resolving Red Giant Winds with the Hubble Space Telescope

Authors: Brian E. Wood, Graham M. Harper, Hans-Reinhard Mueller

Abstract: We describe recent spectroscopic observations of red giant stars made by the Space Telescope Imaging Spectrograph (STIS) instrument on board the Hubble Space Telescope, which have provided spatially resolved observations of the warm chromospheric winds that predominate for early K to mid-M giants. The H I Lyman-alpha lines of a set of 11 red giants observed with the STIS/E140M echelle grating are… ▽ More We describe recent spectroscopic observations of red giant stars made by the Space Telescope Imaging Spectrograph (STIS) instrument on board the Hubble Space Telescope, which have provided spatially resolved observations of the warm chromospheric winds that predominate for early K to mid-M giants. The H I Lyman-alpha lines of a set of 11 red giants observed with the STIS/E140M echelle grating are first analyzed to ascertain wind H I column densities and total wind mass-loss rates. The M giants have estimated mass-loss rates of Mdot=(14-86)e-11 Msun/yr, while the K giants with detected wind absorption have weaker winds with Mdot=(1.5-2.8)e-11 Msun/yr. We use long-slit spectra of H I Lyman-alpha for two particular red giants, Alpha Tau (K5 III) and Gamma Cru (M3.5 III), to study the spatial extent of the Lyman-alpha emission. From these data we estimate limits for the extent of detectable emission, which are r=193 Rstar for Gamma Cru and r=44 Rstar for Alpha Tau. Cross-dispersion emission profiles in the STIS echelle spectra of the larger sample of red giants also show evidence for spatial resolution, not only for H I Lyman-alpha but for other lines with visible wind absorption, such as Fe II, Mg II, Mg I, O I, and C II. We characterize the nature of these spatial signatures. The spatial extent is far more apparent for the M giants than for the K giants, consistent with the stronger winds found for the M giants from the Lyman-alpha analysis. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 31 pages, 15 figures, to appear in The Astrophysical Journal

Discovery of thionylimide, HNSO, in space: the first N-, S- and O-bearing interstellar molecule

Authors: Miguel Sanz-Novo, Víctor M. Rivilla, Holger S. P. Müller, Izaskun Jiménez-Serra, Jesús Martín-Pintado, Laura Colzi, Shaoshan Zeng, Andrés Megías, Álvaro López-Gallifa, Antonio Martínez-Henares, Belén Tercero, Pablo de Vicente, David San Andrés, Sergio Martín, Miguel A. Requena-Torres

Abstract: We present the first detection in space of thionylimide (HNSO) toward the Galactic Center molecular cloud G+0.693-0.027, thanks to the superb sensitivity of an ultradeep molecular line survey carried out with the Yebes 40$\,$m and IRAM 30$\,$m telescopes. This molecule is the first species detected in the interstellar medium containing, simultaneously, N, S and O. We have identified numerous $K$… ▽ More We present the first detection in space of thionylimide (HNSO) toward the Galactic Center molecular cloud G+0.693-0.027, thanks to the superb sensitivity of an ultradeep molecular line survey carried out with the Yebes 40$\,$m and IRAM 30$\,$m telescopes. This molecule is the first species detected in the interstellar medium containing, simultaneously, N, S and O. We have identified numerous $K$$_a$ = 0, 1 and 2 transitions belonging to HNSO covering from $J$$_{\rm up}$ = 2 to $J$$_{\rm up}$ = 10, including several completely unblended features. We derive a molecular column density of $N$ = (8 $\pm$ 1)$\times$10$^{13}$ cm$^{-2}$, yielding a fractional abundance relative to H$_2$ of $\sim$6$\times$10$^{-10}$, which is about $\sim$37 and $\sim$4.8 times less abundant than SO and SO2, respectively. Although there are still many unknowns in the interstellar chemistry of NSO-bearing molecules, we propose that HNSO is likely formed through the reaction of the NSO radical and atomic H on the surface of icy grains, with alternative routes also deserving exploration. Finally, HNSO appears as a promising link between N- , S- and O- interstellar chemistry and its discovery paves the route to the detection of a new family of molecules in space. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Forthcoming paper in The Astrophysical Journal Letters (in press)

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Using multiobjective optimization to reconstruct interferometric data (II): polarimetry and time dynamics

Authors: Alejandro Mus, Hendrik Müller, Ivan Martí-Vidal, Andrei Lobanov

Abstract: In Very Long Baseline Interferometry (VLBI), signals from multiple antennas combine to create a sparsely sampled virtual aperture, its effective diameter determined by the largest antenna separation. The inherent sparsity makes VLBI imaging an ill-posed inverse problem, prompting the use of algorithms like the Multiobjective Evolutionary Algorithm by Decomposition (MOEA/D), as proposed in the firs… ▽ More In Very Long Baseline Interferometry (VLBI), signals from multiple antennas combine to create a sparsely sampled virtual aperture, its effective diameter determined by the largest antenna separation. The inherent sparsity makes VLBI imaging an ill-posed inverse problem, prompting the use of algorithms like the Multiobjective Evolutionary Algorithm by Decomposition (MOEA/D), as proposed in the first paper of this series. This study focuses on extending MOEA/D to polarimetric and time dynamic reconstructions, particularly relevant for the VLBI community and the Event Horizon Telescope Collaboration (EHTC). MOEA/D's success in providing a unique, fast, and largely unsupervised representation of image structure serves as the basis for exploring these extensions. The extension involves incorporating penalty terms specific to total intensity imaging, time-variable, and polarimetric variants within MOEA/D's multiobjective, evolutionary framework. The Pareto front, representing non-dominated solutions, is computed, revealing clusters of proximities. Testing MOEA/D with synthetic datasets representative of EHTC's main targets demonstrates successful recovery of polarimetric and time-dynamic signatures despite sparsity and realistic data corruptions. MOEA/D's extension proves effective in the anticipated EHTC setting, offering an alternative and independent claim to existing methods. It not only explores the problem globally but also eliminates the need for parameter surveys, distinguishing it from Regularized Maximum Likelihood (RML) methods. MOEA/D emerges as a novel and useful tool for robustly characterizing polarimetric and dynamic signatures in VLBI datasets with minimal user-based choices. Future work aims to address the last remaining limitation of MOEA/D, specifically regarding the number of pixels and numerical performance, to establish it within the VLBI data reduction pipeline. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: Both first authors have contributed equally to this work. To appear in A&A

Identifying synergies between VLBI and STIX imaging

Authors: Hendrik Müller, Paolo Massa, Alejandro Mus, Jong-Seo Kim, Emma Perracchione

Abstract: Reconstructing an image from sparsely sampled Fourier data is an ill-posed inverse problem that occurs in a variety of subjects within science, including the data analysis for Very Long Baseline Interferometry (VLBI) and the Spectrometer/Telescope for Imaging X-rays (STIX) for solar observations. Despite ongoing parallel developments of novel imaging algorithms, synergies remain unexplored. We stu… ▽ More Reconstructing an image from sparsely sampled Fourier data is an ill-posed inverse problem that occurs in a variety of subjects within science, including the data analysis for Very Long Baseline Interferometry (VLBI) and the Spectrometer/Telescope for Imaging X-rays (STIX) for solar observations. Despite ongoing parallel developments of novel imaging algorithms, synergies remain unexplored. We study the synergies between the data analysis for the STIX instrument and VLBI, compare the methodologies and evaluate their potential. In this way, we identify key trends in the performance of several algorithmic ideas and draw recommendations for the future. To this end, we organized a semi-blind imaging challenge with data sets and source structures that are typical for sparse VLBI, specifically in the context of the Event Horizon Telescope (EHT), and for STIX observations. 17 different algorithms from both communities, from 6 different imaging frameworks, participated in the challenge, marking this work the largest scale code comparisons for STIX and VLBI to date. Strong synergies between the two communities have been identified, as can be proven by the success of the imaging methods proposed for STIX in imaging VLBI data sets and vice versa. Novel imaging methods outperform the standard CLEAN algorithm significantly in every test-case. Improvements over the performance of CLEAN make deeper updates to the inverse modeling pipeline necessary, or consequently replacing inverse modeling with forward modeling. Entropy-based and Bayesian methods perform best on STIX data. The more complex imaging algorithms utilizing multiple regularization terms (recently proposed for VLBI) add little to no additional improvements for STIX, but outperform the other methods on EHT data. This work demonstrates the great synergy between the STIX and VLBI imaging efforts and the great potential for common developments. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: accepted for publication in A&A

The compact multi-planet system GJ 9827 revisited with ESPRESSO

Authors: V. M. Passegger, A. Suárez Mascareño, R. Allart, J. I. González Hernández, C. Lovis, B. Lavie, A. M. Silva, H. M. Müller, H. M. Tabernero, S. Cristiani, F. Pepe, R. Rebolo, N. C. Santos, V. Adibekyan, Y. Alibert, C. Allende Prieto, S. C. C. Barros, F. Bouchy, A. Castro-González, V. D'Odorico, X. Dumusque, P. Di Marcantonio, D. Ehrenreich, P. Figueira, R. Génova Santos , et al. (14 additional authors not shown)

Abstract: GJ 9827 is a bright, nearby K7V star orbited by two super-Earths and one mini-Neptune on close-in orbits. The system was first discovered using K2 data and then further characterized by other spectroscopic and photometric instruments. Previous literature studies provide several mass measurements for the three planets, however, with large variations and uncertainties. To better constrain the planet… ▽ More GJ 9827 is a bright, nearby K7V star orbited by two super-Earths and one mini-Neptune on close-in orbits. The system was first discovered using K2 data and then further characterized by other spectroscopic and photometric instruments. Previous literature studies provide several mass measurements for the three planets, however, with large variations and uncertainties. To better constrain the planetary masses, we added high-precision radial velocity measurements from ESPRESSO to published datasets from HARPS, HARPS-N, and HIRES and we performed a Gaussian process analysis combining radial velocity and photometric datasets from K2 and TESS. This method allowed us to model the stellar activity signal and derive precise planetary parameters. We determined planetary masses of $M_b = 4.28_{-0.33}^{+0.35}$ M${_\oplus}$, $M_c = 1.86_{-0.39}^{+0.37}$ M${_\oplus}$, and $M_d = 3.02_{-0.57}^{+0.58}$ M${_\oplus}$, and orbital periods of $1.208974 \pm 0.000001$ days for planet b, $3.648103_{-0.000010}^{+0.000013}$ days for planet c, and $6.201812 \pm 0.000009$ days for planet d. We compared our results to literature values and found that our derived uncertainties for the planetary mass, period, and radial velocity amplitude are smaller than the previously determined uncertainties. We modeled the interior composition of the three planets using the machine-learning-based tool ExoMDN and conclude that GJ 9827 b and c have an Earth-like composition, whereas GJ 9827 d has an hydrogen envelope, which, together with its density, places it in the mini-Neptune regime. △ Less

Submitted 16 January, 2024; v1 submitted 11 January, 2024; originally announced January 2024.

Comments: accepted in A&A

A deep search for large complex organic species toward IRAS16293-2422 B at 3 mm with ALMA

Authors: P. Nazari, J. S. Y. Cheung, J. Ferrer Asensio, N. M. Murillo, E. F. van Dishoeck, J. K. Jørgensen, T. L. Bourke, K. -J. Chuang, M. N. Drozdovskaya, G. Fedoseev, R. T. Garrod, S. Ioppolo, H. Linnartz, B. A. McGuire, H. S. P. Müller, D. Qasim, S. F. Wampfler

Abstract: Complex organic molecules (COMs) have been detected ubiquitously in protostellar systems. However, at shorter wavelengths (~0.8mm) it is more difficult to detect larger molecules than at longer wavelengths (~3mm) because of the increase of millimeter dust opacity, line confusion, and unfavorable partition function. We aim to search for large molecules (>8 atoms) in the ALMA Band 3 spectrum of IRAS… ▽ More Complex organic molecules (COMs) have been detected ubiquitously in protostellar systems. However, at shorter wavelengths (~0.8mm) it is more difficult to detect larger molecules than at longer wavelengths (~3mm) because of the increase of millimeter dust opacity, line confusion, and unfavorable partition function. We aim to search for large molecules (>8 atoms) in the ALMA Band 3 spectrum of IRAS 16293-2422 B. We search for more than 70 molecules and identify as many lines as possible in the spectrum. The spectral settings were set to specifically target three-carbon species such as propanol and glycerol. We identify lines of 31 molecules including many oxygen-bearing COMs such as CH3OH and c-C2H4O and a few nitrogen- and sulfur-bearing ones such as HOCH2CN and CH3SH. The largest detected molecules are gGg-(CH2OH)2 and CH3COCH3. We do not detect glycerol or propanol but provide upper limits for them which are in line with previous laboratory and observational studies. The line density in Band 3 is only ~2.5 times lower in frequency space than in Band 7. From the detected lines in Band 3 at a $\gtrsim 6σ$ level, ~25-30% of them could not be identified indicating the need for more laboratory data of rotational spectra. We find similar column densities and column density ratios of COMs (within a factor ~2) between Band 3 and Band 7. The effect of dust optical depth for IRAS 16293-2422 B at an off-source location on column densities and column density ratios is minimal. Moreover, for warm protostars, long wavelength spectra are not only crowded, but also take longer integration times to reach the same sensitivity limit. The 3mm search has not yet resulted in detection of larger and more complex molecules in warm sources. A full deep ALMA Band 2-3 (i.e., 3-4 mm) survey is needed to assess whether low frequency data have the potential to reveal more complex molecules in warm sources. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 686, A59 (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

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Shocking Sgr B2(N1) with its own outflow: A new perspective on segregation between O- and N-bearing molecules

Authors: Laura A. Busch, Arnaud Belloche, Robin T. Garrod, Holger S. P. Mueller, Karl M. Menten

Abstract: We want to investigate the influence of the powerful outflow driven by the hot core Sgr B2(N1) on the gas molecular inventory of the surrounding medium. We used the data taken as part of the 3 mm imaging spectral-line survey ReMoCA (Re-exploring Molecular Complexity with ALMA). Integrated intensity maps of SO and SiO emission reveal a bipolar structure with blue-shifted emission dominantly extendi… ▽ More We want to investigate the influence of the powerful outflow driven by the hot core Sgr B2(N1) on the gas molecular inventory of the surrounding medium. We used the data taken as part of the 3 mm imaging spectral-line survey ReMoCA (Re-exploring Molecular Complexity with ALMA). Integrated intensity maps of SO and SiO emission reveal a bipolar structure with blue-shifted emission dominantly extending to the SE from the centre of the hot core and red-shifted emission to the NW. This is also prominently observed in emission of other S-bearing molecules and species that only contain N as a heavy element, including COMs, but also CH3OH, CH3CHO, HNCO, and NH2CHO. For a selection of COMs and simpler species, spectra were modelled under the assumption of LTE and population diagrams were derived at two positions, one in each outflow lobe. From this analysis, we obtained rotational temperatures, which are in a range of ~100-200K, and column densities. Abundances were subsequently compared to predictions of astrochemical models and to observations of L1157-B1, a position located in the well-studied outflow of the low-mass protostar L1157, and the source G+0.693-0.027, located in the Sgr B2 molecular cloud complex. Given the short distance of the analysed outflow positions to the centre of Sgr B2(N1), we propose a scenario in which a phase of hot-core chemistry (i.e. thermal desorption of ice species and high-temperature gas-phase chemistry) preceded a shock wave. The subsequent compression and further heating of the material resulted in the accelerated destruction of (mainly O-bearing) molecules. Gas-phase formation of cyanides seems to be able to compete with their destruction in the post-shock gas. Abundances of HCnN (n=3,5) are enhanced in the outflow component pointing to (additional) gas-phase formation. To confirm such a scenario, appropriate chemical shock models need to be run. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A

Millimetre and submillimetre spectroscopy of isobutene and its detection in the molecular cloud G+0.693

Authors: Mariyam Fatima, Holger S. P. Müller, Oliver Zingsheim, Frank Lewen, Víctor M. Rivilla, Izaskun Jiménez-Serra, Jesús Martín-Pintado, Stephan Schlemmer

Abstract: Isobutene ((CH$_3$)$_2$C=CH$_2$) is one of the four isomers of butene (C$_4$H$_8$). Given the detection of propene (CH$_3$CH=CH$_2$) toward TMC-1, and also in the warmer environment of the solar-type protostellar system IRAS 16293$-$2422, one of the next alkenes, isobutene, is a promising candidate to be searched for in space. We aim to extend the limited line lists of the main isotopologue of iso… ▽ More Isobutene ((CH$_3$)$_2$C=CH$_2$) is one of the four isomers of butene (C$_4$H$_8$). Given the detection of propene (CH$_3$CH=CH$_2$) toward TMC-1, and also in the warmer environment of the solar-type protostellar system IRAS 16293$-$2422, one of the next alkenes, isobutene, is a promising candidate to be searched for in space. We aim to extend the limited line lists of the main isotopologue of isobutene from the microwave to the millimetre region in order to obtain a highly precise set of rest frequencies and to facilitate its detection in the interstellar medium. We investigated the rotational spectrum of isobutene in the 35$-$370 GHz range using absorption spectroscopy at room temperature. Quantum-chemical calculations were carried out to evaluate vibrational frequencies. We determined new or improved spectroscopic parameters for isobutene up to a sixth-order distortion constant. These new results enabled its detection in the G+0.693 molecular cloud for the first time, where propene was also recently found. The propene to isobutene column density ratio was determined to be about 3:1. The observed spectroscopic parameters for isobutene are sufficiently accurate that calculated transition frequencies should be reliable up to 700 GHz. This will further help in observing this alkene in other, warmer regions of the ISM. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Comments: 11 pages including figures, tables, references, and appendix; Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 680, A25 (2023)

Rotational spectroscopy of the thioformaldehyde isotopologues H$_2$CS and H$_2$C$^{34}$S in four interacting excited vibrational states and an account on the rotational spectrum of thioketene, H$_2$CCS

Authors: Holger S. P. Müller, Atsuko Maeda, Frank Lewen, Stephan Schlemmer, Ivan R. Medvedev, Eric Herbst

Abstract: An investigation of the rotational spectrum of the interstellar molecule thioformaldehyde between 110 and 377 GHz through a pyrolysis reaction revealed a multitude of absorption lines assignable to H$_2$CS and H$_2$C$^{34}$S in their lowest four excited vibrational states besides lines of numerous thioformaldehyde isotopologues in their ground vibrational states reported earlier as well as lines p… ▽ More An investigation of the rotational spectrum of the interstellar molecule thioformaldehyde between 110 and 377 GHz through a pyrolysis reaction revealed a multitude of absorption lines assignable to H$_2$CS and H$_2$C$^{34}$S in their lowest four excited vibrational states besides lines of numerous thioformaldehyde isotopologues in their ground vibrational states reported earlier as well as lines pertaining to several by-products. Additional transitions of H$_2$CS in its lowest four excited vibrational states were recorded in selected regions between 571 and 1386 GHz. Slight to strong Coriolis interactions occur between all four vibrational states with the exception of the two highest lying states because both are totally symmetric vibrations. We present combined analyses of the ground and the four interacting states for our rotational data of H$_2$CS and H$_2$C$^{34}$S. The H$_2$CS data were supplemented with two sets of high-resultion IR data in two separate analyses. The $v_2 = 1$ state has been included in analyses of Coriolis interactions of low-lying fundamental states of H$_2$CS for the first time and this improved the quality of the fits substantially. We extended furthermore assignments in $J$ of transition frequencies of thioketene in its ground vibrational state. △ Less

Submitted 23 September, 2023; v1 submitted 16 September, 2023; originally announced September 2023.

Comments: 23 pages including figures, tables, and references; Mol. Phys., accepted (for the Tim Lee memorial issue)

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Investigation of the rotational spectrum of CD$_3$OD and an astronomical search toward IRAS 16293$-$2422

Authors: V. V. Ilyushin, H. S. P. Müller, J. K. Jørgensen, S. Bauerecker, C. Maul, R. Porohovoi, E. A. Alekseev, O. Dorovskaya, F. Lewen, S. Schlemmer, R. M. Lees

Abstract: Solar-type prestellar cores and protostars display large amounts of deuterated organic molecules. Recent findings on CHD$_2$OH and CD$_3$OH toward IRAS 16293-2422 suggest that even fully deuterated methanol, CD$_3$OD, may be detectable as well. However, searches for CD$_3$OD are hampered in particular by the lack of intensity information from a spectroscopic model. The objective of the present inv… ▽ More Solar-type prestellar cores and protostars display large amounts of deuterated organic molecules. Recent findings on CHD$_2$OH and CD$_3$OH toward IRAS 16293-2422 suggest that even fully deuterated methanol, CD$_3$OD, may be detectable as well. However, searches for CD$_3$OD are hampered in particular by the lack of intensity information from a spectroscopic model. The objective of the present investigation is to develop a spectroscopic model of CD$_3$OD in low-lying torsional states that is sufficiently accurate to facilitate searches for this isotopolog in space. We carried out a new measurement campaign for CD$_3$OD involving two spectroscopic laboratories that covers the 34 GHz-1.1 THz range. A torsion-rotation Hamiltonian model based on the rho-axis method was employed for our analysis. Our resulting model describes the ground and first excited torsional states of CD$_3$OD well up to quantum numbers $J \leq 51$ and $K_a \leq 23$. We derived a line list for radio-astronomical observations from this model that is accurate up to at least 1.1 THz and should be sufficient for all types of radio-astronomical searches for this methanol isotopolog. This line list was used to search for CD$_3$OD in data from the Protostellar Interferometric Line Survey of IRAS 16293$-$2422 obtained with the Atacama Large Millimeter/submillimeter Array. While we found several emission features that can be attributed largely to CD$_3$OD, their number is still not sufficiently high enough to establish a clear detection. Nevertheless, the estimate of 2$\times 10^{15}$ cm$^{-2}$ derived for the CD$_3$OD column density may be viewed as an upper limit that can be compared to column densities of CD$_3$OH, CH$_3$OD, and CH$_3$OH. The comparison indicates that the CD$_3$OD column density toward IRAS 16293-2422 is in line with the enhanced D/H ratios observed for multiply deuterated complex organic molecules. △ Less

Submitted 15 July, 2023; originally announced July 2023.

Comments: 10 pages total, including figures, tables, and references; abstract slightly shortened; Astron. Astrophys., accepted

Journal ref: A&A 677, A49 (2023)

Rotation-tunneling spectrum and astrochemical modeling of dimethylamine, CH$_3$NHCH$_3$, and searches for it in space

Authors: H. S. P. Müller, R. T. Garrod, A. Belloche, V. M. Rivilla, K. M. Menten, I. Jiménez-Serra, J. Martín-Pintado, F. Lewen, S. Schlemmer

Abstract: Methylamine has been the only simple alkylamine detected in the interstellar medium for a long time. With the recent secure and tentative detections of vinylamine and ethylamine, respectively, dimethylamine has become a promising target for searches in space. Its rotational spectrum, however, has been known only up to 45 GHz until now. Here we investigate the rotation-tunneling spectrum of dimethy… ▽ More Methylamine has been the only simple alkylamine detected in the interstellar medium for a long time. With the recent secure and tentative detections of vinylamine and ethylamine, respectively, dimethylamine has become a promising target for searches in space. Its rotational spectrum, however, has been known only up to 45 GHz until now. Here we investigate the rotation-tunneling spectrum of dimethylamine in selected regions between 76 and 1091 GHz using three different spectrometers in order to facilitate its detection in space. The quantum number range is extended to $J = 61$ and $K_a = 21$, yielding an extensive set of accurate spectroscopic parameters. To search for dimethylamine, we refer to the spectral line survey ReMoCA carried out with the Atacama Large Millimeter/submillimeter Array toward the high-mass star-forming region Sagittarius B2(N) and a spectral line survey of the molecular cloud G+0.693$-$0.027 employing the IRAM 30 m and Yebes 40 m radio telescopes. We report nondetections of dimethylamine toward the hot molecular cores Sgr B2(N1S) and Sgr B2(N2b) as well as G+0.693$-$0.027 which imply that dimethylamine is at least 14, 4.5 and 39 times less abundant than methylamine toward these sources, respectively. The observational results are compared to computational results from a gas-grain astrochemical model. The modeled methylamine to dimethylamine ratios are compatible with the observational lower limits. However, the model produces too much ethylamine compared with methylamine which could mean that the already fairly low levels of dimethylamine in the models may also be too high. △ Less

Submitted 31 May, 2023; v1 submitted 19 May, 2023; originally announced May 2023.

Comments: Mon. Not. R. Astron. Soc., accepted. 33 pages including tables, figures, and appendix

Journal ref: Mon. Not. R. Astron. Soc., 523, 2887-2917 (2023)

ATOMIUM: Probing the inner wind of evolved O-rich stars with new, highly excited H$_2$O and OH lines

Authors: A. Baudry, K. T. Wong, S. Etoka, A. M. S. Richards, H. S. P. Muller, F. Herpin, T. Danilovich, M. D. Gray, S. Wallstrom, D. Gobrecht, T. Khouri, L. Decin, C. A. Gottlieb, K. M. Menten, W. Homan, T. J. Millar, M. Montarges, B. Pimpanuwat, J. M. C. Plane, P. Kervella

Abstract: Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolutio… ▽ More Water and the hydroxyl radical are major constituents of the envelope of O-rich late-type stars. Transitions involving energy levels that are highly excited have been observed in both H$_2$O and OH. These and more recently discovered transitions can now be observed at a high sensitivity and angular resolution with the ALMA Array. Spectra and maps of H$_2$O and OH observed with an angular resolution of 20 to $\sim$200 mas were obtained at two epochs with the ALMA array. Observations with the Compact Array were also used to check for time variability of water transitions. Radiative transfer models of water were revisited to characterize masing conditions and up-to-date chemical models were used for comparison with our observations. Ten rotational transitions of H$_2$O with energies up to 9000 K were observed in various vibrational states. All but one are new detections in space, and from these we have derived accurate rest frequencies. Hyperfine split $Λ$-doubling transitions in v = 0, J = 27/2 and 29/2 levels of the $^2Π_{3/2}$ state and, $J = 33/2$ and 35/2 of the $^2Π_{1/2}$ state of OH with excitation energies up to 8900 K were also observed. Four of these transitions are new detections in space. Combining our measurements with earlier observations of OH, the v = 0 and v = 1 $Λ$-doubling frequencies have been improved. Our H$_2$O maps show compact emission and extensions up to twelve stellar radii or more. The 268.149 GHz emission line of water in the v$_2$ = 2 state is time variable, tends to be masing with dominant radiative pumping, and is widely excited. The widespread but weaker 262.898 GHz water line in v$_2$ = 1 also shows signs of maser emission. Emission and absorption of both H$_2$O and OH reveal an infall of matter and complex kinematics influenced by binarity. From our observed column densities, we derived OH/H$_2$O abundance ratios in a few stars. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 41 pages with references and 25 figures in main text. 4 Tables in Appendix A. 10 figures in Appendix B. 2 figures in Appendix C. 14 figures in Appendix D. 5 figures in Appendix E. 2 figures in Appendix F

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

Using multiobjective optimization to reconstruct interferometric data (I)

Authors: Hendrik Müller, Alejandro Mus, Andrei Lobanov

Abstract: Imaging in radioastronomy is an ill-posed inverse problem. Particularly the Event Horizon Telescope (EHT) Collaboration investigated the fidelity of their image reconstructions convincingly by large surveys solving the problem with different optimization parameters. This strategy faces a limitation for the existing methods when imaging the active galactic nuclei (AGN): large and expensive surveys… ▽ More Imaging in radioastronomy is an ill-posed inverse problem. Particularly the Event Horizon Telescope (EHT) Collaboration investigated the fidelity of their image reconstructions convincingly by large surveys solving the problem with different optimization parameters. This strategy faces a limitation for the existing methods when imaging the active galactic nuclei (AGN): large and expensive surveys solving the problem with different optimization parameters are time-consumptive. We present a novel nonconvex, multiobjective optimization modeling approach that gives a different type of claim and may provide a pathway to overcome this limitation. To this end we used a multiobjective version of the genetic algorithm (GA): the Multiobjective Evolutionary Algorithm Based on Decomposition, or MOEA/D. GA strategies explore the objective function by evolutionary operations to find the different local minima, and to avoid getting trapped in saddle points. First, we have tested our algorithm (MOEA/D) using synthetic data based on the 2017 Event Horizon Telescope (EHT) array and a possible EHT + next-generation EHT (ngEHT) configuration. We successfully recover a fully evolved Pareto front of non-dominated solutions for these examples. The Pareto front divides into clusters of image morphologies representing the full set of locally optimal solutions. We discuss approaches to find the most natural guess among these solutions and demonstrate its performance on synthetic data. Finally, we apply MOEA/D to observations of the black hole shadow in Messier 87 (M87) with the EHT data in 2017. MOEA/D is very flexible, faster than any other Bayesian method and explores more solutions than Regularized Maximum Likelihood methods (RML). △ Less

Submitted 13 May, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: to appear in A&A

Journal ref: A&A 675, A60 (2023)

Rotational spectroscopy of oxirane-\textit{2,2}-$d_2$, $c$-CD$_2$CH$_2$O, and its tentative detection toward IRAS 16293$-$2422~B

Authors: Holger S. P. Müller, Jes K. Jørgensen, Jean-Claude Guillemin, Frank Lewen, Stephan Schlemmer

Abstract: We prepared a sample of oxirane doubly deuterated at one C atom and studied its rotational spectrum in the laboratory for the first time between 120~GHz and 1094~GHz. Accurate spectroscopic parameters up to eighth order were determined, and the calculated rest frequencies were used to identify $c$-CD$_2$CH$_2$O tentatively in the interstellar medium in the Atacama Large Millimeter/submillimeter Ar… ▽ More We prepared a sample of oxirane doubly deuterated at one C atom and studied its rotational spectrum in the laboratory for the first time between 120~GHz and 1094~GHz. Accurate spectroscopic parameters up to eighth order were determined, and the calculated rest frequencies were used to identify $c$-CD$_2$CH$_2$O tentatively in the interstellar medium in the Atacama Large Millimeter/submillimeter Array Protostellar Interferometric Line Survey (PILS) of the Class 0 protostellar system IRAS 16293$-$2422. The $c$-CD$_2$CH$_2$O to $c$-C$_2$H$_4$O ratio was estimated to be $\sim$0.054 with $T_{\rm rot} = 125$ K. This value translates to a D-to-H ratio of $\sim$0.16 per H atom which is higher by a factor of 4.5 than the $\sim$0.036 per H atom obtained for $c$-C$_2$H$_3$DO. Such increase in the degree of deuteration referenced to one H atom in multiply deuterated isotopologs compared to their singly deuterated variants have been observed commonly in recent years. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: Journal of Molecular Spectroscopy, in press; Per Jensen special issue. 12 pages here

Journal ref: J. Mol. Spectrosc. 394 (2023) 111777

RadioAstron Space VLBI Imaging of the jet in M87: I. Detection of high brightness temperature at 22 GHz

Authors: Jae-Young Kim, Tuomas Savolainen, Petr Voitsik, Evgeniya V. Kravchenko, Mikhail M. Lisakov, Yuri Y. Kovalev, Hendrik Müller, Andrei P. Lobanov, Kirill V. Sokolovsky, Gabriele Bruni, Philip G. Edwards, Cormac Reynolds, Uwe Bach, Leonid I. Gurvits, Thomas P. Krichbaum, Kazuhiro Hada, Marcello Giroletti, Monica Orienti, James M. Anderson, Sang-Sung Lee, Bong Won Sohn, J. Anton Zensus

Abstract: We present results from the first 22 GHz space very-long-baseline interferometric (VLBI) imaging observations of M87 by RadioAstron. As a part of the Nearby AGN Key Science Program, the source was observed in Feb 2014 at 22 GHz with 21 ground stations, reaching projected $(u,v)$-spacings up to $\sim11\,$G$λ$. The imaging experiment was complemented by snapshot RadioAstron data of M87 obtained duri… ▽ More We present results from the first 22 GHz space very-long-baseline interferometric (VLBI) imaging observations of M87 by RadioAstron. As a part of the Nearby AGN Key Science Program, the source was observed in Feb 2014 at 22 GHz with 21 ground stations, reaching projected $(u,v)$-spacings up to $\sim11\,$G$λ$. The imaging experiment was complemented by snapshot RadioAstron data of M87 obtained during 2013--2016 from the AGN Survey Key Science Program. Their longest baselines extend up to $\sim25\,$G$λ$. For all these measurements, fringes are detected only up to $\sim$2.8 Earth Diameter or $\sim$3 G$λ$ baseline lengths, resulting in a new image with angular resolution of $\sim150\,μ$as or $\sim20$ Schwarzschild radii spatial resolution. The new image not only shows edge-brightened jet and counterjet structures down to submilliarcsecond scales but also clearly resolves the VLBI core region. While the overall size of the core is comparable to those reported in the literature, the ground-space fringe detection and slightly super-resolved RadioAstron image suggest the presence of substructures in the nucleus, whose minimum brightness temperature exceeds $T_{\rm B, min}\sim10^{12}\,$K. It is challenging to explain the origin of this record-high $T_{\rm B, min}$ value for M87 by pure Doppler boosting effect with a simple conical jet geometry and known jet speed. Therefore, this can be evidence for more extreme Doppler boosting due to a blazar-like small jet viewing angle or highly efficient particle acceleration processes occurring already at the base of the outflow. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 27 pages, 13 figures, accepted for publication in ApJ

Journal ref: ApJ 952 (2023) 34

Dynamic and polarimetric VLBI imaging with a multiscalar approach

Authors: Hendrik Müller, Andrei Lobanov

Abstract: Recently multiscale imaging approaches such as DoG-HiT were developed to solve the VLBI imaging problem and showed a promising performance: they are fast, accurate, unbiased and automatic. We extend the multiscalar imaging approach to polarimetric imaging, reconstructions of dynamically evolving sources and finally to dynamic polarimetric reconstructions. These extensions (mr-support imaging) util… ▽ More Recently multiscale imaging approaches such as DoG-HiT were developed to solve the VLBI imaging problem and showed a promising performance: they are fast, accurate, unbiased and automatic. We extend the multiscalar imaging approach to polarimetric imaging, reconstructions of dynamically evolving sources and finally to dynamic polarimetric reconstructions. These extensions (mr-support imaging) utilize a multiscalar approach. The time-averaged Stokes I image is decomposed by a wavelet transform into single subbands. We use the set of statistically significant wavelet coefficients, the multiresolution support, computed by DoG-HiT as a prior in a constrained minimization manner: we fit the single-frame (polarimetric) observables by only varying the coefficients in the multiresolution support. The EHT is a VLBI array imaging supermassive black holes. We demonstrate on synthetic data that mr-support imaging offers ample regularization and is able to recover simple geometric dynamics at the horizon scale in a typical EHT setup. The approach is relatively lightweight, fast and largely automatic and data driven. The ngEHT is a planned extension of the EHT designed to recover movies at the event horizon scales of a supermassive black hole. We benchmark the performance of mr-support imaging for the denser ngEHT configuration demonstrating the major improvements the additional ngEHT antennas will bring to dynamic, polarimetric reconstructions. Current and upcoming instruments offer the observational possibility to do polarimetric imaging of dynamically evolving structural patterns with highest spatial and temporal resolution. State-of-the-art dynamic reconstruction methods can capture this motion with a range of temporal regularizers and priors. With this work, we add an additional, simpler regularizer to the list: constraining the reconstruction to the multiresolution support. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: accepted for publication in A&A

Journal ref: A&A 673, A151 (2023)

Multi-scale and Multi-directional VLBI Imaging with CLEAN

Authors: Hendrik Müller, Andrei Lobanov

Abstract: Very long baseline interferometry (VLBI) is a radio-astronomical technique in which the correlated signal from various baselines is combined into an image of highest angular resolution. Due to sparsity of the measurements, this imaging procedure constitutes an ill-posed inverse problem. For decades the CLEAN algorithm was the standard choice in VLBI studies, although having some serious disadvanta… ▽ More Very long baseline interferometry (VLBI) is a radio-astronomical technique in which the correlated signal from various baselines is combined into an image of highest angular resolution. Due to sparsity of the measurements, this imaging procedure constitutes an ill-posed inverse problem. For decades the CLEAN algorithm was the standard choice in VLBI studies, although having some serious disadvantages and pathologies that are challenged by the requirements of modern frontline VLBI applications. We develop a novel multi-scale CLEAN deconvolution method (DoB-CLEAN) based on continuous wavelet transforms that address several pathologies in CLEAN imaging. We benchmark this novel algorithm against CLEAN reconstructions on synthetic data and reanalyze BL Lac observations of RadioAstron with DoB-CLEAN. DoB-CLEAN approaches the image by multi-scalar and multi-directional wavelet dictionaries. Two different dictionaries are used. Firstly, a difference of elliptical spherical Bessel functions dictionary fitted to the uv-coverage of the observation that is used to sparsely represent the features in the dirty image. Secondly, a difference of elliptical Gaussian wavelet dictionary that is well suited to represent relevant image features cleanly. The deconvolution is performed by switching between the dictionaries. DoB-CLEAN achieves super-resolution compared to CLEAN and remedies the spurious regularization properties of CLEAN. In contrast to CLEAN, the representation by basis functions has a physical meaning. Hence, the computed deconvolved image still fits the observed visibilities, opposed to CLEAN. State-of-the-art multi-scalar imaging approaches seem to outperform single-scalar standard approaches in VLBI and are well suited to maximize the extraction of information in ongoing frontline VLBI applications. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 672, A26 (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)

The ngEHT Analysis Challenges

Authors: Freek Roelofs, Lindy Blackburn, Greg Lindahl, Sheperd S. Doeleman, Michael D. Johnson, Philipp Arras, Koushik Chatterjee, Razieh Emami, Christian Fromm, Antonio Fuentes, Jakob Knollmueller, Nikita Kosogorov, Hendrik Mueller, Nimesh Patel, Alexander Raymond, Paul Tiede, Thalia Traianou, Justin Vega

Abstract: The next-generation Event Horizon Telescope (ngEHT) will be a significant enhancement of the Event Horizon Telescope (EHT) array, with $\sim 10$ new antennas and instrumental upgrades of existing antennas. The increased $uv$-coverage, sensitivity, and frequency coverage allow a wide range of new science opportunities to be explored. The ngEHT Analysis Challenges have been launched to inform develo… ▽ More The next-generation Event Horizon Telescope (ngEHT) will be a significant enhancement of the Event Horizon Telescope (EHT) array, with $\sim 10$ new antennas and instrumental upgrades of existing antennas. The increased $uv$-coverage, sensitivity, and frequency coverage allow a wide range of new science opportunities to be explored. The ngEHT Analysis Challenges have been launched to inform development of the ngEHT array design, science objectives, and analysis pathways. For each challenge, synthetic EHT and ngEHT datasets are generated from theoretical source models and released to the challenge participants, who analyze the datasets using image reconstruction and other methods. The submitted analysis results are evaluated with quantitative metrics. In this work, we report on the first two ngEHT Analysis Challenges. These have focused on static and dynamical models of M87* and Sgr A*, and shown that high-quality movies of the extended jet structure of M87* and near-horizon hourly timescale variability of Sgr A* can be reconstructed by the reference ngEHT array in realistic observing conditions, using current analysis algorithms. We identify areas where there is still room for improvement of these algorithms and analysis strategies. Other science cases and arrays will be explored in future challenges. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 32 pages, 14 figures, accepted for publication in Galaxies

Accretion Flow Morphology in Numerical Simulations of Black Holes from the ngEHT Model Library: The Impact of Radiation Physics

Authors: Koushik Chatterjee, Andrew Chael, Paul Tiede, Yosuke Mizuno, Razieh Emami, Christian Fromm, Angelo Ricarte, Lindy Blackburn, Freek Roelofs, Michael D. Johnson, Sheperd S. Doeleman, Philipp Arras, Antonio Fuentes, Jakob Knollmüller, Nikita Kosogorov, Greg Lindahl, Hendrik Müller, Nimesh Patel, Alexander Raymond, Efthalia Traianou, Justin Vega

Abstract: In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) (M87*) and Sagittarius A$^*$ (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecti… ▽ More In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) (M87*) and Sagittarius A$^*$ (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecting the accretion flow and the jet. The ngEHT Analysis Challenges aim to understand the science extractability from synthetic images and movies to inform the ngEHT array design and analysis algorithm development. In this work, we compare the accretion flow structure and dynamics in numerical fluid simulations that specifically target M87* and Sgr A*, and were used to construct the source models in the challenge set. We consider (1) a steady-state axisymmetric radiatively inefficient accretion flow model with a time-dependent shearing hotspot, (2) two time-dependent single fluid general relativistic magnetohydrodynamic (GRMHD) simulations from the H-AMR code, (3) a two-temperature GRMHD simulation from the BHAC code, and (4) a two-temperature radiative GRMHD simulation from the KORAL code. We find that the different models exhibit remarkably similar temporal and spatial properties, except for the electron temperature, since radiative losses substantially cool down electrons near the BH and the jet sheath, signaling the importance of radiative cooling even for slowly accreting BHs such as M87*. We restrict ourselves to standard torus accretion flows, and leave larger explorations of alternate accretion models to future work. △ Less

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

Comments: Accepted in Galaxies; 23 pages, 7 figures

Searching for dilaton fields in the Ly$α$ forest

Authors: Louis Hamaide, Hendrik Müller, David J. E. Marsh

Abstract: Dilatons (and moduli) couple to the masses and coupling constants of ordinary matter, and these quantities are fixed by the local value of the dilaton field. If, in addition, the dilaton with mass $m_φ$ contributes to the cosmic dark matter density, then such quantities oscillate in time at the dilaton Compton frequency. We show how these oscillations lead to broadening and shifting of the Voigt p… ▽ More Dilatons (and moduli) couple to the masses and coupling constants of ordinary matter, and these quantities are fixed by the local value of the dilaton field. If, in addition, the dilaton with mass $m_φ$ contributes to the cosmic dark matter density, then such quantities oscillate in time at the dilaton Compton frequency. We show how these oscillations lead to broadening and shifting of the Voigt profile of the Ly$α$ forest, in a manner that is correlated with the local dark matter density. We further show how tomographic methods allow the effect to be reconstructed by observing the Ly$α$ forest spectrum of distant quasars. We then simulate a large number of quasar lines of sight using the lognormal density field, and forecast the ability of future astronomical surveys to measure this effect. We find that in the ultra low mass range $10^{-32}\text{ eV}\leq m_φ\leq 10^{-28}\text{ eV}$ upcoming observations can improve over existing limits to the dilaton electron mass and fine structure constant couplings set by fifth force searches by up to five orders of magnitude. Our projected limits apply assuming that the ultralight dilaton makes up a few percent of the dark matter density, consistent with upper limits set by the cosmic microwave background anisotropies. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 22 pages, 9 figures, 1 table

Rotational spectroscopy of mono-deuterated oxirane ($c$-C$_2$H$_3$DO) and its detection towards IRAS 16293$-$2422 B

Authors: Holger S. P. Müller, Jes K. Jørgensen, Jean-Claude Guillemin, Frank Lewen, Stephan Schlemmer

Abstract: We prepared a sample of mono-deuterated oxirane and studied its rotational spectrum in the laboratory between 490 GHz and 1060 GHz in order to improve its spectroscopic parameters and consequently the calculated rest frequencies of its rotational transitions. The updated rest frequencies were employed to detect $c$-C$_2$H$_3$DO for the first time in the interstellar medium in the Atacama Large Mil… ▽ More We prepared a sample of mono-deuterated oxirane and studied its rotational spectrum in the laboratory between 490 GHz and 1060 GHz in order to improve its spectroscopic parameters and consequently the calculated rest frequencies of its rotational transitions. The updated rest frequencies were employed to detect $c$-C$_2$H$_3$DO for the first time in the interstellar medium in the Atacama Large Millimetre/submillimetre Array (ALMA) Protostellar Interferometric Line Survey (PILS) of the Class 0 protostellar system IRAS 16293$-$2422. Fits of the detected lines using the rotation diagrams yield a temperature of $T_{\rm rot} = 103 \pm 19$ K, which in turn agrees well with 125 K derived for the $c$-C$_2$H$_4$O main isotopologue previously. The $c$-C$_2$H$_3$DO to $c$-C$_2$H$_4$O ratio is found to be $\sim$0.15 corresponding to a D-to-H ratio of $\sim$0.036 per H atom which is slightly higher than the D-to-H ratio of species such as methanol, formaldehyde, ketene and but lower than those of the larger complex organic species such as ethanol, methylformate and glycolaldehyde. This may reflect that oxirane is formed fairly early in the evolution of the prestellar cores. The identification of doubly deuterated oxirane isotopomers in the PILS data may be possible judged by the amount of mono-deuterated oxirane and the observed trend that multiply deuterated isotopologues have higher deuteration rates than their mono-deuterated variants. △ Less

Submitted 3 September, 2022; originally announced September 2022.

Comments: 10 pages total; Mon. Not. R. Astron. Soc., accepted

Journal ref: Mon. Not. R. Astron. Soc., 518, 185-193 (2023)

Millimeter wave spectrum and search for vinyl isocyanate toward Sgr B2(N) with ALMA

Authors: K. Vávra, L. Kolesniková, A. Belloche, R. T. Garrod, J. Koucký, T. Uhlíková, K. Luková, J. -C. Guillemin, P. Kania, H. S. P. Müller, K. M. Menten, Š. Urban

Abstract: The interstellar detections of isocyanic acid, methyl isocyanate, and very recently also ethyl isocyanate, open the question of the possible detection of vinyl isocyanate in the interstellar medium. The aim of this study is to extend the laboratory rotational spectrum of vinyl isocyanate into the millimeter wave region and to undertake a check for its presence in the high-mass star forming region… ▽ More The interstellar detections of isocyanic acid, methyl isocyanate, and very recently also ethyl isocyanate, open the question of the possible detection of vinyl isocyanate in the interstellar medium. The aim of this study is to extend the laboratory rotational spectrum of vinyl isocyanate into the millimeter wave region and to undertake a check for its presence in the high-mass star forming region Sgr B2. The rotational spectrum of vinyl isocyanate was recorded in the frequency regions 127.5-218 and 285-330 GHz using the Prague millimeter wave spectrometer. The spectral analysis was supported by high-level quantum-chemical calculations. We assumed local thermodynamic equilibrium to compute synthetic spectra of vinyl isocyanate and to search for it in the ReMoCA survey performed with ALMA toward Sgr B2(N). We also searched for ethyl isocyanate in the same source. Accurate values for the rotational and centrifugal distortion constants are reported for the ground vibrational states of trans and cis vinyl isocyanate. We report nondetections of vinyl and ethyl isocyanate toward the main hot core of Sgr B2(N). We find that vinyl and ethyl isocyanate are at least 11 and 3 times less abundant than methyl isocyanate in this source, respectively. Although the precise formation mechanism of interstellar methyl isocyanate itself remains uncertain, we infer from existing astrochemical models that our observational upper limit for the CH3NCO:C2H5NCO ratio in Sgr B2(N) is consistent with ethyl isocyanate being formed on dust grains via the abstraction or photodissociation of an H atom from methyl isocyanate, followed by the addition of a methyl radical. The dominance of such a process for ethyl isocyanate production, combined with the absence of an analogous mechanism for vinyl isocyanate, would indicate that the ratio C2H3NCO:C2H5NCO should be rather less than unity. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A

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

Resolving desorption of complex organic molecules in a hot core: Transition from non-thermal to thermal desorption or two-step thermal desorption?

Authors: Laura A. Busch, Arnaud Belloche, Robin T. Garrod, Holger S. P. Müller, Karl M. Menten

Abstract: Using the high angular resolution provided by the ALMA interferometre we want to resolve the COM emission in the hot molecular core Sagittarius B2(N1) and thereby shed light on the desorption process of Complex Organic Molecules (COMs) in hot cores. We use data taken as part of the 3 mm spectral line survey Re-exploring Molecular Complexity with ALMA (ReMoCA) to investigate the morphology of COM e… ▽ More Using the high angular resolution provided by the ALMA interferometre we want to resolve the COM emission in the hot molecular core Sagittarius B2(N1) and thereby shed light on the desorption process of Complex Organic Molecules (COMs) in hot cores. We use data taken as part of the 3 mm spectral line survey Re-exploring Molecular Complexity with ALMA (ReMoCA) to investigate the morphology of COM emission in Sagittarius B2(N1). Spectra of ten COMs are modelled under the assumption of LTE and population diagrams are derived for positions at various distances to the south and west from the continuum peak. Based on this analysis, resolved COM rotation temperature and COM abundance profiles are derived. Based on the morphology, a rough separation into O- and N-bearing COMs can be done. Temperature profiles are in agreement with expectations of protostellar heating of an envelope with optically thick dust. Abundance profiles reflect a similar trend as seen in the morphology and, to a great extent, agree with results of astrochemical models that, besides the co-desorption with water, predict that O-bearing COMs are mainly formed on dust grain surfaces at low temperatures while at least some N-bearing COMs and CH$_3$CHO are substantially formed in the gas phase at higher temperatures. Our observational results, in comparison with model predictions, suggest that COMs that are exclusively or to a great extent formed on dust grains desorb thermally at ~100 K from the grain surface likely alongside water. Non-zero abundance values below ~100 K suggest that another desorption process is at work at these low temperatures: either non-thermal desorption or partial thermal desorption related to lower binding energies experienced by COMs in the outer, water-poor ice layers. In either case, this is the first time that the transition between two regimes of COM desorption has been resolved in a hot core. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 665, A96 (2022)

DoG-HiT: A novel VLBI Multiscale Imaging Approach

Authors: Hendrik Müller, Andrei Lobanov

Abstract: Reconstructing images from very long baseline interferometry (VLBI) data with sparse sampling of the Fourier domain (uv-coverage) constitutes an ill-posed deconvolution problem. It requires application of robust algorithms maximizing the information extraction from all of the sampled spatial scales and minimizing the influence of the unsampled scales on image quality. We develop a new multiscale w… ▽ More Reconstructing images from very long baseline interferometry (VLBI) data with sparse sampling of the Fourier domain (uv-coverage) constitutes an ill-posed deconvolution problem. It requires application of robust algorithms maximizing the information extraction from all of the sampled spatial scales and minimizing the influence of the unsampled scales on image quality. We develop a new multiscale wavelet deconvolution algorithm DoG-HiT for imaging sparsely sampled interferometric data which combines the difference of Gaussian (DoG) wavelets and hard image thresholding (HiT). Based on DoG-HiT, we propose a multi-step imaging pipeline for analysis of interferometric data. DoG-HiT applies the compressed sensing approach to imaging by employing a flexible DoG wavelet dictionary which is designed to adapt smoothly to the uv-coverage. It uses closure properties as data fidelity terms only initially and perform non-convex, non-smooth optimization by an amplitude conserving and total flux conserving hard thresholding splitting. DoG-HiT calculates a multiresolution support as a side product. The final reconstruction is refined through self-calibration loops and imaging with amplitude and phase information applied for the multiresolution support only. We demonstrate the stability of DoG-HiT and benchmark its performance against image reconstructions made with CLEAN and Regularized Maximum-Likelihood (RML) methods using synthetic data. The comparison shows that DoG-HiT matches the superresolution achieved by the RML reconstructions and surpasses the sensitivity to extended emission reached by CLEAN. Application of regularized maximum likelihood methods outfitted with flexible multiscale wavelet dictionaries to imaging of interferometric data matches the performance of state-of-the art convex optimization imaging algorithms and requires fewer prior and user defined constraints. △ Less

Submitted 19 June, 2022; originally announced June 2022.

Comments: 20 pages, 11 Figures, accepted for publication in A&A, the imaging software underlying this work will be made publicly available soon under the link https://github.com/hmuellergoe/mrbeam

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

Rotational spectroscopy $n$-propanol: $Aa$ and $Ag$ conformers

Authors: O. Zingsheim, J. Maßen, H. S. P. Müller, B. Heyne, M. Fatima, L. Bonah, A. Belloche, F. Lewen, S. Schlemmer

Abstract: The primary alcohol $n$-propanol (i.e., $normal$-propanol or propan-1-ol; C$_3$H$_7$OH) occurs in five different conformers: $Ga$, $Gg$, $Gg'$, $Aa$, and $Ag$. All rotational spectra of the three conformers of the $G$ family are well described, making astronomical search of their spectroscopic signatures possible, as opposed to those of the $Aa$ and $Ag$ conformers. Our goal is to facilitate the a… ▽ More The primary alcohol $n$-propanol (i.e., $normal$-propanol or propan-1-ol; C$_3$H$_7$OH) occurs in five different conformers: $Ga$, $Gg$, $Gg'$, $Aa$, and $Ag$. All rotational spectra of the three conformers of the $G$ family are well described, making astronomical search of their spectroscopic signatures possible, as opposed to those of the $Aa$ and $Ag$ conformers. Our goal is to facilitate the astronomical detection of $Aa$ and $Ag$ conformers of $n$-propanol by characterizing their rotational spectra. We recorded the rotational spectra of $n$-propanol in the frequency domain of 18$-$505\,GHz.Additional double-modulation double-resonance (DM-DR) measurements were performed, more specifically with the goal to unambiguously assign weak transitions of the $Aa$ conformer and to verify assignments of the $Ag$ conformer. We derived a spectroscopic quantum mechanical model with experimental accuracy (with $J_\textrm{max}=70$ and $K_{a,\textrm{max}}=6$) for $Aa$ $n$-propanol. Furthermore, we unambiguously assigned transitions (with $J_\textrm{max}=69$ and $K_{a,\textrm{max}}=9$) of $Ag$ $n$-propanol; in doing so, we prove the existence of two tunneling states, $Ag^+$ and $Ag^-$. The astronomical search of all five conformers of $n$-propanol is now possible via their rotational signatures. These are applied in a companion article on the detection of $n$-propanol toward the hot molecular core Sgr B2(N2). △ Less

Submitted 7 June, 2022; originally announced June 2022.

Journal ref: A&A 662, A111 (2022)

Interstellar detection and chemical modeling of iso-propanol and its normal isomer

Authors: A. Belloche, R. T. Garrod, O. Zingsheim, H. S. P. Müller, K. M. Menten

Abstract: The detection of a branched alkyl molecule in the high-mass star forming protocluster Sgr B2(N) permitted by the advent of ALMA revealed a new dimension of interstellar chemistry. Astrochemical simulations subsequently predicted that beyond a certain degree of molecular complexity, branched molecules could even dominate over their straight-chain isomers. More generally, we aim at probing further t… ▽ More The detection of a branched alkyl molecule in the high-mass star forming protocluster Sgr B2(N) permitted by the advent of ALMA revealed a new dimension of interstellar chemistry. Astrochemical simulations subsequently predicted that beyond a certain degree of molecular complexity, branched molecules could even dominate over their straight-chain isomers. More generally, we aim at probing further the presence in the ISM of complex organic molecules with the capacity to exhibit both a normal and iso form, via the attachment of a functional group to either a primary or secondary carbon atom. We used the imaging spectral line survey survey ReMoCA performed with ALMA and the results of a recent spectroscopic study of propanol to search for the iso and normal isomers of this molecule in the hot molecular core Sgr B2(N2). We expanded the network of the astrochemical model MAGICKAL to explore the formation routes of propanol. We report the first interstellar detection of iso-propanol toward a position of Sgr B2(N2) that shows narrow linewidths. We also report the first secure detection of normal-propanol in a hot core. i-Propanol is found to be nearly as abundant as n-propanol, with an abundance ratio of 0.6 similar to the ratio of 0.4 that we obtained previously for i- and n-propyl cyanide in Sgr B2(N2). The results are in good agreement with the outcomes of our astrochemical models, which indicate that OH-radical addition to propylene in dust-grain ice mantles, driven by water photodissociation, can produce appropriate quantities of n- and i-propanol. The n-to-i ratio in Sgr B2(N2) may be a direct inheritance of the branching ratio of this reaction process. The detection of n- and i-propanol and their ratio indicate that the modest preference toward the normal form of propyl cyanide determined previously may be a more general feature among similarly sized interstellar molecules. [abridged] △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 662, A110 (2022)

The ALMA-PILS survey: First tentative detection of 3-hydroxypropenal (HOCHCHCHO) in the interstellar medium and chemical modeling of the C$_3$H$_4$O$_2$ isomers

Authors: A. Coutens, J. -C. Loison, A. Boulanger, E. Caux, H. S. P. Müller, V. Wakelam, S. Manigand, J. K. Jørgensen

Abstract: Characterizing the molecular composition of solar-type protostars is useful for improving our understanding of the physico-chemical conditions under which the Sun and its planets formed. In this work, we analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) data of the Protostellar Interferometric Line Survey (PILS), an unbiased spectral survey of the solar-type protostar IRAS~16293--24… ▽ More Characterizing the molecular composition of solar-type protostars is useful for improving our understanding of the physico-chemical conditions under which the Sun and its planets formed. In this work, we analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) data of the Protostellar Interferometric Line Survey (PILS), an unbiased spectral survey of the solar-type protostar IRAS~16293--2422, and we tentatively detected 3-hydroxypropenal (HOCHCHCHO) for the first time in the interstellar medium towards source B. Based on the observed line intensities and assuming local thermodynamic equilibrium, its column density is constrained to be $\sim$10$^{15}$ cm$^{-2}$, corresponding to an abundance of 10$^{-4}$ relative to methanol, CH$_3$OH. Additional spectroscopic studies are needed to constrain the excitation temperature of this molecule. We included HOCHCHCHO and five of its isomers in the chemical network presented in Manigand et al. (2021) and we predicted their chemical evolution with the Nautilus code. The model reproduces the abundance of HOCHCHCHO within the uncertainties. This species is mainly formed through the grain surface reaction CH$_2$CHO + HCO $\rightarrow$ HCOCH$_2$CHO, followed by the tautomerization of HCOCH$_2$CHO into HOCHCHCHO. Two isomers, CH$_3$COCHO and CH$_2$COHCHO, are predicted to be even more abundant than HOCHCHCHO. Spectroscopic studies of these molecules are essential in searching for them in IRAS~16293--2422 and other astrophysical sources. △ Less

Submitted 26 March, 2022; originally announced March 2022.

Comments: Accepted in A&A Letters

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

Toward the limits of complexity of interstellar chemistry: Rotational spectroscopy and astronomical search for n- and i-butanal

Authors: M. Sanz-Novo, A. Belloche, V. M. Rivilla, R. T. Garrod, J. L. Alonso, P. Redondo, C. Barrientos, L. Kolesniková, J. C. Valle, L. Rodríguez-Almeida, I. Jiménez-Serra, J. Martín-Pintado, H. S. P. Muller, K. Menten

Abstract: In recent times, large organic molecules of exceptional complexity have been found in diverse regions of the interstellar medium. In this context, we aim to provide accurate frequencies of the ground vibrational state of two key aliphatic aldehydes, n-butanal and its branched-chain isomer, i-butanal, to enable their eventual detection in the interstellar medium. We employ a frequency modulation mi… ▽ More In recent times, large organic molecules of exceptional complexity have been found in diverse regions of the interstellar medium. In this context, we aim to provide accurate frequencies of the ground vibrational state of two key aliphatic aldehydes, n-butanal and its branched-chain isomer, i-butanal, to enable their eventual detection in the interstellar medium. We employ a frequency modulation millimeter-wave absorption spectrometer to measure the rotational features of n- and i-butanal. We use the spectral line survey ReMoCA performed with the Atacama Large Millimeter/submillimeter Array to search for n- and i-butanal toward the star-forming region Sgr B2(N). We also search for both aldehydes toward the molecular cloud G+0.693-0.027 with IRAM 30 m and Yebes 40 m observations. Several thousand rotational transitions belonging to the lowest-energy conformers have been assigned in the laboratory spectra up to 325 GHz. A precise set of the relevant rotational spectroscopic constants has been determined for each structure. We report non-detections of n- and i-butanal toward both sources, Sgr B2(N1S) and G+0.693-0.027. We find that n- and i-butanal are at least 2-6 and 6-18 times less abundant than acetaldehyde toward Sgr B2(N1S), respectively, and that n-butanal is at least 63 times less abundant than acetaldehyde toward G+0.693-0.027. Comparison with astrochemical models indicates good agreement between observed and simulated abundances (where available). Grain-surface chemistry appears sufficient to reproduce aldehyde ratios in G+0.693-0.027; gas-phase production may play a more active role in Sgr B2(N1S). Our astronomical results indicate that the family of interstellar aldehydes in the Galactic center region is characterized by a drop of one order of magnitude in abundance at each incrementation in the level of molecular complexity. △ Less

Submitted 14 March, 2022; originally announced March 2022.

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

Recoil imaging for directional detection of dark matter, neutrinos, and physics beyond the Standard Model

Authors: C. A. J. O'Hare, D. Loomba, K. Altenmüller, H. Álvarez-Pol, F. D. Amaro, H. M. Araújo, D. Aristizabal Sierra, J. Asaadi, D. Attié, S. Aune, C. Awe, Y. Ayyad, E. Baracchini, P. Barbeau, J. B. R. Battat, N. F. Bell, B. Biasuzzi, L. J. Bignell, C. Boehm, I. Bolognino, F. M. Brunbauer, M. Caamaño, C. Cabo, D. Caratelli, J. M. Carmona , et al. (142 additional authors not shown)

Abstract: Recoil imaging entails the detection of spatially resolved ionization tracks generated by particle interactions. This is a highly sought-after capability in many classes of detector, with broad applications across particle and astroparticle physics. However, at low energies, where ionization signatures are small in size, recoil imaging only seems to be a practical goal for micro-pattern gas detect… ▽ More Recoil imaging entails the detection of spatially resolved ionization tracks generated by particle interactions. This is a highly sought-after capability in many classes of detector, with broad applications across particle and astroparticle physics. However, at low energies, where ionization signatures are small in size, recoil imaging only seems to be a practical goal for micro-pattern gas detectors. This white paper outlines the physics case for recoil imaging, and puts forward a decadal plan to advance towards the directional detection of low-energy recoils with sensitivity and resolution close to fundamental performance limits. The science case covered includes: the discovery of dark matter into the neutrino fog, directional detection of sub-MeV solar neutrinos, the precision study of coherent-elastic neutrino-nucleus scattering, the detection of solar axions, the measurement of the Migdal effect, X-ray polarimetry, and several other applied physics goals. We also outline the R&D programs necessary to test concepts that are crucial to advance detector performance towards their fundamental limit: single primary electron sensitivity with full 3D spatial resolution at the $\sim$100 micron-scale. These advancements include: the use of negative ion drift, electron counting with high-definition electronic readout, time projection chambers with optical readout, and the possibility for nuclear recoil tracking in high-density gases such as argon. We also discuss the readout and electronics systems needed to scale-up such detectors to the ton-scale and beyond. △ Less

Submitted 17 July, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: 77 pages, 20 figures. Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

Rotational spectroscopy of isotopic oxirane, $c$-C$_2$H$_4$O

Authors: Holger S. P. Müller, Jean-Claude Guillemin, Frank Lewen, Stephan Schlemmer

Abstract: We studied the rotational spectrum of oxirane in a sample of natural isotopic composition in selected regions between 158 GHz and 1093 GHz. Investigations of the isotopologs with one $^{13}$C or one $^{18}$O were the primary focus in order to facilitate searches for them in space. We also examined the main isotopic species mainly to look into the performance of Watson's A and S reductions both in… ▽ More We studied the rotational spectrum of oxirane in a sample of natural isotopic composition in selected regions between 158 GHz and 1093 GHz. Investigations of the isotopologs with one $^{13}$C or one $^{18}$O were the primary focus in order to facilitate searches for them in space. We also examined the main isotopic species mainly to look into the performance of Watson's A and S reductions both in an oblate and in a prolate representation. Even though oxirane is a rather asymmetric oblate rotor, the A reduction in the III$^l$ representation did not yield a satisfactory fit, as was observed frequently earlier for other molecules. The other three combinations yielded satisfactory fits of similar quality among each other; the A reduction in the I$^r$ representation required two parameters less than both S reduction fits. △ Less

Submitted 16 February, 2022; v1 submitted 23 January, 2022; originally announced January 2022.

Comments: 11 pages, including tables and figure; minor update; J. Mol. Spectrosc., published

Journal ref: J. Mol. Spectrosc. 384 (2022) Art. No. 111584

(Sub-)millimeter-wave spectroscopy of gauche-propanal

Authors: Oliver Zingsheim, Holger S. P. Müller, Luis Bonah, Frank Lewen, Sven Thorwirth, Stephan Schlemmer

Abstract: A detailed analysis of (sub-)millimeter-wave spectra of the vibrational ground state ($\upsilon=0$) combined with the energetically lowest excited vibrational state ($\upsilon_{24}=1$; aldehyde torsion) of gauche-propanal (g-C$_2$H$_5$CHO) up to 500 GHz is presented. Both vibrational states, $\upsilon=0$ and $\upsilon_{24}=1$, are treated with tunneling rotation interactions between their two resp… ▽ More A detailed analysis of (sub-)millimeter-wave spectra of the vibrational ground state ($\upsilon=0$) combined with the energetically lowest excited vibrational state ($\upsilon_{24}=1$; aldehyde torsion) of gauche-propanal (g-C$_2$H$_5$CHO) up to 500 GHz is presented. Both vibrational states, $\upsilon=0$ and $\upsilon_{24}=1$, are treated with tunneling rotation interactions between their two respective tunneling states, which originate from two stable degenerate gauche-conformers; left- and right-handed configurations separated by a small potential barrier. Thanks to double-modulation double-resonance (DM-DR) measurements, important but weak $c$-type transitions connecting the tunneling states could be unambiguously assigned. In addition, Coriolis interaction as well as Fermi resonance between the two vibrational states needed to be taken into account to derive fits with experimental accuracy using Pickett's SPFIT program in a reduced axis system (RAS). Based on the rotational analysis, the fundamental vibrational frequency $ν_{24}$ of gauche-propanal is redetermined to 68.75037(30) cm$^{-1}$. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: accepted manuscript by Journal of Molecular Spectroscopy; 16 pages (11 main paper + 5 SM); 12 Figures (8+4)

Journal ref: J. Mol. Spectrosc. 384 (2022) Art. No. 111565

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)

Laboratory rotational spectroscopy of acrylamide and search for acrylamide and propionamide toward Sgr B2(N) with ALMA

Authors: L. Kolesniková, A. Belloche, J. Koucký, E. R. Alonso, R. T. Garrod, K. Luková, K. M. Menten, H. S. P. Müller, P. Kania, Š. Urban

Abstract: Numerous complex organic molecules have been detected in the universe among which amides are considered as models for species containing the peptide linkage. Acrylamide bears in its backbone not only the peptide bond, but also the vinyl functional group which is a common motif in many interstellar compounds. This makes acrylamide an interesting target for a search in space. In addition, a tentativ… ▽ More Numerous complex organic molecules have been detected in the universe among which amides are considered as models for species containing the peptide linkage. Acrylamide bears in its backbone not only the peptide bond, but also the vinyl functional group which is a common motif in many interstellar compounds. This makes acrylamide an interesting target for a search in space. In addition, a tentative detection of the related molecule propionamide has been recently claimed toward Sgr B2(N). We report accurate laboratory measurements and analyses of thousands rotational transitions for the syn and skew forms of acrylamide between 75 and 480 GHz. Tunneling through a low energy barrier between two symmetrically equivalent configurations has been revealed for the higher-energy skew species. We searched for emission of acrylamide in the imaging spectral line survey ReMoCA performed with ALMA toward Sgr B2(N). We also searched for propionamide in the same source. Neither acrylamide nor propionamide were detected toward the two main hot molecular cores of Sgr B2(N). We did not detect propionamide either toward a position located to the east of the main hot core, thereby not confirming the recent claim of its interstellar detection toward this position. We find that acrylamide and propionamide are at least 26 and 14 times, respectively, less abundant than acetamide toward the main hot core Sgr B2(N1S), and at least 6 and 3 times, respectively, less abundant than acetamide toward the secondary hot core Sgr B2(N2). A comparison with results of astrochemical kinetics model for related species suggests that acrylamide may be a few hundred times less abundant than acetamide, corresponding to a value at least an order of magnitude lower than the observational upper limits. Propionamide may be as little as only a factor of two less abundant than the upper limit derived toward Sgr B2(N1S). △ Less

Submitted 5 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A

Journal ref: Astron. Astrophys., 659, A111 (2022)

Rotational and rovibrational spectroscopy of CD$_3$OH with an account of CD$_3$OH toward IRAS 16293$-$2422

Authors: V. V. Ilyushin, H. S. P. Müller, J. K. Jørgensen, S. Bauerecker, C. Maul, Y. Bakhmat, E. A. Alekseev, O. Dorovskaya, S. Vlasenko, F. Lewen, S. Schlemmer, K. Berezkin, R. M. Lees

Abstract: Solar-type protostars harbor highly deuterated complex organics. While this degree of deuteration may provide important clues in studying the formation of these species, spectroscopic information on multiply deuterated isotopologs is often insufficient. In particular, searches for triply deuterated methanol, CD$_3$OH, are hampered by the lack of intensity information from a spectroscopic model. Th… ▽ More Solar-type protostars harbor highly deuterated complex organics. While this degree of deuteration may provide important clues in studying the formation of these species, spectroscopic information on multiply deuterated isotopologs is often insufficient. In particular, searches for triply deuterated methanol, CD$_3$OH, are hampered by the lack of intensity information from a spectroscopic model. The aim of the present study is to develop such a model of CD$_3$OH in low-lying torsional states that is sufficiently accurate to facilitate further searches for CD$_3$OH in space. We performed a new measurement campaign for CD$_3$OH involving three spectroscopic laboratories that covers the 34 GHz-1.1 THz and the 20-900 cm$^{-1}$ ranges. The analysis was performed using the torsion-rotation Hamiltonian model based on the rho-axis method. We determined a model that describes the ground and first excited torsional states of CD$_3$OH, up to quantum numbers $J \leqslant 55$ and $K_a \leqslant 23$, and we derived a line list for radio-astronomical observations. This list is accurate up to at least 1.1 THz and should be sufficient for all types of radio-astronomical searches for this methanol isotopolog. It was used to search for CD$_3$OH in data from the Protostellar Interferometric Line Survey of IRAS 16293-2422 using ALMA. CD$_3$OH is securely detected in the data, with a large number of clearly separated and well-reproduced lines. We detected lines belonging to the ground and the first excited torsional states. The derived abundance of CD$_3$OH relative to non-deuterated isotopolog confirm the significant enhancement of this multiply deuterated variant. This finding is in line with other observations of multiply deuterated complex organic molecules and may serve as an important constraint on their formation models. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: Astronomy and Astrophysics, in press. 21 pages, including figures, tables, and appendices; abstract shortened

Journal ref: A&A 658, A127 (2022)

Millimetre-wave laboratory study of glycinamide and search for it with ALMA toward Sagittarius B2(N)

Authors: Z. Kisiel, L. Kolesniková, A. Belloche, J. -C. Guillemin, L. Pszczółkowski, E. R. Alonso, R. T. Garrod, E. Białkowska-Jaworska, I. León, H. S. P. Müller, K. M. Menten, J. L. Alonso

Abstract: Glycinamide is considered to be one of the possible precursors of the simplest amino acid glycine. Its only rotational spectrum reported so far has been in the cm-wave region. The aim of this work is to extend its laboratory spectrum into the mm wave region to support its searches in the ISM. Glycinamide was synthesised chemically and was studied with broadband rotational spectroscopy in the 90-32… ▽ More Glycinamide is considered to be one of the possible precursors of the simplest amino acid glycine. Its only rotational spectrum reported so far has been in the cm-wave region. The aim of this work is to extend its laboratory spectrum into the mm wave region to support its searches in the ISM. Glycinamide was synthesised chemically and was studied with broadband rotational spectroscopy in the 90-329 GHz region. Tunneling across a low energy barrier between two symmetry equivalent configurations of the molecule resulted in splitting of each vibrational state and many perturbations in associated rotational energy levels, requiring careful coupled state fits for each vibrational doublet. We searched for emission of glycinamide in the imaging spectral line survey ReMoCA performed with ALMA toward Sgr B2(N). We report the first analysis of the mm-wave rotational spectrum of glycinamide, resulting in fitting to experimental measurement accuracy of over 1200 transition frequencies for the ground state tunneling doublet, of many lines for tunneling doublets for two singly excited vibrational states, and determination of precise vibrational separation in each doublet. We did not detect emission from glycinamide in the hot core Sgr B2(N1S). We found that glycinamide is at least seven times less abundant than aminoacetonitrile and 1.8 times less abundant than urea in this source. A comparison with results of astrochemical kinetics models for species related to glycinamide suggests that its abundance may be at least one order of magnitude below the upper limit obtained toward Sgr B2(N1S). This means that glycinamide emission in this source likely lies well below the spectral confusion limit in the frequency range covered by the ReMoCA survey. Targetting sources with a lower level of spectral confusion, such as the Galactic Center shocked region G+0.693-0.027, may be a promising avenue. [abridged] △ Less

Submitted 22 October, 2021; originally announced October 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 657, A99 (2022)

Orbital obliquity sampling in the Kepler-20 system using the 3D animation software Blender

Authors: Holger M. Müller, Panagiotis Ioannidis, Jürgen H. M. M. Schmitt

Abstract: The mutual orbital alignment in multiple planetary systems is an important parameter for understanding their formation. There are a number of elaborate techniques to determine the alignment parameters using photometric or spectroscopic data. Planet--planet occultations (PPOs), which can occur in multiple transiting systems, are one intuitive example. While the presence of PPOs constrains the orbit… ▽ More The mutual orbital alignment in multiple planetary systems is an important parameter for understanding their formation. There are a number of elaborate techniques to determine the alignment parameters using photometric or spectroscopic data. Planet--planet occultations (PPOs), which can occur in multiple transiting systems, are one intuitive example. While the presence of PPOs constrains the orbital alignment, the absence at first glance does not. Planetary systems, for which the measurement of orbital obliquities with conventional techniques remains elusive, call for new methods whereby at least some information on the alignments can be obtained. Here we develop a method that uses photometric data to gain this kind of information from multi-transit events. In our approach we synthesize multi-transit light curves of the exoplanets in question via the construction of a grid of projected orbital tilt angles $α$, while keeping all transit parameters constant. These model light curves contain PPOs for some values of $α$. To compute the model light curves, we use the 3D animation software Blender for our transit simulations, which allows the use of arbitrary surface brightness distributions of the star, such as limb darkening from model atmospheres. The resulting model light curves are then compared to actual measurements. We present a detailed study of the multi-transiting planetary system Kepler-20, including parameter fits of the transiting planets and an analysis of the stellar activity. We apply our method to Kepler-20 b and c, where we are able to exclude some orbital geometries, and find a tendency of these planets to eclipse in front of different stellar hemispheres in a prograde direction. [...] △ Less

Submitted 18 October, 2021; originally announced October 2021.

Comments: 14 pages, 8 figures, 6 tables. Accepted for publication in Astronomy & Astrophysics (A&A)

Journal ref: A&A 657, A37 (2022)

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)

New Observational Constraints on the Winds of M Dwarf Stars

Authors: Brian E. Wood, Hans-Reinhard Mueller, Seth Redfield, Fallon Konow, Hunter Vannier, Jeffrey L. Linsky, Allison Youngblood, Aline A. Vidotto, Moira Jardine, Julian D. Alvarado-Gomez, Jeremy J. Drake

Abstract: High resolution UV spectra of stellar H I Lyman-alpha lines from the Hubble Space Telescope (HST) provide observational constraints on the winds of coronal main sequence stars, thanks to an astrospheric absorption signature created by the interaction between the stellar winds and the interstellar medium. We report the results of a new HST survey of M dwarf stars, yielding six new detections of ast… ▽ More High resolution UV spectra of stellar H I Lyman-alpha lines from the Hubble Space Telescope (HST) provide observational constraints on the winds of coronal main sequence stars, thanks to an astrospheric absorption signature created by the interaction between the stellar winds and the interstellar medium. We report the results of a new HST survey of M dwarf stars, yielding six new detections of astrospheric absorption. We estimate mass-loss rates for these detections, and upper limits for nondetections. These new constraints allow us to characterize the nature of M dwarf winds and their dependence on coronal activity for the first time. For a clear majority of the M dwarfs, we find winds that are weaker or comparable in strength to that of the Sun, i.e. Mdot<=1 Mdot_sun. However, two of the M dwarfs have much stronger winds: YZ CMi (M4 Ve; Mdot=30 Mdot_sun) and GJ 15AB (M2 V+M3.5 V; Mdot=10 Mdot_sun). Even these winds are much weaker than expectations if the solar relation between flare energy and coronal mass ejection (CME) mass extended to M dwarfs. Thus, the solar flare/CME relation does not appear to apply to M dwarfs, with important ramifications for the habitability of exoplanets around M dwarfs. There is evidence for some increase in Mdot with coronal activity as quantified by X-ray flux, but with much scatter. One or more other factors must be involved in determining wind strength besides spectral type and coronal activity, with magnetic topology being one clear possibility. △ Less

Submitted 30 April, 2021; originally announced May 2021.

Comments: 40 pages, 11 figures, to appear in The Astrophysical Journal