Abstract visibility
Abstract text size

Papers for Friday, Sep 23 2022

Papers with local authors

Yu-Peng Chen, Shu Zhang, Long Ji, Shuang-Nan Zhang, Ling-Da Kong, Peng-Ju Wang, Zhi Chang, Jing-Qiang Peng, Jian Li, Jin-Lu Qu, Zhao-Sheng Li, Lian Tao, Ming-Yu Ge

arXiv admin note: substantial text overlap with arXiv:2208.13556 ; text overlap with arXiv:2208.12124

0 votes
Paper 24 — arXiv:2209.10721
0 votes
Paper 24 — arXiv:2209.10721

After in quiescence for 49 years, 4U~1730--22 became active and had two outbursts in 2021 \& 2022; ten thermonuclear X-ray bursts were detected with Insight-HXMT. Among them, the faintest burst showed a double-peaked profile, placing the source as the 5th accreting neutron star (NS) exhibiting double/triple-peaked type-I X-ray bursts; the other bursts showed photospheric radius expansion (PRE). The properties of double-peaked non-PRE burst indicate that it could be related to a stalled burning front. For the five bright PRE bursts, apart from the emission from the neutron star (NS) surface, we find the residuals both in the soft ($<$3 keV) and hard ($>$10 keV) X-ray band. Time-resolved spectroscopy reveals that the excess can be attributed to an enhanced pre-burst/persistent emission or the Comptonization of the burst emission by the corona/boundary-layer. We find, the burst emission shows a rise until the photosphere touches down to the NS surface rather than the theoretical predicted constant Eddington luminosity. The shortage of the burst emission in the early rising phase is beyond the occlusion by the disk. We speculate that the findings above correspond to that the obscured part (not only the lower part) of the NS surface is exposed to the line of sight due to the evaporation of the obscured material by the burst emission, or the burst emission is anisotropic ($\xi>1$) in the burst early phase. In addition, based on the average flux of PRE bursts at their touch-down time, we derive a distance estimation as 10.4 kpc.

Xin-Ying Song, Shuang-Nan Zhang, Ming-Yu Ge, Shu Zhang

15 pages, 7 figures, accepted for publication in MNRAS on 2022 September 20

0 votes
Paper 31 — arXiv:2209.10832
0 votes
Paper 31 — arXiv:2209.10832

GRB 220426A is a bright gamma-ray burst (GRB) dominated by the photospheric emission. We perform several tests to speculate the origin of this photospheric emission. The dimensionless entropy $\eta$ is large, which is not usual if we assume that it is a pure hot fireball launched by neutrino-antineutrino annihilation mechanism only. Moreover, the outflow has larger $\eta$ with lower luminosity $L$ in the first few seconds, so that the trend of time-resolved $\eta-L$ can not be described as a monotonically positive correlation between $\eta$ and $L$. A hybrid outflow with almost completely thermalized Poynting flux could account for the quasi-thermal spectrum as well as large $\eta$. More importantly, the existence of magnetic field could affect the proton density and neutron-proton coupling effect, so that it could account for the observed trend of time-resolved $\eta-L$. The other origins for the photospheric emission, such as non-dissipative hybrid outflow or magnetic reconnection, are not supported because their radiation efficiencies are low, which is not consistent with non-detection of the afterglow for GRB 220426A. Therefore, we think the hybrid outflow may be the most likely origin.

Shubham Kanodia, Suvrath Mahadevan, Jessica Libby-Roberts, Gudmundur Stefansson, Caleb I. Canas, Anjali A. A. Piette, Alan Boss, Johanna Teske, John Chambers, Greg Zeimann, Andrew Monson, Paul Robertson, Joe P. Ninan, Andrea S.J. Lin, Chad F. Bender, William D. Cochran, Scott A. Diddams, Arvind F. Gupta, Samuel Halverson, Suzanne Hawley, Henry A. Kobulnicky, Andrew J. Metcalf, Brock A. Parker, Luke Powers, Lawrence W. Ramsey, Arpita Roy, Christian Schwab, Tera N. Swaby, Ryan C. Terrien, John Wisniewski

Submitted to ApJ. Comments are welcome. arXiv admin note: text overlap with arXiv:2203.07178

0 votes
Paper 53 — arXiv:2209.11160
0 votes
Paper 53 — arXiv:2209.11160

