Abstract visibility
Abstract text size

Papers for Friday, Dec 23 2022

Ryota Ikeda, Takahiro Morishita, Takafumi Tsukui, Benedetta Vulcani, Michele Trenti, Benjamin Metha, Ana Acebron, Pietro Bergamini, Claudio Grillo, Daisuke Iono, Amata Mercurio, Piero Rosati, Eros Vanzella

Submitted to MNRAS, 10 pages, 9 figures, 2 tables

We present a search and characterization of ultra-diffuse galaxies (UDGs) in the Frontier Fields cluster Abell 2744 at $z=0.308$. We use JWST/NIRISS F200W observations, acquired as part of the GLASS-JWST Early Release Science Program, aiming to characterize morphologies of cluster UDGs and their diffuse stellar components. A total number of 22 UDGs are identified by our selection criteria using morphological parameters, down to stellar mass of $\sim10^{7}M_{\odot}$. The selected UDGs are systematically larger in effective radius in F200W than in HST/ACS F814W images, which implies that some of them would not have been identified as UDGs when selected at rest-frame optical wavelengths. In fact, we find that about one third of the UDGs were not previously identified based on the F814W data. We observe a flat distribution of the UDGs in the stellar mass-size plane, similar to what is found for cluster quiescent galaxies at comparable mass. We also find 10 potential candidates with disturbed morphologies, a previously overlooked but important population as a possible progenitor of the local UDGs. Our pilot study using the new JWST F200W filter showcases the efficiency of searching UDGs at cosmological distances, with 1/30 of the exposure time of the previous deep observing campaign with HST. Further studies with JWST focusing on spatially-resolved properties of individual sources will provide insight into their origin.

Aleksey Generozov, Hagai B. Perets

10 pages, 5 figures, submitted to MNRAS

The majority of massive black holes (MBHs) likely hosted gas discs during their lifetimes. These could either be long-lived active galactic nuclei (AGN) discs, or shorter-lived discs formed following singular gas infall events, as was likely the case in our own Galactic Center. Stars and compact objects in such environments are therefore expected to interact with the gaseous disc as they go through it, and potentially become aligned and fully embedded within it. The interactions of embedded stars with the gas could give rise to a plethora of physical processes affecting the stars, including growth through accretion of gas, migration in the disc, stellar captures, and mergers with other stars. The impact of such processes strongly depends on the population of stars that eventually align with the disc and become embedded in it. Here we make use of analytic tools to analyze the alignment process, accounting for both geometric drag and gas dynamical friction. We find that up to $\sim$10% of stars and stellar mass black holes can align with AGN disc in the Galactic Center. The orbits of aligned stars are typically circularized and are prograde with respect to the AGN disc. Furthermore, alignment and accretion are intimately linked, and the capture of stars by an AGN disc can potentially explain the origin of the young stellar disc in the Galactic Center with a top-heavy mass function, even without the need for a star-formation event.

Arvind K. Dattatrey, R.K.S. Yadav, Sharmila Rani, Annapurni Subramaniam, Gaurav Singh, Snehalata Sahu, Ravi S. Singh

16 pages, 10 figures, accepted for publication in ApJ

We report the discovery of extremely low-mass white dwarfs (ELM WDs) as a companion of blue straggler stars (BSSs) in the Galactic globular cluster NGC 362 using images from AstroSats Ultra Violet Imaging Telescope (UVIT). Spectral Energy Distributions (SEDs) for 26 FUV bright member BSSs are created using data from the UVIT, UVOT, Gaia EDR3, and the 2.2 m ESO/MPI telescope. A single SED is fitted to 14 BSSs, whereas double-SED fits revealed ELM WDs as binary companions in 12 of the 26 BSSs studied. The effective temperature, radius, luminosity and mass of the 12 ELM WDs are found to have a range (Teff = 9750-18000 K, R = 0.1-0.4 Rsun, L = 0.4-3.3 Lsun, and M =0.16-0.20 Msun). These suggest that 12 BSSs are post-mass-transfer systems formed through Case A/B mass transfer pathway. To the best of our knowledge, this is the first finding of ELM WDs as companions to BSS in globular clusters. This cluster is known to have a binary BSS sequence and the 12 binary and 14 single BSSs (as classified by the SEDs) follow the mass transfer and collisional sequence of BSS in the colour-magnitude diagram. The cooling ages of 9 of the ELM WDs are found to be younger than 500 Myr. Though the binary BSSs may have formed during the core-collapse (approx 200 Myr) or as part of the dynamical evolution of the cluster, they provide new insights on the dynamics of this cluster.

Emília Valença Ferreira de Aragão, Luca Mancini, Noelia Faginas-Lago, Marzio Rosi, Nadia Balucani, Fernando Pirani

13 pages, 5 figures. Preprint version submitted to LNCS (Springer) ICCSA2021

In this work we characterize an initial van der Waals adduct in the potential energy surface of reaction between cyanoacetylene HC3N and the cyano radical. The geometry of the CN-HC3N adduct has been optimized through calculations employing ab initio methods. Results show that the energy of the adduct lays below the reactants. Additionally, a saddle point that connects the adduct to an important intermediate of the PES has been localized, with energy below the reactants. Calculations of the intermolecular potential have been performed and results show that the energy of the van der Waals adduct is higher than estimated with the ab initio methods.

Prathamesh D. Tamhane, Brian R. McNamara, Helen R. Russell, Francoise Combes, Yu Qiu, Alastair C. Edge, Roberto Maiolino, Andrew C. Fabian, Paul E.J. Nulsen, R. Johnstone, Stefano Carniani

Accepted for publication in MNRAS. 22 pages, 13 figures 1 table

We present XSHOOTER observations with previous ALMA, MUSE and $HST$ observations to study the nature of radio-jet triggered star formation and the interaction of radio jets with the interstellar medium in the brightest cluster galaxy (BCG) in the Abell 1795 cluster. Using $HST$ UV data we determined an ongoing star formation rate of 9.3 M$_\odot$ yr$^{-1}$. The star formation follows the global Kennicutt-Schmidt law, however, it has a low efficiency compared to circumnuclear starbursts in nearby galaxies with an average depletion time of $\sim$1 Gyr. The star formation and molecular gas are offset by $\sim1$ kpc indicating that stars have decoupled from the gas. We detected an arc of high linewidth in ionized gas where electron densities are elevated by a factor of $\sim$4 suggesting a shock front driven by radio jets or peculiar motion of the BCG. An analysis of nebular emission line flux ratios suggests that the gas is predominantly ionized by star formation with a small contribution from shocks. We also calculated the velocity structure function (VSF) of the ionized and molecular gases using velocity maps to characterize turbulent motion in the gas. The ionized gas VSF suggests that the radio jets are driving supersonic turbulence in the gas. Thus radio jets can not only heat the atmosphere on large scales and may quench star formation on longer timescales while triggering star formation in positive feedback on short timescales of a few million years.

S. Samuroff, R. Mandelbaum, J. Blazek, A. Campos, N. MacCrann, G. Zacharegkas, A. Amon, J. Prat, S. Singh, J. Elvin-Poole, A. J. Ross, A. Alarcon, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, C. Chang, R. Chen, A. Choi, M. Crocce, C. Davis, J. DeRose, S. Dodelson, C. Doux, A. Drlica-Wagner, K. Eckert, S. Everett, A. Ferté, M. Gatti, G. Giannini, D. Gruen, R. A. Gruendl, I. Harrison, K. Herner, E. M. Huff, M. Jarvis, N. Kuropatkin, P.-F. Leget, P. Lemos, J. McCullough, J. Myles, A. Navarro-Alsina, S. Pandey, A. Porredon, M. Raveri, M. Rodriguez-Monroy, R. P. Rollins, A. Roodman, G. Rossi, E. S. Rykoff, C. Sánchez, L. F. Secco, I. Sevilla-Noarbe, E. Sheldon, T. Shin, M. A. Troxel, I. Tutusaus, N. Weaverdyck, B. Yanny, B. Yin, et al. (44 additional authors not shown)

Submitted to MNRAS; 22 pages, 11 figures

We present direct constraints on galaxy intrinsic alignments using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift $z\sim0.2-1.0$, luminous red galaxies (LRGs) from eBOSS at $z\sim0.8$, and also a SDSS-III BOSS CMASS sample at $z\sim0.5$. We measure two point intrinsic alignment correlations, which we fit using a model that includes lensing, magnification and photometric redshift error. Fitting on scales $6<r_{\rm p} < 70$ Mpc$/h$, we make a detection of intrinsic alignments in each sample, at $5\sigma-22\sigma$ (assuming a simple one parameter model for IAs). Using these red samples, we measure the IA-luminosity relation. Our results are statistically consistent with previous results, but offer a significant improvement in constraining power, particularly at low luminosity. With this improved precision, we see detectable dependence on colour between broadly defined red samples. It is likely that a more sophisticated approach than a binary red/blue split, which jointly considers colour and luminosity dependence in the IA signal, will be needed in future. We also compare the various signal components at the best fitting point in parameter space for each sample, and find that magnification and lensing contribute $\sim2-18\%$ of the total signal. As precision continues to improve, it will certainly be necessary to account for these effects in future direct IA measurements. Finally, we make equivalent measurements on a sample of Emission Line Galaxies (ELGs) from eBOSS at $z\sim 0.8$. We report a null detection, constraining the IA amplitude (assuming the nonlinear alignment model) to be $A_1=0.07^{+0.32}_{-0.42}$ ($|A_1|<0.78$ at $95\%$ CL).

