Abstract visibility
Abstract text size

Papers for Friday, Jul 22 2022

Papers with local authors

Philipp Kempski, Eliot Quataert, Jonathan Squire

Submitted to MNRAS

0 votes
Paper 6 — arXiv:2207.10107
0 votes
Paper 6 — arXiv:2207.10107

Active Galactic Nuclei (AGN) are believed to provide the energy that prevents runaway cooling of gas in the cores of galaxy clusters. However, how this energy is transported and thermalized throughout the Intracluster Medium (ICM) remains unclear. In recent work we showed that streaming cosmic rays (CRs) destabilise sound waves in dilute ICM plasmas. Here we show that CR streaming in the presence of gravity also destabilises a pressure-balanced wave. We term this new instability the CR buoyancy instability (CRBI). In stark contrast to standard results without CRs, the pressure-balanced mode is highly compressible at short wavelengths due to CR streaming. Maximal growth rates are of order $(p_c / p_g) \beta^{1/2} \omega_{\rm ff}$, where $p_c/p_g$ is the ratio of CR pressure to thermal gas pressure, $\beta$ is the ratio of thermal to magnetic pressure and $\omega_{\rm ff}$ is the free-fall frequency. The CRBI operates alongside buoyancy instabilities driven by background heat fluxes, i.e. the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI). When the thermal mean free path $l_{\rm mfp}$ is $\ll$ the gas scale height $H$, the HBI/MTI set the growth rate on large scales, while the CRBI sets the growth rate on small scales. Conversely, when $l_{\rm mfp} \sim H$ and $(p_c/p_g) \beta^{1/2} \gtrsim 1$, CRBI growth rates exceed HBI/MTI growth rates even on large scales. Our results suggest that CR-driven instabilities may be partially responsible for the sound waves/weak shocks and turbulence observed in galaxy clusters. CR-driven instabilities generated near radio bubbles may also play an important role redistributing AGN energy throughout clusters.

Ruchi Pandey, Gesesew R. Habtie, Rahul Bandyopadhyay, Ramkrishna Das, François Teyssier, Joan Guarro Flo

16 Pages, 7 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS) on 19 July 2022

0 votes
Paper 26 — arXiv:2207.10473
0 votes
Paper 26 — arXiv:2207.10473

We present the evolution of the optical spectra of the 2021 outburst of RS Ophiuchi (RS Oph) over about a month after the outburst. The spectral evolution is similar to the previous outbursts. Early spectra show prominent P Cygni profiles of hydrogen Balmer, \ion{Fe}{ii}, and \ion{He}{i} lines. The emission lines were very broad during the initial days, which later became narrower and sharper as the nova evolved. This is interpreted as the expanding shocked material into the winds of the red giant companion. We find that the nova ejecta expanded freely for $\sim 4$ days, and afterward, the shock velocity decreased monotonically with time as $v\propto t^{-0.6}$. The physical and chemical parameters associated with the system are derived using the photoionization code \textsc{cloudy}. The best-fit \textsc{cloudy} model shows the presence of a hot central white dwarf source with a roughly constant luminosity of $\sim$1.00 $\times$ 10$^{37}$ erg s$^{-1}$. The best-fit photoionization models yield absolute abundance values by number, relative to solar of He/H $\sim 1.4 - 1.9$, N/H = $70 - 95$, O/H = $0.60 - 2.60$, and Fe/H $\sim 1.0 - 1.9$ for the ejecta during the first month after the outburst. Nitrogen is found to be heavily overabundant in the ejecta. The ejected hydrogen shell mass of the system is estimated to be in the range of $3.54 - 3.83 \times 10^{-6} M_{\odot}$. The 3D morpho-kinematic modelling shows a bipolar morphology and an inclination angle of $i=30^{\circ}$ for the RS Oph binary system.

All other papers

Bin Liu, Daniel J. D'Orazio, Alejandro Vigna-Gómez, Johan Samsing

10 pages, 8 figures, 1 supplemental material. To be submitted

We study the dynamics of a solar-type star orbiting around a black hole binary (BHB) in a nearly coplanar system. We present a novel effect that can prompt a growth and significant oscillations of the eccentricity of the stellar orbit when the system encounters an "apsidal precession resonance", where the apsidal precession rate of the outer stellar orbit matches that of the inner BHB. The eccentricity excitation requires the inner binary to have a non-zero eccentricity and unequal masses, and can be created even in non-coplanar triples. We show that the secular variability of the stellar orbit's apocenter induced by the changing eccentricity could be potentially detectable with \textit{Gaia} if the inner BHB is a LISA source, thereby providing a distinctive probe on the existence of BHBs.

F. Mernier, N. Werner, J. Bagchi, M.-L. Gendron-Marsolais, Gopal-Krishna, M. Guainazzi, A. Richard-Laferrière, T. W. Shimwell, A. Simionescu

5 pages, 2 figures, submitted to MNRAS Letters. Comments are welcome

Observed in a significant fraction of clusters and groups of galaxies, diffuse radio synchrotron emission reveals the presence of relativistic electrons and magnetic fields permeating large-scale systems. Although these non-thermal electrons are expected to upscatter cosmic microwave background photons up to hard X-ray energies, such inverse-Compton (IC) X-ray emission has so far not been unambiguously detected on cluster/group scales. Here we report the first robust ($4.6\,\sigma$) detection of extended IC X-ray emission from MRC$\,$0116+111, a group of galaxies. This unambiguous detection provides the most direct and least model-dependent estimate of the volume-averaged magnetic field within a galaxy group. Such estimates can serve as a fulcrum to theories of magnetic field generation within the largest gravitationally bound systems in the Universe.

J.E. Rodríguez Martín, R. M. González Delgado, G. Martínez-Solaeche, L. A. Díaz-García, A. de Amorim, R. García-Benito, E. Pérez, R. Cid Fernandes, E. R. Carrasco, M. Maturi, A. Finoguenov, P. A. A. Lopes, A. Cortesi, G. Lucatelli, J. M. Diego, A. L. Chies-Santos, R. A. Dupke, Y. Jiménez-Teja, J. M. Vílchez, L. R. Abramo, J. Alcaniz, N. Benítez, S. Bonoli, A. J. Cenarro, D. Cristóbal-Hornillos, A. Ederoclite, A. Hernán-Caballero, C. López-Sanjuan, A. Marín-Franch, C. Mendes de Oliveira, M. Moles, L. Sodré Jr., K. Taylor, J. Varela, H. Vázquez Ramió, I. Márquez

24 pages, 17 figures, 4 tables, accepted in Astronomy & Astrophysics

The miniJPAS is a 1 deg$^2$ survey that uses the Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) filter system (54 narrow-band filters) with the Pathfinder camera. We study mJPC2470-1771, the most massive cluster detected in miniJPAS. We study the stellar population properties of the members, their star formation rates (SFR), star formation histories (SFH), the emission line galaxy (ELG) population, their spatial distribution, and the effect of the environment on them, showing the power of J-PAS to study the role of environment in galaxy evolution. We use a spectral energy distribution (SED) fitting code to derive the stellar population properties of the galaxy members: stellar mass, extinction, metallicity, colours, ages, SFH (a delayed-$\tau$ model), and SFRs. Artificial Neural Networks are used for the identification of the ELG population through the detection of H$\alpha$, [NII], H$\beta$, and [OIII] nebular emission. We use the WHAN and BPT diagrams to separate them into star-forming galaxies and AGNs. We find that the fraction of red galaxies increases with the cluster-centric radius. We select 49 ELG, 65.3\% of the them are probably star forming galaxies, and they are dominated by blue galaxies. 24% are likely to host an AGN (Seyfert or LINER galaxies). The rest are difficult to classify and are most likely composite galaxies. Our results are compatible with an scenario where galaxy members were formed roughly at the same epoch, but blue galaxies have had more recent star formation episodes, and they are quenching from inside-out of the cluster centre. The spatial distribution of red galaxies and their properties suggest that they were quenched prior to the cluster accretion or an earlier cluster accretion epoch. AGN feedback and/or mass might also be intervening in the quenching of these galaxies.

J. W. Petley (1), L. K. Morabito (1 and 2), D. M. Alexander (1), A. L. Rankine (3), V. A. Fawcett (1), D. J. Rosario (4), J. H. Matthews (5), T. M. Shimwell (6 and 7), A. Drabent (8) ((1) Centre for Extragalactic Astronomy, Durham University, (2) Institute for Computational Cosmology, Durham University, (3) Institute for Astronomy, University of Edinburgh, (4) School of Mathematics, Statistics and Physics, Newcastle University, (5) Institute of Astronomy, University of Cambridge, (6) ASTRON, (7) Leiden Observatory, Leiden University, (8) Thuringer Landessternwarte)

16 pages, 14 figures