We present the discovery of TOI-5205b, a transiting Jovian planet orbiting a solar metallicity M4V star, which was discovered using TESS photometry and then confirmed using a combination of precise radial velocities, ground-based photometry, spectra and speckle imaging. The host star TOI-5205 sits near the eponymous `Jao gap', which is the transition region between partially and fully-convective M dwarfs. TOI-5205b has one of the highest mass ratio for M dwarf planets with a mass ratio of almost 0.3$\%$, as it orbits a host star that is just $0.392 \pm 0.015$ $M_{\odot}$. Its planetary radius is $1.03 \pm 0.03~R_J$, while the mass is $1.08 \pm 0.06~M_J$. Additionally, the large size of the planet orbiting a small star results in a transit depth of $\sim 7\%$, making it one of the deepest transits of a confirmed exoplanet orbiting a main-sequence star. The large transit depth makes TOI-5205b a compelling target to probe its atmospheric properties, as a means of tracing the potential formation pathways. While there have been radial velocity-only discoveries of giant planets around mid M dwarfs, this is the first transiting Jupiter with a mass measurement discovered around such a low-mass host star. The high mass of TOI-5205b stretches conventional theories of planet formation and disk scaling relations that cannot easily recreate the conditions required to form such planets.

Papers reserved for later discussion

These are papers reserved by people for discussion at a later date. All reservations are kept for 2 days after the date of the reservation.

Dion Noordhuis, Anirudh Prabhu, Samuel J. Witte, Alexander Y. Chen, Fábio Cruz, Christoph Weniger

5 pages, 2 figures + Supplementary Material

0 votes
09/22/2022: arXiv:2209.09917
0 votes
09/22/2022: arXiv:2209.09917

Axions can be copiously produced in localized regions of neutron star magnetospheres where the ambient plasma is unable to efficiently screen the induced electric field. As these axions stream away from the neutron star they can resonantly transition into photons, generating a large broadband contribution to the neutron star's intrinsic radio flux. In this work, we develop a comprehensive end-to-end framework to model this process from the initial production of axions to the final detection of radio photons, and derive constraints on the axion-photon coupling, $g_{a\gamma\gamma}$, using observations of 27 nearby pulsars. We study the modeling uncertainty in the sourced axion spectrum by comparing predictions from 2.5 dimensional Particle-in-Cell simulations with those derived using a semi-analytic model; these results show remarkable agreement, leading to constraints on the axion-photon coupling that typically differ by a factor of no more than $\sim 3$. The limits presented here are the strongest to date for axion masses $10^{-9} \, {\rm eV} \lesssim m_a \lesssim 10^{-5} \, \rm eV$, and crucially do not rely on the assumption that axions are dark matter.

All other papers

Srimoyee Samaddar, Karthik Venkataramani, Scott. M. Bailey

We use a new Atmospheric Chemistry and Energetics one-dimensional (ACE1D) thermospheric model to show that the energies deposited by the solar soft x-rays in the lower thermosphere at altitudes between 100 -150 km (Bailey et al. 2000), affects the temperature of the entire Earth's thermosphere even at altitudes well above 300 km. By turning off the input solar flux in the different wavelength bins of the model iteratively, we are able to demonstrate that the maximum change in exospheric temperature is due to the changes in the soft x-ray solar bins. We also show, using the thermodynamic heat equation, that the molecular diffusion via non-thermal photoelectrons, is the main source of heat transfer to the upper ionosphere/thermosphere and results in the increase of the temperature of the neutral atmosphere. Moreover, these temperature change and heating effects of the solar soft x-rays are comparable to that of the strong HeII 30.4nm emission. Lastly, we show that the uncertainties in the solar flux irradiance at these soft x-rays wavelengths result in corresponding uncertainties in modeled exospheric temperature and the uncertainties increase with increased solar activity.

Lyra Cao, Marc H. Pinsonneault

25 pages, 19 figures. Accepted for publication in MNRAS

We measure starspot filling fractions for 240 stars in the Pleiades and M67 open star clusters using APOGEE high-resolution H-band spectra. For this work we developed a modified spectroscopic pipeline which solves for starspot filling fraction and starspot temperature contrast. We exclude binary stars, finding that the large majority of binaries in these clusters (80%) can be identified from Gaia DR3 and APOGEE criteria -- important for field star applications. Our data agree well with independent activity proxies, indicating that this technique recovers real starspot signals. In the Pleiades, filling fractions saturate at a mean level of 0.248$\pm$0.005 for active stars with a decline at slower rotation; we present fitting functions as a function of Rossby number. In M67, we recover low mean filling fractions of 0.030$\pm$0.008 and 0.003$\pm$0.002 for main sequence GK stars and evolved red giants respectively, confirming that the technique does not produce spurious spot signals in inactive stars. Starspots also modify the derived spectroscopic effective temperatures and convective overturn timescales. Effective temperatures for active stars are offset from inactive ones by -109$\pm$11 K, in agreement with the Pecaut & Mamajek empirical scale. Starspot filling fractions at the level measured in active stars changes their inferred overturn timescale, which biases the derived threshold for saturation. Finally, we identify a population of stars statistically discrepant from mean activity-Rossby relations and present evidence that these are genuine departures from a Rossby scaling. Our technique is applicable to the full APOGEE catalog, with broad applications to stellar, galactic, and exoplanetary astrophysics.

We use the BayeSN hierarchical probabilistic SED model to analyse the optical-NIR ($BVriYJH$) light curves of 86 Type Ia supernovae (SNe Ia) from the Carnegie Supernova Project to investigate the SN Ia host galaxy dust law distribution and correlations between SN Ia Hubble residuals and host mass. Our Bayesian analysis simultaneously constrains the mass step and dust $R_V$ population distribution by leveraging optical-NIR colour information. We demonstrate how a simplistic analysis where individual $R_V$ values are first estimated for each SN separately, and then the sample variance of these point estimates is computed, overestimates the $R_V$ population variance $\sigma_R^2$. This bias is exacerbated when neglecting residual intrinsic colour variation beyond that due to light curve shape. Instead, Bayesian shrinkage estimates of $\sigma_R$ are more accurate, with fully hierarchical analysis of the light curves being ideal. For the 75 SNe with low-to-moderate reddening (peak apparent $B-V\leq0.3$), we estimate an $R_V$ distribution with population mean $\mu_R=2.59\pm0.14$, and standard deviation $\sigma_R=0.62\pm0.16$. Splitting this subsample at the median host galaxy mass ($10^{10.57}~\mathrm{M}_\odot$) yields consistent estimated $R_V$ distributions between low- and high-mass galaxies, with $\mu_R=2.79\pm0.18$, $\sigma_R=0.42\pm0.24$, and $\mu_R=2.35\pm0.27$, $\sigma_R=0.74\pm0.36$, respectively. When estimating distances from the full optical-NIR light curves while marginalising over various forms of the dust $R_V$ distribution, a mass step of $\gtrsim0.06$ mag persists in the Hubble residuals at the median host mass.

Sumit K. Sarbadhicary, Todd A. Thompson, Laura A. Lopez, Smita Mathur

9 pages, 4 figures. Submitted to MNRAAS

The origin of arcmin-sized Odd Radio Circles (ORCs) found in modern all-sky radio surveys remain uncertain, with explanations ranging from starburst/AGN-driven shocks to supernova remnants (SNRs) in the low-density ambient medium. Using well-calibrated radio light curve models, we assess the possibility that ORCs are radio SNRs evolving in low ambient densities. Our models imply that ORCs 1-5 and J0624-6948 (near the LMC) as SNRs must be within 200 kpc and 100 kpc from the Sun respectively, given their observed flux densities and angular sizes. To be evolving in the circumgalactic medium of the Milky Way, our models require ORCs 1-5 to be ejecta-dominated SNRs within 50 kpc, evolving in ambient densities of $(0.2-1.2) \times 10^{-3}$ cm$^{-3}$. However, this is statistically unlikely because ORCs 1-5 would have ages $<640$ yrs, much smaller than their expected lifetimes of $\gtrsim$10$^5$ yrs at these densities, and because the low SN rate and steep profile of the stellar halo imply a negligible number of ORC-like SNRs within 50 kpc. The circumgalactic medium SNR scenario for J0624-6948 is more likely (though still low probability) compared to ORCs 1-5, as our models allow J0624-6948 to be $\lesssim$3000 yrs. On the other hand, the interpretation of J0624-6948 as a Sedov-Taylor SNR at 50 kpc (LMC) distance is possible for a wide range of ambient densities ($6 \times 10^{-4} - 0.5$ cm$^{-3}$) and ages $\sim$$(0.2-2.6) \times 10^4$ yr, while also being consistent with the local HI environment.

Sam M. Ward, Stephen Thorp, Kaisey S. Mandel, Suhail Dhawan, David O. Jones, Kirsty Taggart, Ryan J. Foley, Gautham Narayan, Kenneth C. Chambers, David A. Coulter, Kyle W. Davis, Thomas de Boer, Kaylee de Soto, Nicholas Earl, Alex Gagliano, Hua Gao, Jens Hjorth, Mark E. Huber, Luca Izzo, Danial Langeroodi, Eugene A. Magnier, Peter McGill, Armin Rest, César Rojas-Bravo, Radosław Wojtak

Submitted to MNRAS; 30 pages, 22 figures

To improve Type Ia supernova (SN Ia) standardisability, the consistency of distance estimates to siblings -- SNe in the same host galaxy -- should be investigated. We present Young Supernova Experiment Pan-STARRS-1 $grizy$ photometry of SN 2021hpr, the third spectroscopically confirmed SN Ia in the high-stellar-mass Cepheid-calibrator galaxy NGC 3147. We analyse NGC 3147's trio of SN Ia siblings: SNe 1997bq, 2008fv and 2021hpr, using a new version of the BayeSN model of SN Ia spectral-energy distributions, retrained simultaneously using optical-NIR $BgVrizYJH$ (0.35--1.8 $\mu$m) data. The distance estimates to each sibling are consistent, with a sample standard deviation $\lesssim$0.01 mag, much smaller than the total intrinsic scatter in the training sample: $\sigma_0\approx0.09$ mag. Fitting normal SN Ia siblings in three additional galaxies, we estimate a $\approx$90% probability that the siblings' intrinsic scatter is smaller than $\sigma_0$. We build a new hierarchical model that fits light curves of siblings in a single galaxy simultaneously; this yields more precise estimates of the common distance and the dust parameters. Fitting the trio for a common dust law shape yields $R_V=2.69\pm0.52$. Our work motivates future hierarchical modelling of more siblings, to tightly constrain their intrinsic scatter, and better understand SN-host correlations. Finally, we estimate the Hubble constant, using a Cepheid distance to NGC 3147, the siblings trio, and 109 Hubble flow ($0.01 < z_{\rm{CMB}} < 0.08$) SNe Ia; marginalising over the siblings' and population's intrinsic scatters, and the peculiar velocity dispersion, yields $H_0=77.9\pm6.5 \text{ km s}^{-1}\text{Mpc}^{-1}$.

James W. Nightingale, Qiuhan He, Xiaoyue Cao, Aristeidis Amvrosiadis, Amy Etherington, Carlos S. Frenk, Richard G. Hayes, Andrew Robertson, Shaun Cole, Samuel Lange, Ran Li, Richard Massey

27 Pages, 15 Figures

The cold dark matter (DM) model predicts that every galaxy contains thousands of DM subhalos; almost all other DM models include a physical process that smooths away the subhalos. The subhalos are invisible, but could be detected via strong gravitational lensing, if they lie on the line of sight to a multiply-imaged background source, and perturb its apparent shape. We present an automated strong lens analysis framework, and scan for DM subhalos in Hubble Space Telescope imaging of 54 strong lenses. We identify two compelling DM subhalo candidates (including one previously found in SLACS0946+1006), where a subhalo is favoured after every systematic test we perform. We find that the detectability of subhalos depends upon the assumed parametric form for the lens galaxy's mass distribution. Comparing fits which assume several more complex mass models reveals $5$ additional (generally lower mass) DM subhalo candidates worthy of further study, and the removal of 11 false positives. We identify 44 non-detections, which are vital to building up enough statistical power to test DM models. Future work will apply even more flexible models to the results of this study, to constrain different DM models. Our full analysis results are available at https://github.com/Jammy2211/autolens_subhalo.

Aniket Bhagwat, Benedetta Ciardi, Erik Zackrisson, Joop Schaye

Accepted in MNRAS. Comments welcome

Overdense, metal-rich regions, shielded from stellar radiation might remain neutral throughout reionization and produce metal as well as 21 cm absorption lines. Simultaneous absorption from metals and 21 cm can complement each other as probes of underlying gas properties. We use Aurora, a suite of high resolution radiation-hydrodynamical simulations of galaxy formation, to investigate the occurrence of such "aligned" absorbers. We calculate absorption spectra for 21 cm, OI, CII, SiII and FeII. We find velocity windows with absorption from at least one metal and 21 cm, and classify the aligned absorbers into two categories: 'aligned and cospatial absorbers' and 'aligned but not cospatial absorbers'. While 'aligned and cospatial absorbers' originate from overdense structures and can be used to trace gas properties, 'aligned but not cospatial absorbers' are due to peculiar velocity effects. The incidence of absorbers is redshift dependent, as it is dictated by the interplay between reionization and metal enrichment, and shows a peak at $z \approx 8$ for the aligned and cospatial absorbers. While aligned but not cospatial absorbers disappear towards the end of reionization because of the lack of an ambient 21 cm forest, aligned and cospatial absorbers are associated with overdense pockets of neutral gas which can be found at lower redshift. We produce mock observations for a set of sightlines for the next generation of telescopes like the ELT and SKA1-LOW, finding that given a sufficiently bright background quasar, these telescopes will be able to detect both types of absorbers throughout reionization.

P. Gandhi (Southampton), T. Kawamuro (UDP), M. Díaz Trigo (ESO), J.A. Paice, P.G. Boorman, M. Cappi, C. Done, A.C. Fabian, K. Fukumura, J.A. Garcia, C.L. Greenwell, M. Guainazzi, K. Makishima, M.S. Tashiro, R. Tomaru, F. Tombesi, Y. Ueda

To appear in Nature Astronomy as a review. Author version, before final editorial and style revisions

Microcalorimeters have demonstrated success in delivering high spectral resolution, and have paved the path to revolutionary new science possibilities in the coming decade of X-ray astronomy. There are several research areas in compact object science that can only be addressed with energy resolution Delta(E)<~5 eV at photon energies of a few keV, corresponding to velocity resolution of <~a few hundred km/s, to be ushered in by microcalorimeters. Here, we review some of these outstanding questions, focusing on how the research landscape is set to be transformed (i) at the interface between accreting supermassive black holes and their host galaxies, (ii) in unravelling the structures of accretion environments, (iii) in resolving long-standing issues on the origins of energy and matter feedback, and (iv) to test mass-scaled unification of accretion and feedback. The need to learn lessons from Hitomi and to make improvements in laboratory atomic data precision as well as plasma modeling are highlighted.

Scott W McIntosh, Robert J. Leamon, R. Egeland

16 pages, 10 figures - Submitted to Frontiers

McIntosh and colleagues identified an event in the solar timeline that appeared to play a role in how Sunspot Cycle 23 (SC23) transitioned into Sunspot Cycle 24 (SC24). The timeframe for this transition was rapid, taking place in as short as time as a solar rotation. M2014 inferred that the transition observed was a critical episode for the Sun's global-scale magnetic field that was being manifest in the spatially and temporally overlapping and magnetic systems belonging to the Sun's 22-year (Hale) magnetic cycle. These events have been dubbed as Hale Cycle terminations, or `terminators' for short. Further exploration revealed a relationship between terminator separation (as a measure of overlap in the Hale Cycles) and the upcoming sunspot cycle amplitude. McIntosh and colleagues extrapolated upon this relationship to identify the termination of the SC24 carrying Hale Cycle band in Mid-2020 and inferred that this would result in a very large Sunspot Cycle 25 (SC25). This paper presents observational analysis of the end of SC24 and the initial months of SC25 growth following a terminator that occurred in mid-December 2021 (approximately 12/13/2021). We use the December 2021 terminator to finalize the forecast of SC25 amplitude 184 ($\pm$17 with 95\% confidence, and $\pm$63 with 68\% confidence). Finally, we use other terminator-related superposed epoch analyses to project the timing of SC25 maxima in late 2023 to mid 2024.