Muriel Zoë Stiefel, Andrea Francesco Battaglia, Krzysztof Barczynski, Hannah Collier, Anna Volpara, Paolo Massa, Conrad Schwanitz, Sofia Tynelius, Louise Harra, Säm Krucker

12 pages, 10 figures, accepted for publication in A&A

Context. We present an analysis of a GOES M1.8 flare with excellent observational coverage in UV, EUV, and X-ray, including observations from the instruments IRIS, SDO with AIA, Hinode/EIS, Hinode/XRT, and Solar Orbiter with the Spectrometer/Telescope for Imaging X-rays (STIX). Hard X-ray emission is often observed at the footpoints of flare loops and is occasionally observed in the corona. In this flare, four nonthermal hard X-ray sources are seen. Aim. Our aim is to understand why we can observe four individual nonthermal sources in this flare and how we can characterize the physical properties of these four sources. Methods. We used the multiwavelength approach to analyze the flare and characterize the four sources. To do this, we combined imaging at different wavelengths and spectroscopic fitting in the EUV and X-ray range. Results. The flare is eruptive with an associated coronal mass ejection, and it shows the classical flare picture of a heated flare loop seen in EUV and X-rays, and two nonthermal hard X-ray footpoints at the loop ends. In addition to the main flare sources, we observed two outer sources in the UV, EUV, and nonthermal X-ray range located away from the main flare loop to the east and west. The two outer sources are clearly correlated in time, and they are only seen during the first two minutes of the impulsive phase, which lasts a total of about four minutes. Conclusions. Based on the analysis, we determine that the outer sources are the anchor points of an erupting filament. The hard X-ray emission is interpreted as flare-accelerated electrons that are injected upward into the filament and then precipitate along the filament toward the chromosphere, producing Bremsstrahlung. While sources like this have been speculated to exist, this is the first report of their detection.

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

32 pages, 14 figures, accepted for publication in Galaxies

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.

Curtis McCully, Matthew Daily, G. Mirek Brandt, Marshall C. Johnson, Mark Bowman, Daniel-Rolf Harbeck

12 pages, 8 figures, SPIE Proceedings 2022

Time domain astronomy has both increased the data volume and the urgency of data reduction in recent years. Spectra provide key insights into astrophysical phenomena but require complex reductions. Las Cumbres Observatory has six spectrographs: two low-dispersion FLOYDS instruments and four NRES high-resolution echelle spectrographs. We present an extension of the data reduction framework, BANZAI, to process spectra automatically, with no human interaction. We also present interactive tools we have developed for human vetting and improvement of the spectroscopic reduction. Tools like those presented here are essential to maximize the scientific yield from current and future time domain astronomy.

Mathieu Vrard, Margarida S. Cunha, Diego Bossini, Pedro P. Avelino, Enrico Corsaro, Benoit Mosser

30 pages, 4 figures, 1 table

Red giants are stars in the late stages of stellar evolution. Because they have exhausted the supply of hydrogen in their core, they burn the hydrogen in the surrounding shell. Once the helium in the core starts fusing, the star enters the clump phase, which is identified as a striking feature in the color-magnitude diagram. Since clump stars share similar observational properties, they are heavily used in astrophysical studies, as probes of distance, extinction through the galaxy, galaxy density, and stellar chemical evolution. In this work, we perform the detailed observational characterization of the deepest layers of clump stars using asteroseismic data from Kepler. We find evidence for large core structural discontinuities in about 6.7% of the stars in our sample, implying that the region of mixing beyond the convective core boundary has a radiative thermal stratification. These stars are otherwise similar to the remaining stars in our sample, which may indicate that the building of the discontinuities is an intermittent phenomenon.

Michail Damoulakis, Rodolfo Barniol Duran, Dimitrios Giannios

13 pages, 14 figures, 1 table. Submitted to MNRAS. Corresponding author: Michail Damoulakis

Despite a generally accepted framework for describing the Gamma-Ray Burst (GRB) afterglows, the nature of the compact object at the central engine and the mechanism behind the prompt emission remain debated. The striped jet model is a promising venue to connect the various GRB stages since it gives a robust prediction for the relation of jet bulk acceleration, magnetization and dissipation profile as a function of distance. Here, we use the constraints of the magnetization and bulk Lorentz of the jet flow at the large scales where the jet starts interacting with the ambient gas in a large sample of bursts to (i) test the striped jet model for the GRB flow and (ii) study its predictions for the prompt emission and the constraints on the nature of the central engine. We find that the peak of the photospheric component of the emission predicted by the model is in agreement with the observed prompt emission spectra in the majority of the bursts in our sample, with a radiative efficiency of about 10 per cent. Furthermore, we adopt two different approaches to correlate the peak energies of the bursts with the type of central engine to find that more bursts are compatible with a neutron star central engine compared to a black hole one. Lastly, we conclude that the model favors broader distribution of stripe length-scales which results in a more gradual dissipation profile in comparison to the case where the jet stripes are characterized by a single length-scale.

Peter F. Bloser (Los Alamos National Laboratory), David Murphy (Centre for Space Research and School of Physics, University College Dublin), Fabrizio Fiore (INAF, Osservatorio Astronomico di Trieste), Jeremy Perkins (NASA Goddard Space Flight Center)

Book chapter for the "Handbook of X-ray and Gamma-ray Astrophysics", Section "Optics and Detectors for Gamma-ray Astrophysics" (Editors in chief: C. Bambi and A. Santangelo, Springer Singapore). 33 pages, 11 figures

After many years of flying in space primarily for educational purposes, CubeSats - tiny satellites with form factors corresponding to arrangements of "1U" units, or cubes, each 10 cm on a side - have come into their own as valuable platforms for technology advancement and scientific investigations. CubeSats offer comparatively rapid, low-cost access to space for payloads that be built, tested, and operated by relatively small teams, with substantial contributions from students and early career researchers. Continuing advances in compact, low-power detectors, readout electronics, and flight computers have now enabled X-ray and gamma-ray sensing payloads that can fit within the constraints of CubeSat missions, permitting in-orbit demonstrations of new techniques and innovative high-energy astronomy observations. Gamma-ray-sensing CubeSats are certain to make an important contribution in the new era of multi-messenger, time-domain astronomy by detecting and localizing bright transients such as gamma-ray bursts, solar flares, and terrestrial gamma-ray flashes; however, other astrophysical science areas requiring long observations in a low-background environment, including gamma-ray polarimetry, studies of nuclear lines, and measurement of diffuse backgrounds, will likely benefit as well. We present the primary benefits of CubeSats for high-energy astronomy, highlight the scientific areas currently or soon to be studied, and review the missions that are currently operating, under development, or proposed. A rich portfolio of CubeSats for gamma-ray astronomy already exists, and the potential for a broad range of creative and scientifically productive missions in the near future is very high.

Sebastian Zieba, Laura Kreidberg

8 pages, 2 figures, Published in JOSS, GitHub: this https URL

Here we present PACMAN, an end-to-end pipeline developed to reduce and analyze HST/WFC3 data. The pipeline includes both spectral extraction and light curve fitting. The foundation of PACMAN has been already used in numerous publications (e.g., Kreidberg et al., 2014; Kreidberg et al., 2018) and these papers have already accumulated hundreds of citations. The Hubble Space Telescope (HST) has become the preeminent workhorse facility for the characterization of extrasolar planets. HST currently has two of the most powerful space-based tools for characterizing exoplanets over a broad spectral range: The Space Telescope Imaging Spectrograph (STIS) in the UV and the Wide Field Camera 3 (WFC3) in the Near Infrared. With the introduction of a spatial scan mode on WFC3 where the star moves perpendicular to the dispersion direction during an exposure, WFC3 observations have become very efficient due to the reduction of overhead time and the possibility of longer exposures without saturation. For exoplanet characterization, WFC3 is used for transit and secondary eclipse spectroscopy, and phase curve observations. The instrument has two different grisms: G102 with a spectral range from 800 nm to up to 1150 nm and G141 encompassing 1075 nm to about 1700 nm. The spectral range of WFC3/G141 is primarily sensitive to molecular absorption from water at approximately 1.4 microns. This led to the successful detection of water in the atmosphere of over a dozen of exoplanets. The bluer part of WFC3, the G102 grism, is also sensitive to water and most notably led to the first detection of a helium exosphere.

We report near-infrared (2.5--5 micron) long-slit (~ 30 arcsec) spectroscopy of a young stellar object in the direction toward the Galactic center with the Infrared Camera on board the AKARI satellite. The present target is suggested to be AFGL 2006 based on its very red color and close location. The spectra show strong absorption features of H$_2$O and CO$_2$ ices, and emission of HI Br alpha recombination line and the 3.3 micron band, the latter of which originates from polycyclic aromatic hydrocarbons (PAHs) or materials containing PAHs. The spectra show a broad, complex absorption feature at 4.65 micron, which is well explained by a combination of absorption features of CO ice, CO gas, and XCN, and HI Pf beta emission. The spectra also indicate excess emission at 4.4 micron. The characteristics of the spectra suggest that the object is a massive young stellar object. The XCN feature shows a good correlation with the Br alpha emission, suggesting that the photolysis by ultraviolet photons plays an important role in the formation of the XCN carriers, part of which are attributed to OCN$^-$. The 4.4 micron emission shows a good correlation with the 3.3 micron PAH emission, providing supporting evidence that it comes from the aromatic C-D stretching vibration. The formation of OCN$^-$ is of importance for the formation process of prebiotic matter in the interstellar medium (ISM), while the detection of aromatic C-D emission provides valuable information on the deuteration process of PAHs in the ISM and implications on the hiding site of the missing deuterium in the ISM.