Broad Absorption Line Quasars (BALQSOs) show strong signatures of powerful outflows, with the potential to alter the cosmic history of their host galaxies. These signatures are only seen in ~10% of optically selected quasars, although the fraction significantly increases in IR and radio selected samples. A proven physical explanation for this observed fraction has yet to be found, along with a determination of why this fraction increases at radio wavelengths. We present the largest sample of radio matched BALQSOs using the LOFAR Two-metre Sky Survey Data Release 2 and employ it to investigate radio properties of BALQSOs. Within the DR2 footprint, there are 3537 BALQSOs from Sloan Digital Sky Survey DR12 with continuum signal to noise >5. We find radio-detections for 1108 BALQSOs, with an important sub-population of 120 LoBALs, an unprecedented sample size for radio matched BALQSOs given the LoTSS sky coverage to date. BALQSOs are a radio-quiet population that show an increase of $\times 1.50$ radio-detection fraction compared to non-BALQSOs. LoBALs show an increase of $\times 2.22$ that of non-BALQSO quasars. We show that this detection fraction correlates with wind-strength, reddening and C_{IV} emission properties of BALQSOs and that these features may be connected, although no single property can fully explain the enhanced radio detection fraction. We create composite spectra for sub-classes of BALQSOs based on wind strength and colour, finding differences in the absorption profiles of radio-detected and radio-undetected sources, particularly for LoBALs. Overall, we favour a wind-ISM interaction explanation for the increased radio-detection fraction of BALQSOs.

Zoe Ansari, Christa Gall, Roger Wesson, Oswin Krause

24 pages,19 figures, 7 tables, submitted to A&A 10/01/2022

Context. Determining properties of dust formed in and around supernovae from observations remains challenging. This may be due to either incomplete coverage of data in wavelength or time but also due to often inconspicuous signatures of dust in the observed data. Aims. Here we address this challenge using modern machine learning methods to determine the amount, composition and temperature of dust from a large set of simulated data. We aim to determine whether such methods are suitable to infer these properties from future observations of supernovae. Methods. We calculate spectral energy distributions (SEDs) of dusty shells around supernovae. We develop a neural network consisting of eight fully connected layers and an output layer with specified activation functions that allow us to predict the dust mass, temperature and composition and their respective uncertainties from each SED. We conduct a feature importance analysis via SHapley Additive exPlanations (SHAP) to find the minimum set of JWST filters required to accurately predict these properties. Results. We find that our neural network predicts dust masses and temperatures with a root-mean-square error (RMSE) of $\sim$ 0.12 dex and $\sim$ 38 K, respectively. Moreover, our neural network can well distinguish between the different dust species included in our work, reaching a classification accuracy of up to 95\% for carbon and 99\% for silicate dust. Conclusions. Our analysis shows that the JWST filters NIRCam F070W, F140M, F356W, F480M and MIRI F560W, F770W, F1000W, F1130W, F1500W, F1800W are likely the most important needed to determine the properties of dust formed in and around supernovae from future observations. We tested this on selected optical to infrared data of SN 1987A at 615 days past explosion and find good agreement with dust masses and temperatures inferred with standard fitting methods in the literature.

S. Schuldt, S. H. Suyu, R. Canameras, Y. Shu, S. Taubenberger, S. Ertl, A. Halkola

17+28 pages, 7+31 figures, 2+5 tables, submitted to A&A

Modeling of strongly gravitationally lensed galaxies is often required in order to use them as astrophysical or cosmological probes. With current and upcoming wide-field imaging surveys, the number of detected lenses is increasing significantly such that automated and fast modeling procedures for ground-based data are urgently needed. This is especially pertinent to short-lived lensed transients in order to plan follow-up observations. Therefore, we present in a companion paper (submitted) a neural network predicting the parameter values with corresponding uncertainties of a Singular Isothermal Ellipsoid (SIE) mass profile with external shear. In this work, we present a newly-developed pipeline glee_auto.py to model consistently any galaxy-scale lensing system. In contrast to previous automated modeling pipelines that require high-resolution images, glee_auto.py is optimized for ground-based images such as those from the Hyper-Suprime-Cam (HSC) or the upcoming Rubin Observatory Legacy Survey of Space and Time. We further present glee_tools.py, a flexible automation code for individual modeling that has no direct decisions and assumptions implemented. Both pipelines, in addition to our modeling network, minimize the user input time drastically and thus are important for future modeling efforts. We apply the network to 31 real galaxy-scale lenses of HSC and compare the results to the traditional models. In the direct comparison, we find a very good match for the Einstein radius especially for systems with $\theta_E \gtrsim 2$". The lens mass center and ellipticity show reasonable agreement. The main discrepancies are on the external shear as expected from our tests on mock systems. In general, our study demonstrates that neural networks are a viable and ultra fast approach for measuring the lens-galaxy masses from ground-based data in the upcoming era with $\sim10^5$ lenses expected.

L. Doyle, H. M. Cegla, E. Bryant, D. Bayliss, M. Lafarga, D. R. Anderson, R. Allart, V. Bourrier, M. Brogi, N. Buchschacher, V. Kunovac, M. Lendl, C. Lovis, M. Moyano, N. Roguet-Kern, J. V. Seidel, D. Sosnowska, P. J. Wheatley, J. S. Acton, M. R. Burleigh, S. L. Casewell, S. Gill, M. R. Goad, B.A. Henderson, J. S. Jenkins, R. H. Tilbrook, R. G. West

18 Pages, 13 Figures, 4 Tables. Accepted for publication in MNRAS

In this paper, we present high-resolution spectroscopic transit observations from ESPRESSO of the super-Neptune WASP-166~b. In addition to spectroscopic ESPRESSO data, we analyse photometric data from {\sl TESS} of six WASP-166~b transits along with simultaneous NGTS observations of the ESPRESSO runs. These observations were used to fit for the planetary parameters as well as assessing the level of stellar activity (e.g. spot crossings, flares) present during the ESPRESSO observations. We utilise the Reloaded Rossiter McLaughlin (RRM) technique to spatially resolve the stellar surface, characterising the centre-to-limb convection-induced variations, and to refine the star-planet obliquity. We find WASP-166~b has a projected obliquity of $\lambda = -15.52^{+2.85}_{-2.76}$$^{\circ}$ and $v\sin(i) = 4.97 \pm 0.09$~kms$^{-1}$ which is consistent with the literature. We were able to characterise centre-to-limb convective variations as a result of granulation on the surface of the star on the order of a few kms$^{-1}$ for the first time. We modelled the centre-to-limb convective variations using a linear, quadratic and cubic model with the cubic being preferred. In addition, by modelling the differential rotation and centre-to-limb convective variations simultaneously we were able to retrieve a potential anti-solar differential rotational shear ($\alpha \sim$ -0.5) and stellar inclination ($i_*$ either 42.03$^{+9.13}_{-9.60}$$^{\circ}$ or 133.64$^{+8.42}_{-7.98}$$^{\circ}$ if the star is pointing towards or away from us). Finally, we investigate how the shape of the cross-correlation functions change as a function of limb angle and compare our results to magnetohydrodynamic simulations.

Steven Ehlert, Chien-Ting Chen, Doug A. Swartz, Ryan C. Hickox, Alexander A. Lutovinov, Andrey N. Semena, Roman Krivonos, Andrey E. Shtykovsky, Alexey Tkachenko

15 pages, 9 figures, Accepted for Publication in MNRAS

We present a new statistical method for constructing background subtracted measurements from event list data gathered by X-ray and gamma ray observatories. This method was initially developed specifically to construct images that account for the high background fraction and low overall count rates observed in survey data from the Mikhail Pavlinsky ART-XC telescope aboard the Spektrum R\"{o}ntgen Gamma (SRG) mission, although the mathematical underpinnings are valid for data taken with other imaging missions and analysis applications. This method fully accounts for the expected Poisson fluctuations in both the sky photon and non X-ray background count rates in a manner that does not result in unphysical negative counts. We derive the formulae for arbitrary confidence intervals for the source counts and show that our new measurement converges exactly to the standard background subtraction calculation in the high signal limit. Utilizing these results, we discuss several variants of images designed to optimize different science goals for both pointed and slewing telescopes. Using realistic simulated data of a galaxy cluster as observed by ART-XC we show that our method provides a more significant and robust detection of the cluster emission as compared to a standard background subtraction. We also demonstrate its advantages using real observations of a point source from the ART-XC telescope. These calculations may have widespread applications for a number of source classes observed with high energy telescopes.

S. Agayeva, V. Aivazyan, S. Alishov, M. Almualla, C. Andrade, S. Antier, J.-M. Bai, A. Baransky, S. Basa, P. Bendjoya, Z. Benkhaldoun, S. Beradze, D. Berezin, U. Bhardwaj, M.-A. Bizouard, M. Blazek, O. Burkhonov, E. Burns, S. Caudill, N. Christensen, F. Colas, A. Coleiro, W. Corradi, M. W. Coughlin, T. Culino, D. Darson, D. Datashvili, T. Dietrich, F. Dolon, D. Dornic, J. Dubouil, J.-G. Ducoin, P.-A. Duverne, A. Esamdin, A. Fouad, F. Guo, V. Godunova, P. Gokuldass, N. Guessoum, E. Gurbanov, R. Hainich, E. Hasanov, P. Hello, T. Hussenot-Desenonges, R. Inasaridze, A. Iskandar, E. E. O. Ishida, N. Ismailov, T. Jegou du Laz, D. A. Kann, G. Kapanadze, S. Karpov, R. W. Kiendrebeogo, A. Klotz, N. Kochiashvili, A. Kaeouach, J.-P. Kneib, W. Kou, K. Kruiswijk, S. Lombardo, M. Lamoureux, N. Leroy, A. Le Van Su, et al. (39 additional authors not shown)