Galactic outflows can be powered either by nuclear starbursts (SB) or active galactic nuclei (AGN). It has been argued that extreme starbursts can power extreme outflows, without the need to invoke AGN feedback. However, contributions from past and/or hidden AGN activity cannot be ruled out. Here, we constrain the potential role of the central black hole in driving powerful outflows in starburst galaxies (with no sign of ongoing AGN activity). We examine whether the galactic outflows can be explained by AGN luminosity evolution in the framework of our AGN `radiative dusty feedback' scenario. We show that the outflow energetics of starburst galaxies in the local Universe can be quantitatively reproduced by power-law and exponential luminosity decays, coupled with radiation trapping. Likewise, a combination of heavy obscuration and mild luminosity decay may account for the energetics of galactic outflows observed in dusty star-forming galaxies in the early Universe. We discuss different physical arguments for SB vs. AGN outflow-driving, and conclude that the latter can have a major impact on the evolution of galaxies.

Vo Hong Minh Phan, Sarah Recchia, Philipp Mertsch, Stefano Gabici

14 pages, 5 figures

Cosmic rays are the only agent able to penetrate into the interior of dense molecular clouds. Depositing (part of) their energy through ionisation, cosmic rays play an essential role in determining the physical and chemical evolution of star-forming regions. To a first approximation their effect can be quantified by the cosmic-ray induced ionization rate. Interestingly, theoretical estimates of the ionization rate assuming the cosmic-ray spectra observed in the local interstellar medium result in an ionization rate that is one to two orders of magnitude below the values inferred from observations. However, due to the discrete nature of sources, the local spectra of MeV cosmic rays are in general not representative for the spectra elsewhere in the Galaxy. Such stochasticity effects have the potential of reconciling modelled ionization rates with measured ones. Here, we model the distribution of low-energy cosmic-ray spectra expected from a statistical population of supernova remnants in the Milky Way. The corresponding distribution for the ionization rate is derived and confronted with data. We find that the stochastic uncertainty helps with explaining the surprisingly high ionization rates observed in many molecular clouds.

Allen North, Timothy Banks

Accepted by JAAVSO, 12 pages, 5 figures, 2 tables

Eleven transit light curves for the exoplanet WASP-140b were studied with the primary objective to investigate the possibility of transit timing variations (TTVs). Previously unstudied MicroObservatory and Las Cumbres Global Telescope Network photometry were analysed using Markov Chain Monte Carlo techniques, including new observations collected by this study of a transit in December 2021. No evidence was found for TTVs. We used two transit models coupled with Bayesian optimization to explore the physical parameters of the system. The radius for WASP-140b was estimated to be $1.38^{+0.18}_{-0.17}$ Jupiter radii, with the planet orbiting its host star in $2.235987 \pm 0.000008$ days at an inclination of $85.75 \pm 0.75$ degrees. The derived parameters are in formal agreement with those in the exoplanet discovery paper of 2016, and somewhat larger than a recent independent study based on photometry by the TESS space telescope.

Melaine Saillenfest, Zeeve Rogoszinski, Giacomo Lari, Kevin Baillié, Gwenaël Boué, Aurélien Crida, Valéry Lainey

Accepted for publication in Astronomy and Astrophysics

Context. The 98{\deg}-obliquity of Uranus is commonly attributed to giant impacts that occurred at the end of the planetary formation. This picture, however, is not devoid of weaknesses. Aims. On a billion-year timescale, the tidal migration of the satellites of Jupiter and Saturn has been shown to strongly affect their spin-axis dynamics. We aim to revisit the scenario of tilting Uranus in light of this mechanism. Methods. We analyse the precession spectrum of Uranus and identify the candidate secular spin-orbit resonances that could be responsible for the tilting. We determine the properties of the hypothetical ancient satellite required for a capture and explore the dynamics numerically. Results. If it migrates over 10 Uranus' radii, a single satellite with minimum mass 4e-4 Uranus' mass is able to tilt Uranus from a small obliquity and make it converge towards 90{\deg}. In order to achieve the tilting in less than the age of the Solar System, the mean drift rate of the satellite must be comparable to the Moon's current orbital expansion. Under these conditions, simulations show that Uranus is readily tilted over 80{\deg}. Beyond this point, the satellite is strongly destabilised and triggers a phase of chaotic motion for the planet's spin axis. The chaotic phase ends when the satellite collides into the planet, ultimately freezing the planet's obliquity in either a prograde, or plainly retrograde state (as Uranus today). Spin states resembling that of Uranus can be obtained with probabilities as large as 80%, but a bigger satellite is favoured, with mass 1.7e-3 Uranus' mass or more. Yet, a smaller ancient satellite is not categorically ruled out, and there is room for improving this basic scenario in future studies. Interactions among several pre-existing satellites is a promising possibility.

G. Maciejewski, M. Fernandez, A. Sota, P. J. Amado, D. Dimitrov, Y. Nikolov, J. Ohlert, M. Mugrauer, R. Bischoff, T. Heyne, F. Hildebrandt, W. Stenglein, A. A. Arevalo, S. Neira, L. A. Riesco, V. Sanchez Martinez, M. M. Verdugo

Accepted for publication in A&A

Hot Jupiters on extremely short-period orbits are expected to be unstable to tidal dissipation and spiral toward their host stars. That is because they transfer the angular momentum of the orbital motion through tidal dissipation into the stellar interior. Although the magnitude of this phenomenon is related to the physical properties of a specific star-planet system, statistical studies show that tidal dissipation might shape the architecture of hot Jupiter systems during the stellar lifetime on the main sequence. The efficiency of tidal dissipation remains poorly constrained in star-planet systems. Stellar interior models show that the dissipation of dynamical tides in radiation zones could be the dominant mechanism driving planetary orbital decay. These theoretical predictions can be verified with the transit timing method. We acquired new precise transit mid-times for five planets. They were previously identified as the best candidates for which orbital decay might be detected. Analysis of the timing data allowed us to place tighter constraints on the orbital decay rate. No statistically significant changes in their orbital periods were detected for all five hot Jupiters in systems HAT-P-23, KELT-1, KELT-16, WASP-18, and WASP-103. For planets HAT-P-23 b, WASP-18 b, and WASP-103 b, observations show that the mechanism of the dynamical tides dissipation probably does not operate in their host stars, preventing them from rapid orbital decay. This finding aligns with the models of stellar interiors of F-type stars, in which dynamical tides are not fully damped due to convective cores. For KELT-16 b, the span of transit timing data was not long enough to verify the theoretical predictions. KELT-1 b was identified as a potential laboratory for studying the dissipative tidal interactions of inertial waves in a convective layer.

F. Espinoza-Arancibia, M. Catelan, G. Hajdu, N. Rodríguez-Segovia, G. Boggiano, K. Joachimi, C. Muñoz-López, C. Ordenes-Huanca, C. Orquera-Rojas, P. Torres, Á. Valenzuela-Navarro

19 pages, 17 figures. Accepted by MNRAS

Pulsating stars, such as Cepheids and RR Lyrae, offer us a window to measure and study changes due to stellar evolution. In this work, we study the former by calculating a set of evolutionary tracks of stars with an initial mass of 4 to 7 $M_\odot$, varying the initial rotation rate and metallicity, using the stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA). Using Radial Stellar Pulsations (RSP), a recently added functionality of MESA, we obtained theoretical instability strip (IS) edges and linear periods for the radial fundamental mode. Period-age, period-age-temperature, period-luminosity, and period-luminosity-temperature relationships were derived for three rotation rates and metallicities, showing a dependence on crossing number, position in the IS, rotation, and metallicity. We calculated period change rates (PCRs) based on the linear periods from RSP. We compared our models with literature results using the Geneva code, and found large differences, as expected due to the different implementations of rotation between codes. In addition, we compared our theoretical PCRs with those measured in our recent work for Large Magellanic Cloud Cepheids. We found good overall agreement, even though our models do not reach the short-period regime exhibited by the empirical data. Implementations of physical processes not yet included in our models, such as pulsation-driven mass loss, an improved treatment of convection that may lead to a better description of the instability strip edges, as well as consideration of a wider initial mass range, could all help improve the agreement with the observed PCRs.

D. Cont, F. Yan, A. Reiners, L. Nortmann, K. Molaverdikhani, E. Pallé, Th. Henning, I. Ribas, A. Quirrenbach, J. A. Caballero, P. J. Amado, S. Czesla, F. Lesjak, M. López-Puertas, P. Mollière, D. Montes, G. Morello, E. Nagel, S. Pedraz, A. Sánchez-López

Accepted for publication in A&A

Ultra-hot Jupiters are highly irradiated gas giant exoplanets on close-in orbits around their host stars. We analyzed high-resolution spectra from CARMENES, HARPS-N, and ESPaDOnS taken over eight observation nights to study the emission spectrum of WASP-33b and draw conclusions about its atmosphere. By applying the cross-correlation technique, we detected the spectral signatures of Ti I, V I, and a tentative signal of Ti II for the first time via emission spectroscopy. These detections are an important finding because of the fundamental role of Ti- and V-bearing species in the planetary energy balance. Moreover, we assessed and confirm the presence of OH, Fe I, and Si I from previous studies. The spectral lines are all detected in emission, which unambiguously proves the presence of an inverted temperature profile in the planetary atmosphere. By performing retrievals on the emission lines of all the detected species, we determined a relatively weak atmospheric thermal inversion extending from approximately 3400 K to 4000 K. We infer a supersolar metallicity close to 1.5 dex in the planetary atmosphere, and find that its emission signature undergoes significant line broadening with a Gaussian FWHM of about 4.5 km/s. Also, we find that the atmospheric temperature profile retrieved at orbital phases far from the secondary eclipse is about 300 K to 700 K cooler than that measured close to the secondary eclipse, which is consistent with different day- and nightside temperatures. Moreover, retrievals performed on the emission lines of the individual chemical species lead to consistent results, which gives additional confidence to our retrieval method. Increasing the number of species included in the retrieval and expanding the set of retrieved atmospheric parameters will further advance our understanding of exoplanet atmospheres.

L. Colzi, D. Romano, F. Fontani, V. M. Rivilla, L. Bizzocchi, M. T. Beltrán, P. Caselli, D. Elia, L. Magrini

19 pages, 8 figures, 6 tables, 4 appendix - Accepted for publication in A&A

Nitrogen isotopic ratios are a key tool for tracing Galactic stellar nucleosynthesis. We present the first study of the $^{14}$N/$^{15}$N abundance ratio in the outer regions of the Milky Way (namely, for galactocentric distances, $R_{\rm GC}$, from 12 kpc up to 19 kpc), with the aim to study the stellar nucleosynthesis effects in the global Galactic trend. We analysed IRAM 30m observations towards a sample of 35 sources in the context of the CHEMical complexity in star-forming regions of the OUTer Galaxy (CHEMOUT) project. We derived the $^{14}$N/$^{15}$N ratios from HCN and HNC for 14 and 3 sources, respectively, using the $J$ = 1-0 rotational transition of HN$^{13}$C, H$^{15}$NC, H$^{13}$CN, and HC$^{15}$N. The results found in the outer Galaxy have been combined with previous measurements obtained in the inner Galaxy. We find an overall linear decreasing H$^{13}$CN/HC$^{15}$N ratio with increasing $R_{\rm GC}$. This translates to a parabolic $^{14}$N/$^{15}$N ratio with a peak at 11 kpc. Updated Galactic chemical evolution models have been taken into account and compared with the observations. The parabolic trend of the $^{14}$N/$^{15}$N ratio with $R_{\rm GC}$ can be naturally explained (i) by a model that assumes novae as the main $^{15}$N producers on long timescales ($\ge$1 Gyr) and (ii) by updated stellar yields for low- and intermediate-mass stars.

Sarvesh Gharat, Bhaskar Bose

Output data will be made available on request after 2 years of publication

Solar Plages are bright chromospheric features observed in Ca II K photographic observations of the sun. These are regions of high magnetic field concentration thus tracer of magnetic activity of the Sun and are one of the most important features to study long-term variability of the Sun as Ca II K spectroheliograms are recorded for more than a century. . However, detection of the plages from century-long databases is a non-trivial task and need significant human resources for doing it manually. Hence, in this study, we propose an image processing algorithm that can identify solar plages from Ca II K photographic observations. The proposed study has been implemented on archival data from Kodaikanal Solar Observatory. To ensure that the algorithm works, irrespective of noise level, brightness, and other image properties, we randomly draw a sample of images from the data archive to test our algorithm.

Kaveh Pahlevan, Laura Schaefer, Linda T. Elkins-Tanton, Steven J. Desch, Peter R. Buseck

5 figures

The deuterium-to-hydrogen (D/H or 2H/1H) ratio of Martian atmospheric water (~6x standard mean ocean water, SMOW) is higher than that of known sources, requiring planetary enrichment. A recent measurement by NASA's Mars Science Laboratory rover Curiosity of >3 Gyr clays yields a D/H ratio ~3x SMOW, demonstrating that most enrichment occurs early in Mars's history. As on Venus, Mars's D/H enrichment is thought to reflect preferential loss to space of 1H (protium) relative to 2H (deuterium), but the global environmental context of large and early hydrogen losses remain to be determined. Here, we apply a recent model of primordial atmosphere evolution to Mars, link the magma ocean of the accretion epoch with a subsequent water-ocean epoch, and calculate the behavior of deuterium for comparison with the observed record. We find that a ~2-3x hydrospheric deuterium-enrichment is produced if the Martian magma ocean is chemically reducing at last equilibration with the primordial atmosphere, making H2-CO the initially dominant species, with minor abundances of H2O-CO2. Reducing gases - in particular H2 - can cause greenhouse warming and prevent a water ocean from freezing immediately after the magma ocean epoch. Moreover, the pressure-temperature conditions are high enough to produce ocean-atmosphere H2O-H2 isotopic equilibrium such that surface H2O strongly concentrates deuterium relative to H2, which preferentially takes up protium and escapes from the primordial atmosphere. The proposed scenario of primordial H2-rich outgassing and escape suggests significant durations (>Myr) of chemical conditions on the Martian surface conducive to prebiotic chemistry immediately following Martian accretion.

Quanzhi Ye, Peter Jenniskens

To appear in Comets III

Earth occasionally crosses the debris streams produced by comets and other active bodies in our solar system. These manifest meteor showers that provide an opportunity to explore these bodies without a need to visit them in-situ. Observations of meteor showers provide unique insights into the physical and dynamical properties of their parent bodies, as well as into the compositions and the structure of near-surface dust. In this chapter, we discuss the development and current state of affairs of meteor science, with a focus on its role as a tool to study comets, and review the established parent body -- meteor shower linkages.

Yingjie Zhu, Judit Szente, Enrico Landi

Submitted to ApJ, 26 pages, 18 figures. Jupyter notebooks are available at this https URL Comments are welcome

Physical quantities, such as ion temperature and nonthermal velocity, provide critical information about the heating mechanism of the million-degree solar corona. We determined the possible ion temperature $T_i$ intervals using extreme ultraviolet (EUV) line widths, only assuming that the plasma nonthermal velocity is the same for all ions. We measured ion temperatures at the polar coronal hole boundary simultaneously observed in 2007 by the EUV Imaging Spectrometer (EIS) on board the Hinode satellite and the Solar Ultraviolet Measurements of Emitted Radiation (SUMER) on board the Solar and Heliospheric Observatory (SOHO). The temperatures of ions with the charge-to-mass ratio ($Z/A$) less than 0.20 or greater than 0.33 are much higher than the local electron temperature. The measured ion temperature decreases with the $Z/A$ to 0.25 and then increases with the charge-to-mass ratio. We ran the Alfv\'en Wave Solar Model-realtime (AWSoM-R) and the SPECTRUM module to validate the ion temperature diagnostic technique and to help interpret the results. We suggest that the widths of hot lines in the coronal hole (e.g., Fe XII, Fe XIII) are also affected by the solar wind bulk motions along the line of sight. We discussed the factors that might affect the line width fitting, including the instrumental width and non-Gaussian wings in some bright SUMER lines that can be fitted by a double-Gaussian or a $\kappa$ distribution. Our study confirms the presence of preferential heating of heavy ions in coronal holes and provides new constraints to coronal heating models.

J. Berthier (1), B. Carry (2), M. Mahlke (2), J. Normand (1) ((1) IMCCE/Observatoire de Paris, (2) OCA/Laboratoire Lagrange)

The sample of Solar system objects has dramatically increased over the last decade. The amount of measured properties (e.g., diameter, taxonomy, rotation period, thermal inertia) has grown even faster. However, this wealth of information is spread over a myriad of articles, under many different designations per object. We provide a solution to the identification of Solar system objects from any of their multiple names or designations. We also compile and rationalize their properties to provide an easy access to them. We aim to continuously update the database as new measurements become available. We built a Web Service, SsODNet, that offers four access points, each corresponding to an identified necessity in the community: name resolution (quaero), compilation of a large corpus of properties (datacloud), determination of the best estimate among compiled values (ssoCard), and statistical description of the population (ssoBFT). The SsODNet interfaces are fully operational and freely accessible to everyone. The name resolver quaero translates any of the ~5.3 million designations of objects into their current official designation. The datacloud compiles about 105 million parameters (osculating and proper elements, pair and family membership, diameter, albedo, mass, density, rotation period, spin coordinates, phase function parameters, colors, taxonomy, thermal inertia, and Yarkovsky drift) from over 3,000 articles (and growing). For each of the known asteroids and dwarf planets (~1.2 million), a ssoCard providing a single best-estimate for each parameter is available. The SsODNet service provides these resources in a fraction of second upon query. Finally, the large ssoBFT table compiles all the best-estimates in a single table for population-wide studies.

Xiang-Lei Chen, A-Li Luo, Jian-Jun Chen, Rui Wang, Xiao-Bin Zhang, Wen Hou, Bo Qiu, Fang Zuo

14 pages, 21 figures, accepted by A&A

Based on luminosity contributions, we develop a spectroscopic modelling method to derive atmospheric parameters of component stars in binary systems. The method is designed for those spectra of binaries which show double-lined features due to the radial velocities differences between the component stars. We first derive the orbital parameters and the stellar radii by solving the light and radial velocity curves. Then the luminosity contributions in different phases can be calculated. The synthesised double-lined spectra model is constructed by superposing theoretical single-star spectra according to the luminosity contributions. Finally, we derive the atmospheric parameters of each component star by the model fitting method. For multi-epoch double-lined spectra observed by the Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) Medium Resolution Survey ($R \sim 7500$), our method gives robust results for detached eclipsing binary systems observed in different orbital phases. Furthermore, this method can also be applied to other spectroscopic data with different resolutions as long as the systems are detached eclipsing binaries with nearly spherical stars.

Jong-Ho Shinn

10 pages, 6 figures, MNRAS in press, online supplementary data are included in the source file

In my previous reanalysis of the local star-forming galaxies observed in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) pilot survey, I reported that the overestimation of E(B - V), hence the star formation rate (SFR), undermined the claim of new galaxy population discovery in the original study. Herein, I re-examine whether the E(B - V) overestimation problem can be alleviated in the Bayesian parameter estimation framework by adopting scientifically motivated new priors. I modelled the emission-line fluxes of galaxies using the strong-line method and four model parameters-metallicity 12 + log (O/H), nebula emission-line colour excess E(B - V), intrinsic [O III] $\lambda$5007 line flux and intrinsic [N II] $\lambda$6584 line flux. Based on mock-data tests, I found that all input values can be recovered within and around the 1-$\sigma$ credible interval by adopting suitable priors for the intrinsic [O III] $\lambda$5007 and [N II] $\lambda$6584 line fluxes: the inverse gamma distribution reflecting the logical constraint that an intrinsic emission-line flux must exceed the observed (reddened) emission-line flux. The mock-data tests were performed for two metallicity calibrations, three colour excess input values [E(B - V) = 0.1, 0.3 and 0.5] and two metallicity input values [12 + log (O/H) = 8.0 and 8.5]. I also found that the new prior can diminish the SFR overestimation eightfold. This study demonstrates how the Bayesian parameter estimation can achieve more accurate estimates with no further observations when the likelihood does not constrain the model parameters correctly.

James R. Beattie, Christoph Federrath, Neco Kriel, Philip Mocz, Amit Seta

18 pages. 14 figures. Submitted to MNRAS. Comments welcome

The small-scale turbulent dynamo (SSD) is likely to be responsible for the magnetisation of the interstellar medium (ISM) that we observe in the Universe today. The SSD efficiently converts kinetic energy $E_{\rm kin}$ into magnetic energy $E_{\rm mag}$, and is often used to explain how an initially weak magnetic field with $E_{\rm mag} \ll E_{\rm kin}$ is amplified, and then maintained at a level $E_{\rm mag} \lesssim E_{\rm kin}$. Usually, this process is studied by initialising a weak seed magnetic field and letting the turbulence grow it to saturation. However, in this study, using three-dimensional, non-ideal magnetohydrodynamical turbulence simulations, we show that the same saturated state can also be achieved if initially $E_{\rm mag} \gg E_{\rm kin}$ or $E_{\rm mag} \sim E_{\rm kin}$. This is realised through a two-stage exponential decay (1. a slow backreaction that converts $E_{\rm mag}$ into $E_{\rm kin}$, and 2. Ohmic dissipation concentrated in anisotropic current sheets) into the saturated state, for which we provide an analytical model. This means that even if there are temporary local enhancements of $E_{\rm mag}$ in the ISM, such that $E_{\rm mag} > E_{\rm kin}$, e.g., through amplifications such as compressions, over a long enough time the field will decay into the saturated state set by the SSD, which is determined by the turbulence and magnetic dissipation. However, we also provide analytical models for the decay times and utilise wait-time statistics from compressive supernova events to show that if the magnetic field is enhanced above the saturated state, it will not have enough time to decay the field before the next supernova event. Hence, unless there exists a mechanism for destroying magnetic fields that is not in our non-ideal MHD models, the amplitudes of interstellar magnetic fields may also be a... (abridged).

Simone R. Hagey, Billy Edwards, Aaron C. Boley

65 pages, 10 figures, accepted to AJ

We analyze a large number of citizen science data and identify eight Hot Jupiter systems that show evidence for deviations from a constant orbital period: HAT-P-19 b, HAT-P-32 b, TrES-1 b, TrES-2 b, TrES-5 b, WASP-4 b, WASP-10 b, and WASP-12 b. The latter system is already well known to exhibit strong evidence for tidal orbital decay and serves as an important control for this study. Several other systems we identify have disputed period drifts in the literature, allowing the results here to serve as an independent analysis. The citizen science data are from the Exoplanet Transit Database (ETD), which is a global project established in 2008 by the Variable Star and Exoplanet Section of the Czech Astronomical Society. With over 400 planets and 12,000 contributed observations spanning 15 years, the ETD is brimming with potential for studying the long-term orbital evolution of close-in Hot Jupiters. We use our results to discuss prioritization of targets for follow up investigations, which will be necessary to confirm the period drifts and their causes.

Yingkang Zhang, Tao An, Sandor Frey, Krisztina Eva Gabanyi, Yulia Sotnikova

28 pages, 4 figures, 8 tables. To appear in ApJ

Up to now, jet kinematic studies of radio quasars have barely reached beyond the redshift range at $z>3.5$. This significantly limits our knowledge of high-redshift jets, which can provide key information for understanding the jet nature and the growth of the black holes in the early Universe. In this paper, we selected 9 radio-loud quasars at $z>3.5$ which display milliarcsec-scale jet morphology. We provided evidence on the source nature by presenting high-resolution very long baseline interferometry (VLBI) images of the sample at 8.4~GHz frequency and making spectral index maps. We also consider Gaia optical positions that are available for 7 out of the 9 quasars, for a better identification of the jet components within the radio structures. We find that 6 sources can be classified as core--jet blazars. The remaining 3 objects are more likely young, jetted radio sources, compact symmetric objects. By including multi-epoch archival VLBI data, we also obtained jet component proper motions of the sample and estimated the jet kinematic and geometric parameters (Doppler factor, Lorentz factor, viewing angle). Our results show that at $z>3.5$, the jet apparent transverse speeds do not exceed 20 times the speed of light ($c$). This is consistent with earlier high-redshift quasar measurements in the literature and the tendency derived from low-redshift blazars that fast jet speeds ($>40\,c$) only occur at low redshifts. The results from this paper contribute to the understanding of the cosmological evolution of radio AGN.

We propose a theoretical explanation of absorption/emission line systems in classical novae based on a fully self-consistent nova explosion model. We found that a reverse shock is formed far outside the photosphere ($\gtrsim 10^{13}$ cm) because later-ejected mass with a faster velocity collides with earlier-ejected matter. Optically thick winds blow continuously at a rate of $\sim 10^{-4} ~M_\odot$ yr$^{-1}$ near optical maximum but its velocity decreases toward optical maximum and increases afterward, so that the shock arises only after optical maximum. The nova ejecta is divided by the shock into three parts, outermost expanding gas (earliest wind before maximum), shocked shell, and inner fast wind, which respectively contribute to premaximum, principal, and diffuse enhanced absorption/emission line systems. A large part of nova ejecta is eventually confined to the shocked shell. The appearance of principal system is consistent with the emergence of shock. This shock is strong enough to explain thermal hard X-ray emissions. The shocked layer has a high temperature of $k T_{\rm sh} \sim 1$ keV $\times ((v_{\rm wind} - v_{\rm shock})/{\rm 1000 ~km~s}^{-1})^2 ={\rm 1 ~keV}\times ((v_{\rm d} - v_{\rm p})/{\rm 1000 ~km~s}^{-1})^2$, where $v_{\rm d} - v_{\rm p}$ is the velocity difference between the diffuse enhanced ($v_{\rm d}$) and principal ($v_{\rm p}$) systems. We compare a 1.3 $M_\odot$ white dwarf model with the observational properties of the GeV gamma-ray detected classical nova V5856 Sgr (ASASSN-16ma) and discuss what kind of novae can produce GeV gamma-ray emissions.

John M. Dickey, Jennifer West, Alec J.M. Thomson, T.L. Landecker, A. Bracco, E. Carretti, J.L. Han, A.S. Hill, Y.K. Ma, S. A. Mao, A. Ordog, Jo-Anne C. Brown, K. A. Douglas, A. Erceg, V. Jelic, R. Kothes, M. Wolleben

37 pages, 26 figures, Ap. J. accepted

Magnetic fields in the ionized medium of the disk and halo of the Milky Way impose Faraday rotation on linearly polarized radio emission. We compare two surveys mapping the Galactic Faraday rotation, one showing the rotation measures of extragalactic sources seen through the Galaxy (from Hutschenreuter et al 2022), and one showing the Faraday depth of the diffuse Galactic synchrotron emission from the Global Magneto-Ionic Medium Survey. Comparing the two data sets in 5deg x 10deg bins shows good agreement at intermediate latitudes, 10 < |b| < 50 deg, and little correlation between them at lower and higher latitudes. Where they agree, both tracers show clear patterns as a function of Galactic longitude: in the Northern Hemisphere a strong sin(2 x longitude) pattern, and in the Southern hemisphere a sin(longitude + pi) pattern. Pulsars with height above or below the plane |z| > 300 pc show similar longitude dependence in their rotation measures. Nearby non-thermal structures show rotation measure shadows as does the Orion-Eridanus superbubble. We describe families of dynamo models that could explain the observed patterns in the two hemispheres. We suggest that a field reversal, known to cross the plane a few hundred pc inside the solar circle, could shift to positive z with increasing Galactic radius to explain the sin(2xlongitude) pattern in the Northern Hemisphere. Correlation shows that rotation measures from extragalactic sources are one to two times the corresponding rotation measure of the diffuse emission, implying Faraday complexity along some lines of sight, especially in the Southern hemisphere.

Ka-Xing Lu, Jin-Ming Bai, Jian-Min Wang, Chen Hu, Yan-Rong Li, Pu Du, Ming Xiao, Hai-Cheng Feng, Sha-Sha Li, Jian-Guo Wang, Zhi-Xiang Zhang, Ying-Ke Huang

23 pages, 15 figures, accepted for publication in ApJS

NGC 5548 is one of the active galactic nuclei (AGN) selected for our long-term spectroscopic monitoring with the Lijiang 2.4~m telescope, aiming at investigating the origin and evolution of the broad-line regions (BLRs), accurately measuring the mass of the supermassive black holes (SMBHs), and understanding structure and evolution of the AGN. We have performed five-season observations for NGC~5548 with the median sampling interval ranging from 1.25 to 3 days. The light curves of the 5100~\AA\ continuum and broad emission lines are measured after subtracting contamination of the host galaxy starlight. The time lags of the broad He~{\sc ii}, He~{\sc i}, H$\gamma$, and H$\beta$ lines with respect to the 5100~\AA\ continuum are obtained for each season and their mean time lags over the five seasons are 0.69, 4.66, 4.60, 8.43 days, respectively. The H$\gamma$ and H$\beta$ velocity-resolved lag profiles in the seasons of 2015, 2018, 2019, and 2021 are constructed, from which an ``M-shaped'' structure is found in 2015 but disappears after 2018. Our five-season reverberation mapping (RM) yields an averaged virial SMBH mass of $M_\bullet/10^7M_\odot=14.22$, with a small standard deviation of $1.89$. By combining the previous 18 RM campaigns and our five-season campaign for NGC~5548, we find that there exists a time lag of 3.5~years between the changes in the BLR size and optical luminosity. In addition, we also construct the BLR radius$-$luminosity relation and the virial relation for NGC~5548.

S.E. Anderson, J.-M. Petit, B. Noyelles, O. Mousis, P. Rousselot

7 pages, 5 figures

Comet C/2016 R2 PanSTARRS (hereafter C/2016 R2) presents an unusually high N2/CO abundance ratio, as well as a heavy depletion in H2O, making it the only known comet of its kind. Understanding its dynamical history is therefore of essential importance as it would allow us to gain a clearer understanding of the evolution of planetesimal formation in our Solar System. Two studies have independently estimated the possible origin of this comet from building blocks formed in a peculiar region of the protoplanetary disk, near the ice line of CO and N2. We intend to investigate the fates of objects formed from the building blocks in these regions. We hope to find a possible explanation for the lack of C/2016 R2-like comets in our Solar System. Using a numerical simulation of the early stages of Solar System formation, we track the dynamics of these objects in the Jumping Neptune scenario based on five different initial conditions for the protosolar disk. We integrate the positions of 250 000 planetesimals over time in order to analyze the evolution of their orbits and create a statistical profile of their expected permanent orbit. Results. We find that objects formed in the region of the CO- and N2- ice lines are highly likely to be sent towards the Oort Cloud or possibly ejected from the Solar System altogether on a relatively short timescale. In all our simulations, over 90% of clones formed in this region evolved into a hyperbolic trajectory, and between 1% and 10% were potentially captured by the Oort Cloud. The handful of comets that remained were either on long-period, highly eccentric orbits like C/2016 R2, or absorbed into the Edgeworth-Kuiper belt. Comets formed <15 au were predominantly ejected early in the formation timeline. As this is the formation zone likely to produce comets of this composition, this process could explain the lack of similar comets observed in the Solar System

J. Lezmy, Y. Copin, M. Rigault, M. Smith, J. D. Neill

14 pages, 17 figures, submitted to Astronomy & Astrophysics

Recent developments in time domain astronomy, like the Zwicky Transient Facility, have made possible a daily scan of the entire visible sky, leading to the discovery of hundreds of new transients every night. Among them, 10 to 15 are supernovae (SNe), which have to be classified prior to cosmological use. The Spectral Energy Distribution machine (SEDm), a low resolution Integral Field Spectrograph, has been designed, built, and operated to spectroscopically classify targets detected by the ZTF main camera. The current Pysedm pipeline is limited by contamination when the transient is too close to its host galaxy core; this can lead to an incorrect typing and ultimately bias the cosmological analyses, and affect the SN sample homogeneity in terms of local environment properties. We present a new scene modeler to extract the transient spectrum from its structured background, aiming at improving the typing efficiency of the SEDm. HyperGal is a fully chromatic scene modeler, which uses pre-transient photometric images to generate a hyperspectral model of the host galaxy; it is based on the CIGALE SED fitter used as a physically-motivated spectral interpolator. The galaxy model, complemented by a point source and a diffuse background component, is projected onto the SEDm spectro-spatial observation space and adjusted to observations. The full procedure is validated on 5000 simulated cubes. We introduce the contrast as the transient-to-total flux ratio at SN location. From estimated contrast distribution of real SEDm observations, we show that HyperGal correctly classifies ~95% of SNe Ia. Compared to the standard extraction method, HyperGal correctly classifies 10% more SNe Ia. The false positive rate is less than 2%, half as much as the standard extraction method. Assuming a similar contrast distribution for core-collapse SNe, HyperGal classifies 14% (11%) additional SNe II (Ibc).

Edoardo Giangrandi, Violetta Sagun, Oleksii Ivanytskyi, Constança Providência, Tim Dietrich

14 pages, 5 figures

We propose a model of asymmetric bosonic dark matter (DM) with self-repulsion mediated by the vector field coupled to the complex scalar particles. By adopting the two-fluid formalism, we study different DM distribution regimes, either, fully condensed inside the core of a star or, otherwise, distributed in a dilute halo around a neutron star (NS). We show that DM condensed in a core leads to a decrease of the total gravitational mass, radius and tidal deformability compared to a pure baryonic star with the same central density, which we will perceive as an effective softening of the equation of state (EoS). On the other hand, the presence of a DM halo increases the tidal deformability and total gravitational mass. As a result, an accumulated DM inside compact stars could mimic an apparent stiffening of strongly interacting matter equation of state and constraints we impose on it at high densities. From the performed analysis of the effect of DM particles in a MeV-GeV mass-scale, interaction strength, and relative DM fractions inside NSs we obtained a rigorous constraint on model parameters. Finally, we discuss several smoking guns of the presence of DM that are free from the above mentioned apparent modification of the strongly interacting matter equation of state. With this we could be probed with the future astrophysical and gravitational wave (GW) surveys.

A recent suggestion that acetamide, \ce{CH3C(O)NH2}, could be readily formed on water-ice grains by the acid induced addition of water across the \ce{CN} bond is now shown to be valid. Computational modelling of the reaction between \ce{R-CN} (R = H, \ce{CH3}) and a cluster of 32 molecules of water and one \ce{H3O+} proceeds auto-catalytically to form firstly a hydroxy imine \ce{R-C(OH)=NH} and secondly an amide \ce{R-C(O)NH2}. Quantum mechanical tunnelling, computed from small-curvature estimates, plays a key role in the rates of these reactions. This work represents the first credible effort to show how amides can be formed from abundant substrates, namely nitriles and water, reacting on a water-ice cluster containing catalytic amounts of hydrons in the interstellar medium with consequential implications towards the origins of life.

Dong Li, Fanpeng Shi, Haisheng Zhao, Shaolin Xiong, Liming Song, Wenxi Peng, Xinqiao Li, Wei Chen, Zongjun Ning

27 pages, accepted by Frontiers in Frontiers in Astronomy and Space Sciences

Quasi-periodic pulsations (QPPs), which carry time features and plasma characteristics of flare emissions, are frequently observed in light curves of solar/stellar flares. In this paper, we investigated non-stationary QPPs associated with recurrent jets during an M1.2 flare on 2022 July 14. A quasi-period of about 45$\pm$10 s, determined by the wavelet transform technique, is simultaneously identified at wavelengths of soft/hard X-ray and microwave emissions, which are recorded by the Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor, Fermi, and the Nobeyama Radio Polarimeters, respectively. A group of recurrent jets with an intermittent cadence of about 45$\pm$10 s are found in Atmospheric Imaging Assembly (AIA) image series at 304 {\AA}, but they are 180-s earlier than the flare QPP. All observational facts suggest that the flare QPP could be excited by recurrent jets, and they should be associated with nonthermal electrons that are periodically accelerated by a repeated energy release process, like repetitive magnetic reconnection. Moreover, the same quasi-period is discovered at double footpoints connected by a hot flare loop in AIA 94 {\AA}, and the phase speed is measured to 1420 km/s. Based on the differential emission measure, the average temperatures, number densities, and magnetic field strengths at the loop top and footpoint are estimated to 7.7/6.7 MK, 7.5/3.6*10^{10} cm ^{-3}, and 143/99 G, respectively. Our measurements indicate that the 45-s QPP is probably modulated by the kink-mode wave of the flare loop.

Here we present the spatially resolved study of six Galactic planetary nebulae (PNe), namely IC 4593, Hen 2-186, Hen 2-429, NGC 3918, NGC 6543 and NGC 6905, from intermediate-resolution spectra of the 2.5 m Isaac Newton Telescope and the 1.54 m Danish telescope. The physical conditions (electron densities, N$_{e}$, and temperatures, T$_{e}$), chemical compositions and dominant excitation mechanisms for the different regions of these objects are derived, in an attempt to go deeper on the knowledge of the low-ionization structures (LISs) hosted by these PNe. We reinforce the previous conclusions that LISs are characterized by lower (or at most equal) N$_{e}$ than their associated rims and shells. As for the T$_{e}$, we point out a \textit{possible} different trend between the N and O diagnostics. T$_e$[NII] does not show significant variations throughout the nebular components, whereas T$_e$[OIII] appears to be slightly higher for LISs. The much larger uncertainties associated with the T$_e$[OIII] of LISs do not allow robust conclusions. Moreover, the chemical abundances show no variation from one to another PN components, not even contrasting LISs with rims and shells, as also found in a number of other works. By discussing the ionization photon flux due to shocks and stellar radiation, we explore the possible mechanisms responsible for the excitation of LISs. We argue that the presence of shocks in LISs is not negligible, although there is a strong dependence on the orientation of the host PNe and LISs.

Liliya Imasheva (1, 2), H.-Thomas Janka (1, 3), Achim Weiss (1, 2) ((1) MPI Astrophysics, Garching, (2) LMU, Munich, (3) TUM, Garching)

24 pages, 12 figures

Thermal bombs are a widely used method to artificially trigger explosions of core-collapse supernovae (CCSNe) to determine their nucleosynthesis or ejecta and remnant properties. Recently, their use in spherically symmetric (1D) hydrodynamic simulations led to the result that {56,57}Ni and 44Ti are massively underproduced compared to observational estimates for Supernova 1987A, if the explosions are slow, i.e., if the explosion mechanism of CCSNe releases the explosion energy on long timescales. It was concluded that rapid explosions are required to match observed abundances, i.e., the explosion mechanism must provide the CCSN energy nearly instantaneously on timescales of some ten to order 100 ms. This result, if valid, would disfavor the neutrino-heating mechanism, which releases the CCSN energy on timescales of seconds. Here, we demonstrate by 1D hydrodynamic simulations and nucleosynthetic post-processing that these conclusions are a consequence of disregarding the initial collapse of the stellar core in the thermal-bomb modelling before the bomb releases the explosion energy. We demonstrate that the anti-correlation of 56Ni yield and energy-injection timescale vanishes when the initial collapse is included and that it can even be reversed, i.e., more 56Ni is made by slower explosions, when the collapse proceeds to small radii similar to those where neutrino heating takes place in CCSNe. We also show that the 56Ni production in thermal-bomb explosions is sensitive to the chosen mass cut and that a fixed mass layer or fixed volume for the energy deposition cause only secondary differences. Moreover, we propose a most appropriate setup for thermal bombs.

Rhys Taylor, Joachim Koppen, Pavel Jachym, Robert Minchin, Jan Palous, Jessica Rosenberg, Steven Schneider, Richard Wunsch, Boris Deshev

16 pages (plus 4 page appendix), 1 table, 7 figures, accepted for publication in AJ

Using data from the Arecibo Galaxy Environment Survey, we report the discovery of five HI clouds in the Leo I group without detected optical counterparts. Three of the clouds are found midway between M96 and M95, one is only 10$^{\prime}$ from the south-east side of the well-known Leo Ring, and the fifth is relatively isolated. HI masses range from 2.6$\times$10$^{6}$ - 9.0$\times$10$^{6}$M$_{\odot}$ and velocity widths (W50) from 16 - 42 km/s. Although a tidal origin is the most obvious explanation, this formation mechanism faces several challenges. For the most isolated cloud, the difficulties are its distance from neighbouring galaxies and the lack of any signs of disturbance in the HI discs of those systems. Some of the clouds also appear to follow the baryonic Tully-Fisher relation between mass and velocity width for normal, stable galaxies which is not expected if they are tidal in origin. Three clouds are found between M96 and M95 which have no optical counterparts, but have otherwise similar properties and location to the optically detected galaxy LeG 13. While overall we favour a tidal debris scenario to explain the clouds, we cannot rule out a primordial origin. If the clouds were produced in the same event that gave rise to the Leo Ring, they may provide important constraints on any model attempting to explain that structure

Jian-He Zheng, Xiang-Yu Wang, Ruo-Yu Liu

9 pages, 4 figures, submitted, comments are welcome

Three tidal disruption event (TDE) candidates (AT2019dsg, AT2019fdr, AT2019aalc) have been found to be coincident with high energy astrophysical neutrinos in multi-messenger follow-ups. Recent studies suggest the presence of a quasi-spherical, optically thick envelope around the supermassive black holes in TDEs, resulted from stellar debris after the disruption. We study whether the neutrino signal can be explained by choked relativistic jets inside the envelope. While powerful jets, such as that in Swift J1644+57, can successfully break out the envelope, those with relatively weak power could be choked by the envelope. Choked jets can still accelerate cosmic rays through internal shocks or reverse shocks deep in the envelope, which further produce high-energy neutrinos via interaction with the thermal photons in the envelope. We explore the parameter space of the jets that can produce detectable neutrino flux while being choked. Under reasonable assumption about the envelope mass, we find that the cumulative neutrino numbers of three TDEs are consistent with the expected range imposed by observations. Compared with other proposed models, the relativistic bulk motion of the jets in our model can magnify the neutrino flux by Lorentz boosting. The neutrino time delay relative to the optical peak time of TDEs can be explained as the jet propagation time in the envelope before being choked. The discovery of TDE-associated neutrino events may suggest that jets might have been commonly formed in TDEs, as expected from super-Eddington accretion, but most of them are too weak to break out of the envelopes.

P. Molaro, L. Izzo, P. Selvelli, P. Bonifacio, E. Aydi, G. Cescutti, E. Guido, E. J. Harvey, M. Hernanz, M. Della Valle

Accepted 2022 September 17. Received 2022 September 16; in original form 2022 April 20. 12 pages, 13 figures

The recurrent nova RS Oph underwent a new outburst on August 8, 2021, reaching a visible brightness of V = 4.8 mag. Observations of the 2021 outburst made with the high resolution UVES spectrograph at the Kueyen-UT2 telescope of ESO-VLT in Paranal enabled detection of the possible presence of 7Be freshly made in the thermonuclear runaway reactions. The 7Be yields can be estimated in N(Be)/N(H) = 5.7 x 10^(-6), which are close to the lowest yields measured in classical novae so far. 7Be is short-lived and decays only into 7Li. By means of a spectrum taken during the nebular phase we estimated an ejected mass of about 1.1 x 10^(-5) Msun, providing an amount of about 4.4x 10^(-10) Msun of 7Li created in the 2021 event. Recurrent novae of the kind of RS Oph may synthesize slightly lower amount of 7Li per event as classical novae, but occur 10^3 times more frequently. The recurrent novae fraction is in the range of 10-30% and they could have contributed to the making of 7Li we observe today. The detection of 7Be in RS Oph provides further support to the recent suggestion that novae are the most effective source of 7Li in the Galaxy.

L. Fronte (1), B. Mazzon (1), F. Metruccio (1), N. Munaretto (1), M. Doro (2, 3), P. Giommi (4, 5), I. Viale (2, 3), U. Barres de Almeida (6). (1) Liceo Scientifico Statale U. Morin, via Asseggiano 39, I-30174, Venezia, Italy (2) University of Padova, Dep. of Physics, Astronomy, via Marzolo 8, I-35131, Padova, Italy (3) INFN sez. Padova, via Marzolo 8, I-35131, Padova, Italy (4) Center for Astro, Particle, Planetary Physics (CAP3), New York University Abu Dhabi, PO Box 129188 Abu Dhabi, United Arab Emirates (5) Institute for Advanced Study, Technische Universitat Muenchen, Lichtenbergstrasse 2a, D-85748 Garching bei Muenchen, Germany, (6) Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brasil

Proceedings of the 12th Cosmic Ray International Symposium (CRIS 2022), 12-16 September 2022, Naples (Italy). Send correspondence to: michele.doro@unipd.it, giommipaolo@gmail.com

Blazars are active galactic nuclei whose ultra-relativistic jets are coaligned with the observer direction. They emit throughout the whole e.m. spectrum, from radio waves to VHE gamma rays. Not all blazars are discovered. In this work, we propose a catalog of 54 new candidates based on the association of HE gamma ray emission and radio, X-ray an optical signatures. The relevance of this work is also that it was performed by four high school students from the Liceo Scientifico Statale Ugo Morin in Venice, Italy using the open-source platform Open Universe, in collaboration with the University of Padova. The framework of the activity is the Italian MIUR PCTO programme. The success of this citizen-science experience and results are hereafter reported and discussed.

Harry Johnston, Dana Sophia Westbeek, Sjoerd Weide, Nora Elisa Chisari, Yohan Dubois, Julien Devriendt, Christophe Pichon

23 pages, 10 figures (2 pages, 4 figures in appendix), submitted to MNRAS

Residuals between measured galactic radii and those predicted by the Fundamental Plane (FP) are possible tracers of weak lensing magnification. However, observations have shown these to be systematically correlated with the large-scale structure. We use the Horizon-AGN hydrodynamical cosmological simulation to analyse these intrinsic size correlations (ISCs) for both elliptical (early-type) and spiral (late-type) galaxies at $z=0.06$. We fit separate FPs to each sample, finding similarly distributed radius residuals, $\lambda$, in each case. We find persistent $\lambda\lambda$ correlations over three-dimensional separations $0.5-17\,h^{-1}{\rm{Mpc}}$ in the case of spiral galaxies, at $>3\sigma$ significance. When relaxing a mass-selection, applied for better agreement with galaxy clustering constraints, the spiral $\lambda\lambda$ detection strengthens to $9\sigma$; we detect a $5\sigma$ density-$\lambda$ correlation; and we observe intrinsically-large spirals to cluster more strongly than small spirals over scales $\lesssim10\,h^{-1}{\rm{Mpc}}$, at $>5\sigma$ significance. Conversely, and in agreement with the literature, we observe lower-mass, intrinsically-small ellipticals to cluster more strongly than their large counterparts over scales $0.5-17\,h^{-1}{\rm{Mpc}}$, at $>5\sigma$ significance. We model $\lambda\lambda$ correlations using a phenomenological non-linear size model, and predict the level of contamination for cosmic convergence analyses. We find the systematic contribution to be of similar order to, or dominant over the cosmological signal. We make a mock measurement of an intrinsic, systematic contribution to the projected surface mass density $\Sigma(r)$ and find statistically significant, low-amplitude, positive (negative) contributions from lower-mass spirals (ellipticals), which may be of concern for large-scale ($\gtrsim\,7\,h^{-1}$ Mpc) measurements.

Johanna Jurcsik, Aron Juhasz

11 pages, 8 figures, accepted for publication in MNRAS

In this paper, we report on the detection of RRab stars with quasi-identical-shape light curves but period differences as large as $0.05-0.21$ d using the Galactic bulge data of the OGLE-IV survey. We have examined stars with shorter periods than the Oosterhoff~I ridge of the bulge. These stars generally have smaller amplitudes and larger Fourier phase-differences than the typical bulge RRab stars have at the same period. Many of these "anomalous" stars have good-quality light curves without any sign of the Blazhko modulation. Examining their Fourier parameters revealed that several of these stars show very similar light-curve to the typical bulge RR Lyrae. We found hundreds of quasi-identical-shape light-curve pairs with different periods between the "anomalous"- and the "normal"-position RRab stars based on the OGLE $I$-band data. The OGLE $V$-band, and the archive VVV and MACHO surveys $K_s$-, $b$- and $r$-band data of these stars were also checked for light-curve-shape similarity. Finally, 149 pairs with identical-shape light curves in each available photometric band were identified. Calculating the physical properties of the variables using empirical formulae, on average, $-0.5$~dex, $-0.13$~mag, 0.67, and 165~K differences between the [Fe/H], $M_V$, $R/R_\odot$, and $T_{\mathrm{eff}}$ values of the members of the pairs are derived, being the short-period stars less metal-poor, fainter, smaller and hotter than the long-period variables. To explain the existence of variables with different physical properties and pulsation periods but with identical-shape light curves is a challenging task for modelling.

Devesh P. Sariya, Ing-Guey Jiang, D. Bisht, R. K. S. Yadav, G. Rangwal

22 pages, 12 figures and 2 tables. Accepted by New Astronomy

In this paper, we have used the Gaia's Data Release-3 (DR3) data to study an intermediate-age open cluster Berkeley 27 (Be 27). A total of 131 most probable cluster members are picked within the cluster's radius based on the membership probability ($> 80\%$). The cluster's radius was estimated as 3.74 arcmin. The mean proper motion (PM) of Be 27 was determined to be ($\mu_{\alpha} cos{\delta}$, $\mu_{\delta}$)= ($-1.076\pm0.008$, $0.152\pm0.007$)~mas~yr$^{-1}$. The blue straggler stars (BSS) of the cluster were found to be located in the central region. Theoretical isochrones of metallicity Z$_{metal}$= 0.008 were compared to the color-magnitude diagram (CMD) of Be 27. As a result, a heliocentric distance of 4.8$\pm$0.2 kpc and log (age) = 9.36$\pm$0.03 were determined for Be 27. The Galactic orbits are derived using the Galactic potential model which demonstrate that Be 27 follows a circular path around the Galactic center. The cluster does not seem to be affected much by the tidal forces from the Galactic thin disk.

Satadru Bag, Wuhyun Sohn, Arman Shafieloo, Kai Liao

20 pages, 8 figures, 1 table

Strong gravitational lensed quasars (QSOs) have emerged as powerful and novel cosmic probes as they can deliver crucial cosmological information, such as a measurement of the Hubble constant, independent of other probes. Although the upcoming LSST survey is expected to discover $10^3-10^4$ lensed QSOs, a large fraction will remain unresolved due to seeing. The stochastic nature of the quasar intrinsic flux makes it challenging to identify lensed ones and measure the time delays using unresolved light curve data only. In this regard, Bag et al (2022) introduced a data-driven technique based on the minimization of the fluctuation in the reconstructed image light curves. In this article, we delve deeper into the mathematical foundation of this approach. We show that the lensing signal in the fluctuation curve is dominated by the auto-correlation function (ACF) of the derivative of the joint light curve. This explains why the fluctuation curve enables the detection of the lensed QSOs only using the joint light curve, without making assumptions about QSO flux variability, nor requiring any additional information. We show that the ACF of the derivative of the joint light curve is more reliable than the ACF of the joint light curve itself because intrinsic quasar flux variability shows significant auto-correlation up to a few hundred days (as they follow a red power spectrum). In addition, we show that the minimization of fluctuation approach provides even better precision and recall as compared to the ACF of the derivative of the joint light curve when the data have significant observational noise.

Marie Karjalainen, Raine Karjalainen, Artie P. Hatzes, Holger Lehmann, Pierre Kervella, Saskia Hekker, Hans Van Winckel, Jakub Überlauer, Michaela Vítková, Marek Skarka, Petr Kabáth, Saskia Prins, Andrew Tkachenko, William D. Cochran, Alain Jorissen

14 pages, 11 figures. Accepted to A&A

Context. Our knowledge of populations and occurrence of planets orbiting evolved intermediate-mass stars is still incomplete. In 2010 we started a planet-search program among 95 giant stars observed by the Kepler mission to increase the sample of giant stars with planets and with reliable estimates of stellar masses and radii. Aims. We present the two systems KIC 3526061 and HD 187878 from our planet-search program for which we could characterise their companions. Methods. We used precise stellar radial velocity measurements taken with four different echelle spectrographs to derive an orbital solution. We used Gaia astrometric measurements to obtain the inclination of the HD 187878 system and Kepler photometric observations to estimate the stellar mass and radius. Results. We report the discovery of a sub-stellar and a stellar companion around two intermediate-mass red giant branch stars. KIC 3526061 b is most likely a brown dwarf with a minimum mass of 18.15 Jupiter masses in a long-period eccentric orbit, with the orbital period 3552 d and orbital eccentricity 0.85. It is the most evolved system found having a sub-stellar companion with such a large eccentricity and wide separation. HD 187878 B has a minimum mass of 78.4 Jupiter masses. Combining the spectroscopic orbital parameters with the astrometric proper motion anomaly we derived an orbital inclination 9.8 deg, which corresponds to the companion's mass in the stellar regime of 0.51 Sun mass. Conclusions. A sub-stellar companion of KIC 3526061 extends the sample of known red giant branch stars with sub-stellar companions on very eccentric wide orbits and might provide a probe of the dynamical evolution of such systems over time.

We investigate two new observational perspectives in the context of torsional gravitational modification of general relativity, i.e., the $f(T)$ gravity: i) We use Pantheon data of type Ia supernovae motivated by a time variation of the Newton's constant on the supernovae distance modulus relation, and find that a joint analysis with Baryon Acoustic Oscillations (BAO) and Big Bang Nucleosynthesis (BBN), i.e., Pantheon+BAO+BBN, provides constraints on the effective free parameter of the theory to be well compatible with the $\Lambda$CDM prediction; ii) We present the framework of $f(T)$ gravity at the level of linear perturbations with the phenomenological functions, namely the effective gravitational coupling $\mu$ and the light deflection parameter $\Sigma$, which are commonly used to parameterize possible modifications of the Poisson equation relating the matter density contrast to the lensing and the Newtonian potentials, respectively. We use the available Cosmic Microwave Background (CMB) data sets from the Planck 2018 release to constrain the free parameters of the $f(T)$ gravity and $\Lambda$CDM models. We find that CMB data, and its joint analyses with Pantheon and BAO data constrain the $f(T)$ gravity scenario to be practically indistinguishable from the $\Lambda$CDM model. We obtain the strongest limits ever reported on $f(T)$ gravity scenario at the cosmological level.

Christopher Thompson (CITA)

Submitted to MNRAS, 23 pages, 4 figures

The emission of intense radio pulses by flaring magnetars is investigated. Small-scale current gradients can be imprinted into a strongly magnetized outflow by the same processes that source fireball radiation in the closed magnetosphere. This structure arises from a combination of crustal yielding, internal tearing, and turbulent cascade. We consider the quasi-linear development of weak, small-scale currents as (i) they are stretched out and frozen by relativistic expansion and then (ii) pass through a shock. In particular, we derive the amplitudes of the ordinary and fast waves that emerge downstream of a relativistic, magnetized shock in response to a mode that is frozen into the upstream flow (a frozen Alfv\'en wave or entropy wave). An upstream mode with comoving wavelength exceeding the skin depth can linearly convert to a secondary mode propagating above the plasma frequency. A simple and accurate treatment of shocks with extreme magnetization is developed, and the formation of internal shocks in the outflow from a bursting, rotating magnetar is outlined. The emission process described here does not require a strong shock or cool $e^\pm$ pairs (in contrast with the electromagnetic maser shock instability). In some cases, a high-frequency wave is reflected back to the observer, but with a minuscule amplitude that makes it subdominant to other emission channels. The dominant secondary electromagnetic mode is superluminal at emission, is subject to weak induced scattering within the outflow, and can reach the observer in the radio band.

Nikita Rawat, J. C. Pandey, Arti Joshi, Simone Scaringi, Umesh Yadava

Accepted for publication in MNRAS

Based on the X-ray observations from XMM-Newton and Swift, and optical observations from Transiting Exoplanet Survey Satellite (TESS) and AAVSO, we present temporal and spectral properties of probable intermediate polar SWIFT J0503.7-2819. The X-ray light curve shows two distinctive features, where possibly the second pole seems to be active during the middle of the XMM-Newton observations. Present analysis confirms and also refines the previously reported orbital period of SWIFT J0503.7-2819 as 81.65$\pm$0.04 min. The X-ray and optical variations of this target have been found to occur at the period of $\sim$ 65 min, which we propose as the spin period of the white dwarf (WD). The energy-dependent modulation at this period, which are due to the photoelectric absorption in the accretion flow, also assures this conjecture. Two temperature thermal plasma model well explains the X-ray spectra with temperatures of $\sim$ 150 eV and $\sim$ 18.5 keV, which is absorbed by a dense material with an average equivalent hydrogen column density of 3.8 $\times$ 10$^{22}$ cm$^{-2}$ that partially covers $\sim$ 27% of the X-ray source. An attempt is made to understand the accretion flow in this system using the present data of SWIFT J0503.7-2819. If the proposed spin period is indeed the actual period, then SWIFT J0503.7-2819 could be the first nearly synchronous intermediate polar below the period gap.

Marlin B. Schäfer, Ondřej Zelenka, Alexander H. Nitz, He Wang, Shichao Wu, Zong-Kuan Guo, Zhoujian Cao, Zhixiang Ren, Paraskevi Nousi, Nikolaos Stergioulas, Panagiotis Iosif, Alexandra E. Koloniari, Anastasios Tefas, Nikolaos Passalis, Francesco Salemi, Gabriele Vedovato, Sergey Klimenko, Tanmaya Mishra, Bernd Brügmann, Elena Cuoco, E. A. Huerta, Chris Messenger, Frank Ohme

25 pages, 6 figures, 4 tables, additional material available at this https URL

We present the results of the first Machine Learning Gravitational-Wave Search Mock Data Challenge (MLGWSC-1). For this challenge, participating groups had to identify gravitational-wave signals from binary black hole mergers of increasing complexity and duration embedded in progressively more realistic noise. The final of the 4 provided datasets contained real noise from the O3a observing run and signals up to a duration of 20 seconds with the inclusion of precession effects and higher order modes. We present the average sensitivity distance and runtime for the 6 entered algorithms derived from 1 month of test data unknown to the participants prior to submission. Of these, 4 are machine learning algorithms. We find that the best machine learning based algorithms are able to achieve up to 95% of the sensitive distance of matched-filtering based production analyses for simulated Gaussian noise at a false-alarm rate (FAR) of one per month. In contrast, for real noise, the leading machine learning search achieved 70%. For higher FARs the differences in sensitive distance shrink to the point where select machine learning submissions outperform traditional search algorithms at FARs $\geq 200$ per month on some datasets. Our results show that current machine learning search algorithms may already be sensitive enough in limited parameter regions to be useful for some production settings. To improve the state-of-the-art, machine learning algorithms need to reduce the false-alarm rates at which they are capable of detecting signals and extend their validity to regions of parameter space where modeled searches are computationally expensive to run. Based on our findings we compile a list of research areas that we believe are the most important to elevate machine learning searches to an invaluable tool in gravitational-wave signal detection.

J. Borovicka, P. Spurny, L. Shrbeny, R. Stork, L. Kotkova, J. Fuchs, J. Keclikova, H. Zichova, J. Manek, P. Vachova, I. Macourkova, J. Svoren, H. Mucke

accepted in Astronomy and Astrophysics

A catalog of 824 fireballs (bright meteors), observed by a dedicated network of all-sky digital photographic cameras in central Europe in the years 2017-2018 is presented. The status of the European Fireball Network, established in 1963, is described. The cameras collect digital images of meteors brighter than an absolute magnitude of about -2 and radiometric light curves with a high temporal resolution of those brighter than a magnitude ~ -4. All meteoroids larger than 5 grams, corresponding to sizes of about 2 cm, are detected regardless of their entry velocity. High-velocity meteoroids are detected down to masses of about 0.1 gram. The largest observed meteoroid in the reported period 2017-2018 had a mass of about 100 kg and a size of about 40 cm. The methods of data analysis are explained and all catalog entries are described in detail. The provided data include the fireball date and time, atmospheric trajectory and velocity, the radiant in various coordinate systems, heliocentric orbital elements, maximum brightness, radiated energy, initial and terminal masses, maximum encountered dynamic pressure, physical classification, and possible shower membership. Basic information on the fireball spectrum is available for some bright fireballs (apparent magnitude < -7). A simple statistical evaluation of the whole sample is provided. The scientific analysis is presented in an accompanying paper.

Ahmed Faris Al-Refaie, Olivia Venot, Quentin Changeat, Billy Edwards

13 pages, 8 figures

We introduce a new chemical kinetic code FRECKLL (Full and Reduced Exoplanet Chemical Kinetics distiLLed) to evolve large chemical networks efficiently. FRECKLL employs `distillation' in computing the reaction rates, which minimizes the error bounds to the minimum allowed by double precision values ($\epsilon \leq 10^{-15}$). FRECKLL requires less than 5 minutes to evolve the full Venot2020 network in a 130 layers atmosphere and 30 seconds to evolve the Venot2020 reduced scheme. Packaged with FRECKLL is a TauREx 3.1 plugin for usage in forward modelling and retrievals. We present TauREx retrievals performed on a simulated HD189733 JWST spectra using the full and reduced Venot2020 chemical networks and demonstrate the viability of total disequilibrium chemistry retrievals and the ability for JWST to detect disequilibrium processes.