Kirsten L. Larson, Janice C. Lee, David A. Thilker, Bradley C. Whitmore, Sinan Deger, James Lilly, Rupali Chandar, Daniel A. Dale, Frank Bigiel, Kathryn Grasha, Brent Groves, Ralf S. Klessen, Kathryn Kreckel, J. M. Diederik Kruijssen, Adam K. Leroy, Hsi-An Pan, Erik Rosolowsky, Eva Schinnerer, Andreas Schruba, Elizabeth J. Watkins, Thomas G. Williams

Submitted to MNRAS. Referee report received with minor comments, and "request to clarify if the smaller associations are always included in the larger ones and how this may affect the photometric fitting of the larger association if the groups have different ages." Revision in progress

We develop a method to identify and determine the physical properties of stellar associations using Hubble Space Telescope (HST) NUV-U-B-V-I imaging of nearby galaxies from the PHANGS-HST survey. We apply a watershed algorithm to density maps constructed from point source catalogues Gaussian smoothed to multiple physical scales from 8 to 64 pc. We develop our method on two galaxies that span the distance range in the PHANGS-HST sample: NGC 3351 (10 Mpc), NGC 1566 (18 Mpc). We test our algorithm with different parameters such as the choice of detection band for the point source catalogue (NUV or V), source density image filtering methods, and absolute magnitude limits. We characterise the properties of the resulting multi-scale associations, including sizes, number of tracer stars, number of associations, photometry, as well as ages, masses, and reddening from Spectral Energy Distribution fitting. Our method successfully identifies structures that occupy loci in the UBVI colour-colour diagram consistent with previously published catalogues of clusters and associations. The median ages of the associations increases from log(age/yr) = 6.6 to log(age/yr) = 6.9 as the spatial scale increases from 8 pc to 64 pc for both galaxies. We find that the youngest stellar associations, with ages < 3 Myr, indeed closely trace H ii regions in H$\alpha$ imaging, and that older associations are increasingly anti-correlated with the H$\alpha$ emission. Owing to our new method, the PHANGS-HST multi-scale associations provide a far more complete census of recent star formation activity than found with previous cluster and compact association catalogues. The method presented here will be applied to the full sample of 38 PHANGS-HST galaxies.

Xi Liu, Xiaolei Guo, Yuliang Xin, Fengrong Zhu, Siming Liu

9 pages, 5 figures

HESS J1356-645 is considered to be a pulsar wind nebula (PWN) associated with the pulsar PSR J1357-6429. We reanalyze the GeV gamma-ray emission in the direction of HESS J1356-645 with more than 13 years of Fermi Large Area Telescope (LAT) data. The extended gamma-ray emission above 5 GeV is found to be spatially coincident with HESS J1356-645. The spectrum in the energy range of 1 GeV-1 TeV can be described by a power law with an index of $\Gamma=1.51\pm0.10$. The broadband spectrum of HESS J1356-645 can be reproduced by a leptonic model with a broken power-law electronic spectrum. In addition, we found evidence that the morphology of the GeV emission from HESS J1356-645 varies with energy, a behavior which is similar to that of the PWN Vela-X. More broadband observations will be helpful to study the energy-dependent characteristics of HESS J1356-645.

Patrícia da Silva, R. B. Menezes, Y. Díaz, A. Rodríguez-Ardila, E. López-Navas, P. Arévalo, L. Hernández-García

20 pages, 16 figures, 7 tables

The central regions of galaxies harbouring active galactic nuclei (AGNs) can be quite complex, especially at high activity, presenting, besides variability, a variety of phenomena related, e.g. to ionization/excitation mechanisms. A detailed study is necessary in order to understand better those objects. For that reason, we performed a multiwavelength analysis of the nuclear region of the nearby Seyfert galaxy NGC 7314, using an optical data cube obtained with the Integral Field Unit from the Gemini Multi-Object Spectrograph, together with Hubble Space Telescope images, X-ray data from the XMM-Newton and the Nuclear Spectroscopic Telescope Array and radio data from Atacama Large Millimeter/Submillimeter Array. The goals were to study the nuclear and circumnuclear emission, the emission of the AGN and the gas kinematics. The optical spectrum shows the emission of a Seyfert nucleus, with broad components in the H$\alpha$ and H$\beta$ emission lines, characterising a type 1 AGN, with a spectrum rich in coronal emission lines. The spatial morphology of the [OIII]$\lambda$5007 suggests the presence of an ionization cone, west of the nucleus, meanwhile the east cone seems to be obscured by dust. An extended [FeVII]$\lambda$6087 emission was also detected, which could be possibly explained by a scenario involving photoionization+shocks mechanisms. X-rays analyses showed that there are variations in the flux; however, we did not detect any variations in the column density along the line of sight. Its variability may be a consequence of changes in the AGN accretion rate.

Leaf Swordy, Heidi Jo Newberg, Thomas Ditto

25 pages, 12 figures, submitted to JATIS

We develop the theoretical foundation for primary objective grating (POG) telescopy. In recent years, a wide range of telescope designs that collect the light over a large grating and focus it with a secondary receiving optic that is placed at grazing exodus have been proposed by Thomas D. Ditto, and are sometimes referred to as Dittoscopes. Applications include discovery and characterization of exoplanets, discovery of near-Earth asteroids, and spectroscopic surveys of the sky. These telescopes would have small aerial mass, and therefore provide a path forward to launch large telescopes into space. Because this series of telescope designs departs from traditional telescope designs, it has been difficult to evaluate which applications are most advantageous for this design. Here, we define a new figure of merit, the "modified etendue," that characterizes the photon collection capability of a POG. It is demonstrated that the diffraction limit for observations is determined by the length of the grating. We evaluate the effects of atmospheric seeing for ground-based applications and the disambiguation of position vs. wavelength in the focal plane using a second dispersing element. Finally, some strategies for fully reaping the benefits of POG optical characteristics are discussed.

This work focuses on providing closed form analytical expressions to define frozen orbits under the effects of the zonal harmonics of an Earth-like planet. Particularly, the perturbation effects from the terms J2, J3, J4, J5, J6, and J7 are considered in this work. This is done using a power series expansion in the small parameter that allow not only to provide an approximate solution to the system, but also to determine the analytical expressions that define the initial osculating conditions that generate frozen orbits. As a result of that, the proposed methodology allows to study the bifurcation of frozen orbits close to the critical inclination by purely analytical methods. Additionally, the derivation of the secular variation of the orbital elements as well as the transformation from osculating to mean elements is provided based on the second order analytical solution proposed in this work. Examples of application are also provided to show the error performance of the results included in this document.

We investigate the cosmic evolutions in the extended Starobinsky model (eSM) obtained by adding one R ab R ab term to the Starobinsky model. We discuss the possibility of various forms of cosmic evolution with a special focus on the radiation-dominated era (RDE). Using simple assumptions, a second-order non-linear differential equation describing the various cosmic evolutions in the eSM is introduced. By solving this non-linear equation numerically, we show that various forms of cosmic evolution, such as the standard cosmic evolution ($a \propto t^{1/2}$) and a unique oscillating cosmic evolution, are feasible due to the effects of higher-order terms introduced beyond Einstein gravity. Furthermore, we consider big bang nucleosynthesis (BBN), which is the most important observational results in the RDE, to constrain the free parameters of the eSM. The primordial abundances of the light elements, such as $^4$He, D, $^3$He, $^7$Li and $^6$Li by the cosmic evolutions are compared with the most recent observation data. It turns out that most cases of non-standard cosmic evolutions can not easily satisfy these BBN constraints, but a free parameter of the viable models with the oscillating cosmic evolution is shown to have an upper limit by the constraints. In particular, we find that the free parameter is most sensitive to deuterium and $^4$He abundances, which are being precisely measured among other elements. Therefore, more accurate measurements in the near future may enable us to distinguish the eSM from the standard model as well as other models.

We investigate the dependence of the angular momentum transport (AMT) on the spatial scales with numerical simulation of solar-like stars. It is thought that turbulence has an essential role in constructing solar differential rotation (DR). In a widely used method to analyse the construction mechanism of DR, the flow is divided into two components, `mean flow' and `turbulence', where `turbulence' includes a broad spectrum of spatial scales. The features of the AMT are expected to depend on the scale. In this study, we decompose the angular momentum flux (AMF) to investigate the dependence of the AMF on the spatial scale. We compare the results with anti-solar- (fast pole) and solar-type (fast equator) DR. Our conclusions are summarized as 1. Radially outward AMT is seen on a large scale (60 Mm < L <120 Mm) in rotationally constrained systems. 2. Even when the scale-integrated AMF is negative, we sometimes observe positive AMF on certain scales. 3. Small-scale turbulence tends to transport the angular momentum radially inward and causes the anti-solar DR, indicating that high-resolution simulation is a negative factor for solar-like DR. Our method to decompose the AMF provides a deep understanding of the angular momentum and construction mechanism of DR.

Maokai Hu, Lei Hu, Ji-an Jiang, Lin Xiao, Lulu Fan, Junjie Wei, Xuefeng Wu

Accepted by Universe