Submitted to the Proceedings of the SPIE, Astronomical Telescopes and Instrumentation 2022

GRANDMA is a world-wide collaboration with the primary scientific goal of studying gravitational-wave sources, discovering their electromagnetic counterparts and characterizing their emission. GRANDMA involves astronomers, astrophysicists, gravitational-wave physicists, and theorists. GRANDMA is now a truly global network of telescopes, with (so far) 30 telescopes in both hemispheres. It incorporates a citizen science programme (Kilonova-Catcher) which constitutes an opportunity to spread the interest in time-domain astronomy. The telescope network is an heterogeneous set of already-existing observing facilities that operate coordinated as a single observatory. Within the network there are wide-field imagers that can observe large areas of the sky to search for optical counterparts, narrow-field instruments that do targeted searches within a predefined list of host-galaxy candidates, and larger telescopes that are devoted to characterization and follow-up of the identified counterparts. Here we present an overview of GRANDMA after the third observing run of the LIGO/VIRGO gravitational-wave observatories in $2019-2020$ and its ongoing preparation for the forthcoming fourth observational campaign (O4). Additionally, we review the potential of GRANDMA for the discovery and follow-up of other types of astronomical transients.

Arsalan Adil, Andreas Albrecht, Lloyd Knox

18 pages, 8 figures, 1 appendix

Several cosmological tensions have emerged in light of recent data, most notably in the inferences of the parameters $H_0$ and $\sigma_8$. We explore the possibility of alleviating both these tensions {\it simultaneously} by means of the Albrecht-Skordis ``quintessence'' potential. The field can reduce the size of the sound horizon $r_s^*$ while concurrently suppressing the power in matter density fluctuations before it comes to dominate the energy density budget today. Interestingly, this rich set of dynamics is governed entirely by one free parameter that is of $\mathcal{O}(10)$ in Planck units. We find that the inferred value of $H_0$ can be increased, while that of $\sigma_8$ can be decreased, both by $\approx 1\sigma$ compared to the $\Lambda$CDM case. However, ultimately the model is disfavored by Planck and BAO data alone, compared to the standard $\Lambda$CDM model, with a $\Delta \chi^2 \approx +6$. When including large scale structure and supernova data $\Delta \chi^2 \approx +1$. We note that historically much attention has been focused on preserving the three angular scales $\theta_D$, $\theta_{EQ}$, and $\theta_s^*$ to their $\Lambda$CDM values. Our work presents an example of how, while doing so indeed maintains a relatively good fit to the CMB data for an increased number of ultra-relativistic species, it is a-priori insufficient in maintaining such a fit in more general model spaces.

Longkang Dai, Xi Zhang, Wencheng D. Shao, Carver J. Bierson, Jun Cui

The current Venus climate is largely regulated by globally-covered concentrated sulfuric acid clouds from binary condensation of sulfuric acid (H2SO4) and water (H2O). To understand this complicated H2SO4-H2O gas-cloud system, previous theoretical studies either adopted complicated microphysical calculations or assumed that both H2SO4 and H2O vapor follow their saturation vapor pressure. In this study, we developed a simple one-dimensional cloud condensation model including condensation, diffusion and sedimentation of H2SO4 and H2O but without detailed microphysics. Our model is able to explain the observed vertical structure of cloud and upper haze mass loading, cloud acidity, H2SO4, and H2O vapor, and the mode-2 particle size on Venus. We found that most H2SO4 is stored in the condensed phase above 48 km, while the partitioning of H2O between the vapor and clouds is complicated. The cloud cycle is mostly driven by evaporation and condensation of H2SO4 rather than H2O and is about seven times stronger than the H2SO4 photochemical cycle. Most of the condensed H2O in the upper clouds is evaporated before the falling particles reach the middle clouds. The cloud acidity is affected by the temperature and the condensation-evaporation cycles of both H2SO4 and H2O. Because of the large chemical production of H2SO4 vapor and relatively inefficient cloud condensation, the simulated H2SO4 vapor above 60 km is largely supersaturated by more than two orders of magnitude, which could be tested by future observations.

The Venusian clouds originate from the binary condensation of H$_{2}$SO$_{4}$ and H$_{2}$O. The two components strongly interact with each other via chemistry and cloud formation. Previous works adopted sophisticated microphysical approaches to understand the clouds. Here we show that the observed vapor and cloud distributions on Venus can be well explained by a semi-analytical model. Our model assumes local thermodynamical equilibrium for water vapor but not for sulfuric acid vapor, and includes the feedback of cloud condensation and acidity to vapor distributions. The model predicts strong supersaturation of the H$_{2}$SO$_{4}$ vapor above 60 km, consistent with our recent cloud condensation model. The semi-analytical model is 100 times faster than the condensation model and 1000 times faster than the microphysical models. This allows us to quickly explore a large parameter space of the sulfuric acid gas-cloud system. We found that the cloud mass loading in the upper clouds has an opposite response of that in the lower clouds to the vapor mixing ratios in the lower atmosphere. The transport of water vapor influences the cloud acidity in all cloud layers while the transport of sulfuric acid vapor only dominates in the lower clouds. This cloud model is fast enough to be coupled with the climate models and chemistry models to understand the cloudy atmospheres of Venus and Venus-like extra-solar planets.

C. B. Oliveira Jr., A.C. Krabbe, J. A. Hernandez-Jimenez, O. L. Dors Jr., I. A. Zinchenko, G. F. Hägele, M. V. Cardaci, A. F. Monteiro

16 pages, 14 figures, 3 tables, accepted for publication in MNRAS

The ionizing source of Low Ionization Nuclear Emission Regions (LINERs) is uncertain. Because of this, an empirical relation to determine the chemical abundances of these objects has not been proposed. In this work, for the first time, we derived two semi-empirical calibrations based on photoionization models to estimate the oxygen abundance of LINERS as a function of the $N2$ and $O3N2$ emission-line intensity ratios. These relations were calibrated using oxygen abundance estimations obtained by comparing the observational emission-line ratios of 43 LINER galaxies (taken from the MaNGA survey) and grids of photoionization models built with the {\sc Cloudy} code assuming post-Asymptotic Giant Branch (post-AGB) stars with different temperatures. We found that the oxygen abundance of LINERs in our sample is in the $\rm 8.48 <~ 12+log(O/H) <~ 8.84$ range, with a mean value of $\rm 12+\log(O/H)=8.65$. We recommend the use of the $N2$ index to estimate the oxygen abundances of LINERs, since the calibration with this index presented a much smaller dispersion than the $O3N2$ index. In addition, the estimated metallicities are in good agreement with those derived by extrapolating the disk oxygen abundance gradients to the centre of the galaxies showing that the assumptions of the models are suitable for LINERs. We also obtained a calibration between the logarithm of the ionization parameter and the [OIII]/[OII] emission-line ratio.

Ryosuke T. Tominaga, Hiroshi Kobayashi, Shu-ichiro Inutsuka

19 pages, 12 figures, accepted for publication in ApJ

Our previous linear analysis presents a new instability driven by dust coagulation in protoplanetary disks. The coagulation instability has the potential to concentrate dust grains into rings and assist dust coagulation and planetesimal formation. In this series of papers, we perform numerical simulations and investigate nonlinear outcome of coagulation instability. In this paper (Paper I), we first conduct local simulations to demonstrate the existence of coagulation instability. Linear growth observed in the simulations is in good agreement with the previous linear analysis. We next conduct radially global simulations to demonstrate that coagulation instability develops during the inside-out disk evolution due to dust growth. To isolate the various effects on dust concentration and growth, we neglect effects of backreaction to a gas disk and dust fragmentation in Paper I. This simplified simulation shows that either of backreaction or fragmentation is not prerequisite for local dust concentration via the instability. In most runs with weak turbulence, dust concentration via coagulation instability overcomes dust depletion due to radial drift, leading to the formation of multiple dust rings. The nonlinear development of coagulation instability also accelerates dust growth, and the dimensionless stopping time $\tau_{\mathrm{s}}$ reaches unity even at outer radii (>10 au). Therefore, coagulation instability is one promising process to retain dust grains and to accelerate dust growth beyond the drift barrier.

Gao-Lei Hou, Olga V. Lushchikova, Joost M. Bakker, Peter Lievens, Leen Decin, Ewald Janssens

45 pages (23 pages for main text and 22 pages for supplementary information) 21 figures (10 figures for main text and 11 figures for supplementary information) 6 tables (1 tables for main text and 5 tables for supplementary information)

