Abstract visibility
Abstract text size

Papers for Wednesday, May 04 2022

Papers with local authors

David Rubin, G. Aldering, P. Antilogus, C. Aragon, S. Bailey, C. Baltay, S. Bongard, K. Boone, C. Buton, Y. Copin, S. Dixon, D. Fouchez, E. Gangler, R. Gupta, B. Hayden, W. Hillebrandt, A. G. Kim, M. Kowalski, D. Kuesters, P.-F. Leget, F. Mondon, J. Nordin, R. Pain, E. Pecontal, R. Pereira, S. Perlmutter, K. A. Ponder, D. Rabinowitz, M. Rigault, K. Runge, C. Saunders, G. Smadja, N. Suzuki, C. Tao, S. Taubenberger, R. C. Thomas, M. Vincenzi, The Nearby Supernova Factory

Under review at ApJS

0 votes
Paper 8 — arXiv:2205.01116
0 votes
Paper 8 — arXiv:2205.01116

We calibrate spectrophotometric optical spectra of 32 stars commonly used as standard stars, referenced to 14 stars already on the HST-based CALSPEC flux system. Observations of CALSPEC and non-CALSPEC stars were obtained with the SuperNova Integral Field Spectrograph over the wavelength range 3300 A to 9400 A as calibration for the Nearby Supernova Factory cosmology experiment. In total, this analysis used 4289 standard-star spectra taken on photometric nights. As a modern cosmology analysis, all pre-submission methodological decisions were made with the flux scale and external comparison results blinded. The large number of spectra per star allows us to treat the wavelength-by-wavelength calibration for all nights simultaneously with a Bayesian hierarchical model, thereby enabling a consistent treatment of the Type Ia supernova cosmology analysis and the calibration on which it critically relies. We determine the typical per-observation repeatability (median 14 mmag for exposures >~ 5 s), the Maunakea atmospheric transmission distribution (median dispersion of 7 mmag with uncertainty 1 mmag), and the scatter internal to our CALSPEC reference stars (median of 8 mmag). We also check our standards against literature filter photometry, finding generally good agreement over the full 12-magnitude range. Overall, the mean of our system is calibrated to the mean of CALSPEC at the level of ~ 3 mmag. With our large number of observations, careful crosschecks, and 14 reference stars, our results are the best calibration yet achieved with an integral-field spectrograph, and among the best calibrated surveys.

All other papers

David Leiser (1), Stefan Loehle (1), Stefanos Fasoulas (2) ((1) High Enthalpy Flow Diagnostics Group, Institute of Space Systems, University of Stuttgart, Germany, (2) Institute of Space Systems, University of Stuttgart, Germany)

submitted to Journal of Spacecraft and Rockets

The fragmentation of two aerospace aluminum alloys is investigated in a ground testing facility including mechanical loads as occurring due to aerodynamic forces in a real atmospheric entry event at three trajectory points. The emission spectroscopic analysis shows that these materials fail after distinct alkali metal features are observed in the spectra. The two alloys feature characteristic emissions of the different alkali metals. The presence of lithium lines that have previously been exclusively attributed to battery failure in observation campaigns may be considered as a marker for aluminum breakup. This is particularly interesting for future entry observations because it allows a new insight into the structural failure processes of the demising spacecraft. The lack of emission of alloying elements points to these spectra being a candidate for the determination of spacecraft demise. The identification of such features in ground testing will allow a more certain identification of specific break-up events

Aldo G. Sepulveda, Daniel Huber, Zhoujian Zhang, Gang Li, Michael C. Liu, Timothy R. Bedding

9 pages, 4 figures, 1 table. Submitted to ApJ Letters

51 Eri is well known for hosting a directly-imaged giant planet and for its membership to the $\beta$ Pictoris moving group. Using two-minute cadence photometry from the Transiting Exoplanet Survey Satellite (TESS), we detect multi-periodic variability in 51 Eri that is consistent with pulsations of Gamma Doradus ($\gamma$ Dor) stars. We identify the most significant pulsation modes (with frequencies between $\sim$0.5-3.9 cycles/day and amplitudes ranging between $\sim$1-2 mmag) as dipole and quadrupole gravity-modes, as well as Rossby modes, as previously observed in Kepler $\gamma$ Dor stars. Our results demonstrate that previously reported variability attributed to stellar rotation is instead likely due to $\gamma$ Dor pulsations. Using the mean frequency of the $l=1$ gravity-modes, together with empirical trends of the Kepler $\gamma$ Dor population, we estimate a plausible stellar core rotation period of 0.9$^{+0.3}_{-0.1}$ days for 51 Eri. We find no significant evidence for transiting companions around 51 Eri in the residual light curve. The detection of $\gamma$ Dor pulsations presented here, together with follow-up observations and modeling, may enable the determination of an asteroseismic age for this benchmark system. Future TESS observations would allow a constraint on the stellar core rotation rate, which in turn traces the surface rotation rate, and thus would help clarify whether or not the stellar equatorial plane and orbit of 51 Eri b are coplanar.

Enrique Lopez-Rodriguez, Sui Ann Mao, Rainer Beck, Alejandro S. Borlaff, Evangelia Ntormousi, Konstantinos Tassis, Daniel A. Dale, Julia Roman-Duval, Kandaswamy Subramanian, Sergio Martin-Alvarez, Pamela M. Marcum, Susan E. Clark, William T. Reach, Doyal A. Harper, Ellen G. Zweibel