Type Ia Supernovae (SNe Ia) are the thermonuclear explosion of a carbon-oxygen white dwarf (WD) and are well-known as a distance indicator. However, it is still unclear how WDs increase their mass near the Chandrasekhar limit and how the thermonuclear runaway happens. The observational clues associated with these open questions, such as the photometric data within hours to days since the explosion, are scarce. Thus, an essential way is to discover SNe Ia at specific epochs with optimal surveys. The 2.5-m Wide Field Survey Telescope (WFST) is an upcoming survey facility deployed in western China. In this paper, we assess the detecability of SNe Ia with mock observations of WFST. Followed by the volumetric rate, we generate a spectral series of SNe Ia based on a data-based model and introduce the line-of-sight extinction to calculate the brightness from the observer. By comparing with the detection limit of WFST, which is affected by the observing conditions, we can count the number of SNe Ia discovered by mock WFST observations. We expect that WFST can find more than $3.0\times10^{4}$ pre-maximum SNe Ia within one-year running. In particular, WFST could discover about 45 bright SNe Ia, 99 early-phase SNe Ia, or $1.1\times10^{4}$ well-observed SNe Ia with the hypothesized Wide, Deep, or Medium mode, respectively, suggesting WFST will be an influential facility in time-domain astronomy.

Solar wind back-mapping is a combination of ballistic mapping and magnetic mapping. By examining the different model ingredients that can affect the derived back-mapped position, we aim to provide a more precise estimate of the source location and a measure of confidence in the mapping procedure. This can be used to improve the connection of remote sensing with in situ measurements. For the ballistic mapping we created custom velocity profiles. These profiles are constrained by observations of the fast solar wind close to the Sun and are used to examine the mapping uncertainty. The coronal magnetic field topology from the solar surface up to the source surface is modeled with a PFSS extrapolation. The sensitivity of the extrapolated field is examined by adding noise to the input magnetogram and performing a Monte Carlo simulation, where for multiple noise realizations we calculate the source position of the solar wind. Next, the effect of free parameters, like the height of the source surface, is examined and statistical estimates are derived. We used Gaussian Mixture clustering to group the back-mapped points, due to different sources of uncertainty, and provide a confidence area for the source location of the solar wind. Furthermore, we computed a number of metrics to evaluate the back-mapping results and assessed their statistical significance by examining 3 high speed stream events. Lastly, we explored the effect of corotation, close to the Sun, on the source region of the solar wind. Our results show that the height of the source surface produces the largest uncertainty in the source region of the fast solar wind, followed by the choice of the velocity profile and the noise in the input magnetogram. Additionally, we display the ability to derive a confidence area on the solar surface that represents the potential source region of the in-situ measured fast solar wind.

Huanbin Chi, Shoulin Wei, Feng Wang, Zhongmu Li

16 Pages, 14 figures, 4 tables, accepted by APJs

Open clusters (OCs) are regarded as tracers to understand stellar evolution theory and validate stellar models. In this study, we presented a robust approach to identifying OCs. A hybrid method of pyUPMASK and RF is first used to remove field stars and determine more reliable members. An identification model based on the RF algorithm built based on 3714 OC samples from Gaia DR2 and EDR3 is then applied to identify OC candidates. The OC candidates are obtained after isochrone fitting, the advanced stellar population synthesis (ASPS) model fitting, and visual inspection. Using the proposed approach, we revisited 868 candidates and preliminarily clustered them by the friends-of-friends algorithm in Gaia EDR3. Excluding the open clusters that have already been reported, we focused on the remaining 300 unknown candidates. From high to low fitting quality, these unrevealed candidates were further classified into Class A (59), Class B (21), and Class C (220), respectively. As a result, 46 new reliable open cluster candidates among classes A and B are identified after visual inspection.

Ka-Yui Au, Jay Strader, Samuel J. Swihart, Lupin C. C. Lin, Albert K. H. Kong, Jumpei Takata, Chung-Yue Hui, Teresa Panurach, Isabella Molina, Elias Aydi, Kirill Sokolovsky, Kwan-Lok Li

Accepted for publication in ApJ

We present the study of multi-wavelength observations of an unidentified Fermi Large Area Telescope (LAT) source, 4FGL J1910.7-5320, a new candidate redback millisecond pulsar binary. In the 4FGL 95% error region of 4FGL J1910.7-5320, we find a possible binary with a 8.36-hr orbital period from the Catalina Real-Time Transient Survey (CRTS), confirmed by optical spectroscopy using the SOAR telescope. This optical source was recently independently discovered as a redback pulsar by the TRAPUM project, confirming our prediction. We fit the optical spectral energy distributions of 4FGL J1910.7-5320 with a blackbody model, inferring a maximum distance of 4.1 kpc by assuming that the companion fills its Roche-lobe with a radius of R = 0.7R_sun. Using a 12.6 ks Chandra X-ray observation, we identified an X-ray counterpart for 4FGL J1910.7-5320, with a spectrum that can be described by an absorbed power-law with a photon index of 1.0+/-0.4. The spectrally hard X-ray emission shows tentative evidence for orbital variability. Using more than 12 years of Fermi-LAT data, we refined the position of the {\gamma}-ray source, and the optical candidate still lies within the 68% positional error circle. In addition to 4FGL J1910.7-5320, we find a variable optical source with a periodic signal of 4.28-hr inside the 4FGL catalog 95% error region of another unidentified Fermi source, 4FGL J2029.5-4237. However, the {\gamma}-ray source does not have a significant X-ray counterpart in a 11.7 ks Chandra observation, with a 3-{\sigma} flux upper limit of 2.4*10^-14 erg cm^-2 s^-1 (0.3-7 keV). Moreover, the optical source is outside our updated Fermi-LAT 95% error circle. These observational facts all suggest that this new redback millisecond pulsar powers the {\gamma}-ray source 4FGL J1910.7-5320 while 4FGL J2029.5-4237 is unlikely the {\gamma}-ray counterpart to the 4.28-hr variable.

Themiya Nanayakkara, Karl Glazebrook, Colin Jacobs, Corentin Schreiber, Gabriel Brammer, James Esdaile, Glenn G. Kacprzak, Ivo Labbe, Claudia Lagos, Danilo Marchesini, Z. Cemile Marsan, Hasti Nateghi, Pascal A. Oesch, Casey Papovich, Rhea-Silvia Remus, Kim-Vy H. Tran

Submitted

The confirmation of the presence of very massive quiescent galaxies at epochs only 1-2 Gyr after the Big Bang [1-8] has challenged models of cosmology and galaxy formation [9]. Producing sufficient numbers of these requires abundant numbers of the host dark matter halos to have been assembled and sufficient time for star formation to proceed extremely quickly and then cease just as rapidly. Ground-based spectroscopy has suggested ages of 200-300 Myr[3] at redshifts $3<z<4$. The true number and ages of these objects have however been highly uncertain as ground-based spectra has been limited to the brightest of them [e.g. 3, 5], at wavelengths ~2$\mu$m, which introduces a signficant potential bias towards younger objects [7]. The launch of the James Webb Space Telescope (JWST) enables dramatically more sensitive and constraining spectroscopic observations due to the very low sky background, sharp image quality, and access to wavelengths beyond 2{\mu}m. Here we report JWST NIRSpec [10] (0.6-5.3$\mu$m) observations of five new quiescent galaxy candidates that were beyond the limit of previous ground-based spectroscopy. The high signal:noise spectra of galaxies with continuum significantly fainter than earlier confirmations show that they are also at redshifts 3<z<4, and that they have substantial stellar masses of ~0.5-1.2x1011 M$\odot$ comparable to massive galaxies in the nearby Universe. One of the galaxies has been quenched for >~ 1 billion years pointing to a presence of substantially older and fainter galaxies than those revealed so far by ground-based spectroscopy. This suggests that some of the massive galaxies have very early formation epochs (during the epoch of reionization, z >~ 6) pointing to a need for high conversion rates of baryons to stars in the first massive galaxy halos in the early Universe [11, 12].

Tadafumi Matsuno, Else Starkenburg, Eduardo Balbinot, Amina Helmi

The catalog will be publicly made available when accepted, or before upon reasonable request

Context: In the latest Gaia data release (DR3), the GSP-Spec module has provided stellar parameters and chemical abundances measured from the RVS spectra alone. However, very metal-poor stars (VMP stars; $[\mathrm{Fe/H}]<-2$) suffer from parameter degeneracy due to a lack of information in their spectra, making it difficult to obtain reliable stellar parameters and metallicities for many of them. Aim: We aim to improve metallicity estimates for VMP stars analysed by the GSP-Spec module. Methods: We compute the Ca triplet equivalent widths from the published set of GSP-Spec stellar parameters. We then convert these obtained equivalent widths to metallicities adopting photometric temperatures and surface gravities that we derive based on Gaia and 2MASS catalogs. Results: Comparison to high-resolution studies shows that our approach drastically reduces the cases where the estimated metallicities are far off for VMP stars. Now only $23\%$ of VMP stars have a metallicity different by more than $0.5\,\mathrm{dex}$ from a high-resolution value, while this fraction is $76\%$ with the original metallicity estimates by the GSP-Spec module. We explore possible ways to remove stars with poor metallicity estimates while keeping as many stars with reliable metallicity as possible. Our improved metallicity estimates and new quality cut result in producing a catalog of bright ($G\lesssim 13$) metal-poor stars, containing 57 stars at $[\mathrm{Fe/H}]<-3$ and 2202 VMP stars. These numbers increase to 174 and 2794 if one allows a low level of metal-rich contaminants. Conclusion: The inclusion of photometric information greatly contributes to breaking parameter degeneracy, enabling precise metallicity estimates for VMP stars from Gaia RVS spectra. We produce a publicly-available catalog of bright metal-poor stars suitable for high-resolution follow-up.

Shubham Yadav, M. Mishra, Tapomoy Guha Sarkar, Captain R. Singh

12 pages, 27 figures