Debanjan Sengupta, Orkan M. Umurhan

55 pages, 27 figures, accepted for publication in ApJ

We examine the settled particle layers of planet forming disks in which the streaming instability (SI) is thought to be either weak or inactive. A suite of low-to-moderate resolution three-dimensional simulations in a $0.2H$ sized box, where $H$ is the pressure scale height, are performed using PENCIL for two Stokes numbers, \St$=0.04$ and $0.2$, at 1\% disk metallicity. We find a complex of Ekman-layer jet-flows emerge subject to three co-acting linearly growing processes: (1) the Kelvin-Helmholtz instability (KHI), (2) the planet-forming disk analog of the baroclinic Symmetric Instability (SymI), and (3) a later-time weakly acting secondary transition process, possibly a manifestation of the SI, producing a radially propagating pattern state. For \St$=0.2$, KHI is dominant and manifests as off-midplane axisymmetric rolls, while for \St$=0.04$ the axisymmetric SymI mainly drives turbulence. SymI is analytically developed in a model disk flow, predicting that it becomes strongly active when the Richardson number (Ri) of the particle-gas midplane layer transitions below 1, exhibiting growth rates $\le\sqrt{2/\Ri - 2}\cdot\Omega$, where $\Omega$ is local disk rotation rate. For fairly general situations absent external sources of turbulence it is conjectured that the SI, when and if initiated, emerges out of a turbulent state primarily driven and shaped by at least SymI and/or KHI. We also find that turbulence produced in $256^3$ resolution simulations are not statistically converged and that corresponding $512^3$ simulations may be converged for \St$=0.2$. Furthermore, we report that our numerical simulations significantly dissipate turbulent kinetic energy on scales less than 6-8 grid points.