42 pages, 28 figures. Submitted to ApJ (Comments are welcome)

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $\mu$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{d} = [19,48]$ K, and dense, $\log_{10}(N_{HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154\mu m} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $\mu$m. The falling $53-154$ $\mu$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy.

D. Eckert, S. Ettori, E. Pointecouteau, R. F. J. van der Burg, S. I. Loubser

Accepted for publication in A&A

The mass profiles of massive dark matter halos are highly sensitive to the nature of dark matter and potential modifications of the theory of gravity on large scales. The $\Lambda$CDM paradigm makes strong predictions on the shape of dark matter halos and on the dependence of the shape parameters on halo mass, such that any deviation from the predicted universal shape would have important implications for the fundamental properties of dark matter. Here we use a set of 12 galaxy clusters with available deep X-ray and Sunyaev-Zeldovich data to constrain the shape of the gravitational field with an unprecedented level of precision over two decades in radius. On average, we find that the NFW profile provides an excellent description of the recovered mass profiles, with deviations of less than 10% over a wide radial range. However, there appears to be more diversity in the shape of individual profiles than can be captured by the NFW model. The average NFW concentration and its scatter agree very well with the prediction of the $\Lambda$CDM framework. For a subset of systems, we disentangle the gravitational field into the contribution of baryonic components (gas, brightest cluster galaxy, and satellite galaxies) and that of dark matter. The stellar content dominates the gravitational field inside $\sim0.02R_{500}$ but is responsible for only 1-2% of the total gravitational field inside $R_{200}$. The total baryon fraction reaches the cosmic value at $R_{200}$ and slightly exceeds it beyond this point, possibly indicating a mild level of nonthermal pressure support ($10-20\%$) in cluster outskirts. Finally, the relation between observed and baryonic acceleration exhibits a complex shape that strongly departs from the radial acceleration relation in spiral galaxies, which shows that the aforementioned relation does not hold at the galaxy-cluster scale.

L. G. Bouma, R. Kerr, J. L. Curtis, H. Isaacson, L. A. Hillenbrand, A. W. Howard, A. L. Kraus, A. Bieryla, D. W. Latham, E. A Petigura, D. Huber

AAS journals submitted. 22 pages, 6 figures, 3 appendices, 4 tables. Machine readable versions of tables 2 and 3 uploaded

Stellar positions and velocities from Gaia are yielding a new view of open cluster dispersal. Here we present an analysis of a group of stars spanning Cepheus to Hercules, hereafter the Cep-Her complex. The group includes four Kepler Objects of Interest: Kepler-1643 b ($2.32 \pm 0.14$ Earth-radii, 5.3 day orbital period), KOI-7368 b ($2.22 \pm 0.12$ Earth-radii, 6.8 days), KOI-7913 Ab ($2.34 \pm 0.18$ Earth-radii, 24.2 days), and Kepler-1627 Ab ($3.85 \pm 0.11$ Earth-radii, 7.2 days). The latter Neptune-sized planet is in part of the Cep-Her complex called the $\delta$ Lyr cluster (Bouma et al. 2022). Here we focus on the former three systems, which are in other regions of the association. Based on kinematic evidence from Gaia, stellar rotation periods from TESS, and spectroscopy, these three objects are also approximately 40 million years (Myr) old. More specifically, we find that Kepler-1643 is $46^{+9}_{-7}$ Myr old, based on its membership in a dense sub-cluster of the complex called RSG-5. KOI-7368 and KOI-7913 are $36^{+10}_{-8}$ Myr old, and are in a diffuse region that we call CH-2. Based on the transit shapes and high resolution imaging, all three objects are most likely planets, with false positive probabilities of $6 \times 10^{-9}$, $4 \times 10^{-3}$, and $1 \times 10^{-4}$ for Kepler-1643, KOI-7368, and KOI-7913 respectively. These planets demonstrate that mini-Neptunes with sizes of approximately 2 Earth radii exist at ages of 40 million years.

E. Bertola, C. Vignali, G. Lanzuisi, M. Dadina, M. Cappi, R. Gilli, G. A. Matzeu, G. Chartas, E. Piconcelli, A. Comastri

Accepted for publication in A&A. 14 pages, 4 figures, 6 tables

We present new joint XMM-Newton and NuSTAR observations of APM 08279+5255, a gravitationally-lensed, broad-absorption line quasar ($z=3.91$). After showing a fairly stable flux ($f_{\rm2-10}\simeq4-5.5\times10^{-13}\rm~erg~s^{-1}$) from 2000 to 2008, APM 08279+5255 was found in a fainter state in the latest X-ray exposures ($f_{\rm2-10}\simeq2.7\times10^{-13}\rm~erg~s^{-1}$), which can likely be ascribed to a lower X-ray activity. Moreover, the 2019 data present a prominent Fe K$\alpha$ emission line and do not show any significant absorption line. This fainter state, coupled to the first hard X-ray sampling of APM 08279+5255, allowed us to measure X-ray reflection and the high-energy cutoff in this source for the first time. From the analysis of previous XMM-Newton and Chandra observations, X-ray reflection is demonstrated to be a long-lasting feature of this source, but less prominent prior to 2008, possibly due to a stronger primary emission. The estimated high-energy cutoff ($E_{\rm cut}=99_{-35}^{+91}$ keV) sets a new redshift record for the farthest ever measured and places APM 08279+5255 in the allowed region of the compactness-temperature diagram of X-ray coronae, in agreement with previous results on high-$z$ quasars.

A. J. Barger, L. L. Cowie, A. H. Blair, L. H. Jones

24 pages, 15 figures, accepted for publication in ApJ

We present deep SCUBA-2 450 micron imaging of the two GOODS fields, achieving a central rms of 1.14 mJy for the GOODS-N and 1.86 mJy for the GOODS-S. For each field, we give a catalog of >4-sigma detections (79 and 16 sources, respectively). We construct the 450 micron number counts, finding excellent agreement with others from the literature. We match the 450 micron sources to 20 cm data (both fields) and ALMA 870 micron data (GOODS-S) to gauge the accuracy of the 450 micron positions. We use the extensive redshift information available on the fields to test how well redshift can be estimated from simple flux ratios (450 micron/850 micron and 20 cm/850 micron), finding tight correlations. We provide a catalog of candidate high-redshift submillimeter galaxies. We look for evolution in dust temperature with redshift by fitting the spectral energy distributions of the sources, but we do not find any significant redshift evolution after accounting for the far-infrared luminosity. We do not find evidence for 450 micron selection picking out warmer sources than 850 micron selection. We find that a 450 micron selected sample only adds low-redshift (z<1.5) galaxies beyond an 850 micron sample.

D. Eckert, S. Ettori, A. Robertson, R. Massey, E. Pointecouteau, D. Harvey, I. G. McCarthy

Re-submitted to A&A after revision

TThe fundamental properties of the postulated dark matter (DM) affect the internal structure of gravitationally-bound structures. In the cold dark matter paradigm, DM particles interact only via gravity. Their distribution is well represented by an Einasto profile with shape parameter $\alpha\approx0.18$, in the smallest dwarf galaxies or the most massive galaxy clusters alike. Conversely, if dark matter particles self-interact via additional forces, we expect the mass density profiles of DM halos to flatten in their central regions, thereby increasing the Einasto shape parameter. We measure the structural properties of the 12 massive X-COP galaxy clusters from observations of their hot gaseous atmosphere using the X-ray observatory XMM-Newton, and of the Sunyaev-Zeldovich effect using the Planck all-sky survey. After removing morphologically disturbed systems, we measure Einasto shape parameters with mean $\langle\alpha\rangle = 0.19 \pm 0.03$ and intrinsic scatter $\sigma_{\alpha}=0.06$, in close agreement with the prediction of the cold dark matter paradigm. We use cosmological hydrodynamical simulations of cluster formation with self-interacting DM (BAHAMAS-SIDM) to determine how the Einasto shape parameter depends on the self-interaction cross section. We use the fitted relation to turn our measurements of $\alpha$ into constraints on the self-interaction cross section, which imply $\sigma/m < 0.19$ cm$^2$/g (95% confidence level) at collision velocity $v_\mathrm{DM-DM}\sim1,000$ km/s. This is lower than the interaction cross-section required for dark matter self-interactions to solve the core-cusp problem in dwarf spheroidal galaxies, unless the cross section is a strong function of velocity.

Genevieve Schroeder (Northwestern/CIERA), Tanmoy Laskar, Wen-fai Fong, Anya E. Nugent, Edo Berger, Ryan Chornock, Kate D. Alexander, Jennifer Andrews, R. Shane Bussmann, Alberto J. Castro-Tirado, Armaan V. Goyal, Charles D. Kilpatrick, Maura Lally, Adam Miller, Peter Milne, Kerry Paterson, Alicia Rouco Escorial, Michael C. Stroh, Giacomo Terreran, Bevin Ashley ZaudererCente

40 pages, 15 figures, submitted to ApJ

We present cm-band and mm-band afterglow observations of five long-duration $\gamma$-ray bursts (GRBs; GRB 130131A, 130420B, 130609A, 131229A, 140713A) with dust-obscured optical afterglow emission, known as "dark" GRBs. We detect the radio afterglow of two of the dark GRBs (GRB 130131A and 140713A), along with a tentative detection of a third (GRB 131229A) with the Karl G. Jansky Very Large Array (VLA). Supplemented by three additional VLA-detected dark GRBs from the literature, we present uniform modeling of their broadband afterglows. We derive high line-of-sight dust extinctions of $A_{V, \rm GRB} \gtrsim 2.2 - 10.6~{\rm mag}$. Additionally, we model the host galaxies of the six bursts in our sample, and derive host galaxy dust extinctions of $A_{V, \rm Host} \approx 0.3-4.7~{\rm mag}$. Across all tested $\gamma$-ray (fluence and duration) and afterglow properties (energy scales, geometries and circumburst densities), we find dark GRBs to be representative of more typical unobscured long GRBs, except in fluence, for which observational biases and inconsistent classification may influence the dark GRB distribution. Additionally, we find that $A_{V, \rm GRB}$ is not related to a uniform distribution of dust throughout the host, nor to the extremely local environment of the burst, indicating that a larger scale patchy dust distribution is the cause of the high line-of-sight extinction. Since radio observations are invaluable to revealing heavily dust-obscured GRBs, we make predictions for the detection of radio emission from host star formation with the next generation VLA.

E. S. Kammoun, D. Barret, P. Peille, R. Willingale, T. Dauser, J. Wilms, M. Guainazzi, J. M. Miller

Accepted for publication in A&A

The X-ray Integral Field Unit (X-IFU) is the high resolution X-ray spectrometer of ESA's Athena X-ray observatory. It will deliver X-ray data in the 0.2-12 keV band with an unprecedented spectral resolution of 2.5 eV up to 7 keV. During the observation of very bright X-ray sources, the X-IFU detectors will receive high photon rates. The count rate capability of the X-IFU will be improved by using the defocusing option, which will enable the observations of extremely bright sources with fluxes up to $\simeq 1$ Crab. In the defocused mode, the point spread function (PSF) of the telescope will be spread over a large number of pixels. In this case, each pixel receives a small fraction of the overall flux. Due to the energy dependence of the PSF, this mode will generate energy dependent artefacts increasing with count rate if not analysed properly. To account for the degradation of the energy resolution with pulse separation in a pixel, a grading scheme (here four grades) will be defined to affect the proper energy response to each event. This will create selection effects preventing the use of the nominal Auxiliary Response File (ARF) for all events. We present a new method for the reconstruction of the spectra obtained from observations performed with a PSF that varies as a function of energy. We apply our method to the case of the X-IFU spectra obtained during the defocused observations. We use the end-to-end SIXTE simulator to model defocused X-IFU observations. Then we estimate new ARF for each of the grades by calculating the effective area at the level of each pixel. Our method allows us to successfully reconstruct the spectra of bright sources when employed in the defocused mode, without any bias. Finally, we address how various sources of uncertainty related to our knowledge of the PSF as a function of energy affect our results.

R. Scott Barrows, Julia M. Comerford, Daniel Stern, Marianne Heida

21 pages, 17 figures. Accepted for publication in the Astrophysical Journal

Ultraluminous X-ray sources (ULXs) are thought to be powerful X-ray binaries (XRBs) and may contribute significantly to the redshift-dependent X-ray emission from star forming galaxies. We have assembled a uniform sample of 259 ULXs over the redshift range z=0.002-0.51 to constrain their physical nature and their contribution to the Cosmic X-Ray Background (CXB). The sample is constructed by crossmatching galaxies from the Sloan Digital Sky Survey with the Chandra Source Catalog and selecting off-nuclear X-ray sources after applying astrometric corrections. The fraction of contaminants is ~30% and shows no evolution with redshift. The host galaxy star formation rates (SFRs) are systematically elevated relative to the parent sample when matched in host stellar mass. The specific SFRs suggest a slight preference for high-mass XRBs, and the X-ray luminosity scaling relations with host galaxy stellar mass and SFR indicate that the highest redshift sources represent relatively luminous XRB populations that dominate their host galaxy X-ray emission. The fraction of galaxies hosting at least one ULX of a given luminosity increases with redshift over the full range of our sample, as expected if ULXs are preferentially found in galaxies with high SFRs and low metallicities. At z~0.5, the ULX X-ray flux is consistent with the X-ray emission from star-forming galaxies. Moreover, ULXs may account for up to ~40% of the integrated flux from XRBs in the normal galaxy population out to z~0.5, suggesting they may contribute significantly to the overall ionizing radiation from galaxies.

Matthew R Buckley, Sung Hak Lim, Eric Putney, David Shih

23 pages, 9 figures

Measuring the density profile of dark matter in the Solar neighborhood has important implications for both dark matter theory and experiment. In this work, we apply autoregressive flows to stars from a realistic simulation of a Milky Way-type galaxy to learn -- in an unsupervised way -- the stellar phase space density and its derivatives. With these as inputs, and under the assumption of dynamic equilibrium, the gravitational acceleration field and mass density can be calculated directly from the Boltzmann Equation without the need to assume either cylindrical symmetry or specific functional forms for the galaxy's mass density. We demonstrate our approach can accurately reconstruct the mass density and acceleration profiles of the simulated galaxy, even in the presence of Gaia-like errors in the kinematic measurements.

Sam F. Sweere, Ivan Valtchanov, Maggie Lieu, Antonia Vojtekova, Eva Verdugo, Maria Santos-Lleo, Florian Pacaud, Alexia Briassouli, Daniel Cámpora Pérez

The field of artificial intelligence based image enhancement has been rapidly evolving over the last few years and is able to produce impressive results on non-astronomical images. In this work we present the first application of Machine Learning based super-resolution (SR) and de-noising (DN) to enhance X-ray images from the European Space Agency's XMM-Newton telescope. Using XMM-Newton images in band [0.5,2] keV from the European Photon Imaging Camera pn detector (EPIC-pn), we develop XMM-SuperRes and XMM-DeNoise deep learning-based models that can generate enhanced SR and DN images from real observations. The models are trained on realistic XMM-Newton simulations such that XMM-SuperRes will output images with two times smaller point-spread function and with improved noise characteristics. The XMM-DeNoise model is trained to produce images with 2.5x the input exposure time from 20 to 50 ks. When tested on real images, DN improves the image quality by 8.2%, as quantified by the global peak-signal-to-noise ratio. These enhanced images allow identification of features that are otherwise hard or impossible to perceive in the original or in filtered/smoothed images with traditional methods. We demonstrate the feasibility of using our deep learning models to enhance XMM-Newton X-ray images to increase their scientific value in a way that could benefit the legacy of the XMM-Newton archive.

In this article, we do a thorough investigation of the competency of the forthcoming Cosmic Microwave Background (CMB) and Galaxy surveys in probing the features in the primordial power spectrum. Primordial features are specific model-dependent corrections on top of the standard power-law inflationary power spectrum; the functional form being given by different inflationary scenarios. Signature of any significant departure from the feature-less power spectrum will enable us to decipher the intricacies of the inflationary Universe. Here, we delve into three major yet distinct features, namely, Bump feature, Sharp feature signal, and Resonance feature signal. To analyse the features, we adopt a specific template for each feature model. We estimate the possible constraints on the feature parameters by employing Fisher matrix forecast analysis for the upcoming CMB missions such as CMB-S4, CORE-M5, LiteBIRD, PICO conjointly with DESI, and EUCLID galaxy surveys. To this end, we make use of four distinct observations to forecast on the bounds on the model parameters, namely, CMB, Baryon Acoustic Oscillations (BAO), Galaxy Clustering and Gravitational Weak Lensing or Cosmic Shear and their permissible synergy. For large scale structure (LSS) information, we consider different upper limits of scale for different redshifts for the purpose of circumventing the propagation of the errors stemming from the uncertainties on nonlinear scales into the constraints on the feature parameters. A comparative analysis of all three features has been done to estimate relative capabilities of these upcoming observations in shedding light on this crucial aspect of precision cosmology.

Most amino acids and sugars molecules occur in mirror, or chiral, images of each other, knowns as enantiomers. However, life on Earth is mostly homochiral: proteins contain almost exclusively L-amino acids, while only D-sugars appear in RNA and DNA. The mechanism behind this fundamental asymmetry of life remains unknown, despite much progress in the theoretical and experimental understanding of homochirality in the past decades. We review three potential mechanisms for the emergence of biological homochirality on primal Earth and explore their implications for astrobiology: the first, that biological homochirality is a stochastic process driven by local environmental fluctuations; the second, that it is driven by circularly-polarized ultraviolet radiation in star-forming regions; and the third, that it is driven by parity violation at the elementary particle level. We argue that each of these mechanisms leads to different observational consequences for the existence of enantiomeric excesses in our solar system and in exoplanets, pointing to the possibility that the search for life elsewhere will help elucidate the origins of homochirality on Earth.

Jaime E. Pineda, Jorma Harju, Paola Caselli, Olli Sipilä, Mika Juvela, Charlotte Vastel, Erik Rosolowsky, Andreas Burkert, Rachel K. Friesen, Yancy Shirley, María José Maureira, Spandan Choudhury, Dominique M. Segura-Cox, Rolf Güsten, Anna Punanova, Luca Bizzocchi, Alyssa A. Goodman

Accepted for publication in AJ. 18 Pages, 16 Figures, 3 Tables

Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores.

S. Duarte Puertas, J. M. Vilchez, J. Iglesias-Páramo, M. Mollá, E. Pérez-Montero, C. Kehrig, L. S. Pilyugin, I. A. Zinchenko

Accepted for publication in A&A, 18 pages, 13 figures, and 1 table

In this work we study the stellar mass -- metallicity relation (MZR) of an extended sample of star-forming galaxies in the local Universe and its possible dependence with the star formation rate (SFR). A sample of $\sim$195000 Sloan Digital Sky Survey (SDSS) star-forming galaxies has been selected up to z=0.22 with the aim of analysing the behaviour of the relation of MZR with respect to SFR and taking into account the age of their stellar populations. For this sample we have obtained, for the first time, aperture corrected oxygen and nitrogen-to-oxygen abundances (O/H and N/O, respectively) and SFR using the empirical prescriptions from the Calar Alto Legacy Integral Field Area (CALIFA) survey. To perform this study we make use also of the stellar mass of the galaxies and the parameter Dn(4000) as a proxy of the age of the stellar population. We derive a robust MZR locus, which is found to be fully consistent with the "anchoring" points of a selected set of well studied nearby galaxies with a direct derivation of the chemical abundance. A complex relation between MZR and SFR across the whole range of galaxy mass and metallicity has been observed, where the slope changes seen in the O/H -- SFR plane present a pattern which seems to be tuned to the galaxies' stellar age, and therefore, stellar age has to be taken into account in the stellar mass -- metallicity -- SFR relation. In order to provide an answer to the question of whether or not the MZR depends on the SFR it is essential to take into account the age of the stellar populations of galaxies. A strong dependence between the MZR and SFR is observed mainly for star-forming galaxies with strong SFR values and low Dn(4000). The youngest galaxies of our SDSS sample show the highest SFR measured for their stellar mass.

Lister Staveley-Smith, Emma Bond, Kenji Bekki, Tobias Westmeier

11 pages, 4 figures, 3 tables; accepted by Publications of the Astronomical Society of Australia

A plausible formation scenario for the Galactic globular clusters 47 Tuc and Omega Cen is that they are tidally stripped remnants of dwarf galaxies, in which case they are likely to have retained a fraction of their dark matter cores. In this study, we have used the ultra-wide band receiver on the Parkes telescope (Murriyang) to place upper limits on the annihilation rate of exotic Light Dark Matter particles (X) via the XX/$e^+e^-$ channel using measurements of the recombination rate of positronium (Ps). This is an extension of a technique previously used to search for Ps in the Galactic Centre. However, by stacking of spectral data at multiple line frequencies, we have been able to improve sensitivity. Our measurements have resulted in 3-sigma flux density (recombination rate) upper limits of 1.7 mJy (1.4x10^43/s) and 0.8 mJy (1.1x10^43/s) for 47 Tuc and Omega Cen, respectively. Within the Parkes beam at the cluster distances, which varies from 10 to 23 pc depending on the frequency of the recombination line, we calculate upper limits to the dark matter mass and rms dark matter density of <1.2-1.3x10$^5$ fn$^{-0.5}$ (mX/MeV c^-2) solar masses and <48-54 fn$^{-0.5}$ (mX/MeV c$^{-2}$) solar masses pc$^{-3}$ for the clusters, where fn=Rn/Rp is the ratio of Ps recombination transitions to annihilations, estimated to be 0.01. The radio limits for Omega Cen suggest that, for a fiducial dark/luminous mass ratio of 0.05, any contribution from Light Dark Matter is small unless the cross section <7.9x10^-28 (mX/MeV c<$^{-2}$)$^2$ cm$^3$/s. Owing to the compactness and proximity of the clusters, archival 511-keV measurements suggest even tighter limits than permitted by CMB anisotropies, <8.6x10$^{-31}$ (mX/MeV c$^{-2}$)$^2$ cm$^3$/s. Due to the very low synchrotron radiation background, our recombination rate limits substantially improve on previous radio limits for the Milky Way.

D. Alina, J.-Ph. Bernard, K.H. Yuen, A. Lazarian, A. Hughes, M. Iskakova, A. Akimkhan, A. Mukanova

14 pages, 22 figures

Polarization of interstellar dust emission is a powerful probe of dust properties and magnetic field structure. Yet studies of external galaxies are hampered by foreground dust contribution. The aim of this study is to separate the polarised signal from the Large Magellanic Cloud (LMC) from that of the Milky Way (MW) in order to construct a wide-field, spatially complete map of dust polarization using the Planck 353 GHz data. To estimate the foreground polarization direction, we used velocity gradients in HI spectral line data and assessed the performance of the output by comparing to starlight extinction polarization. We estimate the foreground intensity using dust properties derived previously from the Planck data and we assume the foreground polarization to be uniform and equal to the average of the MW around the galaxy. After foreground removal, the geometry of the plane-of-the-sky magnetic field in the LMC tends to follow the structure of the atomic gas in the LMC. This is notably the case along the molecular ridges extending south and south-east of the 30 Doradus star-forming complex, and along more diffuse southern arm extending towards the Small Magellanic Cloud. There is also an alignment between the magnetic field and the outer arm in the western part. The median polarization fraction in the LMC is slightly lower than that observed for the MW while the anti-correlation between the polarization angle dispersion function and the polarization fraction is slightly larger. Overall, polarization fraction distribution is similar to that observed in the MW.

M. S. Rumenskikh, I. F. Shaikhislamov, M. L. Khodachenko, H. Lammer, I. B. Miroshnichenko, A. G. Berezutsky, L. Fossati

A 3D fully self-consistent multi-fluid hydrodynamic aeronomy model is applied to simulate the hydrogen-helium expanding upper atmosphere of the hot Jupiter HD189733b, and related absorption in the Lya line and the 10830 A line of metastable helium. We studied the influence of a high-energy stellar flux, stellar wind, and Lya cooling to reproduce the available observations. We found that to fit the width of the absorption profile in 10830 A line the escaping upper atmosphere of planet should be close to the energy limited escape achieved with a significantly reduced Lya cooling at the altitudes with HI density higher than 3*10^6 cm^-3. Based on the preformed simulations, we constrain the helium abundance in the upper atmosphere of HD189733b by a rather low value of He/H~0.005. We show that under conditions of a moderate stellar wind similar to that of the Sun the absorption of Lya line takes place mostly within the Roche lobe due to thermal broadening at a level of about 7%. At an order of magnitude stronger wind, a significant absorption of about 15% at high blue shifted velocities of up to 100 km/s is generated in the bowshock region, due to Doppler broadening. These blue shifted velocities are still lower than those (~200 km/s) detected in one of the observations. We explain the differences between performed observations, though not in all the details, by the stellar activity and the related fluctuations of the ionizing radiation (in case of 10830 A line), and stellar wind (in case of Lya line).

John D. Ilee, Catherine Walsh, Jeff Jennings, Richard A. Booth, Giovanni P. Rosotti, Richard Teague, Takashi Tsukagoshi, Hideko Nomura

6 pages, 4 figures, accepted for publication in MNRAS Letters

The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $\mu$Jy beam$^{-1}$) moderate resolution (0.37") ALMA observations at 1 mm of the nearest protoplanetary disc, TW Hya. Using the visibility analysis tool `frank', we reveal a structured millimetre intensity distribution out to 100 au, well beyond previous estimates of 60-70 au. Our analysis suggests the presence of a new millimetre continuum gap at 82 au, coincident with similar features seen in optical/near-infrared scattered light and millimetre molecular line observations. Examination of the fit residuals confirms the presence of the previously reported au-scale continuum excess at 52 au (P.A. = 242.5 degrees). Our results demonstrate the utility of combining deep, moderate resolution observations with super-resolution analysis techniques to probe the faintest regions of protoplanetary discs.

Nathalie Ysard, Marc-Antoine Miville-Deschênes, Laurent Verstraete, Anthony Peter Jones

accepted for publication by A&A

Context. Excess microwave emission commonly known as `anomalous microwave emission' (AME) is now routinely detected in the Milky Way. Although its link with the rotation of interstellar (carbonaceous) nano-grains seems to be relatively well established at cloud scales, large-scale observations show a lack of correlation between the different tracers of nano-carbons and AME, which has led the community to question the viability of this link. Aims. Using ancillary data and spinning dust models for nano-carbons and nano-silicates, we explore the extent to which the out of Galactic Plane AME could come from one carrier or the other. Methods. Contrary to previous large-scale studies, our method is not built on comparing the correlations of the different dust tracers with each other, but rather on comparing the poor correlations predicted by the models with the correlations actually observed. This is based on as realistic as possible estimates of the gas ionisation state and grain charge as a function of the local radiation field and gas density. Results. First, nano-carbon dust can explain all the observations for medium properties in agreement with the latest findings about the separation of cold and warm neutral medium in the diffuse ISM. The dispersion in the observations can be accounted for with little variations in the dust size distribution, abundance or electric dipole moment. Second, whatever the properties and the abundance of the nano-silicate dust considered, spinning nano-silicates are excluded as being the sole source of the AME. Third, the best agreement with the observations is obtained if the emission of spinning nano-carbons alone is taken into account. However, a marginal participation of nano-silicates to the AME production cannot be excluded as long as their abundance does not exceed $Y_{\rm Si} \sim 1\%$.

Miguel Cárcamo, Anna M. M. Scaife, Emma L. Alexander, J. Patrick Leahy

Submitted to MNRAS

The reconstruction of Faraday depth structure from incomplete spectral polarization radio measurements using the RM Synthesis technique is an under-constrained problem requiring additional regularisation. In this paper we present cs-romer: a novel object-oriented compressed sensing framework to reconstruct Faraday depth signals from spectro-polarization radio data. Unlike previous compressed sensing applications, this framework is designed to work directly with data that are irregularly sampled in wavelength-squared space and to incorporate multiple forms of compressed sensing regularisation. We demonstrate the framework using simulated data for the JVLA telescope under a variety of observing conditions, and we introduce a methodology for identifying the optimal basis function for reconstruction of these data, using an approach that can also be applied to datasets from other telescopes and over different frequency ranges. In this work we show that the delta basis function provides optimal reconstruction for JVLA L-band data and we use this basis with observations of the low-mass galaxy cluster Abell 1314 in order to reconstruct the Faraday depth of its constituent cluster galaxies. We use the cs-romer framework to de-rotate the Galactic Faraday depth contribution directly from the wavelength-squared data and to handle the spectral behaviour of different radio sources in a direction-dependent manner. The results of this analysis show that individual galaxies within Abell 1314 deviate from the behaviour expected for a Faraday-thin screen such as the intra-cluster medium and instead suggest that the Faraday rotation exhibited by these galaxies is dominated by their local environments.

J van Dongen, L Prokhorov, S J Cooper, M A Barton, E Bonilla, K L Dooley, J C Driggers, A Effler, N A Holland, A Huddart, M Kasprzack, J S Kissel, B Lantz, A L Mitchell, J O'Dell, A Pele, C Robertson, C M Mow-Lowry

13 pages, 16 figures

Control noise is a limiting factor in the low-frequency performance of the LIGO gravitational wave detectors. In this paper we model the effects of using new sensors called HoQIs to control the suspension resonances. We show if we were to use HoQIs, instead of the standard shadow sensors, we can suppress resonance peaks up to tenfold more while simultaneously reducing the noise injected by the damping system. Through a cascade of effects this will reduce the resonant cross-coupling, allow for improved stability for feed-forward control, and result in improved sensitivity of the detector in the 10-20 Hz band. This analysis shows that local sensors such as HoQIs should be used in current and future detectors to improve low-frequency performance.

G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Accepted for publication in A&A. 14 pages, 7 figures, 2 tables

We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with $\sim 17000$ sources. Both samples have the same photometric coverage, from optical to far-infrared. We construct the spectral energy distributions (SEDs) of all sources and fit them using the CIGALE code, utilizing the same templates and parametric grid for both samples. We account for the mass incompleteness and exclude quiescent systems from both datasets. These allow us to compare the SFR of the two populations in a uniform manner, minimising systematic effects. Based on our analysis, AGN at low and moderate L$_X$ ($\rm L_{X,2-10keV} < 10^{44}\,ergs^{-1}$), have SFR that is lower, or at most, equal to that of star-forming galaxies, in agreement with previous studies. The large number of luminous X-ray AGN available in the eFEDS dataset, enable us to expand our investigations at higher L$_X$ to test previous, tentative results. At $\rm L_{X,2-10keV} > 10^{44.2}\,ergs^{-1}$, the SFR of AGN appears enhanced, by $\sim 30\%$, compared to that of star-forming sources, for systems with stellar mass, $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$, confirming indications found in previous studies. The most massive sources $\rm log\,[M_*(M_\odot)] > 11.5$, present a flat SFR$_{norm}$-L$_X$ relation up to $\rm L_{X,2-10keV} \sim 10^{44.5}\,ergs^{-1}$, with SFR similar to that of star-forming galaxies. However, at higher L$_X$ ($\rm L_{X,2-10keV} \sim 10^{45}\,ergs^{-1}$), we find indications that the SFR of these massive AGN hosts may be enhanced compared to that of non-AGN systems.

Chengdong Li, James Binney

17 pages, 20 figures, submitted to MNRAS

We investigate the structure of our Galaxy's young stellar disc by fitting the distribution functions (DFs) of a new family to five-dimensional Gaia data for a sample of $47\,000$ OB stars. Tests of the fitting procedure show that the young disc's DF would be strongly constrained by Gaia data if the distribution of Galactic dust were accurately known. The DF that best fits the real data accurately predicts the kinematics of stars at their observed locations, but it predicts the spatial distribution of stars poorly, almost certainly on account of errors in the best-available dust map. We argue that dust models could be greatly improved by modifying the dust model until the spatial distribution of stars predicted by a DF agreed with the data. The surface density of OB stars is predicted to peak at $R\simeq5.5\mbox{kpc}$, slightly outside the reported peak in the surface density of molecular gas; we suggest that the latter radius may have been under-estimated through the use of poor kinematic distances. The velocity distributions predicted by the best-fit DF for stars with measured line-of-sight velocities $v_\parallel$ reveal that the outer disc is disturbed at the level of 10 $\mbox{km}~\mbox{s}^{-1}$ in agreement with earlier studies, and that the measured values of $v_\parallel$ have significant contributions from the orbital velocities of binaries. Hence the outer disc is colder than it is sometimes reported to be.

P. F. L. Maxted, N. J. Miller, S. Hoyer, V. Adibekyan, S. G. Sousa, N. Billot, A. Fortier, A. E. Simon, A. Collier Cameron, M. I. Sawyne, P. Gutermann, A. H. M. J. Triaud, J. Southworth, Y. Alibert, R. Alonso, G. Anglada, T. Bárczy, D. Barrado y Navascues, S. C. C. Barros, W. Baumjohann, M. Beck, T. Beck, W. Benz, X. Bonfils, A. Brandeker, C. Broeg, M. Buder, J. Cabrera, S. Charnoz, C. Corral van Damme, Sz. Csizmadia, M. B. Davies, M. Deleuil, L. Delrez, O. Demangeon, B.-O. Demory, D. Ehrenreich, A. Erikson, L. Fossati, M. Fridlund, D. Gandolfi, M. Gillon, M. Güdel, K. Heng, J. E. Hernández Leon, K. G. Isaak, L. L. Kiss, J. Laskar, A. Lecavelier des Etangs, M. Lendl, C. Lovis, D. Magrin, M. Munari, V. Nascimbeni, G. Olofsson, R. Ottensamer, I. Pagano, E. Pallé, G. Peter, et al. (18 additional authors not shown)

Accepted for publication in MNRAS. 15 pages, 9 figures

EBLM J0113+31 is moderately bright (V=10.1), metal-poor ([Fe/H]$\approx-0.3$) G0V star with a much fainter M dwarf companion on a wide, eccentric orbit (=14.3 d). We have used near-infrared spectroscopy obtained with the SPIRou spectrograph to measure the semi-amplitude of the M dwarf's spectroscopic orbit, and high-precision photometry of the eclipse and transit from the CHEOPS and TESS space missions to measure the geometry of this binary system. From the combined analysis of these data together with previously published observations we obtain the following model-independent masses and radii: $M_1 = 1.029 \pm 0.025 M_{\odot}$, $M_2 = 0.197 \pm 0.003 M_{\odot}$, $R_1 = 1.417 \pm 0.014 R_{\odot}$, $R_2 = 0.215 \pm 0.002 R_{\odot}$. Using $R_1$ and the parallax from Gaia EDR3 we find that this star's angular diameter is $\theta = 0.0745 \pm 0.0007$ mas. The apparent bolometric flux of the G0V star corrected for both extinction and the contribution from the M dwarf ($<0.2$ per cent) is ${\mathcal F}_{\oplus,0} = (2.62\pm 0.05)\times10^{-9}$ erg.cm$^{-2}$.s$^{-1}$. Hence, this G0V star has an effective temperature $T_{\rm eff,1} = 6124{\rm\,K} \pm 40{\rm \,K\,(rnd.)} \pm 10 {\rm \,K\,(sys.)}$. EBLM J0113+31 is an ideal benchmark star that can be used for "end-to-end" tests of the stellar parameters measured by large-scale spectroscopic surveys, or stellar parameters derived from asteroseismology with PLATO. The techniques developed here can be applied to many other eclipsing binaries in order to create a network of such benchmark stars.

Qinghui Sun, Sharon Xuesong Wang, Tianjun Gan, Andrew W. Mann

3 figures, 3 tables, accepted for publication in Research in Astronomy and Astrophysics

We report the results of our search of planet candidates in Open Clusters and Young Stellar Associations based on the TESS Objects of Interest Catalog. We find one confirmed planet, one promising candidate, one brown dwarf, and three unverified planet candidates in a sample of 1229 Open Clusters from the second Gaia data release. We discuss individual planet-star systems based on their basic parameters, membership probability, and the observation notes from the ExoFOP-TESS website. We also find ten planet candidates (P $>$ 95%) in Young Stellar Associations by using the BANYAN $\Sigma$ Bayesian Algorithm. Among the ten candidates, five are known planet systems. We estimate the rotation periods of the host stars using the TESS light curves and estimate their ages based on gyrochronology. Two candidates with periodic variations are likely to be young planets, but their exact memberships to Young Stellar Associations remain unknown.

Anatole Storck, Ross Church

9 pages, 9 Figures, comments welcome

The upcoming space-based gravitational wave observatory, LISA, will detect binaries containing pairs of neutron stars. LISA observations of double neutron stars, particularly in combination with electromagnetic observations, have the potential to improve our understanding of neutron star physics and binary evolution. In this work we assess the effect of changing the model of the Milky Way on predictions of the population of double neutron star binaries that will be detected and resolved by LISA. We conclude that the spatial and gravitational-wave frequency distributions of these binaries are insensitive to the choice of galactic models, compared to the stochastic variation induced by the small sample size. In particular, the time-consuming computation of the binaries' Galactic orbits is not necessary. However, the distribution of eccentricities is strongly affected by the choice of star-formation history. Binaries with eccentricities e>0.1 that can be measured by LISA observations are all younger than 100 Myr. We caution that comparisons between different predictions for LISA observations need to use consistent star formation histories, and that the Galactic star formation history should be taken into account in the analysis of the observations themselves.

In recent years new types of coordinate transformations have appeared in cosmology on top of the standard gauge transformations, such as the dilatations and special conformal transformations, or the ones leading to (conformal) Fermi coordinates. Some of these can remove effects that are invariant under the standard gauge transformations and also affect asymptotic boundary conditions, thus introducing a non-trivial ambiguity in our cosmological modeling. In this short note we point out that this ambiguity is irrelevant for the quantities we use to compare our model with observations -- the cosmological observable relations -- as they are invariant under all of these transformations. Importantly, this invariance holds only if one takes into account all the relativistic contributions to an observable, which is not the case in the literature in general. We finally also show that the practically-relevant property of conformal Fermi coordinates (a FLRW metric up to second order in distance) can be achieved through a globally-defined standard gauge transformation.

Marco Palla, Pablo Santos-Peral, Alejandra Recio-Blanco, Francesca Matteucci

16 pages, 16 figures. Accepted for publication in Astronomy and Astrophysics (A&A)

Context. The [Mg/Fe] abundance ratios are a fundamental fossil signature to trace the chemical evolution of the disc. Despite of the huge observational and theoretical effort, discrepancies between models and data are still present and several explanations have been put forward to explain the [$\alpha$/Fe] bimodality. Aims. In this work, we take advantage of a new AMBRE:HARPS dataset, which provides new and more precise [Mg/Fe] estimations, as well as reliable stellar ages for a subsample of stars, to study the evolution of the solar neighbourhood. Methods. The above data are compared with detailed chemical evolution models for the Milky Way, exploring the most used prescriptions for stellar yields and different formation scenarios for the Galactic disc, i.e. the delayed two-infall and the parallel model, also including prescriptions for stellar radial migration. Results. We see that most of the stellar yields struggle to reproduce the observed trend of the data and that semi-empirical yields are still the best to describe the [Mg/Fe] evolution in the thick and thin discs. In particular, most of the yields still predict a steeper decrease of the [Mg/Fe] ratio at high metallicity than what is shown by the data. The bulk of the data are well reproduced by the parallel and two-infall scenarios, but both scenarios have problems in explaining the most metal-rich and metal-poor tails of the low-$\alpha$ data. These tails can be explained in light of radial migration from inner and outer disc regions, respectively. Conclusions. Despite of the evidence of stellar migration, it is difficult to estimate the actual contribution of stars from other parts of the disc to the solar vicinity. However, the comparison between data and models suggests that peculiar histories of star formation, such as that of the two-infall model, are still needed to reproduce the observed distribution of stars.

Jacob W. Isbell (1), Klaus Meisenheimer (1), Jörg-Uwe Pott (1), Marko Stalevski (2, 3), Konrad R. W. Tristram (4), Joel Sanchez-Bermudez (1, 5), Karl-Heinz Hofmann (6), Violeta Gámez Rosas (7), Walter Jaffe (7), Leonard Burtscher (7), James Leftley (8), Romain Petrov (8), Bruno Lopez (8), Thomas Henning (1), Gerd Weigelt (6), Fatme Allouche (8), Philippe Berio (8), Felix Bettonvil (7), Pierre Cruzalebes (8), Carsten Dominik (9), Matthias Heininger (6), Michiel Hogerheijde (7, 9), Stéphane Lagarde (8), Michael Lehmitz (1), Alexis Matter (8), Florentin Millour (8), Sylvie Robbe-Dubois (8), Dieter Schertl (6), Roy van Boekel (1), Josef Varga (7), Julien Woillez (10) ((1) Max Planck Institute for Astronomy, (2) Astronomical Observatory of Belgrade, (3) Sterrenkundig Observatorium, Universiteit Gent, (4) European Southern Observatory, Santiago, Chile, (5) Instituto de Astronomía, Universidad Nacional Autónoma de México, (6) Max Planck Institute for Radio Astronomy, (7) Sterrewacht Leiden, (8) Laboratoire Lagrange, Université Côte d'Azur, (9) Anton Pannekoek Institute for Astronomy, University of Amsterdam, (10) European Southern Observatory, Garching, Germany)

Accepted for publication in A&A

Active galactic nuclei play a key role in the evolution of galaxies, but their inner workings and physical connection to the host are poorly understood due to a lack of angular resolution. Infrared interferometry makes it possible to resolve the circumnuclear dust in the nearby Seyfert 2 galaxy, Circinus. Previous observations have revealed complex structures and polar dust emission but interpretation was limited to simple models. MATISSE makes it possible to image these structures for the first time. We observed the Circinus Galaxy with VLTI/MATISSE, producing 150 correlated flux spectra and 100 closure phase spectra. We reconstructed images in the N-band at ~10 mas resolution. We fit blackbody functions with dust extinction to several aperture-extracted fluxes from the images to produce a temperature distribution of central dusty structures. We find significant substructure in the circumnuclear dust: central unresolved flux of ~0.5 Jy, a thin disk 1.9 pc in diameter oriented along ~45 deg,and a ~4x1.5 pc polar emission extending orthogonal to the disk. The polar emission exhibits patchiness, which we attribute to clumpy dust. Flux enhancements to the east and west of the disk are seen for the first time. We distinguish the temperature profiles of the disk and of the polar emission: the disk shows a steep temperature gradient indicative of denser material; the polar profile is flatter, indicating clumpiness and/or lower dust density. The unresolved flux is fitted with a high temperature, ~370 K. The polar dust remains warm (~200 K) out to 1.5 pc from the disk. The recovered morphology and temperature distribution resembles modeling of accretion disks with radiation-driven winds at large scales, but we placed new constraints on the subparsec dust. The subparsec features imaged here place new constraints on the physical modeling of circumnuclear dust in active galaxies.

Adam B. Langeveld, Nikku Madhusudhan, Samuel H. C. Cabot

Submitted to MNRAS after addressing reviewer comments. 22 pages, 10 figures

The alkali metal sodium (Na) is one of the most commonly detected chemical species in the upper atmospheres of giant exoplanets. In this work we conducted a homogeneous survey of Na in a diverse sample of ten highly irradiated giant exoplanets using high-resolution transmission spectroscopy. Our sample includes nine planets with previous Na detections and one new detection. We confirm previous detections and assess multiple approaches for deriving Na line properties from high-resolution transmission spectra. The homogeneously measured sodium line depths were used to constrain the atmospheric heights ($H_{\text{Na}}$) with respect to the planetary radii ($R_{\text{p}}$). We assess an empirical trend describing the relative atmospheric height ($H_{\text{Na}}/R_{\text{p}}$) as a function of planetary equilibrium temperature ($T_{\text{eq}}$) and surface gravity ($g$), in which $H_{\text{Na}}/R_{\text{p}}$ decreases exponentially with $\xi \propto gT_{\text{eq}}$, approaching a constant at large $\xi$. We also report the sodium D2/D1 line ratios across our sample and find that seven targets have line ratios which are consistent with unity. Finally, we measured net blueshifted offsets of the sodium absorption lines from their rest frame wavelengths for all ten planets, corresponding to day-night wind velocities of a few km s$^{-1}$. This suggests that the broad sample of exoplanets share common underlying processes which govern atmospheric dynamics. Our study highlights a promising avenue for using high-resolution transmission spectroscopy to further our understanding of how atmospheric characteristics vary over a diverse sample of exoplanets.

Luca Caloni, Martina Gerbino, Massimiliano Lattanzi, Luca Visinelli

26 pages, 8 figures

We constrain the coupling of thermally-produced axion-like particles (here axions) with photons and gluons, using data from the cosmic microwave background (CMB) spectra and baryon acoustic oscillations. The axion possesses an explicit soft breaking mass term and it is produced thermally in the early Universe from either axion-photon or axion-gluon processes, accounting for the recent progresses in the field. We derive the most stringent bounds on the axion-gluon coupling to date on the mass range considered $10^{-4} \lesssim m_a/{\rm eV} \lesssim 100$, superseding the current bounds from SN1987A. The bounds on the axion-photon coupling are competitive with the results from the CAST collaboration for the axion mass $m_a \gtrsim 3\,$eV. We comment on the forecast reaches that will be available given the sensitivity of future CMB-S4 missions.

A. Migliorini, Z. Kanuchova, S. Ioppolo, M. Barbieri, N.C. Jones, S.V. Hoffmann, G. Strazzulla, F. Tosi, G. Piccioni

Accepted for publication in Icarus

Since its first identification on the surface of Ganymede in 1995, molecular oxygen (O2) ice has been at the center of a scientific debate as the surface temperature of the Jovian moon is on average well above the freezing point of O2. Laboratory evidence suggested that solid O2 may either exist in a cold (<50 K) subsurface layer of the icy surface of Ganymede, or it is in an atmospheric haze of the moon. Alternatively, O2 is constantly replenished at the surface through ion irradiation of water-containing ices. A conclusive answer on the existence of solid O2 on the surface of Ganymede is hampered by the lack of detailed, extensive observational datasets. We present new ground-based, high-resolution spectroscopic observations of Ganymede's surface obtained at the Telescopio Nazionale Galileo. These are combined with dedicated laboratory measurements of ultraviolet-visible (UV-vis) photoabsorption spectra of O2 ice, both pure and mixed with other species of potential interest for the Galilean satellites. Our study confirms that the two bands identified in the visible spectra of Ganymede's surface are due to the (1,0) and (0,0) transition bands of O2 ice. Oxygen-rich ice mixtures including water (H2O) and carbon dioxide (CO2) can reproduce observational reflectance data of the Ganymede's surface better than pure O2 ice in the temperature range 20-35 K. Solid H2O and CO2 also provide an environment where O2 ice can be trapped at higher temperatures than its pure ice desorption under vacuum space conditions. Our experiments at different temperatures show also that the (1,0)/(0,0) ratio in case of the CO2:O2=1:2 ice mixture at 35 K has the closest value to observations, while at 30 K the (1,0)/(0,0) ratio seems to be mixture independent with the exception of the N2:O2=1:2 ice mixture. The present work will support the ESA/JUICE mission to the Jovian system.

Robert D. Kavanagh, Aline A. Vidotto, Harish K. Vedantham, Moira M. Jardine, Joseph R. Callingham, Julien Morin

15 pages, 12 figures. Accepted for publication in MNRAS

The nearby M dwarf WX UMa has recently been detected at radio wavelengths with LOFAR. The combination of its observed brightness temperature and circular polarisation fraction suggests that the emission is generated via the electron-cyclotron maser instability. Two distinct mechanisms have been proposed to power such emission from low-mass stars: either a sub-Alfv\'enic interaction between the stellar magnetic field and an orbiting planet, or reconnection at the edge of the stellar magnetosphere. In this paper, we investigate the feasibility of both mechanisms, utilising the information about the star's surrounding plasma environment obtained from modelling its stellar wind. Using this information, we show that a Neptune-sized exoplanet with a magnetic field strength of 10-100 G orbiting at ~0.034 au can accurately reproduce the observed radio emission from the star, with corresponding orbital periods of 7.4 days. Due to the stellar inclination, a planet in an equatorial orbit is unlikely to transit the star. While such a planet could induce radial velocity semi-amplitudes from 7 to 396 m s$^{-1}$, it is unlikely that this signal could be detected with current techniques due to the activity of the host star. The application of our planet-induced radio emission model here illustrates its exciting potential as a new tool for identifying planet-hosting candidates from long-term radio monitoring. We also develop a model to investigate the reconnection-powered emission scenario. While this approach produces less favourable results than the planet-induced scenario, it nevertheless serves as a potential alternative emission mechanism which is worth exploring further.

André de Gouvêa, Ivan Martinez-Soler, Yuber F. Perez-Gonzalez, Manibrata Sen

15 pages, 9 figures

The relic neutrinos from old supernova explosions are among the most ancient neutrino fluxes within experimental reach. Thus, the diffuse supernova neutrino background (DSNB) could teach us if neutrino masses were different in the past (redshifts $z\lesssim 5$). Oscillations inside the supernova depend strongly on the neutrino mass-squared differences and the values of the mixing angles, rendering the DSNB energy spectrum sensitive to variations of these parameters. Considering a purely phenomenological parameterization of the neutrino masses as a function of redshift, we compute the expected local DSNB spectrum here on Earth. Given the current knowledge of neutrino oscillation parameters, specially the fact that $|U_{e3}|^2$ is small, we find that the $\nu_e$ spectrum could be significantly different from standard expectations if neutrinos were effectively massless at $z\gtrsim1$ as long as the neutrino mass ordering is normal. On the other hand, the $\overline{\nu}_e$ flux is not expected to be significantly impacted. Hence, a measurement of both the neutrino and antineutrino components of the DSNB should allow one to test the possibility of recent neutrino mass generation.

Features in the inflationary landscape can inject extra energies to inflation models and produce on-shell particles with masses much larger than the Hubble scale of inflation. This possibility extends the energy reach of the program of cosmological collider physics, in which signals associated with these particles are generically Boltzmann-suppressed. We study the mechanisms of this classical cosmological collider in two categories of primordial features. In the first category, the primordial feature is classical oscillation, which includes the case of coherent oscillation of a massive field and the case of oscillatory features in the inflationary potential. The second category includes any sharp feature in the inflation model. All these classical features can excite unsuppressed quantum modes of other heavy fields which leave observational signatures in primordial non-Gaussianities, including the information about the particle spectra of these heavy degrees of freedom.

Anirban Biswas, Debasish Borah, Nayan Das, Dibyendu Nanda

38 pages, 7 captioned figures

We propose a scenario where dark matter (DM) and dark radiation ($\Delta {\rm N}_{\rm eff}$) can be generated non-thermally due to the presence of a light Dirac neutrino portal between the standard model (SM) and dark sector particles. The SM is minimally extended by three right handed neutrinos ($\nu_R$), a Dirac fermion DM candidate ($\psi$) and a complex scalar ($\phi$), transforming non-trivially under an unbroken $\mathbb{Z}_4$ symmetry while being singlets under the SM gauge group. While DM and $\nu_R$ couplings are considered to be tiny in order to be in the non-thermal or freeze-in regime, $\phi$ can be produced either thermally or non-thermally depending upon the strength of its Higgs portal coupling. We consider both these possibilities and find out the resulting DM abundance and $\Delta {\rm N}_{\rm eff}$ via freeze-in mechanism to constrain the model parameters in the light of Planck 2018 data. We find that the scenario where $\phi$ remains out of equilibrium throughout or after a certain epoch allows more parameter space consistent with DM phenomenology and $\Delta {\rm N}_{\rm eff}$ in view of Planck 2018 data. The next generation experiments like CMB-S4, SPT-3G etc. will have the required sensitivities to probe a major portion of the entire model parameter space, offering a promising way of probing such non-thermal DM scenario which typical direct detection experiments are not much sensitive to.

Emanuel V. Chimanski, Ronaldo V. Lobato, Andre R. Goncalves, Carlos A. Bertulani

The description of stellar interiors remains as a big challenge for the nuclear astrophysics community. The consolidated knowledge is restricted to density regions around the saturation of hadronic matter $\rho _{0} = 2.8\times 10^{14} {\rm\ g\ cm^{-3}}$, regimes where our nuclear models are successfully applied. As one moves towards higher densities and extreme conditions up to five to twenty times $\rho_{0}$, little can be said about the microphysics of such objects. Here, we employ a Markov Chain Monte Carlo (MCMC) strategy in order to access the variability of polytropic three-pircewised models for neutron star equation of states. With a fixed description of the hadronic matter, we explore a variety of models for the high density regimes leading to stellar masses up to 2.5 $M_{\odot}$. In addition, we also discuss the use of a Bayesian power regression model with heteroscedastic error. The set of EoS from LIGO was used as inputs and treated as data set for testing case.

Third-generation gravitational wave (GW) observatories such as Einstein Telescope (ET) and Cosmic Explorer (CE) will be ideal instruments to probe the structure of neutron stars through the GWs they emit when undergoing binary coalescence. In this work we make predictions about how well ET in particular will enable us to reconstruct the neutron star equation of state through observations of tens of binary neutron star coalescences with signal-to-noise ratios in the hundreds. We restrict ourselves to information that can be extracted from the inspiral, which includes tidal effects and possibly r-mode resonances. In treating the latter we go beyond the Newtonian approximation, introducing and utilizing new universal relations. We find that the ability to observe resonant r-modes would have a noticeable impact on neutron star equation of state measurements with ET.

B. Tripathi, A.E. Fraser, P.W. Terry, E.G. Zweibel, M.J. Pueschel

Letter submitted to Physics of Plasmas; A simulation movie embedded in Fig. 5; Supplementary text available at this https URL

Straining of magnetic fields by large-scale shear flow, generally assumed to lead to intensification and generation of small scales, is re-examined in light of the persistent observation of large-scale magnetic fields in astrophysics. It is shown that in magnetohydrodynamic turbulence, unstable shear flows have the unexpected effect of sequestering magnetic energy at large scales, due to counteracting straining motion of nonlinearly excited large-scale stable eigenmodes. This effect is quantified via dissipation rates, energy transfer rates, and visualizations of magnetic field evolution by artificially removing the stable modes.

H. M. Tedila, R. Yuen, N. Wang, J. P. Yuan, Z. G. Wen, W. M. Yan, S. Q. Wang, S. J. Dang, D. Li, P. Wang, W. W. Zhu, J. R. Niu, C. C. Miao, M. Y. Xue, L. Zhang, Z. Y. Tu, R. Rejep, J. T. Xie, FAST Collaboration

We report on the variation in the single-pulse emission from PSR J1900+4221 (CRAFTS 19C10) observed at frequency centered at 1.25 GHz using the Five-hundred-meter Aperture Spherical radio Telescope. The integrated pulse profile shows two distinct components, referred to here as the leading and trailing components, with the latter component also containing a third weak component. The single-pulse sequence reveals different emissions demonstrating as nulling, regular, and bright pulses, each with a particular abundance and duration distribution. There also exists pulses that follow a log-normal distribution suggesting the possibility of another emission, in which the pulsar is radiating weakly. Changes in the profile shape are seen across different emissions. We examine the emission variations in the leading and trailing components collectively and separately, and find moderate correlation between the two components. The inclination angle is estimated to be about 7{\deg} based on pulse-width, and we discuss that nulling in this pulsar does not seem to show correlation with age and rotation period.

We extend the scaling relations of stably stratified turbulence from the geophysical regime of unity Prandtl number to the astrophysical regime of extremely small Prandtl number applicable to stably stratified regions of stars and gas giants. A transition to a new turbulent regime is found to occur when the Prandtl number drops below the inverse of the buoyancy Reynolds number, i.e. $PrRb<1$, which signals a shift of the dominant balance in the buoyancy equation. Application of critical balance arguments then derives new predictions for the anisotropic energy spectrum and dominant balance of the Boussinesq equations in the $PrRb\ll1$ regime. We find that all the standard scaling relations from the unity $Pr$ limit of stably stratified turbulence simply carry over if the Froude number, $Fr$, is replaced by a modified Froude number, $Fr_M\equiv Fr/(PrRb)^{1/4}$. The geophysical and astrophysical regimes are thus smoothly connected across the $PrRb=1$ transition. Applications to vertical transport in stellar radiative zones and modification to the instability criterion for the small-scale dynamo are discussed.

F. Acernese, M. Agathos, A. Ain, S. Albanesi, A. Allocca, A. Amato, T. Andrade, N. Andres, M. Andrés-Carcasona, T. Andrić, S. Ansoldi, S. Antier, T. Apostolatos, E. Z. Appavuravther, M. Arène, N. Arnaud, M. Assiduo, S. Assis de Souza Melo, P. Astone, F. Aubin, S. Babak, F. Badaracco, M. K. M. Bader, S. Bagnasco, J. Baird, T. Baka, G. Ballardin, G. Baltus, B. Banerjee, C. Barbieri, P. Barneo, F. Barone, M. Barsuglia, D. Barta, A. Basti, M. Bawaj, M. Bazzan, F. Beirnaert, M. Bejger, I. Belahcene, V. Benedetto, M. Berbel, S. Bernuzzi, D. Bersanetti, A. Bertolini, U. Bhardwaj, A. Bianchi, S. Bini, M. Bischi, M. Bitossi, M.-A. Bizouard, F. Bobba, M. Boër, G. Bogaert, M. Boldrini, L. D. Bonavena, F. Bondu, R. Bonnand, B. A. Boom, V. Boschi, V. Boudart, Y. Bouffanais, A. Bozzi, et al. (431 additional authors not shown)

86 pages, 33 figures. To be submitted to Class. and Quantum Grav

The Advanced Virgo detector has contributed with its data to the rapid growth of the number of detected gravitational-wave signals in the past few years, alongside the two LIGO instruments. First, during the last month of the Observation Run 2 (O2) in August 2017 (with, most notably, the compact binary mergers GW170814 and GW170817) and then during the full Observation Run 3 (O3): an 11 months data taking period, between April 2019 and March 2020, that led to the addition of about 80 events to the catalog of transient gravitational-wave sources maintained by LIGO, Virgo and KAGRA. These discoveries and the manifold exploitation of the detected waveforms require an accurate characterization of the quality of the data, such as continuous study and monitoring of the detector noise. These activities, collectively named {\em detector characterization} or {\em DetChar}, span the whole workflow of the Virgo data, from the instrument front-end to the final analysis. They are described in details in the following article, with a focus on the associated tools, the results achieved by the Virgo DetChar group during the O3 run and the main prospects for future data-taking periods with an improved detector.

We consider the motion of a point particle in a stationary spacetime under the influence of a scalar, electromagnetic or gravitational self-force. We show that the conservative piece of the first-order self-force gives rise to Hamiltonian dynamics, and we derive an explicit expression for the Hamiltonian on phase space. Specialized to the Kerr spacetime, our result generalizes the Hamiltonian function previously obtained by Fujita et. al., which is valid only for non-resonant orbits. We discuss implications for the first law of binary black hole mechanics.

Shyam Balaji, P. S. Bhupal Dev, Joseph Silk, Yongchao Zhang

20 pages+references, 10 figures, 2 tables

We derive improved stellar luminosity limits on a generic light CP-even scalar field $S$ mixing with the Standard Model (SM) Higgs boson from the supernova SN1987A, the Sun, red giants (RGs) and white dwarfs (WDs). For the first time, we include the geometric effects for the decay and absorption of $S$ particles in the stellar interior. For SN1987A and the Sun, we also take into account the detailed stellar profiles. We find that a broad range of the scalar mass and mixing angle can be excluded by our updated astrophysical constraints. For instance, SN1987A excludes $1.0\times10^{-7} \lesssim \sin\theta \lesssim 4.0\times 10^{-5}$ and scalar mass up to 290 MeV, which covers the cosmological blind spot with a high reheating temperature. The updated solar limit excludes the mixing angle in the range of $3.8\times 10^{-11} < \sin\theta < 0.27$, with scalar mass up to 39 keV. The RG and WD limits are updated to $4.3\times 10^{-13} < \sin \theta < 8.6\times 10^{-5}$ and $2.3\times 10^{-18} < \sin \theta < 3.5\times 10^{-8}$, with scalar mass up to 398 keV and 292 keV, respectively.