Emission properties of compact astrophysical objects such as Neutron Stars have been found to be associated with crucial astronomical observables. In the current work, we obtain the mass, pressure and baryon number density profiles of the non-rotating neutron stars using the modified Tolman Oppenheimer Volkoff (TOV) system of equations in the presence of an intense radially distributed magnetic field. Employing the above profiles, we have determined the cooling rates of spherically symmetric neutron stars as a function of time with and without including the magnetic field using the NSCool code. We used a specific distance-dependent magnetic field in the modified TOV equations to obtain the profiles. We employ three different equation of states to solve the TOV equations by assuming the core of Neutron Stars to be made up of a hadronic matter. Using the above profiles, the cooling rate of neutron stars is obtained by employing NSCool code. Furthermore, based on the cooling rate, we determine the luminosity of Neutrinos, Axions and Photons emitting from the neutron stars in the presence and absence of a magnetic field for different axion masses and three equations of states. Our comparative study indicates that the colling rate and luminosities of axions, photons and neutrinos changes significantly due to the impact of magnetic field.

Valentina Cesare, Ugo Becciani, Alberto Vecchiato, Mario Gilberto Lattanzi, Fabio Pitari, Mario Raciti, Giuseppe Tudisco, Marco Aldinucci, Beatrice Bucciarelli

18 pages, 8 figures, 3 pseudocodes, published in Astronomy and Computing, Volume 41, October 2022, 100660, accepted for publication on 4th October 2022

The Gaia Astrometric Verification Unit-Global Sphere Reconstruction (AVU-GSR) Parallel Solver aims to find the astrometric parameters for $\sim$10$^8$ stars in the Milky Way, the attitude and the instrumental specifications of the Gaia satellite, and the global parameter $\gamma$ of the post Newtonian formalism. The code iteratively solves a system of linear equations, $\mathbf{A} \times \vec{x} = \vec{b}$, where the coefficient matrix $\mathbf{A}$ is large ($\sim$$10^{11} \times 10^8$ elements) and sparse. To solve this system of equations, the code exploits a hybrid implementation of the iterative PC-LSQR algorithm, where the computation related to different horizontal portions of the coefficient matrix is assigned to separate MPI processes. In the original code, each matrix portion is further parallelized over the OpenMP threads. To further improve the code performance, we ported the application to the GPU, replacing the OpenMP parallelization language with OpenACC. In this port, $\sim$95% of the data is copied from the host to the device at the beginning of the entire cycle of iterations, making the code $compute$ $bound$ rather than $data$$-$$transfer$ $bound$. The OpenACC code presents a speedup of $\sim$1.5 over the OpenMP version but further optimizations are in progress to obtain higher gains. The code runs on multiple GPUs and it was tested on the CINECA supercomputer Marconi100, in anticipation of a port to the pre-exascale system Leonardo, that will be installed at CINECA in 2022.

Cosmic rays are relativistic particles that come to the Earth from outer space. Despite a great effort made in both experimental and theoretical research, their origin is still unknown. One of the main keys to understand their nature is the determination of its chemical composition as a function of primary energy. In this paper, we review the measurements of the mass composition above $10^{15}$ eV. We first summarize the main aspects of air shower physics that are relevant in composition analyses. We discuss the composition measurements made by using optical, radio, and surface detectors and the limitations imposed by current high-energy hadronic interaction models that are used to interpret the experimental data. We also review the photons and neutrinos searches conducted in different experiments, which, in addition to being important to understand the nature of cosmic rays, can provide relevant information related to the abundance of heavy or light elements in the flux at the highest energies. Finally, we summarize the future composition measurements that are currently being planned or under development.

Andrew Ridden-Harper, Stevanus Nugroho, Laura Flagg, Ray Jayawardhana, Jake D. Turner, Ernst de Mooij, Ryan MacDonald, Emily Deibert, Motohide Tamura, Takayuki Kotani, Teruyuki Hirano, Masayuki Kuzuhara, Masashi Omiya, Nobuhiko Kusakabe

Submitted to AJ

Terrestrial exoplanets orbiting M-dwarf stars are promising targets for transmission spectroscopy with existing or near-future instrumentation. The atmospheric composition of such rocky planets remains an open question, especially given the high X-ray and ultraviolet flux from their host M-dwarfs that can drive atmospheric escape. The 1.3R$_\oplus$ exoplanet GJ 486b ($T_{\rm{eq}} \sim$ 700 K), orbiting an M3.5 star, is expected to have one of the strongest transmission spectroscopy signals among known terrestrial exoplanets. We observed three transits of GJ 486b using three different high-resolution spectrographs: IRD on Subaru, IGRINS on Gemini South, and SPIRou on the Canada-France-Hawai'i Telescope. We searched for atmospheric absorption from a wide variety of molecular species via the cross-correlation method, but did not detect any robust atmospheric signals. Nevertheless, our observations are sufficiently sensitive to rule out several clear atmospheric scenarios via injection and recovery tests, and extend comparative exoplanetology into the terrestrial regime. Our results suggest that GJ 486b does not possess a clear H$_2$/He dominated atmosphere, nor a clear 100% water-vapor atmosphere. Other high mean-molecular-weight secondary atmospheres or H$_2$/He dominated atmospheres with clouds remain possible. Our findings provide further evidence suggesting that terrestrial planets orbiting M-dwarf stars may experience significant atmospheric loss.

Using $\lambda$-21-cm galactic neutral atomic hydrogen data from the HI4PI survey of Bekhti et al. (2016) and 0.75-30 MeV $\gamma$-ray emission from the Imaging Compton Telescope, we have searched for the origin event that accelerated high velocity cloud Complex M. Radio plots of $l-b$, $l-v$, and $b-v$ show a cavity centered at ($l$, $b$) $\sim$ (150$^{\circ}$, 50.$^{\circ}$) and extending about $\pm$33$^{\circ}$. The best view of the cavity is at a velocity of -25 km s$^{-1}$, which shows a circular cross section on the back (receding) face. Complex M, at -85 km s$^{-1}$, is on the front (approaching) face. The $\gamma$-ray emission reveals several minima, the largest centered at ($l$, $b$) $\sim$ (150$^{\circ}$, 50.$^{\circ}$) and coincident with the position and extent of the cavity seen in the radio data. Using the know distance to Complex M and assuming that the cavity is spherical, we can bootstrap the distance to the original, explosive source of the cavity, D = 307 pc, calculate the radius of the cavity, R = 166 pc, and approximate the expansion velocity, V$_E$ $\approx$ 40 km $s^{-1}$, of the cavity. The total energy of the expanding cavity is 3.0 $\pm$ 1.0 ${\times}$ 10$^{50}$ ergs, well within the range of a single supernova. These results indicate that this explosion took place about four million years ago. As the blast wave from this supernova propagated outwards, it began to sweep up interstellar gas and carved out the Local Chimney, a low-density extension of the Local Bubble that reaches into the galactic halo.

The stability of eight nominal fictitious Uranus Trojan orbits over the age of the Solar system has been measured. The initial inclinations, i0, were 0 deg., 5 deg., 15 deg., and 30 deg. relative to the ecliptic plane. Initial eccentricities ranged from 0 to 0.1 for i0 = 0 deg., 5 deg., and 0 to 0.2 for i0 = 15 deg., 30 deg. Half of the orbits were in the L4 swarm, and half were in the L5 swarm. Orbits in the L4 swarm had mean longitudes 8.8 deg. from the nominal L4 Lagrange point, and orbits in the L5 swarm had mean longitudes 18.2 deg. from the nominal L5 point. I integrated 10,000 massless clones per nominal orbit in the six-body problem (Sun, test particle, four giant planets) for 4.5 Gyr and calculated the half-life for each orbit. A total of 1291 test particles survived for the entire integration time. Of these survivors, 99% were associated with the nominal orbit with i0 = 0 deg. in the L4 swarm. These surviving test particles had initial eccentricities in the range e0 < 0.07. The half-lives associated with L4 orbits were 1258 Myr, 286 Myr, 56 Myr, and 237 Myr for nominal orbits with i0 = 0 deg., 5 deg., 15 deg., and 30 deg., respectively. The half-lives associated with L5 orbits were 103 Myr, 281 Myr, 25 Myr, and 46 Myr, respectively. The overall results showed that the ecliptic plane is one good place to search for primordial Uranus Trojans.

A. Albert (1 and 2), S. Alves (3), M. André (4), M. Ardid (5), S. Ardid (5), J.-J. Aubert (6), J. Aublin (7), B. Baret (7), S. Basa (8), Y. Becherini (7), B. Belhorma (9), M. Bendahman (7 and 10), F. Benfenati (11 and 12), V. Bertin (6), S. Biagi (13), M. Bissinger (14), J. Boumaaza (10), M. Bouta (15), M.C. Bouwhuis (16), H. Brânzaş (17), R. Bruijn (16 and 18), J. Brunner (6), J. Busto (6), B. Caiffi (19), D. Calvo (3), S. Campion (20 and 21), A. Capone (20 and 21), L. Caramete (17), F. Carenini (11 and 12), J. Carr (6), V. Carretero (3), S. Celli (20 and 21), L. Cerisy (6), M. Chabab (22), T. N. Chau (7), R. Cherkaoui El Moursli (10), T. Chiarusi (11), M. Circella (23), J.A.B. Coelho (7), A. Coleiro (7), R. Coniglione (13), P. Coyle (6), A. Creusot (7), A. F. Díaz (24), et al. (110 additional authors not shown)

6 pages, 5 figures