Infrared emission bands with wavelengths between 3-20 {\mu}m are observed in a variety of astrophysical environments [1,2]. They were discovered in the 1970s and are generally attributed to organic compounds [3,4]. However, over 40 years of research efforts still leave the source of these emission bands largely unidentified [5-7]. Here, we report the first laboratory infrared (6-25 {\mu}m) spectra of gas-phase fullerene-metal complexes, [C60-Metal]+ (Metal = Fe and V), and show with density functional theory calculations that complexes of C60 with cosmically abundant metals, including Li, Na, K, Mg, Ca, Al, V, and Fe, all have similar infrared spectral patterns. Comparison with observational infrared spectra from several fullerene-rich planetary nebulae demonstrates a strong positive linear cross-correlation. The infrared features of [C60-Metal]+ coincide with four bands attributed earlier to neutral C60 bands, and in addition also with several to date unexplained bands. Abundance and collision theory estimates furthermore indicate that [C60-Metal]+ could plausibly form and survive in astrophysical environments. Hence, [C60-Metal]+ are proposed as promising carriers, in supplement to C60, of astronomical infrared emission bands, potentially representing the largest molecular species in space other than the bare fullerenes C60, C60+, and C70. This work opens a new chapter for studying cosmic fullerene species and carbon chemistry in the Universe.

K. Kreckel, O. Egorov, F. Belfiore, B. Groves, S. C. O. Glover, R. S. Klessen, K. Sandstrom, F. Bigiel, D. A. Dale, K. Grasha, F. Scheuermann, T. G. Williams

21 pages, 22 figures, accepted by A&A

Aims: Temperature uncertainties plague our understanding of abundance variations within the ISM. Using the PHANGS-MUSE large program, we develop and apply a new technique to model the strong emission lines arising from HII regions in 19 nearby spiral galaxies at ~50 pc resolution and infer electron temperatures for the nebulae. Methods: Due to the charge-exchange coupling of the ionization fraction of the atomic oxygen to that of hydrogen, the emissivity of the observed [OI]6300/Ha line ratio can be modeled as a function of gas phase oxygen abundance (O/H), ionization fraction (f_ion) and electron temperature (T_e). We measure (O/H) using a strong line metallicity calibration, and identify a correlation between f_ion and [SIII]9069/[SII]6716,6730, tracing ionization parameter variations. Results: We solve for T_e, and test the method by reproducing direct measurements of T_e([NII]5755) based on auroral line detections to within ~600 K. We apply this charge-exchange method of calculating T_e to 4,129 HII regions across 19 PHANGS-MUSE galaxies. We uncover radial temperature gradients, increased homogeneity on small scales, and azimuthal temperature variations in the disks that correspond to established abundance patterns. This new technique for measuring electron temperatures leverages the growing availability of optical integral field unit spectroscopic maps across galaxy samples, increasing the statistics available compared to direct auroral line detections.

Shun'ichi Horigome, Kohei Hayashi, Shin'ichiro Ando

31 pages, 9 figures

Dark matter halos of dwarf spheroidal galaxies (dSphs) play important roles in dark matter detection. Generally we estimate the halo profile using a kinematical equation of dSphs but the halo profile has a large uncertainty because we have only a limited number of kinematical dataset. In this paper, we utilize cosmological models of dark matter subhalos to obtain better constraints on halo profile of dSphs. The constraints are realized as two cosmological priors: satellite prior, based on a semi-analytic model of the accretion history of subhalos and their tidal stripping effect, and stellar-to-halo mass relation prior, which estimates halo mass of a galaxy from its stellar mass using empirical correlations. In addition, we adopt a radial dependent likelihood function by considering velocity dispersion profile, which allows us to mitigate the parameter degeneracy in the previous analysis using a radial independent likelihood function with averaged dispersion. Using these priors, we estimate the squared dark matter density integrated over the region-of-interest (so-called $J$-factor) of 8 classical and 27 ultra-faint dSphs. Our method significantly decreases the uncertainty of $J$-factors (upto about $20\%$) compared to the previous radial independent analysis. We confirm the model dependence of $J$-factor estimates by evaluating Bayes factors of different model setups and find that the estimates are still stable even when assuming different cosmological models.

S.R. Bhavanam, Sumohana S. Channappayya, P.K. Srijith, Shantanu Desai

21 pages, 23 figures. Accepted in Astronomy and Computing

Cosmic Ray (CR) hits are the major contaminants in astronomical imaging and spectroscopic observations involving solid-state detectors. Correctly identifying and masking them is a crucial part of the image processing pipeline, since it may otherwise lead to spurious detections. For this purpose, we have developed and tested a novel Deep Learning based framework for the automatic detection of CR hits from astronomical imaging data from two different imagers: Dark Energy Camera (DECam) and Las Cumbres Observatory Global Telescope (LCOGT). We considered two baseline models namely deepCR and Cosmic-CoNN, which are the current state-of-the-art learning based algorithms that were trained using Hubble Space Telescope (HST) ACS/WFC and LCOGT Network images respectively. We have experimented with the idea of augmenting the baseline models using Attention Gates (AGs) to improve the CR detection performance. We have trained our models on DECam data and demonstrate a consistent marginal improvement by adding AGs in True Positive Rate (TPR) at 0.01% False Positive Rate (FPR) and Precision at 95% TPR over the aforementioned baseline models for the DECam dataset. We demonstrate that the proposed AG augmented models provide significant gain in TPR at 0.01% FPR when tested on previously unseen LCO test data having images from three distinct telescope classes. Furthermore, we demonstrate that the proposed baseline models with and without attention augmentation outperform state-of-the-art models such as Astro-SCRAPPY, Maximask (that is trained natively on DECam data) and pre-trained ground-based Cosmic-CoNN. This study demonstrates that the AG module augmentation enables us to get a better deepCR and Cosmic-CoNN models and to improve their generalization capability on unseen data.

Alessandro Maraio, Alex Hall, Andy Taylor

16 pages, 19 figures, submitted to MNRAS

Headline constraints on cosmological parameters from current weak lensing surveys are derived from two-point statistics that are known to be statistically sub-optimal, even in the case of Gaussian fields. We study the performance of a new fast implementation of the Quadratic Maximum Likelihood (QML) estimator, optimal for Gaussian fields, to test the performance of Pseudo-Cl estimators for upcoming weak lensing surveys and quantify the gain from a more optimal method. Through the use of realistic survey geometries, noise levels, and power spectra, we find that there is a decrease in the errors in the statistics of the recovered E-mode spectra to the level of ~20% when using the optimal QML estimator over the Pseudo-Cl estimator on the largest angular scales, while we find significant decreases in the errors associated with the B-modes for the QML estimator. This raises the prospects of being able to constrain new physics through the enhanced sensitivity of B-modes for forthcoming surveys that our implementation of the QML estimator provides. We test the QML method with a new implementation that uses conjugate-gradient and finite-differences differentiation methods resulting in the most efficient implementation of the full-sky QML estimator yet, allowing us to process maps at resolutions that are prohibitively expensive using existing codes. In addition, we investigate the effects of apodisation, B-mode purification, and the use of non-Gaussian maps on the statistical properties of the estimators. Our QML implementation is publicly available and can be accessed from GitHub.

Zhixing Li (Tsinghua), Hong Guo (SHAO), Yi Mao (Tsinghua)

14 pages, 8 figures, 1 table. Submitted to ApJ. Comments welcome

Atomic hydrogen (H I) gas, mostly residing in dark matter halos after cosmic reionization, is the fuel for star formation. Its relation with properties of host halo is the key to understand the cosmic H I distribution. In this work, we propose a flexible, empirical model of H I-halo relation. In this model, while the H I mass depends primarily on the mass of host halo, there is also secondary dependence on other halo properties. We apply our model to the observation data of the Arecibo Fast Legacy ALFA Survey (ALFALFA), and find it can successfully fit to the cosmic H I abundance ($\Omega_{\rm HI}$), average H I-halo mass relation $\langle M_{\rm HI}|M_{\rm h}\rangle$, and the H I clustering. The bestfit of the ALFALFA data rejects with high confidence level the model with no secondary halo dependence of H I mass and the model with secondary dependence on halo spin parameter ($\lambda$), and shows strong dependence on halo formation time ($a_{1/2}$) and halo concentration ($c_{\rm vir}$). In attempt to explain these findings from the perspective of hydrodynamical simulations, the IllustrisTNG simulation confirms the dependence of H I mass on secondary halo parameters. However, the IllustrisTNG results show strong dependence on $\lambda$ and weak dependence on $c_{\rm vir}$ and $a_{1/2}$, and also predict a much larger value of H I clustering on large scales than observations. This discrepancy between the simulation and observation calls for improvements in understanding the H I-halo relation from both theoretical and observational sides.

Amit Singh Ubhi, Leonid Prokhorov, Sam Cooper, Chiara Di Fronzo, John Bryant, David Hoyland, Alexandra Mitchell, Jesse van Dongen, Conor Mow-Lowry, Alan Cumming, Giles Hammond, Denis Martynov