Tim D. Pearce, Florian Kirchschlager, Gaël Rouillé, Steve Ertel, Alexander Bensberg, Alexander V. Krivov, Mark Booth, Sebastian Wolf, Jean-Charles Augereau

18 pages, 9 figures, accepted for publication in MNRAS

Excess near-infrared emission is detected around one fifth of main-sequence stars, but its nature is a mystery. These excesses are interpreted as thermal emission from populations of small, hot dust very close to their stars (`hot exozodis'), but such grains should rapidly sublimate or be blown out of the system. To date, no model has fully explained this phenomenon. One mechanism commonly suggested in the literature is cometary supply, where star-grazing comets deposit dust close to the star, replenishing losses from grain sublimation and blowout. However, we show that this mechanism alone is very unlikely to be responsible for hot exozodis. We model the trajectory and size evolution of dust grains released by star-grazing comets, to establish the dust and comet properties required to reproduce hot-exozodi observations. We find that cometary supply alone can only reproduce observations if dust ejecta has an extremely steep size distribution upon release, and the dust-deposition rate is extraordinarily high. These requirements strongly contradict our current understanding of cometary dust and planetary systems. Cometary supply is therefore unlikely to be solely responsible for hot exozodis, so may need to be combined with some dust-trapping mechanism (such as gas or magnetic trapping) if it is to reproduce observations.

H. Adhikary, K.K. Allison, N. Amin, E.V. Andronov, T. Antićić, I.-C. Arsene, Y. Balkova, M. Baszczyk, D. Battaglia, S. Bhosale, A. Blondel, M. Bogomilov, Y. Bondar, N. Bostan, A. Brandin, A. Bravar, W. Bryliński, J. Brzychczyk, M. Buryakov, M. Ćirković, M. Csanad, J. Cybowska, T. Czopowicz, A. Damyanova, N. Davis, H. Dembinski, A. Dmitriev, W. Dominik, P. Dorosz, J. Dumarchez, R. Engel, G.A. Feofilov, L. Fields, Z. Fodor, M. Friend, A. Garibov, M. Gaździcki, O. Golosov, V. Golovatyuk, M. Golubeva, K. Grebieszkow, F. Guber, A. Haesler, M. Haug, S.N. Igolkin, S. Ilieva, A. Ivashkin, A. Izvestnyy, S.R. Johnson, K. Kadija, N. Kargin, N. Karpushkin, E. Kashirin, M. Kiełbowicz, V.A. Kireyeu, H. Kitagawa, R. Kolesnikov, D. Kolev, et al. (99 additional authors not shown)

We present a measurement of the momentum spectra of $\pi^\pm$, K$^\pm$, p$^\pm$, $\Lambda$, $\bar{\Lambda}$ and K$^{0}_{S}$ produced in interactions of negatively charged pions with carbon nuclei at beam momenta of 158 and 350 GeV/c. The total production cross sections are measured as well. The data were collected with the large-acceptance spectrometer of the fixed target experiment NA61/SHINE at the CERN SPS. The obtained double-differential $p$-$p_T$ spectra provide a unique reference data set with unprecedented precision and large phase-space coverage to tune models used for the simulation of particle production in extensive air showers in which pions are the most numerous projectiles.

This work includes two new results - principally two new constants of motion for the linearised restricted 3-body problem (e.g. for the Trojan asteroids) and an important isosceles triangle generalisation of Lagrange's equilateral triangle solution of the restricted case leading to hidden constants for Hildans as well as Trojans. Both of these results are classical, but we also have included new results on Newtonian quantum gravity emanating from the asymptotics relevant for WIMPish particles, explaining the origin of systems like that of the Trojans. The latter result uses a generalisation of our semi-classical mechanics for Schr\"odinger equations involving vector as well as scalar potentials, presented here for the first time, thereby providing an acid test of our ideas in predicting the quantum curvature and torsion of WIMPish trajectories for our astronomical elliptic states. The combined effect is to give a new celestial mechanics for WIMPs in gravitational systems as well as new results for classical problems. As we shall explain, we believe these results could help to see how spiral galaxies evolve into elliptical ones. A simple classical consequence of our isosceles triangle result gives a Keplerian type $4^{\textrm{th}}$ Law for 3-body problems. This is confined to the Appendix.

Kinematically forbidden channels can set the freeze-out dark matter (DM) relic abundance. These channels are described by DM annihilations into heavier states, which vanish at zero temperature limit, but occur at finite temperatures in the early Universe. For the case that the final state of the forbidden channel is scalar mediators that couple to Standard Model (SM) matter through mixing with the SM Higgs, the signals from DM-nucleon interactions and from mediator-related missing energy or displaced vertices could be detected by direct detections and particle physics experiments, respectively. We thus present a study on the simplest secluded vector dark matter model that can exhibit this scenario in the mass range from sub-GeV to TeV. The dark matter resides in the hidden sector, which is in thermal equilibrium with the SM before freeze-out. During freeze-out, the depletion of its density results from its annihilation into two heavier but metastable scalars, where the coupling can be determined by having the correct relic density and constrained by the perturbative unitarity bound. However, much of the allowed parameter space is insensitive to the mixing angle between the hidden scalar and SM Higgs. We find that a more significant mass splitting between DM and the mediator can be allowed only in the sub-GeV region. This model of the forbidden DM interacting with SM particles through the scalar portal is testable in experiments.

Quentin Baghi, John G. Baker, Jacob Slutsky, James Ira Thorpe

11 pages, 7 figures

We recently introduced the basic concepts of an approach to filtering strongly laser-noise dominated space-based gravitational-wave data, like LISA's phase comparison data streams, which does not rely on independent knowledge of a temporal delays pattern in the dominant noise that generates the data. Instead, our automated Principal Component Interferometry (aPCI) approach only assumes that one can produce some linear combinations of the temporally nearby regularly spaced phase measurements, which cancel the laser noise. Then we let the data reveal those combinations, thus providing us with a set of laser-noise-free data channels. Our basic approach relied on the simplifying additional assumption that laser-noise-cancelling data combinations or the filters which lead to the laser-noise-free data streams are time-independent. In LISA, however, these filters will vary as the constellation armlengths evolve. Here, we discuss a generalization of the basic aPCI concept compatible with data dominated by a still unmodeled but slowly varying dominant noise covariance. We find that despite its independence on any model, the aPCI processing successfully mitigates laser frequency noise below the other noise sources level, and that its sensitivity to gravitational waves is the same as the state-of-the-art second-generation time-delay interferometry, up to a 2% error.

Dark matter (DM)-electron scattering is a prime target of a number of direct DM detection experiments and constitutes a promising avenue for exploring interactions of DM in the sub-GeV mass-range, challenging to probe with nuclear recoils. We extend the recently proposed halo-independent analysis method for DM-electron scattering, which allows to infer the local DM halo properties without any additional assumptions about them, to include in-medium effects through dielectric functions of the target material. We show that in-medium effects could significantly affect halo-independent analysis response functions for germanium and silicon and thus are essential for proper inference of local DM halo characteristics from direct DM detection data.

We study production of dark relics during reheating after the end of inflation in a system consisting of a non-minimally coupled spectator scalar field and the inflaton. We derive a set of renormalized quantum transport equations for the one-point function and the two-point function of the spectator field and solve them numerically. We find that our system can embody both tachyonic and parametric instabilities. The former is an expected result due to the non-minimal coupling, but the latter displays new features driven by a novel interplay of the two-point function with the Ricci scalar. We find that when the parametric instability driven by the two-point function takes place, it dominates the total particle production. The quantitative results are also found to be highly sensitive to the model parameters.

We discuss strange stars admixed with fermionic dark matter in the presence of a strong magnetic field using the two-fluid Tolman-Oppenheimer-Volkov equations. We describe strange quark matter using the MIT bag model and consider magnetic fields in the range $\sim 10^{17}-10^{18}$ G. For the fermionic dark matter, we consider the cases of free particles and strongly self-interacting particles, with dark fermion masses $m=5, 100, 500$ GeV. Even though strong magnetic fields contribute to decreasing the total mass of the star, they attenuate the rate of decrease in the maximum mass brought about by increasing the dark matter fraction in the admixed system.

Eduardo S. Fraga, Rodrigo da Mata, Savvas Pitsinigkos, Andreas Schmitt

6 pages, 4 figures, contribution to proceedings of QCD@Work, 27-30 June 2022

We review the results of a phenomenological model for cold and dense nuclear matter exhibiting a chiral phase transition. The idea is to model the quark-hadron phase transition under neutron star conditions within a single model, but without adding quark degrees of freedom by hand. To this end, strangeness is included in the form of hyperonic degrees of freedom, whose light counterparts provide the strangeness in the chirally restored phase. In the future, the model can be used for instance to compute the surface tension at the (first-order) chiral phase transition and to study the possible existence of inhomogeneous phases.

It was shown a few years back that for the Brans-Dicke theory with a positive cosmological constant $\Lambda$, and a de Sitter or cosmological event horizon in the asymptotic region, not only there exist no non-trivial field configurations, but also the inverse Brans-Dicke parameter $\omega^{-1}$ must be vanishing, thereby essentially reducing the theory to Einstein's General Relativity. The assumption of the existence of the cosmological horizon was crucial for this proof. However, since the Brans-Dicke field $\phi$, couples directly to the $\Lambda$-term in the energy-momentum tensor as well as in its equation of motion, perhaps it is reasonable to ask : can $\phi$ become strong instead and screen the effect of $\Lambda$ at very large scales, so that the asymptotic de Sitter structure is replaced by some physically acceptable alternative non-singular boundary condition? In this work we analytically argue that under the assumption of any generic asymptotic stationary spacetime structure in the absence of a cosmological event horizon, similar non-existence results hold, as long as the spacetime is free of any naked curvature singularity. We further support this result by providing explicit numerical computations. Thus we conclude that in the presence of a positive $\Lambda$, and for any generic asymptotic spacetime structure free from curvature singularity, a stationary black hole or even a star solution in the Brans-Dicke theory always necessitates $\omega^{-1}=0$, and thereby reducing the theory to General Relativity.

Minxi He, Yusuke Mikura, Yuichiro Tada

17 pages, 1 figure

We propose an extension of the Higgs inflation to the hybrid metric-Palatini gravity, where we introduce non-minimal couplings between Higgs and both the metric-type and the Palatini-type Ricci scalars. We study the inflationary phenomenology of our model and find that slow-roll inflation can be realized in the large-field regime, giving the observationally favored predictions. In particular, the scalar spectral index exhibits an attractor behavior to $n_{\mathrm{s}}\sim 0.964$, while the tensor-to-scalar ratio can take an arbitrary value depending on the non-minimal coupling parameters, with the metric-Higgs limit $r\sim10^{-3}$ being the maximum. We also investigate the unitarity property of our model. As the ultraviolet (UV) cutoff as a low-energy effective field theory (EFT) of this model is significantly lower than the Planck scale due to a strong curvature of field-space, we consider a possible candidate of UV-extended theories with an additional scalar field introduced so as to flatten the field-space in five-dimension. While the field-space can be flatten completely and this approach can lead to a weakly-coupled EFT, we gain an implication that Planck-scale EFT can be only realized in the limit of metric-Higgs inflation. We also discuss generalizations of the model up to mass-dimension four.

Florent Leclercq

6 pages, 2 figures. Accepted for publication as proceedings of MaxEnt'22 (18-22 July 2022, IHP, Paris, France, this https URL ). The pySELFI code is publicly available at this http URL and on GitHub ( this https URL )

This paper presents recent methodological advances to perform simulation-based inference (SBI) of a general class of Bayesian hierarchical models (BHMs), while checking for model misspecification. Our approach is based on a two-step framework. First, the latent function that appears as second layer of the BHM is inferred and used to diagnose possible model misspecification. Second, target parameters of the trusted model are inferred via SBI. Simulations used in the first step are recycled for score compression, which is necessary to the second step. As a proof of concept, we apply our framework to a prey-predator model built upon the Lotka-Volterra equations and involving complex observational processes.

Carson Ezell, Alexandre Lazarian, Abraham Loeb

10 pages, 2 figures; submitted for publication

The risk of a catastrophic or existential disaster for our civilization is increasing this century. A significant motivation for near-term space settlement is the opportunity to safeguard civilization in the event of a planetary-scale disaster. While a catastrophic event could destroy significant cultural, scientific, and technological progress on earth, early space settlements can maintain a backup data storage system of human activity, including the events leading to the disaster. The system would improve the ability of early space settlers to recover our civilization after collapse. We show that advances in laser communications and data storage enable the development of a data storage system on the lunar surface with a sufficient uplink data rate and storage capacity to preserve valuable information about the achievements of our civilization and the chronology of the disaster.

Soft x-ray radiation from the sun is responsible for the production of high energy photoelectrons in the D and E regions of the ionosphere, where they deposit most of their ionization energy. The photoelectrons created by this process are the main drivers for dissociation of Nitrogen molecule ($N_2$) below 200 km. Furthermore, the dissociation of $N_2$ is one of main mechanisms of the production of Nitric Oxide (NO) at these altitudes. In order to estimate the dissociation rate of $N_2$ we need its dissociation cross-sections. The dissociation cross-sections for $N_2$ by photoelectrons are primarily estimated from the cross-sections of its excitation states using predissociation factors and dissociative ionization channels. Unfortunately, the lack of cross-sections data, particularly at high electron energies and of higher excited states of $N_2$ and $N_2^+$, introduces a lot of uncertainty in the dissociation rate calculation, which subsequently leads to uncertainties in the NO production rate from this source.