Interactions of cosmic ray protons, atomic nuclei, and electrons in the interstellar medium in the inner part of the Milky Way produce a $\gamma$-ray flux from the Galactic Ridge. If the $\gamma$-ray emission is dominated by proton and nuclei interactions, a neutrino flux comparable to the $\gamma$-ray flux is expected from the same sky region. Data collected by the ANTARES neutrino telescope are used to constrain the neutrino flux from the Galactic Ridge in the 1-100 TeV energy range. Neutrino events reconstructed both as tracks and showers are considered in the analysis and the selection is optimized for the search of an excess in the region $|l| < 30\deg$, $|b| < 2\deg$. The expected background in the search region is estimated using an off region with similar sky coverage. Neutrino signal originating from a power-law spectrum with slope ranging from $\Gamma_\nu=1$ to $4$ is simulated in both channels. The observed energy distributions are fitted to constrain the neutrino emission from the Ridge. The energy distributions in the signal region are inconsistent with the background expectation at $\sim 96\%$ confidence level. The mild excess over the background is consistent with a neutrino flux with a power law with a slope $2.45^{+0.22}_{-0.34}$ and a flux normalization $dN_\nu/dE_\nu = 4.0^{+2.7}_{-2.0} \times 10^{-16} \text{GeV}^{-1} \text{cm}^{-2} \text{s}^{-1} \text{sr}^{-1}$ at 40 TeV reference energy. Such flux is consistent with the expected neutrino signal if the bulk of the observed $\gamma$-ray flux from the Galactic Ridge originates from interactions of cosmic ray protons and nuclei with a power-law spectrum extending well into the PeV energy range.

A. Simaz Bunzel, F. García, J. A. Combi, S. Chaty

8 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

The formation of neutron stars is associated with powerful astrophysical transients such as supernovae. In many cases, asymmetries in the supernova explosions are thought to be responsible for the large observed velocities of neutron stars. We aim to study the complete evolutionary history of one particular eccentric high-mass X-ray binary containing a neutron star, GX 301-2, and characterize the natal kick at the time of neutron star formation. We used the publicly available stellar-evolution code MESA to evolve binaries from their initial stages until the core-collapse scenario. We incorporated a natal kick distribution based on observations to continue the evolution during the X-ray binary phase and search for candidates matching current observations of GX 301-2. We find that the range of initial masses is constrained to be less than around $30$ M$_\odot$ depending on the initial mass ratio, as higher initial masses will most likely end up producing a black hole. In the completely conservative mass-transfer scenario under study, only is an interaction between the stars when the donor is still burning Hydrogen in its core, the so-called Case A of mass transfer, able to produce progenitors for GX 301-2. The natal kick study favours kicks of variable strength, which in turn increases the tilt angle between the orbital angular momentum and the spin of the neutron star. We conclude that only a narrow initial progenitor parameter space is able to produce a binary such as GX 301-2. Additionally, the strength of the natal kick can span a wide range of values, but it can be constrained when considering new data concerning the systemic velocity of the binary. Finally, we derive the fraction of the expected number of binaries such as GX 301-2 in the Galaxy to be $\sim 6 \times 10^{-5}$, implying a really low chance of finding a binary similar to GX 301-2.

G. Del Zanna, J. Samra, A. Monaghan, C. Madsen, P. Bryans, E. DeLuca, H. Mason, B. Berkey, A. de Wijn, Y. Rivera

Accepted for publication

The Airborne Infrared Spectrometer (AIR-Spec) offers an unprecedented opportunity to explore the Near Infra-Red (NIR) wavelength range. It has been flown at two total solar eclipses, in 2017 and 2019. The wavelength range of the much improved instrument on the second flight (July 2, 2019) was shifted to cover two density sensitive lines from S XI. In this paper we study detailed diagnostics for temperature, electron density and elemental abundances by comparing results from AIR-Spec slit positions above the east and the west limb with those from Hinode/EIS, the PolarCam detector and SDO/AIA. We find very good agreement in the electron densities obtained from the EIS EUV line ratios, those from the NIR S XI ratio and those obtained from the polarized brightness PolarCam measurements. Electron densities ranged from Log Ne [cm$^{-3}]$ = 8.4 near the limb, falling to 7.2 at $R_0=1.3$. EIS spectra indicate that the temperature distribution above the west limb is near-isothermal at around 1.3 MK, while that on the east has an additional higher-T component. The AIR-Spec radiances in Si X and S XI as well as the AIA data in the 171, 193, and 211 Angstroms bands are consistent with the EIS results. EIS and AIR-Spec data indicate that the sulphur abundance (relative to silicon) is photospheric in both regions, confirming our previous results of the 2017 eclipse. The AIA data also indicate that the absolute iron abundance is photospheric. Our analysis confirms the importance of the diagnostic potential of the NIR wavelength range, and that this important wavelength range can be used reliably and independently to determine coronal plasma parameters.

Valentin J. M. Le Gouellec, Anaëlle J. Maury, Charles L. H. Hull

High angular resolution observations of Class 0 protostars have produced detailed maps of the polarized dust emission in the envelopes of these young embedded objects. Interestingly, the improved sensitivity brought by ALMA has revealed wide dynamic ranges of polarization fractions, with specific locations harboring surprisingly large amounts of polarized dust emission. Our aim is to characterize the grain alignment conditions and dust properties responsible for the observed polarized dust emission in the inner envelopes (~1000 au) of Class 0 protostars. We analyzed the polarized dust emission maps obtained with ALMA and compared them to molecular line emission maps of specific molecular tracers, mainly CCH, which allowed us to probe one of the key components in dust grain alignment theories: the irradiation field. We show that CCH peaks toward outflow cavity walls, where the polarized dust emission is also enhanced. Our analysis provides a tentative correlation between the morphology of the polarized intensity and CCH emission, suggesting that the radiation field impinging on the cavity walls favors both the grain alignment and the warm carbon chain chemistry in these regions. We propose that shocks happening along outflow cavity walls could potentially represent an additional source of photons contributing to dust grain alignment. However, some parts of the cores, such as the equatorial planes, exhibit enhanced polarized flux, although no radiation driven chemistry is observed, for example where radiative torques are theoretically not efficient enough. This suggests that additional physical conditions, such as source geometry and dust grain evolution, may play a role in grain alignment.

Zdenek Sekanina

16 pages, 3 figures, 11 tables

Discovery of nine populations in a set of 193 select SOHO Kreutz sungrazers (Sekanina 2021) is confirmed for the first time via a histogram of the true longitudes of the ascending node, constructed for a revised set of 220 select sungrazers imaged exclusively by the SOHO's C2 coronagraph. Marsden's orbits are approximately corrected for effects of the out-of-plane nongravitational force. Population I displays two peaks in the histogram, one presumably belonging to a side branch alike to Population Pe, but with no related naked-eye sungrazer known. Swarms/clusters of objects are commonplace, providing evidence on cascading fragmentation proceeding throughout the orbit. Augmentation to all C2-only SOHO Kreutz comets, aimed at removing deliberate bias against Populations I and Pe, reduces the appearance of Populations Ia and Pre-I to bulges along the slope of the histogram because of the swollen wings of Populations I and Pe, respectively. Populations II through IV change very little or not at all. The high Population I-to-II abundance ratio, of 14:1, may be a product of temporal limitations in fragment release. A drop in the number of fragments toward the ends of the nodal-longitude distribution, especially from Population II to IV, is in line with the contact-binary model.

Eleonora Di Valentino, Stefano Gariazzo, William Giarè, Alessandro Melchiorri, Olga Mena, Fabrizio Renzi

14 pages, 3 figures, 9 tables

We present strong model-marginalized limits on mixed hot dark matter scenarios, which consider both thermal neutrinos and thermal QCD axions. A novel aspect of our analyses is the inclusion of small-scale Cosmic Microwave Background (CMB) observations from the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT), together with those from the Planck satellite and Baryon Acoustic Oscillation (BAO) data. After marginalizing over a number of well-motivated non-minimal background cosmologies, the tightest $95\%$ CL upper bound we obtain is $0.21$ eV, both for $\sum m_\nu$ and $m_{\rm a}$, from the combination of ACT, Planck and BAO measurements. Restricting the analyses to the standard $\Lambda$CDM picture, we find $\sum m_\nu<0.16$ eV and $m_{\rm a}<0.18$ eV, both at $95\%$ CL. Interestingly, the best background cosmology is never found within the minimal $\Lambda$CDM plus hot relics, regardless of the data sets exploited in the analyses. The combination of Planck with either BAO, SPT or ACT prefers a universe with a non-zero value of the running in the primordial power spectrum with strong evidence. Small-scale CMB probes, both alone and combined with BAO, either prefer, with substantial evidence, non-flat universes (as in the case of SPT) or a model with a time varying dark energy component (as in the case of ACT).

Matteo Biagetti, Juan Calles, Lina Castiblanco, Katherine González, Jorge Noreña

19 pages, 5 figures

We present a model for the squeezed dark matter bispectrum, where the short modes are deep in the non-linear regime. We exploit the consistency relations for large-scale structures combined with a response function approach to write the squeezed bispectrum in terms of a few unknown functions of the short modes. We provide an ansatz for a fitting function for these response functions, checking that the resulting model is reliable when compared to the one-loop squeezed bispectrum. We then test the model against measured bispectra from numerical simulations for short modes ranging between $k \sim 0.1 \, h/$Mpc, and $k \sim 0.7 \, h/$Mpc at redshift $z=0$. To evaluate the goodness of the fit of our model we implement a non-Gaussian covariance and find agreement within $1$-$\sigma$ standard deviation of the simulated data.