We demonstrate the control scheme of an active platform with a six degree of freedom (6D) seismometer. The inertial sensor simultaneously measures translational and tilt degrees of freedom of the platform and does not require any additional sensors for the stabilisation. We show that a feedforward cancellation scheme can efficiently decouple tilt-to-horizontal coupling of the seismometer in the digital control scheme. We stabilise the platform in the frequency band from 250 mHz up to 10 Hz in the horizontal degrees of freedom and achieve a suppression factor of 100 around 1 Hz. Further suppression of ground vibrations was limited by the non-linear response of the piezo actuators of the platform and by its limited range (5 {\mu}m). In this paper we discuss the 6D seismometer, its control scheme, and the limitations of the test bed.

Duncan V. Mifsud, Zuzana Kaňuchová, Sergio Ioppolo, Péter Herczku, Alejandra Traspas Muiña, Béla Sulik, K.K. Rahul, Sándor T.S. Kovács, Perry A. Hailey, Robert W. McCullough, Nigel J. Mason, Zoltán Juhász

Accepted for publication in Phys. Chem. Chem. Phys

The detection of ozone (O3) in the surface ices of Ganymede, Jupiters largest moon, and of the Saturnian moons Rhea and Dione, has motivated several studies on the route of formation of this species. Previous studies have successfully quantified trends in the production of O3 as a result of the irradiation of pure molecular ices using ultraviolet photons and charged particles (i.e., ions and electrons), such as the abundances of O3 formed after irradiation at different temperatures or using different charged particles. In this study, we extend such results by quantifying the abundance of O3 as a result of the 1 keV electron irradiation of a series of 14 stoichiometrically distinct CO2:O2 astrophysical ice analogues at 20 K. By using mid-infrared spectroscopy as our primary analytical tool, we have also been able to perform a spectral analysis of the asymmetric stretching mode of solid O3 and the variation in its observed shape and profile among the investigated ice mixtures. Our results are important in the context of better understanding the surface composition and chemistry of icy outer Solar System objects, and may thus be of use to future interplanetary space missions such as the ESA Jupiter Icy Moons Explorer and the NASA Europa Clipper missions, as well as the recently launched NASA James Webb Space Telescope.

Jin-Ping Zhu, Xiangyu Ivy Wang, Hui Sun, Yuan-Pei Yang, Zhuo Li, Rui-Chong Hu, Ying Qin, Shichao Wu

Submitted to ApJL on June 23, revised manuscript following the referee's comments

Here we collect three unique bursts, GRBs\,060614, 211211A and 211227A, all characterized by a long-duration main emission (ME) phase and a rebrightening extended emission (EE) phase, to study their observed properties and the potential origin as neutron star-black hole (NSBH) mergers. NS-first-born (BH-first-born) NSBH mergers tend to contain fast-spinning (non-spinning) BHs that more easily (hardly) allow tidal disruption to happen with (without) forming electromagnetic signals. We find that NS-first-born NSBH mergers can well interpret the origins of these three GRBs, supported by that: (1) Their X-ray MEs and EEs show unambiguous fall-back accretion signatures, decreasing as $\propto{t}^{-5/3}$, which might account for their long duration. The EEs can result from the fall-back accretion of $r$-process heating materials, predicted to occur after NSBH mergers. (2) The beaming-corrected local event rate density for this type of merger-origin long-duration GRBs is $\mathcal{R}_0\sim2.4^{+2.3}_{-1.3}\,{\rm{Gpc}}^{-3}\,{\rm{yr}}^{-1}$, consistent with that of NS-first-born NSBH mergers. (3) Our detailed analysis on the EE, afterglow and kilonova of the recently high-impact event GRB\,211211A reveals it could be a merger between a $\sim1.17^{+0.13}_{-0.04}\,M_\odot$ NS and a $\sim9.3^{+1.6}_{-1.6}\,M_\odot$ BH with an aligned-spin of $\chi_{\rm{BH}}\sim0.67^{+0.08}_{-0.10}$, supporting an NS-first-born NSBH formation channel. Long-duration burst with rebrightening fall-back accretion signature after ME, and bright kilonova might be commonly observed features for on-axis NSBHs. We estimate the multimessenger detection rate between gravitational waves, GRBs and kilonovae from NSBH mergers in O4 (O5) is $\sim0.1\,{\rm{yr}}^{-1}$ ($\sim1\,{\rm{yr}}^{-1}$).

H.E.S.S. Collaboration, H. Abdalla, F. Aharonian, F. Ait Benkhali, E.O. Anguner, C. Armand, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlohr, B. Bi, M. Bottcher, J. Bolmont, M. de Bony de Lavergne, R. Brose, F. Brun, F. Cangemi, S. Caroff, M. Cerruti, T. Chand, A. Chen, G. Cotter, J. Damascene Mbarubucyeye, J. Devin, A. Djannati-Ataı, A. Dmytriiev, V. Doroshenko, K. Egberts, A. Fiasson, G. Fichet de Clairfontaine, G. Fontaine, S. Funk, S. Gabici, G. Giavitto, D. Glawion, J.F. Glicenstein, M.-H. Grondin, J.A. Hinton, W. Hofmann, T. L. Holch, M. Holler, D. Horns, Zhiqiu Huang, M. Jamrozy, F. Jankowsky, E. Kasai, K. Katarzynski, U. Katz, B. Khelifi, W. Kluzniak, Nu. Komin, K. Kosack, D. Kostunin, G. Lamanna, M. Lemoine-Goumard, J.-P. Lenain, et al. (81 additional authors not shown)

Accepted in Phys. Rev. Lett., includes Supplemental materials. 5+10 pages, 2+8 figures, 4 tables

The central region of the Milky Way is one of the foremost locations to look for dark matter (DM) signatures. We report the first results on a search for DM particle annihilation signals using new observations from an unprecedented gamma-ray survey of the Galactic Center (GC) region, ${\it i.e.}$, the Inner Galaxy Survey, at very high energies ($\gtrsim$ 100 GeV) performed with the H.E.S.S. array of five ground-based Cherenkov telescopes. No significant gamma-ray excess is found in the search region of the 2014-2020 dataset and a profile likelihood ratio analysis is carried out to set exclusion limits on the annihilation cross section $\langle \sigma v\rangle$. Assuming Einasto and Navarro-Frenk-White (NFW) DM density profiles at the GC, these constraints are the strongest obtained so far in the TeV DM mass range. For the Einasto profile, the constraints reach $\langle \sigma v\rangle$ values of $\rm 3.7\times10^{-26} cm^3s^{-1}$ for 1.5 TeV DM mass in the $W^+W^-$ annihilation channel, and $\rm 1.2 \times 10^{-26} cm^3s^{-1}$ for 0.7 TeV DM mass in the $\tau^+\tau^-$ annihilation channel. With the H.E.S.S. Inner Galaxy Survey, ground-based $\gamma$-ray observations thus probe $\langle \sigma v\rangle$ values expected from thermal-relic annihilating TeV DM particles.

M. Moscibrodzka

9 pages, 4 figures, submitted to ApJ

We describe upgrades to a numerical code which computes synchrotron and inverse-Compton emission from relativistic plasma including full polarization. The introduced upgrades concern scattering kernel which is now capable of scattering the polarized and unpolarized photons on non-thermal population of electrons. We describe the scheme to approach this problem and we test the numerical code against known analytic solution. Finally, using the upgraded code, we predict polarization of light that is scattered off sub-relativistic thermal or relativistic thermal and non-thermal free electrons. The upgraded code enables more realistic simulations of emissions from plasma jets associated with accreting compact objects.

Sowmya G. M., Rajani G., U. Paniveni, R. Srikanth

15 pages; 5 figures; Accepted for publication in Research in Astronomy and Astrophysics

We present findings from an analysis of the fractal dimension of solar supergranulation as a function of latitude, supergranular cell size and solar rotation, employing spectroheliographic data in the Ca II K line of solar cycle no. 23. We find that the fractal dimension tends to decrease from about 1.37 at the equator to about 1 at 20 degree latitude in either hemisphere, suggesting that solar rotation rate has the effect of augmenting the irregularity of supergranular boundaries. Considering that supergranular cell size is directly correlated with fractal dimension, we conclude that the mechanism behind our observation is that solar rotation influences the cell outflow strength, and thereby cell size, with the latitude dependence of the supergranular fractal dimension being a consequence thereof.

Joseph F. V. Allingham, Mathilde Jauzac, David J. Lagattuta, Guillaume Mahler, Céline Bœhm, Geraint F. Lewis, Dominique Eckert, Alastair Edge, Stefano Ettori

20 pages, 14 figures, submitted to MNRAS

We present the strong lensing analysis of two galaxy clusters: MACS J0242.5-2132 (MACS J0242, $z=0.313$) and MACS J0949.8+1708 (MACS J0949, $z=0.383$). Their total matter distributions are constrained thanks to the powerful combination of observations with the Hubble Space Telescope and the MUSE instrument. Using these observations, we precisely measure the redshift of six multiple image systems in MACS J0242, and two in MACS J0949. We also include four multiple image systems in the latter cluster identified in HST imaging without MUSE redshift measurements. For each cluster, our best-fit mass model consists of a single cluster-scale halo, and 57 (170) galaxy-scale halos for MACS J0242 (MACS J0949). Multiple images positions are predicted with a $rms$ 0.39 arcsec and 0.15 arcsec for MACS J0242 and MACS J0949 models respectively. From these mass models, we derive aperture masses of $M(R<$200 kpc$) = 1.67_{-0.05}^{+0.03}\times 10^{14} M_{\odot}$, and $M(R<$200 kpc$) = 2.00_{-0.20}^{+0.05}\times 10^{14} M_{\odot}$. Combining our analysis with X-ray observations from the XMM-Newton Observatory, we show that MACS J0242 appears to be a relatively relaxed cluster, while conversely, MACS J0949 shows a relaxing post-merger state. At 200 kpc, X-ray observations suggest the baryon fraction to be respectively $f_b = 0.115^{+0.003}_{-0.004}$ and $0.053^{+0.007}_{-0.006}$ for MACS J0242 and MACS J0949. MACS J0242 being relaxed, its density profile is very well fit by a NFW distribution, in agreement with X-ray observations. Finally, the strong lensing analysis of MACS J0949 suggests a flat dark matter density distribution in the core, between 10 and 100 kpc, and not following a NFW profile. This appears consistent with X-ray observations.

Javier R. Goicoechea, Sara Cuadrado, Franck Le Petit

To appear in "The European Physical Journal Web of Conferences". Invited paper at "Multi-line Diagnostics of the Interstellar Medium", Nice 2022

Far-ultraviolet photons from OB-type massive stars regulate the heating, ionization, and chemistry of much of the neutral interstellar gas in star-forming galaxies. The interaction of FUV radiation and interstellar matter takes place in environments broadly known as photodissociation regions (PDRs). PDR line diagnostics are the smoking gun of the radiative feedback from massive stars. Improving our understanding of stellar feedback in the ISM requires quantifying the energy budget, gas dynamics, and chemical composition of PDR environments. This goal demands astronomical instrumentation able to deliver multi-line spectroscopic images of the ISM (of the Milky Way and nearby galaxies). It also requires interdisciplinary collaborations to obtain the rate coefficients and cross sections of the many microphysical processes that occur in the ISM and that are included in models such as the Meudon PDR code.

H. Nguyen, M. R. Rugel, C. Murugeshan, K. M. Menten, A. Brunthaler, J. S. Urquhart, R. Dokara, S. A. Dzib, Y. Gong, S. Khan, S-N. X. Medina, G. N. Ortiz-Leon, W. Reich, F. Wyrowski, A. Y. Yang, H. Beuther, W. D. Cotton, J. D. Pandian

Accepted in A&A July 18, 2022

Class II methanol (CH$_{3}$OH) masers are amongst the clearest signposts of recent high-mass star formation (HMSF). A complete catalogue outlines the distribution of star formation in the Galaxy, the number of young star-forming cores, and the physical conditions of their environment. The Global View on Star Formation (GLOSTAR) survey, which is a blind survey in the radio regime of 4$-$8 GHz, maps the Galactic mid-plane in the radio continuum, 6.7 GHz methanol line, the 4.8 GHz formaldehyde line, and several radio recombination lines. We present the analysis of the observations of the 6.7 GHz CH$_{3}$OH maser transition using data from the D-configuration of the Very Large Array (VLA). We analyse the data covering Galactic longitudes from $-2^{\circ}< l <60^{\circ}$ and Galactic latitudes of $|\textit{b}|<1^{\circ}$. We detect a total of 554 methanol masers, out of which 84 are new, and catalogue their positions, velocity components, and integrated fluxes. With a typical noise level of $\sim$18 mJy beam$^{-1}$, this is the most sensitive unbiased methanol survey for methanol masers to date. We search for dust continuum and radio continuum associations, and find that 97% of the sources are associated with dust, and 12% are associated with radio continuum emission.

Ozgur Akarsu, Eoin O. Colgain, Emre Ozulker, Somyadip Thakur, Lu Yin

14 pages, 3 figures

Using the fact that the comoving angular diameter distance to last scattering is strictly constrained almost model-independently, we show that, for any model agreeing with $\Lambda$CDM on its background dynamics at $z\sim0$ and size of the comoving sound horizon at last scattering, the deviations of the Hubble radius from the one of $\Lambda$CDM, should be a member of the set of admissible wavelets. The family of models characterized by this framework also offers non-trivial oscillatory behaviours in various functions that define the kinematics of the universe, even when the wavelets themselves are very simple. We discuss the consequences of attributing these kinematics to, first, dark energy, second, varying gravitational coupling strength. Utilizing some simplest wavelets, we demonstrate the power and flexibility of this framework in describing the BAO data without any modifications to the agreement with CMB measurements. This framework also provides a natural explanation for the bumps found in non-parametric observational reconstructions of the Hubble parameter and dark energy density as compensations of the dips required by the BAO data, and questions the physical reality of their existence. We note that utilizing this framework on top of the models that agree with both the CMB and local $H_0$ measurements but are held back by BAO data, one may resurrect these models through the wiggly nature of wavelets that can naturally accommodate the BAO data. Finally, we also suggest narrowing the plausible set of admissible wavelets to further improve our framework by imposing conditions from expected kinematics of a viable cosmological model or first principle fundamental physics such as energy conditions.

Astrid Anker, Steven W. Barwick, Dave Z. Besson, Sjoerd Bouma, Maddalena Cataldo, Geoffrey Gaswint, Christian Glaser, Steffen Hallmann, Jordan C. Hanson, Jakob Henrichs, Stuart A. Kleinfelder, Robert Lahmann, Zachary S. Meyers, Anna Nelles, Alexander Novikov, Manuel P. Paul, Lilly Pyras, Christopher Persichilli, Ilse Plaisier, Ryan Rice-Smith, Mohammad F.H. Seikh, Joulien Tatar, Christoph Welling, Leshan Zhao

41 pages, 23 figures, to be submitted to JCAP

TAROGE-M is a self-triggered radio antenna array atop the 2700 m high Mt. Melbourne in Antarctica, designed to detect impulsive geomagnetic emission from extensive air showers induced by ultra-high energy (UHE) particles beyond 0.1 EeV, including cosmic rays (CRs), Earth-skimming tau neutrinos, and particularly, the "ANITA anomalous events" (AAEs) from near and below the horizon, which origin remains uncertain and requires more experimental inputs for clarification. The detection concept of TAROGE-M takes advantage of a high altitude with synoptic view toward the horizon as an efficient signal collector, and the radio quietness as well as strong and near vertical geomagnetic field in Antarctica. This approach has a low energy threshold, high duty cycle, and is easy to extend for quickly enlarging statistics. Here we report experimental results from the first TAROGE-M station deployed in 2020, corresponding to $25.3$-days of livetime. The station consists of six receiving antennas operating at 180-450 MHz, and can reconstruct source directions with $\sim0.3^\circ$ angular resolution. To demonstrate its ability to detect UHE air showers, a search for CR signals in the data was conducted, resulting in seven identified events. These events have a mean reconstructed energy of $0.95_{-0.31}^{+0.46}$ EeV and zenith angles between $25^\circ-82^\circ$, with both distributions agreeing with simulations. The estimated CR flux is also consistent with results of other experiments. The TAROGE-M sensitivity to AAEs is approximated by the tau neutrino exposure with simulations, suggesting comparable sensitivity as ANITA's at $~1$ EeV energy with a few station-years of operation. These first results verified the station design and performance in a polar and high-altitude environment, and are promising for further discovery of tau neutrinos and AAEs after an extension in the near future.

Duoming Fan, Donald E. Willcox, Christopher DeGrendele, Michael Zingale, Andrew Nonaka

We demonstrate the use of neural networks to accelerate the reaction steps in the MAESTROeX stellar hydrodynamics code. A traditional MAESTROeX simulation uses a stiff ODE integrator for the reactions; here we employ a ResNet architecture and describe details relating to the architecture, training, and validation of our networks. Our customized approach includes options for the form of the loss functions, a demonstration that the use of parallel neural networks leads to increased accuracy, and a description of a perturbational approach in the training step that robustifies the model. We test our approach on millimeter-scale flames using a single-step, 3-isotope network describing the first stages of carbon fusion occurring in Type Ia supernovae. We train the neural networks using simulation data from a standard MAESTROeX simulation, and show that the resulting model can be effectively applied to different flame configurations. This work lays the groundwork for more complex networks, and iterative time-integration strategies that can leverage the efficiency of the neural networks.

Susmita Adhikari, Arka Banerjee, Kimberly K. Boddy, Francis-Yan Cyr-Racine, Harry Desmond, Cora Dvorkin, Bhuvnesh Jain, Felix Kahlhoefer, Manoj Kaplinghat, Anna Nierenberg, Annika H. G. Peter, Andrew Robertson, Jeremy Sakstein, Jesús Zavala