Daniele Michilli, Mohit Bhardwaj, Charanjot Brar, Chitrang Patel, B.M. Gaensler, Victoria M. Kaspi, Aida Kirichenko, Kiyoshi W. Masui, Ketan R. Sand, Paul Scholz, Kaitlyn Shin, Ingrid Stairs, Tomas Cassanelli, Amanda M. Cook, Matt Dobbs, Fengqiu Adam Dong, Emmanuel Fonseca, Adaeze Ibik, Jane Kaczmarek, Calvin Leung, Aaron B. Pearlman, Emily Petroff, Ziggy Pleunis, Masoud Rafiei-Ravandi, Pranav Sanghavi, Shriharsh P. Tendulkar

We report on improved sky localizations of thirteen repeating fast radio bursts (FRBs) discovered by CHIME/FRB via the use of interferometric techniques on channelized voltages from the telescope. These so-called 'baseband localizations' improve the localization uncertainty area presented in past studies by more than three orders of magnitude. The improved localization regions are provided for the full sample of FRBs to enable follow-up studies. The localization uncertainties, together with limits on the source distances from their dispersion measures (DMs), allow us to identify likely host galaxies for two of the FRB sources. FRB 20180814A lives in a massive passive red spiral at z~0.068 with very little indication of star formation, while FRB 20190303A resides in a merging pair of spiral galaxies at z~0.064 undergoing significant star formation. These galaxies show very different characteristics, further confirming the presence of FRB progenitors in a variety of environments even among the repeating sub-class.

Giovanni Tambalo, Miguel Zumalacárregui, Liang Dai, Mark Ho-Yeuk Cheung

41 pages, 27 figures

Just like light, gravitational waves (GWs) are deflected and magnified by gravitational fields as they propagate through the Universe. However, their low frequency, phase coherence and feeble coupling to matter allow for distinct lensing phenomena, such as diffraction and central images, that are challenging to observe through electromagnetic sources. Here we explore how these phenomena can be used to probe features of gravitational lenses. We focus on two variants of the singular isothermal sphere, with 1) a variable slope of the matter density and 2) a central core. We describe the imprints of these features in the wave- and geometric-optics regimes, including the prospect of detecting central images. We forecast the capacity of LISA and advanced LIGO to study strongly lensed signals and measure the projected lens mass, impact parameter and slope or core size. A broad range of lens masses allows all parameters to be measured with precision up to $\sim 1/{\rm SNR}$, despite large degeneracies. Thanks to wave-optics corrections, all parameters can be measured, even when no central image forms. Although GWs are sensitive to projected quantities, we compute the probability distribution of lens redshift, virial mass and projection scale given a cosmology. As an application, we consider the prospect of constraining self-interacting and ultra-light dark matter, showing the regions of parameter space accessible to strongly-lensed GWs. The distinct GW signatures will enable novel probes of fundamental physics and astrophysics, including the properties of dark matter and the central regions of galactic halos.

B. Holl, C. Fabricius, J. Portell, L. Lindegren, P. Panuzzo, M. Bernet, J. Castañeda, G. Jevardat de Fombelle, M. Audard, C. Ducourant, D.L. Harrison, D.W. Evans, G. Busso, A. Sozzetti, E. Gosset, F. Arenou, F. De Angeli, M. Riello, L. Eyer, L. Rimoldini, P. Gavras, N. Mowlavi, K. Nienartowicz, I. Lecoeur-Taïbi, P. García-Lario, D. Pourbaix

60 Figures, 2 Tables, submitted to A&A (v1: first partial revision incorporating various, but not yet all comments from the referee)

Context: Gaia DR3 time series data may contain spurious signals related to the time-dependent scan angle. Aims: We aim to explain the origin of scan-angle dependent signals and how they can lead to spurious periods, provide statistics to identify them in the data, and suggest how to deal with them in Gaia DR3 data and in future releases. Methods: Using real Gaia data, alongside numerical and analytical models, we visualise and explain the features observed in the data. Results: We demonstrated with Gaia data that source structure (multiplicity or extendedness) or pollution from close-by bright objects can cause biases in the image parameter determination from which photometric, astrometric and (indirectly) radial velocity time series are derived. These biases are a function of the time-dependent scan direction of the instrument and thus can introduce scan-angle dependent signals, which in turn can result in specific spurious periodic signals. Numerical simulations qualitatively reproduce the general structure observed in the spurious period and spatial distribution of photometry and astrometry. A variety of statistics allows for identification of affected sources. Conclusions: The origin of the scan-angle dependent signals and subsequent spurious periods is well-understood and is in majority caused by fixed-orientation optical pairs with separation <0.5" (amongst which binaries with P>>5y) and (cores of) distant galaxies. Though the majority of sources with affected derived parameters have been filtered out from the Gaia archive, there remain Gaia DR3 data that should be treated with care (e.g. gaia_source was untouched). Finally, the various statistics discussed in the paper can not only be used to identify and filter affected sources, but alternatively reveal new information about them not available through other means, especially in terms of binarity on sub-arcsecond scale.

We present differential double-copy relations between gluon and graviton three-point functions in (A)dS$_{d+1}$. We introduce a set of differential operators in (A)dS that naturally generalize on-shell kinematics of scattering amplitudes in flat space. This provides a way to construct (A)dS correlators by replacing the kinematic variables of amplitudes with the corresponding differential operators and suitably ordering them. By construction, the resulting correlators are manifestly conformally invariant, with the correct flat-space limit, and exhibit a differential double-copy structure.

Rouzbeh Allahverdi, Jacek K. Osiński, Ngo Phuc Duc Loc

12 pages, 4 figures

We present a minimal extension of the standard model that includes a long-lived fermion with weak-scale mass and an ${\cal O}({\rm GeV})$ fermionic dark matter candidate both of which are coupled to quarks. Decays of a TeV-scale colored scalar in a radiation-dominated phase bring the former to a thermal abundance while also producing dark matter. The long-lived fermion then dominates the energy density of the Universe and drives a period of early matter domination. It decays to reheat the Universe, mainly through baryon-number-violating interactions that also generate a baryon asymmetry, with a small branching fraction to dark matter. We find the allowed parameter space of the model and show that it can be probed by proposed long-lived particle searches as well as next-generation neutron-antineutron oscillation experiments. This model provides a robust explanation of dark matter and baryogenesis as long as the Universe is in a radiation-dominated phase at $T \gtrsim {\cal O}({\rm TeV})$.

We propose a novel and simple scenario to explain baryon asymmetry and dark matter (DM) by utilizing an early matter-dominated era (EMDE) caused by a heavy metastable particle. Within the EMDE, lack of pressure enhances the formation of primordial black holes (PBHs) which can then contribute to the relic abundance of DM. The eventual decay of heavy metastable particle that has baryon number and $CP$ violating interactions reheats the Universe and gives rise to baryon asymmetry. Since in this setup, PBH serves as a DM candidate, the particle physics model may not require new stable degrees of freedom which leads to more freedom in the model-building side. As an example, we show that a modulus field which dominates the energy density of the Universe prior to its decay, may explain both DM and baryon asymmetry in the Universe in the context of the Minimal Supersymmetric Standard Model (MSSM) while the lightest superparticle is not stable and cannot be a DM candidate due to the R-parity violating interactions needed for baryogenesis.

Sudhanshu Pandey, Maarten van Nistelrooij, Joannes D. Maasakkers, Pratik Sutar, Sander Houweling, Daniel J. Varon, Paul Tol, David Gains, John Worden, Ilse Aben

The twin Sentinel-3 satellites have multi-band radiometers which observe at methane-sensitive shortwave infrared bands with daily global coverage and 500 m ground pixel resolution. We investigate the methane observation capability of Sentinel-3 and how its coverage-resolution combination fits between Sentinel-5p and Sentinel-2. We show that methane plume enhancements can be retrieved from the shortwave infrared bands of Sentinel-3. We report a lowest emission detection by Sentinel-3 of 9 t/h under favorable detection conditions of low wind speeds and high surface albedo. We demonstrate Sentinel-3-based identification and monitoring of methane leaks using two case studies. Near Moscow, Sentinel-3 shows that two major short-term leaks separated by 30 km occurred simultaneously at a gas pipeline and appear as a single methane plume in Sentinel-5p data. For a major Sentinel-5p leak detection near the Hassi Messaoud oil/gas field in Algeria, Sentinel-3 identifies the leaking facility emitting continuously for 6 days, and Sentinel-2 pinpoints the source of the leak at an oil/gas well. Sentinel-2 and Sentinel-3 also show the 6-day leak was followed by a four-month period of burning of the leaking gas, suggesting a gas well blowout to be the cause of the leak. We find similar source rate quantifications from plume detections by Sentinel-3 and Sentinel-2 for these leaks, demonstrating utility of Sentinel-3 for emission quantification. These case studies show that zooming in with Sentinel-3 and Sentinel-2 in synergy allows precise identification and quantification as well as monitoring of the sources corresponding to methane anomalies observed in global scans of Sentinel-5p.

We study the non-linear stability of fixed-point solutions to the $\alpha'$-exact equations from O$(d,d)$ invariant cosmology, with and without matter perturbations. Previous non-linear analysis in the literature is revisited, and its compatibility with known linear perturbation results is shown. Some formal aspects of cosmological perturbations in duality invariant cosmology are discussed, and we show the existence of time-reparameterization invariant variables for perturbations.