87 pages, 20 figures; to be submitted to Reviews of Modern Physics

Self-interacting dark matter (SIDM) arises generically in scenarios for physics beyond the Standard Model that have dark sectors with light mediators or strong dynamics. The self-interactions allow energy and momentum transport through halos, altering their structure and dynamics relative to those produced by collisionless dark matter. SIDM models provide a promising way to explain the diversity of galactic rotation curves, and they form a predictive and versatile framework for interpreting astrophysical phenomena related to dark matter. This review provides a comprehensive explanation of the physical effects of dark matter self-interactions in objects ranging from galactic satellites (dark and luminous) to clusters of galaxies and the large-scale structure. The second major part describes the methods used to constrain SIDM models including current constraints, with the aim of advancing tests with upcoming galaxy surveys. This part also provides a detailed review of the unresolved small-scale structure formation issues and concrete ways to test simple SIDM models. The review is rounded off by a discussion of the theoretical motivation for self-interactions, degeneracies with baryonic and gravitational effects, extensions to the single-component elastic-interactions SIDM framework, and future observational and theoretical prospects.

Katherine A. Suess, Rachel Bezanson, Erica J. Nelson, David J. Setton, Sedona H. Price, Pieter van Dokkum, Gabriel Brammer, Ivo Labbe, Joel Leja, Tim B. Miller, Brant Robertson, John R. Weaver, Katherine E. Whitaker

8 pages, 4 figures, submitted to ApJL. Comments welcome!

Galaxy sizes and their evolution over cosmic time have been studied for decades and serve as key tests of galaxy formation models. However, at $z\gtrsim1$ these studies have been limited by a lack of deep, high-resolution rest-frame infrared imaging that accurately traces galaxy stellar mass distributions. Here, we leverage the new capabilities of the James Webb Space Telescope to measure the 4.4$\mu$m sizes of ${\sim}1000$ galaxies with $\log{\rm{M}_*/\rm{M}_\odot}\ge9$ and $1.0\le z \le 2.5$ from public CEERS imaging in the EGS deep field. We compare the sizes of galaxies measured from NIRCam imaging at 4.4$\mu$m ($\lambda_{\mathrm{rest}}\sim1.6\mu $m) with sizes measured at $1.5\mu$m ($\lambda_{\mathrm{rest}}\sim5500$A). We find that, on average, galaxy half-light radii are $\sim8$\% smaller at 4.4$\mu$m than 1.5$\mu$m in this sample. This size difference is markedly stronger at higher stellar masses and redder rest-frame $V-J$ colors: galaxies with ${\rm M}_* \sim 10^{11}\,{\rm M}_\odot$ have 4.4$\mu$m sizes that are $\sim 25$\,\% smaller than their 1.5$\mu$m sizes. Our results indicate that galaxy mass profiles are significantly more compact than their rest-frame optical light profiles at cosmic noon, and demonstrate that spatial variations in age and attenuation are important, particularly for massive galaxies. The trend that we find here impacts our understanding of the size growth and evolution of galaxies, and suggests that previous studies based on rest-frame optical light may not have captured the mass-weighted structural evolution of galaxies. This paper represents a first step towards a new understanding of the morphologies of early massive galaxies enabled by JWST's infrared window into the distant universe.

Lorenzo Speri, Andrea Antonelli, Laura Sberna, Stanislav Babak, Enrico Barausse, Jonathan R. Gair, Michael L. Katz

Gravitational wave observations of extreme mass-ratio inspirals (EMRIs) offer the opportunity to probe the environments of active galactic nuclei (AGN) through the torques that accretion disks induce on the binary. Within a Bayesian framework, we study how well such environmental effects can be measured using gravitational wave observations from the Laser Interferometer Space Antenna (LISA). We focus on the torque induced by planetary-type migration on quasi-circular inspirals, and use different prescriptions for geometrically-thin and radiatively-efficient disks. We find that LISA could detect migration for a wide range of disk viscosities and accretion rates, for both $\alpha$ and $\beta$ disk prescriptions. For a typical EMRI with masses $50M_\odot+10^6M_\odot$, we find that LISA could distinguish between migration in $\alpha$ and $\beta$ disks and measure the torque amplitude with $\sim 20\%$ relative precision. Provided an accurate torque model, we also show how to turn gravitational-wave measurements of the torque into constraints on the disk properties. Furthermore, if an electromagnetic counterpart is identified, the multimessenger observations of the AGN-EMRI system could yield direct measurements of the disk viscosity. Finally, we investigate the impact of neglecting environmental effects in the analysis of the gravitational-wave signal of our reference system, finding 3-sigma biases in the primary mass and spin. Our analysis can be easily generalized to any environmental effect, provided that the torque has a simple power law-like dependence on the orbital separation.

Christopher Dessert, Joshua W. Foster, Yujin Park, Benjamin R. Safdi, Weishuang Linda Xu

7+20 pages, 3+23 figures

Thermal higgsino dark matter (DM), with mass around 1 TeV, is a well-motivated, minimal DM scenario that arises in supersymmetric extensions of the Standard Model. Higgsinos may naturally be the lightest superpartners in Split-supersymmetry models that decouple the scalar superpartners while keeping higgsinos and gauginos close to the TeV scale. Higgsino DM may annihilate today to give continuum gamma-ray emission at energies less than a TeV in addition to a line-like signature at energies equal to the mass. Previous searches for higgsino DM, for example with the H.E.S.S. gamma-ray telescope, have not reached the necessary sensitivity to probe the higgsino annihilation cross-section. In this work we make use of 14 years of $\textit{Fermi}$ gamma-ray data at energies above $\sim$10 GeV to search for the continuum emission near the Galactic Center from higgsino annihilation. We interpret our results using DM profiles from Milky Way analogue galaxies in the FIRE-2 hydrodynamic cosmological simulations. We set the strongest constraints to-date on higgsino-like DM. Our results show a mild, $\sim$2$\sigma$ preference for higgsino DM with a mass near the thermal higgsino mass and, depending on the DM density profile, the expected cross-section.

Sungwoo Hong, Gowri Kurup, Maxim Perelstein

50 pages incl. appendix

We consider theories in which a dark sector is described by a Conformal Field Theory (CFT) over a broad range of energy scales. A coupling of the dark sector to the Standard Model breaks conformal invariance. While weak at high energies, the breaking grows in the infrared, and at a certain energy scale the theory enters a confined (hadronic) phase. One of the hadronic excitations can play the role of dark matter. We study a "Conformal Freeze-In" cosmological scenario, in which the dark sector is populated through its interactions with the SM at temperatures when it is conformal. In this scenario, the dark matter relic density is determined by the CFT data, such as the dimension of the CFT operator coupled to the Standard Model. We show that this simple and highly predictive model of dark matter is phenomenologically viable. The observed relic density is reproduced for a variety of SM operators ("portals") coupled to the CFT, and the resulting models are consistent with observational constraints. The mass of the COFI dark matter candidate is predicted to be in the keV-MeV range.

We show that very compact axion mini-clusters can form in models where axion-like-particle (ALP) dark matter is produced via the kinetic misalignment mechanism, which is well-motivated in pre-inflationary $U(1)$ symmetry breaking scenarios. This is due to ALP fragmentation. We predict denser halos than what has been obtained so far in the literature from standard misalignment in post-inflationary $U(1)$ breaking scenarios or from large misalignment. The main reason is that adiabatic fluctuations are significant at early times; therefore, even if amplification from parametric resonance effects is moderate, the final size of ALP fluctuations is larger in kinetic misalignment. We compare halo mass functions and halo spectra obtained in kinetic misalignment, large misalignment, and standard misalignment, respectively. Our analysis does not depend on the specific model realization of the kinetic misalignment mechanism. We present our results generally as a function of the ALP mass and the ALP decay constant only. We show that a sizable region of this ALP parameter space can be tested by future experiments that probe small-scale structures.

The weak proton-proton fusion into a deuteron ($^2$H) is the driving reaction in the energy production in the Sun, as well as similar main sequence stars. Its reaction rate in the solar interior is determined only theoretically. Here, we provide a new determination of the rate of this reaction in solar conditions $S^{11}(0)$, and analyze theoretical and experimental sources for uncertainties, using effective field theory of quantum chromo-dynamics without explicit pions at next-to-leading order. We find an enhancement of $S^{11}$ by $1-4\%$ over the previously recommended value. This change reduces the calculated fluxes of neutrinos originating in $^8$B and $^7$Be nuclear reactions in the Sun, in a way that increases the tension between Solar models and observables, known as the "Solar Composition Problem".

Tetsuya Shiromizu, Keisuke Izumi, Kangjae Lee, Diego Soligon

8 pages

In accordance with current models of the accelerating universe as a spacetime with a positive cosmological constant, new results about a cosmological upper bound for the area of stable marginally outer trapped surfaces are found taking into account angular momentum, gravitational waves and matter. Compared to previous results which take into account only some of the aforementioned variables, the bound is found to be tighter, giving a concrete limit to the size of black holes especially relevant in the early universe.