The current design of space-based gravitational wave detectors utilizes heterodyne laser interferometry in inter-satellite science measurements. Frequency variations of the heterodyne beatnotes are predominantly caused by the Doppler effect from relative satellite motion along lines of sight. Generally considered to be outside the measurement band, the Doppler effect appears to have been largely overlooked in literature on numerical simulations of time-delay interferometry (TDI). However, the potential impact on the effectiveness of TDI should be assessed. The issue is particularly relevant to TianQin that features geocentric orbits, because of strong gravity disturbances from the Earth-Moon system at $<1\times 10^{-4}$ Hz. In this paper, based on high-precision orbital data obtained from detailed gravity field modeling, we incorporate the Doppler shift in the generation of TianQin's beatnote phase signals. To remove the large-scale Doppler phase drift at $<1\times 10^{-4}$ Hz, we develop a high-performance high-pass filter and consider two possible processing sequences, i.e., applying the filter before or after TDI combinations. Our simulation results favor the former and demonstrate successful removal of the low-frequency gravity disturbances for TianQin without degrading the TDI performance, assuming 10 m pseudo-ranging uncertainty. The filtering scheme can be used in developing the initial noise-reduction pipeline for TianQin.

TianQin is a proposed space-based gravitational-wave detector mission to be deployed and operated in high Earth orbits. As a sequel to [Zhang et al. Phys. Rev. D 103, 062001 (2021)], we investigate a type of ``orbital noise'' in TianQin's range acceleration that is caused by gravitational perturbation associated with solar free oscillations. Frequencies of such oscillations are typically within TianQin's measurement band of 0.1 mHz--1 Hz, and the disturbance level needs careful assessment. By using high-precision orbit propagation and adding the Sun's time-variable oblateness $J_2$ to detailed gravity-field models, we examine the effect in the frequency domain and show that the solar free oscillation noise is expected to be two orders of magnitude lower than the noise requirement on single links and hence has little impact on the mission.

Natalia Buzulukova (1 and 2), Bruce Tsurutani (3) ((1) NASA GSFC, Heliophysics Division, Geospace Physics Laboratory, Greenbelt, MD, USA, (2) University of Maryland, Department of Astronomy, College Park, MD, USA, (3) Heliospheric Physics and Astrophysics Section, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA)

Review article, 4 figures

Space Weather is the portion of space physics that has a direct effect on humankind. Space Weather is an old branch of space physics that originates back to 1808 with the publication of a paper by the great naturalist Alexander von Humboldt (von Humboldt, 1808). Space Weather is currently experiencing explosive growth, because its effects on human technologies have become more and more diverse. Space Weather is due to the variability of solar processes that cause interplanetary, magnetospheric, ionospheric, atmospheric and ground level effects. Space Weather can at times have strong impacts on technological systems and human health. The threats and risks are not hypothetical, and in the event of extreme Space Weather events the consequences could be quite severe for humankind. The purpose of the review is to give a brief overall view of the full chain of physical processes responsible for Space Weather risks and hazards, tracing them from solar origins to effects and impacts in interplanetary space, in the Earth's magnetosphere and ionosphere and at the ground. The paper shows that the risks associated with Space Weather have not been constant over time; they have evolved as our society becomes more and more technologically advanced. The paper begins with a brief introduction to the Carrington event. Next, the descriptions of the strongest known Space Weather processes are reviewed. The concepts of geomagnetic storms and substorms are briefly introduced. The main effects/impacts of Space Weather are also considered, including geomagnetically induced currents (GICs) which are thought to cause power outages. The effects of radiation on avionics and human health, ionospheric effects and impacts, and thermosphere effects and satellite drag will also be discussed. Finally, we will discuss the current challenges of Space Weather forecasting and examine some of the worst-case scenarios.

Yuichiro Nakai, Ryo Namba, Ippei Obata

27 pages, 6 figures

We explore a mechanism to produce a light dark photon dark matter through a coupling between the dark photon field and a spectator scalar field which plays no role in the inflationary expansion of the Universe while rolling down its potential during the inflation. The motion of the spectator field efficiently produces dark photons with large wavelengths which become non-relativistic before the time of matter-radiation equality. The spectrum of the wavelengths is peaky so that the constraint from the isocurvature perturbation can be evaded. The correct relic abundance is then achieved over a wide range of the dark photon mass down to $10^{-13} \ \text{eV}$. Our mechanism favors high-scale inflation models which can be tested in future observations. Furthermore, fluctuations of the dark photon field during inflation produce gravitational waves detectable at future space-based interferometers and/or pulsar timing array experiments.

A general formalism to find the density profile of a slowly rotating stellar object in modified gravity is presented. We derive a generic Lane-Emden equation and its analytical solution for a wide class of modified theories of gravity.

Luca Mancini, Marco Trinari, Emília Valença Ferreira de Aragão, Marzio Rosi, Nadia Balucani

15 pages, 7 figures. Preprint version submitted to LNCS (Springer) ICCSA2022

We have performed a theoretical investigation of the S$^+$($^4$S) + SiH$_{2}$($^1$A$_1$) reaction, a possible formation route of the HSiS$^+$ and SiSH$^+$ cations that are alleged to be precursors of interstellar silicon sulfide, SiS. Electronic structure calculations allowed us to identify the main reaction pathways for the systems. The reaction has two exothermic channels leading to the isomeric species $^3$HSiS$^{+}$ and $^3$SiSH$^{+}$ formed in conjunction with H atoms. The reaction is not characterized by an entrance barrier and, therefore, it is believed to be fast also under the very low temperature conditions of interstellar clouds. The two ions are formed in their first electronically excited state because of the spin multiplicity of the overall potential energy surface. In addition, following the suggestion that neutral species are formed by proton transfer of protonated cations to ammonia, we have derived the potential energy surface for the reactions $^3$HSiS$^{+}$/$^3$SiSH$^{+}$ + NH$_{3}$($^{1}$A$_1$).

M. R. Pelicer, M. Antonelli, D. P. Menezes, F. Gulminelli

15 pages, 12 figures

The presence of nuclear pasta is expected to modify the transport properties in the mantle of neutron stars. The non-spherical geometry of the pasta nuclear clusters leads to anisotropies in the collision frequencies, impacting the thermal and electrical conductivity. We derive analytical expressions for the anisotropic collision frequencies using the Boltzmann equation in the relaxation time approximation. The average parallel, perpendicular and Hall electrical conductivities are computed in the high-temperature regime above crustal melting, considering incoherent elastic electron-pasta scattering and randomly oriented pasta structures. Numerical values are obtained at different densities and temperatures by using the IUFSU parametrization of the non-linear Walecka model to determine the crustal structure. We find that the anisotropy of the collision frequencies grows with the length of the pasta structures and, independently of the magnetic field, the presence of rod and slab phases decreases the conductivity by more than one order of magnitude. Our numerical results indicate that, even if the pasta structures might survive above the crustal melting point, no strong anisotropies are to be expected in the conduction properties in this temperature regime, even in the presence of a very high magnetic field.

Christian Dioguardi, Antonio Racioppi, Eemeli Tomberg

20 pages, 6 figures

We study single-field slow-roll inflation embedded in a Palatini quadratic $F(R)$ gravity, where the Einstein-Hilbert term has the wrong sign, apparently leading to repulsive gravity. This can be avoided as long as $F'(R)$ and $F''(R)$ stay positive. Surprisingly, consistency of the theory requires the Jordan frame inflaton potential to be unbounded from below. Even more surprisingly, this corresponds to an Einstein frame inflaton potential bounded from below and positive definite. We prove that such a quadratic gravity is an attractor configuration for all the Palatini $F(R)$ that, for infinite curvature, diverge faster than $R^2$.

Yongsoo Jho, Tae-Geun Kim, Jong-Chul Park, Seong Chan Park, Yeji Park

8 pages, 10 figures

Primordial black holes (PBHs) can be significant sources of axions and axion-like particles (ALPs) in the Universe as the Hawking radiation of the PBH includes light particles when the Hawking temperature exceeds the particle's mass. Once produced, as axions predominantly decay into photons, we may detect the enhanced photon spectrum using sensitive detectors. We introduce a new methodology by defining the time-varying decay process for particles to fly and decay over time on a cosmological scale. This paper provides the estimated photon spectrum and the flux under some simplified assumptions about PBH, 1) monochromatic mass spectrum and 2) isotropic distribution at cosmological scales. Future detectors, such as e-ASTROGAM, have great chances of detecting the signal.

Kaustubh Agashe, Jae Hyeok Chang, Steven J. Clark, Bhaskar Dutta, Yuhsin Tsai, Tao Xu

8 pages + references, 5 figures

Future gamma-ray experiments, such as the e-ASTROGAM and AMEGO telescopes, can detect the Hawking radiation of photons from primordial black holes (PBHs) if they make up a fraction or all of dark matter. PBHs can analogously also Hawking radiate new particles, which is especially interesting if these particles are mostly secluded from the Standard Model (SM) sector, since they might therefore be less accessible otherwise. A well-motivated example of this type is axion-like particles (ALPs) with a tiny coupling to photons. We assume that the ALPs produced by PBHs decay into photons well before reaching the earth, so these will augment the photons directly radiated by the PBHs. Remarkably, we find that the peaks in the energy distributions of ALPs produced from PBHs are different than the corresponding ones for Hawking radiated photons due to the spin-dependent greybody factor. Therefore, we demonstrate that this process will in fact distinctively modify the PBHs' gamma-ray spectrum relative to the SM prediction. We use monochromatic asteroid-mass PBHs as an example to show that e-ASTROGAM can observe the PBH-produced ALP gamma-ray signal (for masses up to ~60 MeV) and further distinguish it from Hawking radiation without ALPs. By measuring the gamma-ray signals, e-ASTROGAM can thereby probe yet unexplored parameters in the ALP mass and photon coupling.