Petri J. Käpylä (Göttingen University), Nishant K. Singh (IUCAA)

10 pages, 5 figures, submitted to JFM Rapids

Turbulent motions enhance the diffusion of large-scale flows and temperature gradients. Such diffusion is often parameterized by coefficients of turbulent viscosity ($\nu_{\rm t}$) and turbulent thermal diffusivity ($\chi_{\rm t}$) that are analogous to their microscopic counterparts. We compute the turbulent diffusion coefficients by imposing large-scale velocity and temperature gradients on a turbulent flow and measuring the response of the system. We also confirm our results using experiments where the imposed gradients are allowed to decay. To achieve this, we use weakly compressible three-dimensional hydrodynamic simulations of isotropically forced homogeneous turbulence. We find that the turbulent viscosity and thermal diffusion, as well as their ratio the turbulent Prandtl number, ${\rm Pr}_{\rm t} = \nu_{\rm t}/\chi_{\rm t}$, approach asymptotic values at sufficiently high Reynolds and Pecl\'et numbers. We also do not find a significant dependence of ${\rm Pr}_{\rm t}$ on the microscopic Prandtl number ${\rm Pr} = \nu/\chi$. These findings are in stark contrast to results from the $k-\epsilon$ model which suggests that ${\rm Pr}_{\rm t}$ increases monotonically with decreasing ${\rm Pr}$. The current results are relevant for the ongoing debate of, for example, the nature of the turbulent flows in the very low ${\rm Pr}$ regimes of stellar convection zones.

Kaine A. Bunting, Huw Morgan

This paper has been accepted for publication in the Journal of Space Weather and Space Climate (JSWSC)

Accurate forecasting of the solar wind has grown in importance as society becomes increasingly dependent on technology that is susceptible to space weather events. This work describes an inner boundary condition for ambient solar wind models based on tomography maps of the coronal plasma density gained from coronagraph observations, providing a novel alternative to magnetic extrapolations. The tomographical density maps provide a direct constraint of the coronal structure at heliocentric distances of 4 to 8Rs, thus avoiding the need to model the complex non-radial lower corona. An empirical inverse relationship converts densities to solar wind velocities which are used as an inner boundary condition by the Heliospheric Upwind Extrapolation (HUXt) model to give ambient solar wind velocity at Earth. The dynamic time warping (DTW) algorithm is used to quantify the agreement between tomography/HUXt output and in situ data. An exhaustive search method is then used to adjust the lower boundary velocity range in order to optimize the model. Early results show up to a 32% decrease in mean absolute error between the modelled and observed solar wind velocities compared to that of the coupled MAS/HUXt model. The use of density maps gained from tomography as an inner boundary constraint is thus a valid alternative to coronal magnetic models, and offers a significant advancement in the field given the availability of routine space-based coronagraph observations.

In big bang nucleosynthesis (BBN), the light matter abundance is dictated by the neutron-to-proton ($n/p$) ratio which is controlled by the standard weak processes in the early universe. Here, we study the effect of an extra particle species ($\chi$) which co-annihilates with neutron, thereby potentially changing the ($n/p$) ratio in addition to the former processes. We find a novel interplay between the co-annihilation and the weak interaction in deciding the ($n/p$) ratio and the yield of $\chi$. At the initial stage of BBN for the large co-annihilation strength ($G_D$) in comparison to the weak coupling ($G_F$), more neutrons are removed from the thermal bath modifying the ($n/p$) ratio from its standard evolution. We find that the standard BBN prediction is restored for $G_D/G_F \lesssim 10^{-1}$, while the mass of $\chi$ being much smaller than the neutron mass. When the mass of $\chi$ is comparable to the neutron mass, we can allow large $G_D/G_F$ values, as the thermal abundance of $\chi$ becomes Boltzmann-suppressed. Therefore, the ($n/p$) ratio is restored to its standard value via dominant weak processes in later epochs. We also discuss the viability of this new particle to be a dark matter candidate.

Lavinia Heisenberg, Shayan Hemmatyar, Stefan Zentarra

59 pages, 4 figures

We give a detailed exposition of the formalism of Kinetic Field Theory (KFT) with emphasis on the perturbative determination of observables. KFT is a statistical non-equilibrium classical field theory based on the path integral formulation of classical mechanics, employing the powerful techniques developed in the context of quantum field theory to describe classical systems. Unlike previous work on KFT, we perform the integration over the probability distribution of initial conditions in the very last step. This significantly improves the clarity of the perturbative treatment and allows for physical interpretation of intermediate results. We give an introduction to the general framework, but focus on the application to interacting $N$-body systems. Specializing the results to cosmic structure formation, we reproduce the linear growth of the cosmic density fluctuation power spectrum on all scales from microscopic, Newtonian particle dynamics alone.

V. Kumar, F. Pizzi, A. Giesecke, J. Simkanin, Th. Gundrum, M. Ratajczak, F. Stefani

9 pages, 8 figures

The effect of the nutation angle on the flow inside a precessing cylinder is experimentally explored and compared with numerical simulations. The focus is laid on the typical breakdown of the directly forced m=1 Kelvin mode for increasing precession ratio (Poincar\'e number), and the accompanying transition between a laminar and turbulent flow. Compared to the reference case with a 90{\deg} nutation angle, prograde rotation leads to an earlier breakdown, while in the retrograde case the forced mode continues to exist also for higher Poincar\'e numbers. Depending largely on the occurrence and intensity of an axisymmetric double-roll mode, a kinematic dynamo study reveals a sensitive dependency of the self-excitation condition on the nutation angle and the Poincar\'e number. Optimal dynamo conditions are found for 90{\deg} angle which, however, might shift to slightly retrograde precession for higher Reynolds numbers.

The Hubble tension is one of the most exciting problems that Cosmology faces today. A lot of possible solutions for it have already been proposed in the last few years, with a lot of them using a lot of new and exotic physics ideas to deal with the problem. But it was shown recently that the H0 tension might not require new physics but only a more accurate discussion of measurements and interestingly it was the Heisenberg uncertainty principle which was pivotal for that revelation. Accordingly, if one observed the photon mass beyond the indeterminacy through uncertainty then one could in principle reconcile the tension. We examine this in a greater detail in this work by taking into account modifications of the Compton wavelength, which was pivotal in the initial discussion on the interconnection between uncertainty and the H0 tension. We mainly discuss two types of modifications, one based on generalized uncertainty principles (GUP) and the other on higher dimensional physics considerations. We firstly show that both minimal length and maximal momentum GUP based modifications do not provide photon mass values on the required scales and hence we cannot address the tension in this sense at all in this case. We then show that one can get the photon mass to be on the required scales even after incorporating higher dimensional effects, one cannot reconcile the Hubble Tension in the same sense as in the (3+1) space-time case. In this way, we show that certain new physics considerations cannot address the H0 tension in the same sense as one can using the original Heisenberg uncertainty principle in a (3+1) space-time.

Nahid Ahmadi, Mahdiyar Noorbala, Niloufar Feyzabadi, Fatemeh Eghbalpoor, Zahra Ahmadi

31 pages, 7 figures

Transitions between different inflationary slow-roll scenarios are known to provide short non-slow-roll periods with non-trivial consequences. We consider the effect of quantum diffusion on the inflationary dynamics in a transition process. Using the stochastic {\delta}N formalism, we follow the detailed evolution of noises through a sharp transition modeled by the Starobinsky potential, although some of our results apply to any sharp transition. We find how the stochastic noise induced by the transition affects the coarse-grained fields. We then consider the special case that the potential is flat after the transition. It is found that the particular noise we obtain cannot drive the inflaton past the classically unreachable field values. By deriving the characteristic function, we also study the tail behavior for the distribution of curvature perturbations {\zeta}, which we find to decay faster than e^(-3{\zeta}).

E.A. Kuznetsov, E.A. Mikhailov

20 pages, 7 figures

The process of breaking of inviscid incompressible flows along a rigid body with slipping boundary conditions is studied. Such slipping flows are compressible, which is the main reason for the formation of a singularity for the gradient of the velocity component parallel to rigid border. Slipping flows are studied analytically in the framework of two- and three-dimensional inviscid Prandtl equations. Criteria for a gradient catastrophe are found in both cases. For 2D Prandtl equations breaking takes place both for the parallel velocity along the boundary and for the vorticity gradient. For three-dimensional Prandtl flows, breaking, i.e. the formation of a fold in a finite time, occurs for the symmetric part of the velocity gradient tensor, as well as for the antisymmetric part - vorticity. The problem of the formation of velocity gradients for flows between two parallel plates is studied numerically in the framework of two-dimensional Euler equations. It is shown that the maximum velocity gradient grows exponentially with time on a rigid boundary with a simultaneous increase in the vorticity gradient according to a double exponential law. Careful analysis shows that this process is nothing more than the folding, with a power-law relationship between the maximum velocity gradient and its width: $% \max|u_x|\propto \ell^{-2/3}$.