Abstract visibility
Abstract text size

Papers for Wednesday, Sep 20 2023

Papers with local authors

Alisa Galishnikova, Alexander Philippov, Eliot Quataert

31 pages, 17 figures, 1 table. Accepted for publication in ApJ

0 votes
Paper 2 — arXiv:2309.10029
0 votes
Paper 2 — arXiv:2309.10029

Low-collisionality plasma in a magnetic field generically develops anisotropy in its distribution function with respect to the magnetic field direction. Motivated by the application to radiation from accretion flows and jets, we explore the effect of temperature anisotropy on synchrotron emission. We derive analytically and provide numerical fits for the polarized synchrotron emission and absorption coefficients for a relativistic bi-Maxwellian plasma (we do not consider Faraday conversion/rotation). Temperature anisotropy can significantly change how the synchrotron emission and absorption coefficients depend on observing angle with respect to the magnetic field. The emitted linear polarization fraction does not depend strongly on anisotropy, while the emitted circular polarization does. We apply our results to black hole imaging of Sgr A* and M87* by ray-tracing a GRMHD simulation and assuming that the plasma temperature anisotropy is set by the thresholds of kinetic-scale anisotropy-driven instabilities. We find that the azimuthal asymmetry of the 230 GHz images can change by up to a factor of 3, accentuating ($T_\perp > T_\parallel$) or counteracting ($T_\perp < T_\parallel$) the image asymmetry produced by Doppler beaming. This can change the physical inferences from observations relative to models with an isotropic distribution function, e.g., by allowing for larger inclination between the line of sight and spin direction in Sgr A*. The observed image diameter and the size of the black hole shadow can also vary significantly due to plasma temperature anisotropy. We describe how the anisotropy of the plasma can affect future multi-frequency and photon ring observations. In Appendices we calculate kinetic anisotropy-driven instabilities (mirror, whistler, and firehose) for relativistically hot plasmas.

Kyoung-Soo Lee, Eric Gawiser, Changbom Park, Yujin Yang, Francisco Valdes, Dustin Lang, Vandana Ramakrishnan, Byeongha Moon, Nicole Firestone, Stephen Appleby, Maria Celeste Artale, Moira Andrews, Franz E. Bauer, Barbara Benda, Adam Broussard, Yi-Kuan Chiang, Robin Ciardullo, Arjun Dey, Rameen Farooq, Caryl Gronwall, Lucia Guaita, Yun Huang, Ho Seong Hwang, Sanghyeok Im, Woong-Seob Jeong, Shreya Karthikeyan, Hwihyun Kim, Seongjae Kim, Gautam Nagaraj, Julie Nantais, Nelson Padilla, Jaehong Park, Alexandra Pope, Roxana Popescu, David J. Schlegel, Eunsuk Seo, Akriti Singh, Hyunmi Song, Paulina Troncoso, A. Katherina Vivas, Ann Zabludoff, Alfredo Zenteno

21 pages, 8 figures, submitted to Astrophysical Journal

0 votes
Paper 23 — arXiv:2309.10191
0 votes
Paper 23 — arXiv:2309.10191

We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.2, 7.4, and 9.8 nm, corresponding to Lya at z=2.4, 3.1, and 4.5 and cosmic times of 2.8, 2.1, and 1.4 Gyr, respectively. When combined with even deeper, public broad-band data from Hyper Suprime-Cam, DECam, and in the future, LSST, the ODIN narrow-band images will enable the selection of over 100,000 Lya-emitting (LAE) galaxies at these epochs. ODIN-selected LAEs will identify protoclusters as galaxy overdensities, and the deep narrow-band images enable detection of highly extended Lya blobs (LABs). Primary science goals include measuring the clustering strength and dark matter halo connection of LAEs, LABs, and protoclusters, and their respective relationship to filaments in the cosmic web. The three epochs allow the redshift evolution of these properties to be determined during the period known as Cosmic Noon, where star formation was at its peak. The two narrow-band filter wavelengths are designed to enable interloper rejection and further scientific studies by revealing [O II] and [O III] at z=0.34, Lya and He II 1640 at z=3.1, and Lyman continuum plus Lya at z=4.5. Ancillary science includes similar studies of the lower-redshift emission-line galaxy samples and investigations of nearby star-forming galaxies resolved into numerous [O III] and [S II] emitting regions.

Papers with votes

A. M. Nierenberg, R. E. Keeley, D. Sluse, D. Gilman, S. Birrer, T. Treu, K. N. Abazajian, T. Anguita, A. J. Benson, V. N. Bennert, S. G. Djorgovski, X. Du, C. D. Fassnacht, S. F. Hoenig, A. Kusenko, C. Lemon, M. Malkan, V. Motta, L. A. Moustakas, D. Stern, R. H. Wechsler
1 vote
Paper 14 — arXiv:2309.10101
1 vote
Paper 14 — arXiv:2309.10101

The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$\mu$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models.

Jiří Svoboda, Michal Dovčiak, James F. Steiner, Fabio Muleri, Adam Ingram, Anastasiya Yilmaz, Nicole Rodriguez Cavero, Lorenzo Marra, Juri Poutanen, Alexandra Veledina, Mehrnoosh Rahbardar Mojaver, Stefano Bianchi, Javier Garcia, Philip Kaaret, Henric Krawczynski, Giorgio Matt, Jakub Podgorný, Martin C. Weisskopf, Fabian Kislat, Pierre-Olivier Petrucci, Maimouna Brigitte, Michal Bursa, Sergio Fabiani, Kun Hu, Sohee Chun, Guglielmo Mastroserio, Romana Mikušincová, Ajay Ratheesh, Roger W. Romani, Paolo Soffitta, Francesco Ursini, Silvia Zane, Iván Agudo, Lucio A. Antonelli, Matteo Bachetti, Luca Baldini, Wayne H. Baumgartner, Ronaldo Bellazzini, Stephen D. Bongiorno, Raffaella Bonino, Alessandro Brez, Niccolo Bucciantini, Fiamma Capitanio, Simone Castellano, et al. (68 additional authors not shown)

14 pages, 8 figures, submitted to ApJ

1 vote
Paper 74 — arXiv:2309.10813
1 vote
Paper 74 — arXiv:2309.10813

X-ray polarization is a powerful tool to investigate the geometry of accreting material around black holes, allowing independent measurements of the black hole spin and orientation of the innermost parts of the accretion disk. We perform the X-ray spectro-polarimetric analysis of an X-ray binary system in the Large Magellanic Cloud, LMC X-3, that hosts a stellar-mass black hole, known to be persistently accreting since its discovery. We report the first detection of the X-ray polarization in LMC X-3 with the Imaging X-ray Polarimetry Explorer, and find the average polarization degree of 3.2% +- 0.6% and a constant polarization angle -42 deg +- 6 deg over the 2-8 keV range. Using accompanying spectroscopic observations by NICER, NuSTAR, and the Neil Gehrels Swift observatories, we confirm previous measurements of the black hole spin via the X-ray continuum method, a ~ 0.2. From polarization analysis only, we found consistent results with low black-hole spin, with an upper limit of a < 0.7 at a 90% confidence level. A slight increase of the polarization degree with energy, similar to other black-hole X-ray binaries in the soft state, is suggested from the data but with a low statistical significance.

All other papers

Yu-Jing Qin, Keming Zhang, Joshua Bloom, Jesper Sollerman, Erez A. Zimmerman, Ido Irani, Steve Schulze, Avishay Gal-Yam, Mansi Kasliwal, Michael W. Coughlin, Daniel A. Perley, Christoffer Fremling, Shrinivas Kulkarni

submitted

We identify the progenitor star of SN 2023ixf in the nearby galaxy Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion HST/ACS images. The supernova position, localized with diffraction-spike pattern and high precision relative astrometry, unambiguously coincides with a single progenitor candidate of m_F814W=24.96(-0.04)(+0.05). Forced photometry further recovers 2-sigma detections in the F673N and F675W bands and imposes robust flux limits on the bluer bands. Given the reported infrared excess and semi-regular variability of the progenitor, we fit a time-dependent spectral energy distribution (SED) model of a dusty red supergiant (RSG) to a combined dataset of HST photometry, as well as ground-based near-infrared and Spitzer/IRAC [3.6], [4.5] photometry from the literature. The progenitor closely resembles a RSG of T_eff=3343+/-27 K and logL=5.10+/-0.02, with a 0.11+/-0.01 dex (25.2+/-1.7 per cent) variation over the mean luminosity at a period of P=1128.3+/-6.5 days, heavily obscured by a dust envelope with an optical depth of tau=2.83+/-0.03 at 1 micron (or A_V=10.28+/-0.11 mag). Such observed signatures match a post-main sequence star of 18.1(-1.2)(+0.7) Msun, close to the most massive SN II progenitor, with a pulsation-enhanced mass-loss rate of M_dot=(3.58+/-0.15) x 10^(-4) Msun/yr. The dense and confined circumstellar material is likely ejected during the last episode of radial pulsation before the explosion. Notably, we find strong evidence for periodic variation of tau (or both T_eff and tau) along with luminosity, a necessary assumption to reproduce the wavelength dependence of the variability, which implies dust sublimation and condensation during radial pulsations. Given the observed SED, partial dust obscuration remains a possible scenario, but any unobstructed binary companion over 7.1 Msun can be ruled out.

M. Tristram, A.J. Banday, M. Douspis, X. Garrido, K.M. Górski, S. Henrot-Versillé, S. Ilić, R. Keskitalo, G. Lagache, C.R. Lawrence, B. Partridge, D. Scott

17 pages, 21 figures, 2 appendix pages. submitted to A&A

We present constraints on cosmological parameters using maps from the last Planck data release (PR4). In particular, we detail an upgraded version of the cosmic microwave background likelihood, HiLLiPoP, based on angular power spectra and relying on a physical modelling of the foreground residuals in the spectral domain. This new version of the likelihood retains a larger sky fraction (up to 75%) and uses an extended multipole range. Using this likelihood, along with low-l measurements from LoLLiPoP, we derive constraints on $\Lambda$CDM parameters that are in good agreement with previous Planck 2018 results, but with 10% to 20% smaller uncertainties. We demonstrate that the foregrounds can be accurately described in spectra domain with only negligible impact on $\Lambda$CDM parameters. We also derive constraints on single-parameter extensions to $\Lambda$CDM including $A_L$, $\Omega_K$, $N_{eff}$, and $\sum m_{\nu}$. Noteworthy results from this updated analysis include a lensing amplitude value of $A_L = 1.036 \pm 0.051$, which aligns more closely with theoretical expectations within the $\Lambda$CDM framework. Additionally, our curvature measurement, $\Omega_K = -0.012 \pm 0.010$, now demonstrates complete consistency with a flat universe, and our measurement of $S_8$ is closer to the measurements derived from large-scale structure surveys (at the 1.6$\sigma$ level).

S. D. T. Grant, D. B. Jess, C. J. Dillon, M. Mathioudakis, C. A. Watson, J. A. G. Jackman, D. G. Jackson, P. J. Wheatley, M. R. Goad, S. L. Casewell, D. R. Anderson, M. R. Burleigh, R. G. West, J. I. Vines

16 pages, 5 figures, 5 tables. Accepted to ApJ

Previous examinations of fully-convective M-dwarf stars have highlighted enhanced rates of nanoflare activity on these distant stellar sources. However, the specific role the convective boundary, which is believed to be present for spectral types earlier than M2.5V, plays on the observed nanoflare rates is not yet known. Here, we utilize a combination of statistical and Fourier techniques to examine M-dwarf stellar lightcurves that lie on either side of the convective boundary. We find that fully convective M2.5V (and later sub-types) stars have greatly enhanced nanoflare rates compared with their pre-dynamo mode transition counterparts. Specifically, we derive a flaring power-law index in the region of $3.00 \pm 0.20$, alongside a decay timescale of $200 \pm 100$~s for M2.5V and M3V stars, matching those seen in prior observations of similar stellar sub-types. Interestingly, M4V stars exhibit longer decay timescales of $450 \pm 50$~s, along with an increased power-law index of $3.10 \pm 0.18$, suggesting an interplay between the rate of nanoflare occurrence and the intrinsic plasma parameters, for example, the underlying Lundquist number. In contrast, partially convective (i.e., earlier sub-types from M0V to M2V) M-dwarf stars exhibit very weak nanoflare activity, which is not easily identifiable using statistical or Fourier techniques. This suggests that fully convective stellar atmospheres favor small-scale magnetic reconnection, leading to implications for the flare-energy budgets of these stars. Understanding why small-scale reconnection is enhanced in fully convective atmospheres may help solve questions relating to the dynamo behavior of these stellar sources.

Zoe A. Le Conte, Dimitri A. Gadotti, Leonardo Ferreira, Christopher J. Conselice, Camila de Sá-Freitas, Taehyun Kim, Justus Neumann, Francesca Fragkoudi, E. Athanassoula, Nathan J. Adams

Submitted to MNRAS. 15 pages, 10 figures. Figure 7 summarises the main results

The presence of a stellar bar in a disc galaxy indicates that the galaxy hosts a dynamically settled disc and that bar-driven processes are taking place in shaping the evolution of the galaxy. Studying the cosmic evolution of the bar fraction in disc galaxies is therefore essential to understand galaxy evolution in general. Previous studies have found, using the Hubble Space Telescope (HST), that the bar fraction significantly declines from the local Universe to redshifts near one. Using the first four pointings from the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science Survey (CEERS) and the initial public observations for the Public Release Imaging for Extragalactic Research (PRIMER), we extend the studies on the bar fraction in disc galaxies to redshifts $1 \leq z \leq 3$, i.e., for the first time beyond redshift two. We only use galaxies that are also present in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) on the Extended Groth Strip (EGS) and Ultra Deep Survey (UDS) HST observations. An optimised sample of 768 close-to-face-on galaxies is visually classified to find the fraction of bars in disc galaxies in two redshift bins: $1 \leq z \leq 2$ and $2 < z \leq 3$. The bar fraction decreases from $\sim 18.9^{+ 9.7}_{- 9.4}$ per cent to $\sim 6.6^{+ 7.1}_{- 5.9}$ per cent (from the lower to the higher redshift bin), but is $\sim 3 - 4$ times greater than the bar fraction found in previous studies using bluer HST filters. Our results show that bar-driven evolution commences at early cosmic times and that dynamically settled discs are already present at a lookback time of $\sim 11$ Gyrs.

Anna C. Wright, Jason Tumlinson, Molly S. Peeples, Brian W. O'Shea, Cassandra Lochhaas, Lauren Corlies, Britton D. Smith, Nguyen Binh, Ramona Augustin, Raymond C. Simons

32 pages, 15 figures, submitted to ApJ

Over the next decade, the astronomical community will be commissioning multiple wide-field observatories well-suited for studying stellar halos in both integrated light and resolved stars. In preparation for this, we use five high-resolution cosmological simulations of Milky Way-like galaxies from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) suite to explore the properties and components of stellar halos. At $z=0$, we find that the FOGGIE stellar halos have masses, metallicity gradients, and surface brightness profiles that are consistent with observations. In agreement with other simulations, the FOGGIE stellar halos receive 30-40% of their mass from stars that formed in situ. However, this population tends to be more centrally concentrated in the FOGGIE simulations and therefore does not contribute excess light or mass to the outskirts of the halos. The rest of the stars in each stellar halo are accreted from $\sim$10-50 other galaxies, with the majority of the accreted mass originating in 2-4 of these contributors. While the phase-mixed inner halo ($r<$50 kpc) of each FOGGIE galaxy includes stars from a large number of contributors, the halo outskirts of three of the five galaxies are primarily made up of stars from only a few contributors. We predict that upcoming wide-field observatories, like the Nancy Grace Roman Space Telescope, will probe stellar halos around Milky Way-like galaxies out to $\sim$100 kpc in integrated light and will be able to distinguish the debris of dwarf galaxies with extended star formation histories from the underlying halo with resolved color-magnitude diagrams.

Diego Calderón, Ondřej Pejcha, Brian D. Metzger, Paul C. Duffell

18 pages, 13 figures, 1 table. Submitted to MNRAS

The disruption of a star by the tidal forces of a spinning black hole causes the stellar stream to precess, which affects the conditions for triggering the tidal disruption event (TDE). In this work, we study the effect of the precession of the tidal stream on TDE light curves, as a result of the interaction of the wind and luminosity of the event with the stellar stream wrapped around the black hole. We perform two-dimensional radiation-hydrodynamic simulations using the moving-mesh hydrodynamic code JET with its recently developed radiation treatment module. The model considers an isotropic wind and irradiation following the analytical fallback rate of a stellar polytrope, and their interaction with the remnant of the precessed stream. We study the effect of the black hole mass, the accretion efficiency, and the inclination between the orbital and spin planes on the observed light curves along different lines of sight. From our results, we were able to identify two extreme behaviours depending on the properties of the event: i) models with low-mass black holes ($\sim10^6~M_{\odot}$), low inclination ($\sim0$), and low accretion efficiency ($\sim0.01$) show light curves with a short early peak caused by the interaction of the wind with the inner edge of the stream. The line of sight has little effect on the light curves, since the stream covers only a small fraction of the solid angle, due to the precession occurring in the orbital plane; ii) models with high-mass black holes ($>10^7~M_{\odot}$), high inclination ($\sim90^{\circ}$), and high accretion efficiency ($\sim0.1$) produce light curves with luminosity peaks that can be delayed by up to 50-100~d as the stream blocks the radiation in the early phase of the event. Our results show that black hole spin and misalignment do not imprint a recognisable feature on the light curves but rather can add complications to their analysis.

Funda Bostanci, Tugba Boztepe, Tolga Guver, Tod E. Strohmayer, Yuri Cavecchi, Ersin Gogus, Diego Altamirano, Peter Bult, Deepto Chakrabarty, Sebastien Guillot, Gaurava K. Jaisawal, Christian Malacaria, Giulio C. Mancuso, Andrea Sanna, Jean H. Swank

Accepted for publication in the Astrophysical Journal

We present temporal and time-resolved spectral analyses of all the thermonuclear X-ray bursts observed from the neutron star low-mass X-ray binary (LMXB) 4U 1728-34 with NICER from June 2017 to September 2019. In total, we detected 11 X-ray bursts from the source and performed time-resolved spectroscopy. Unlike some of the earlier results for other bursting sources from NICER, our spectral results indicate that the use of a scaling factor for the persistent emission is not statistically necessary. This is primarily a result of the strong interstellar absorption in the line of sight towards 4U 1728-34, which causes the count rates to be significantly lower at low energies. We also searched for burst oscillations and detected modulations in six different bursts at around the previously known burst oscillation frequency of 363 Hz. Finally, we report the detection of oscillations prior to two bursts at 356 and 359 Hz, respectively. This is the first time in the literature where burst oscillations are detected before the rapid rise in X-ray flux, from any known burster. These oscillations disappear as soon as the burst rise starts and occur at a somewhat lower frequency than the oscillations we detect during the bursts.

Y. Bi, J. J. Yang, Y. Qin, Z. P. Qiang, J. C. Hong, B. Yang, Z. Xu, H. Liu, K. F. Ji

Accepted for publication in Astronomy & Astrophysics

Nanoflares are impulsive energy releases by magnetic reconnection in the braided coronal magnetic field, which is a potential mechanism for heating the corona. However, there are still sporadic observations of the interchange of braiding structure segments and footpoints inside coronal loops, which is predicted to be the morphological evolution of the reconnecting magnetic bundles in the nanoflare picture. This work aims to detect the evolutions of the pairs of braiding strands within the apparent single coronal loops observed in Atmospheric Imaging Assembly (AIA) images. The loop strands are detected on two kinds of upsampled AIA 193 \AA\ images, which are obtained by upscaling the Point Spread Function matched AIA images via Bicubic interpolation and are generated using a super-resolution convolutional neural network, respectively. The architecture of the network is designed to map the AIA images to unprecedentedly high spatial resolution coronal images taken by High-resolution Coronal Imager (Hi-C) during its brief flight. At times, pairs of separate strands that appear braided together later evolved into pairs of almost parallel strands with completely exchanged parts. These evolutions offer morphological evidence that magnetic reconnections between the braiding strands have taken place, which is further supported by the appearance of transient hot emissions containing significant high-temperature components (T > 5MK) at the footpoints of the braiding structures. The brief appearances of the two rearranging strands support that magnetic reconnections have occurred within what appears to be a single AIA loop.

F. H. Vincent, M. Wielgus, N. Aimar, T. Paumard, G. Perrin

20 pages, 14 figures, submitted to A&A

[Abridged] Context. The Galactic Center supermassive black hole is well known to exhibit transient peaks of flux density on a daily basis across the spectrum. Recent infrared and millimeter observations have strengthened the case for the association between these flares and circular orbital motion in the vicinity of the event horizon. The strongly polarized synchrotron radiation associated with these events leads to specific observables called QU loops, that is, looping motion in the Stokes QU plane of linear polarization. Aims. We want to deepen the understanding of the QU loops associated with orbiting hot spots. We compute such loops in Minkowski and Schwarzschild spacetimes in order to determine which aspects of the observed patterns are due to special- or general-relativistic phenomena. Results. We show that QU loops in Minkowski spacetime at low or moderate inclination i < 45 deg share all qualitative features of Schwarzschild QU loops: there exist QU loops for all setups considered (including face-on view and vertical magnetic field), there may be one or two QU loops per orbital period for a vertical magnetic field configuration, there are always two QU loops in case of a toroidal magnetic field. We provide analytical formulas in Minkowski spacetime to explain the details of this behavior. Moreover, we analyze the flux variation of the hot spot and show that it is dictated either by the angular dependence of the radiative transfer coefficients, or by relativistic beaming. In the former case, this can lead to extreme flux ratios even at moderate inclination. Finally, we highlight the increasing mirror asymmetry of the Schwarzschild QU track with increasing inclination and show that this behavior is a specific Schwarzschild feature caused by light bending.

Jeniveve Pearson, David J. Sand, Peter Lundqvist, Lluís Galbany, Jennifer E. Andrews, K. Azalee Bostroem, Yize Dong, Emily Hoang, Griffin Hosseinzadeh, Daryl Janzen, Jacob E. Jencson, Michael J. Lundquist, Darshana Mehta, Nicolás Meza Retamal, Manisha Shrestha, Stefano Valenti, Samuel Wyatt, Joseph P. Anderson, Chris Ashall, Katie Auchettl, Eddie Baron, Stéphane Blondin, Christopher R. Burns, Yongzhi Cai, Ting-Wan Chen, Laura Chomiuk, David A. Coulter, Dane Cross, Kyle W. Davis, Thomas de Jaeger, James M. DerKacy, Dhvanil D. Desai, Georgios Dimitriadis, Aaron Do, Joseph R. Farah, Ryan J. Foley, Mariusz Gromadzki, Claudia P. Gutiérrez, Joshua Haislip, Jonay I. González Hernández, Jason T. Hinkle, Willem B. Hoogendam, D. Andrew Howell, Peter Hoeflich, Eric Hsiao, et al. (43 additional authors not shown)

We present optical, infrared, ultraviolet, and radio observations of SN 2022xkq, an underluminous fast-declining type Ia supernova (SN Ia) in NGC 1784 ($\mathrm{D}\approx31$ Mpc), from $<1$ to 180 days after explosion. The high-cadence observations of SN 2022xkq, a photometrically transitional and spectroscopically 91bg-like SN Ia, cover the first days and weeks following explosion which are critical to distinguishing between explosion scenarios. The early light curve of SN 2022xkq has a red early color and exhibits a flux excess which is more prominent in redder bands; this is the first time such a feature has been seen in a transitional/91bg-like SN Ia. We also present 92 optical and 19 near-infrared (NIR) spectra, beginning 0.4 days after explosion in the optical and 2.6 days after explosion in the NIR. SN 2022xkq exhibits a long-lived C I 1.0693 $\mu$m feature which persists until 5 days post-maximum. We also detect C II $\lambda$6580 in the pre-maximum optical spectra. These lines are evidence for unburnt carbon that is difficult to reconcile with the double detonation of a sub-Chandrasekhar mass white dwarf. No existing explosion model can fully explain the photometric and spectroscopic dataset of SN 2022xkq, but the considerable breadth of the observations is ideal for furthering our understanding of the processes which produce faint SNe Ia.

Zhe-Yu Daniel Lin, Zhi-Yun Li, Ian W. Stephens, Manuel Fernández-López, Carlos Carrasco-González, Claire J. Chandler, Alice Pasetto, Leslie W. Looney, Haifeng Yang, Rachel E. Harrison, Sarah I. Sadavoy, Thomas Henning, A. Meredith Hughes, Akimasa Kataoka, Woojin Kwon, Takayuki Muto, Dominique Segura-Cox

22 pages, 13 figures, submitted to MNRAS

Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength $\lambda$. The polarization morphology at Band 7 is predominantly parallel to the disk minor axis but appears azimuthally oriented at Band 3, with the morphology at Band 6 in between the two. We present new ~0.2" (29 au) polarization observations at Q-Band (7.0 mm) using the Karl G. Jansky Very Large Array (VLA) and at Bands 4 (2.1 mm), 5 (1.5 mm), and 7 using ALMA, consolidating HL Tau's position as the protoplanetary disk with the most complete wavelength coverage in dust polarization. The polarization patterns at Bands 4 and 5 continue to follow the morphological transition with wavelength previously identified in Bands 3, 6, and 7. Based on the azimuthal variation, we decompose the polarization into contributions from scattering ($s$) and thermal emission ($t$). We find that $s$ decreases slowly with increasing $\lambda$, and $t$ increases more rapidly with $\lambda$ which are expected from optical depth effects of toroidally aligned, scattering prolate grains. The relatively weak $\lambda$ dependence of $s$ is consistent with large, porous grains. The sparse polarization detections from the Q-band image are also consistent with toroidally aligned prolate grains.

B. Bucciarelli, L. Stanghellini

10 pages, 8 Figures

We study the population of Galactic planetary Nebulae (PNe) and their central stars (CSs) through the analysis of their heliocentric distances and Galactic distribution. Distances are obtained by means of a revised statistical scale, based on an astrometrically-defined sample of CSs parallaxes from Gaia DR3 as calibrators. The statistical scale is applied to infer distances of a significant number (~850) of Galactic PNe, for which we deliver a new catalog of PN distances. By adopting a circular velocity curve of the Galaxy, we also derive 3D peculiar velocities from DR3 proper motions and published radial velocities of a large sample (~300) of PN CSs. We date PN progenitors based both on the best chemical abundances culled from the literature and on CS kinematic properties, finding a confirmation of the first method with the second. The slope of the radial oxygen gradient of the Galactic Disk traced by the complete PNe sample amounts to -0.0144 +/- 0.00385 [dex/kpc]. Furthermore, by distinguishing between PNe with old (> 7.5 Gyr) and young (< 1 Gyr) progenitors, we estimate the gradient to be respectively -0.0121 +/- 0.00465 and -0.022 +/- 0.00758 [dex/kpc], thus disclosing a mild steepening since Galaxy formation, with a slope change of 0.01 dex. These results are in broad agreement with previous PN studies, but now based on DR3 Gaia analysis, and also in agreement with what traced by most other Galactic probes.

Micael Andrade, Aion Viana (for the SWGO collaboration)

Proc. ICRC 2023, 7 pages, 3 figures

Dark matter is thought to make up most of the matter density of the Universe, yet its true nature remains uncertain. Among dark matter theories, Weakly Interacting Massive Particles (WIMPs) are a prominent candidate for dark matter because they can reproduce the observed abundance of dark matter in the universe. There are various methods for searching for WIMPs, one of which is indirect detection, which involves looking for the Standard Model particles produced by the decay or self-annihilation of dark matter particles. Within the mass range of GeV to PeV for the dark matter particle, this type of search can be conducted by detecting $\gamma$-rays in astrophysical objects with high concentrations of dark matter. Dwarf galaxies, although not the most dense, are excellent targets for this type of observation since they are dominated by dark matter, are relatively close to Earth, and have a low astrophysical background. In this work, the detectability of dark matter annihilation or decay signals from dwarf galaxies is predicted using the Southern Wide-field Gamma-ray Observatory (SWGO), a future $\gamma$-ray observatory that will be built in South America. This wide field-of-view survey instrument will be able to study many important dark matter targets in the Southern Hemisphere, and the combined observation of all targets will provide competitive, if not the best, limits for dark matter with masses in the range of hundreds of GeV to PeV.

Camila Angulo-Valdez, Rosa L. Becerra, Margarita Pereyra, Felipe Vargas, Keneth Garcia-Cifuentes, Alan M. Watson, Fabio De Colle, Nissim Fraija, Nathaniel R. Butler, Maria G. Dainotti, Simone Dichiara, William H. Lee, Eleonora Troja, Joshua S. Bloom, J. Jesús González, Alexander S. Kutyrev, J. Xavier Prochaska, Enrico Ramirez-Ruiz, Michael G. Richer

Submitted to MNRAS, 11 pages, 6 figures

We present analytical and numerical models of the bright long GRB 210822A at $z=1.736$. The intrinsic extreme brightness exhibited in the optical, which is very similar to other bright GRBs (e.g., GRBs 080319B, 130427A, 160625A 190114C, and 221009A), makes GRB 210822A an ideal case for studying the evolution of this particular kind of GRB. We use optical data from the RATIR instrument starting at $T+315.9$ s, with publicly available optical data from other ground-based observatories, as well as X-ray data from the Swift/X-ray Telescope (XRT) and data from the Swift/Ultraviolet/Optical Telescope (UVOT). The temporal profiles and spectral properties during the late stages align consistently with the conventional forward shock model, complemented by a reverse shock element that dominates optical emissions during the initial phases ($T<300$ s). Furthermore, we observe a break at $T=80000$ s that we interpreted as evidence of a jet break, which constrains the opening angle to be about $\theta_\mathrm{j}=(3-5)$ degrees. Finally, we apply a machine-learning technique to model the multi-wavelength light curve of GRB 210822A using the AFTERGLOWPY library. We estimate the angle of sight $\theta_{obs}=(6.4 \pm 0.1) \times 10^{-1}$ degrees, the energy $E_0= (7.9 \pm 1.6)\times 10^{53}$ ergs, the electron index $p=2.52 \pm 0.12$, the thermal energy fraction in electrons $\epsilon_e=(4.63 \pm 0.91) \times 10^{-5}$ and in the magnetic field $\epsilon_B= (8.66 \pm 1.01) \times 10^{-6}$, the efficiency $\chi = 0.89 \pm 0.01$, and the density of the surrounding medium $n_\mathrm{0} = 0.85 \pm 0.01$.

Tobias Cornelius Hinse, Bertil F. Dorch, Lars V. T. Occhionero, Jakob P. Holck

20 pages, 13 figures, accepted, Frontiers in Astronomy and Space Sciences

The 450th anniversary of the discovery of the SN 1572 supernova event was celebrated in 2022. A closer look at the historical development of the field of supernova astronomy reveals the scientific importance of Tycho Brahe's 1572 observations of this "new star". In their quest to learn more about the new type of stellar explosion and subsequent evolution, the initial protagonists in this field (Baader and Zwicky among others) gradually turned their attention to the final remnant state of these supernova events. Since the remnant object thought to be associated with the extragalactic supernova event was found to be very dim, the focus quickly shifted toward nearby galactic events. It is at this point where Tycho Brahe's observations played an important and often overlooked role in the context of the development of stellar evolution as a scientific field. Tycho Brahe's meticulous and detailed recordings of the change in brightness of the new star not only allowed modern astronomers to classify SN 1572 as a supernova event but also helped them pinpoint the exact astrometric location of SN 1572. These findings helped to empirically link extragalactic supernova events to nearby past supernova remnants in the Milky Way. This enabled subsequent observations allowing further characterization. Transforming the historical recordings to a standardized photometric system also allowed the classification of SN 1572 as a type I supernova event.

We derive an effective field theory (EFT) for cosmological Lyman alpha forest fluctuations valid for the power spectrum at the one-loop order. The ``bottom-up'' EFT expansion at the level of the transmitted flux is identical to the line-of-sight dependent bias model first derived by Desjacques et al. We confirm this result by a ``top-down'' derivation based on the exponential map of the optical depth field. Specifically, we show that the combination of the exponential map and conditions of renormalizability generates the same EFT expansion as the ``bottom-up'' approach. In passing, we point out inconsistencies of the tree-level perturbative expansion of the exponential map without counterterms. To facilitate practical applications, we generalize the FFTLog method for efficient calculations of one-loop integrals from new line-of-sight dependent operators. Finally, we compare the one-loop EFT model against data from the Sherwood hydrodynamic simulation. The theory fits the data with sub percent accuracy up to $k= $ 3 $h$Mpc$^{-1}$ at $z= 2.8$ for both 3D and 1D correlations. Our model can be readily used for cosmological full-shape analyses of Lyman alpha forest data.

P. Tisserand, C.L. Crawford, J. Soon, G.C. Clayton, A.J. Ruiter, I.R. Seitenzahl

9 pages, 6 figures, submitted to A&A

The Gaia DR3 release includes heliocentric radial velocity measurements and velocity variability indices for tens of millions of stars observed over 34 months. In this study, we utilise these indices to investigate the intrinsic radial velocity variations of Hydrogen-deficient carbon (HdC) stars and Extreme Helium (EHe) stars across their large ranges of temperature and brightness. Taking advantage of the newly defined HdC temperature classes, we examine the evolution of the total velocity amplitude with effective temperature. Additionally, we analyse the variation in dust production rate in R Coronae Borealis (RCB) stars with temperature, using two different proxies for RCB stars' photometric state: one from Gaia and another from the 2MASS survey. Furthermore, we provide a list of heliocentric radial velocities and associated errors for targets not observed by Gaia DR3. Some of these velocities were obtained from our 2.3m/WiFeS survey, while others were retrieved from the literature. Our observations reveal an interesting trend in the evolution of the maximum radial velocity amplitude across each HdC temperature class, indicating possibly that the helium shell-burning giant stage starts with strong atmospheric motions that decrease in strength, up to $\sim$6000 K, before picking up again as the HdC star atmosphere shrink further in size and reach warmer temperature. We also observe a correlation between stellar temperature and the dust production rate. The dust formation rate appears to be much higher in colder RCB stars compared to warmer ones.

P. Tisserand, C.L. Crawford, J. Soon, G.C. Clayton, A.J. Ruiter, I.R. Seitenzahl

16 pages, 17 figures, submitted to A&A

Upon its release the Gaia DR3 catalogue has led to tremendous progress in multiple fields of astronomy by providing the complete astrometric solution for nearly 1.5 billion sources. We analysed the photometric and astrometric results for Hydrogen-deficient (HdC), Extreme Helium (EHe), and DYPer type stars to identify any potential biases. This analysis aimed to select stars suitable for kinematic and spatial distribution studies. We investigated the information obtained from the Gaia IPD process, which was cross-matched with Gaia light curves. One main objective was to understand the impact of photometric declines in R Coronae Borealis (RCB) stars on Gaia astrometry. Based on the evidence gathered, we have reached the conclusion that the astrometric fits for numerous RCB stars are not valid due to the Gaia PSF chromaticity effect in both shape and centroid. The astrometric results of all stars with significant time-dependent colour variation should be similarly affected. RCB stars might thus be promising sources to correct this effect in future Gaia releases. Furthermore, after validating the Gaia astrometric results for 92 stars, we observed that the majority of HdC and EHe stars are distributed across the three old stellar structures, the thick disk, the bulge and the halo. However, we have also uncovered evidence indicating that some HdC and EHe stars exhibit orbits characteristic of the thin disk. This is also particularly true for all DYPer type stars under study. Finally, we have provided a list of star memberships for each Galactic substructure. We are beginning to observe a relationship between kinematics, stellar population, and metallicity in RCB and EHe stars. That relation can be explained, within the double degenerate scenario, by the large range in the delay time distribution expected from population synthesis simulations, particularly through the HybCO merger channel.

Anastasia Lavrukhina, Konstantin Malanchev, Matwey V. Kornilov

4 pages, 2 figures, Extended version of a note published in RNAAS

Modern astronomical surveys produce millions of light curves of variable sources. These massive data sets challenge the community to create automatic light-curve processing methods for detection, classification, and characterisation of variable stars. In this paper, we present a novel method for extracting the variable components of a light curve based on Otsu's thresholding method. To validate the effectiveness of this method, we apply it to the light curves of detached eclipsing binaries and dwarf novae, sourced from OGLE catalogues.

Arthur D. Adams, Michael R. Meyer, Alex R. Howe, Ben Burningham, Sebastian Daemgen, Jonathan Fortney, Mike Line, Mark Marley, Sascha P. Quanz, Kamen Todorov

39 pages, 21 figures. Accepted for publication in the Astronomical Journal

We present model constraints on the atmospheric structure of HD 106906 b, a planetary-mass companion orbiting at a ~700 AU projected separation around a 15 Myr-old stellar binary, using the APOLLO retrieval code on spectral data spanning 1.1-2.5 $\mu$m. C/O ratios can provide evidence for companion formation pathways, as such pathways are ambiguous both at wide separations and at star-to-companion mass ratios in the overlap between the distributions of planets and brown dwarfs. We benchmark our code against an existing retrieval of the field L dwarf 2M2224-0158, returning a C/O ratio consistent with previous fits to the same JHKs data, but disagreeing in the thermal structure, cloud properties, and atmospheric scale height. For HD 106906 b, we retrieve C/O $=0.53^{+0.15}_{-0.25}$, consistent with the C/O ratios expected for HD 106906's stellar association and therefore consistent with a stellar-like formation for the companion. We find abundances of H$_2$O and CO near chemical equilibrium values for a solar metallicity, but a surface gravity lower than expected, as well as a thermal profile with sharp transitions in the temperature gradient. Despite high signal-to-noise and spectral resolution, more accurate constraints necessitate data across a broader wavelength range. This work serves as preparation for subsequent retrievals in the era of JWST, as JWST's spectral range provides a promising opportunity to resolve difficulties in fitting low-gravity L dwarfs, and also underscores the need for simultaneous comparative retrievals on L dwarf companions with multiple retrieval codes.

The soft excess is a significant emission component in the Soft (<1 keV) X-ray spectra of many AGN. It has been explained by disk reflection, a warm corona and other models. Understanding its origin is crucial for the energy budget of AGN emission, and for using it to study the inner accretion disk. Here, we track the weeks-to-months variability of several AGN that show different levels of soft excess strength with NICER. We use the variability time scales to compare the relative size of the soft excess emission region to the corona producing the hard X-ray emission above 1 keV. We find that the size of the soft excess emission region relative to the corona is not the same for the three sources studied. For TON S180, the soft excess region is comparable in size to the hard corona. While for MRK 335 and 1H0707-495, the soft excess region is larger than the corona by a factor of 2-4. This is the first time the relative sizes are quantified independently of the assumptions of the spectral models.

Michael Rodruck, Jane Charlton, Sanchayeeta Borthakur, Aparna Chitre, Patrick R. Durrell, Debra Elmegreen, Jayanne English, Sarah C. Gallagher, Caryl Gronwall, Karen Knierman, Iraklis Konstantopoulos, Yuexing Li, Moupiya Maji, Brendan Mullan, Gelys Trancho, William Vacca

26 pages, 27 figures. Accepted at MNRAS

We present results of a Hubble Space Telescope (HST) UBVI-band study of star clusters in tidal tails, using new WFC3 and ACS imaging to complement existing WFPC2 data. We survey 12 tidal tails across seven merging systems, deriving ages and masses for 425 star cluster candidates (SCCs). The stacked mass distribution across all systems follows a power law of the form $dN/dM \propto M^{\beta}$, with $\beta = -2.02 \pm 0.15$, consistent with what is seen in other star forming environments. GALEX and Swift UV imaging provide star formation rates (SFRs) for our tidal tails, which when compared with ages and masses of our SCCs, allows for a determination of the cluster formation efficiency (CFE). We find the CFE increases with increasing SFR surface density, matching the theoretical model. We confirm this fit down at SFR densities lower than previously measured (log $\Sigma_\text{SFR} \: (\text{M}_\odot \: \text{yr}^{-1} \: \text{kpc}^{-2}) \approx -4.2$), as related to the CFE. We determine the half-light radii for a refined sample of 57 SCCs with our HST WFC3 and ACS imaging, and calculate their dynamical age, finding the majority of them to be gravitationally bound. We also provide evidence of only low-mass ($< 10^4 \: \text{M}_\odot$) cluster formation in our nearest galaxy, NGC 1487, consistent with the theory that this system is a dwarf merger.

Chirag Modi, Oliver H.E. Philcox

6 pages, 3 figures

We propose a new framework for the analysis of current and future cosmological surveys, which combines perturbative methods (PT) on large scales with conditional simulation-based implicit inference (SBI) on small scales. This enables modeling of a wide range of statistics across all scales using only small-volume simulations, drastically reducing computational costs, and avoids the assumption of an explicit small-scale likelihood. As a proof-of-principle for this hybrid simulation-based inference (HySBI) approach, we apply it to dark matter density fields and constrain cosmological parameters using both the power spectrum and wavelet coefficients, finding promising results that significantly outperform classical PT methods. We additionally lay out a roadmap for the next steps necessary to implement HySBI on actual survey data, including consideration of bias, systematics, and customized simulations. Our approach provides a realistic way to scale SBI to future survey volumes, avoiding prohibitive computational costs.

This paper describes a repository for ontologies of astronomy, astronautics, and other space-related topics. It may be called AstroPortal (or SpacePortal), AstroHub (or SpaceHub), etc. The creation of this repository will be applicable to academic, research and other data-intensive sectors. It is relevant for space sciences (including astronomy), Earth science, and astronautics (spaceflight), among other data-intensive disciplines. The repository should provide a centralized platform to search, review and create ontologies for astro-related topics. It thereby can decrease research time, while also providing a user-friendly means to study and compare knowledge organization systems or semantic resources of the target domains. With no apparent repository available on the target domain, this paper also expresses a novel concept.

The rotational lightcurves of the Pluto-Charon system were previously believed to be solely attributed to their surfaces. However, a proposed scenario of haze cooling \citep{2017Natur.551..352Z} suggests that the atmospheric haze of Pluto could significantly contribute to mid-infrared emission, which calls for a revisit of previous analyses. In this study, we employ a Bayesian retrieval approach to constrain the haze emission from the rotational lightcurves of the Pluto-Charon system. The lightcurves were observed by the Spitzer and Herschel telescopes at 24 and 70 $\mu$m, and were combined with the latest surface albedo maps of Pluto and Charon from the New Horizons spacecraft. Our results show that including the haze emission is consistent with all current observations, with the best-fit haze flux around 1.63 mJy. This is in agreement with the composition of Titan-like tholins. However, the ``surface only" scenario, which excludes the haze contribution, can still explain the observations. We conclude that the current data at 24 $\mu$m cannot constrain Pluto's haze emission due to the degeneracy with Charon's surface emission. Regardless, some surface properties of Pluto are well constrained by the shape of the lightcurves, with a thermal inertia of approximately 8--10 MKS and a relatively low CH$_4$ emissivity of 0.3--0.5. We suggest that observations by the JWST telescope at 18 $\mu$m, which can resolve Pluto from Charon, could directly probe the haze emission of Pluto due to the low surface emission at that wavelength.

James J. Buchanan, Michael D. Schneider, Kerianne Pruett, Robert E. Armstrong

34 pages, 16 figures

We apply Markov Chain Monte Carlo (MCMC) to the problem of parametric galaxy modeling, estimating posterior distributions of galaxy properties such as ellipticity and brightness for more than 100,000 images of galaxies taken from DC2, a simulated telescope survey resembling the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We use a physically informed prior and apply selection corrections to the likelihood. The resulting posterior samples enable rigorous probabilistic inference of galaxy model parameters and their uncertainties. These posteriors are one key ingredient in a fully probabilistic description of galaxy catalogs, which can ultimately enable a refined Bayesian estimate of cosmological parameters. We systematically examine the reliability of the posterior mean as a point estimator of galaxy parameters, and of the posterior width as a measure of uncertainty, under some common modeling approximations. We implement the probabilistic modeling and MCMC inference using the JIF (Joint Image Framework) tool, which we make freely available online.

Xilong Liang, Suk-Jin Yoon, Jingkun Zhao, Zhaoyu Li, Jiajun Zhang, Yaqian Wu

15 pages, 5 figures. Accepted by ApJ

Using Gaia DR3 data and the wavelet transformation technique, we study the substructures of the Hercules moving group (HMG): Hercules 1 (H1) and Hercules 2 (H2). Spectroscopic survey data from LAMOST, APOGEE, and GALAH are used to obtain metallicities and ages of stars belonging to the HMG. Our analysis leads to several key findings as follows: ($a$) the HMG is on average richer in metallicity than the Galactic disk, with H2 being metal richer than H1; ($b$) the HMG likely has a radial metallicity gradient distinct from that of the disk; ($c$) the HMG is on average older than the disk, with H2 being older than H1; ($d$) the HMG likely has a radial age gradient distinct from that of the disk; and ($e$) the metallicity and age distributions of the HMG depend mainly on the Galactic radius but show no dependence on the azimuthal velocity. Taken all together, we conclude that the HMG is composed primarily of stars undergoing radial migration. We suggest that the HMG is associated with a higher-order dynamical resonance of the bar of the Galaxy.

Keith R. Dienes, Lucien Heurtier, Fei Huang, Tim M.P. Tait, Brooks Thomas

61 pages, LaTeX, 21 figures

Many theories of BSM physics predict the existence of large or infinite towers of decaying states. In a previous paper (arXiv:2111.04753) we pointed out that this can give rise to a surprising cosmological phenomenon that we dubbed "stasis" during which the relative abundances of matter and radiation remain constant across extended cosmological eras even though the universe is expanding. Indeed, such stasis epochs are universal attractors, with the universe necessarily entering (and later exiting) such epochs for a wide variety of initial conditions. Matter/radiation stasis is therefore an important and potentially unavoidable feature of many BSM cosmologies. In this paper we extend our arguments to universes containing significant amounts of vacuum energy, and demonstrate that such universes also give rise to various forms of stasis between vacuum energy and either matter or radiation. We also demonstrate the existence of several forms of "triple stasis" during which the abundances of matter, radiation, and vacuum energy all simultaneously remain fixed despite cosmological expansion. We further describe several close variants of stasis which we call "quasi-stasis" and "oscillatory stasis" and discuss the circumstances under which each of these can arise. Finally, we develop a general formalism for understanding the emergence of stasis within BSM cosmologies irrespective of the number or type of different energy components involved. Taken together, these results greatly expand the range of theoretical and phenomenological possibilities for the physics of the early universe, introducing new types of cosmological eras which may play an intrinsic and potentially inevitable role within numerous BSM cosmologies.

A. Waszewski, J.S. Morgan, R. Chhetri, R. Ekers, M.C.M. Cheung, N.D.R Bhat, M. Johnston-Hollitt

14 pages, 5 figures

We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations being taken at sunspot minimum, this search has revealed several interesting transitory features characterised by elevated scintillation levels. One solar wind enhancement is captured in two observations several hours apart, allowing its radial movement away from the Sun to be measured. We present here a methodology for measuring the plane-of-sky velocity for the moving heliospheric structure. The plane-of-sky velocity was inferred as $0.66\pm0.147\,^{\text{o}}\text{hr}^{-1}$, or $480\pm106\,\text{km}\,\text{s}^{-1}$ assuming a distance of 1AU. After cross-referencing our observed structure with multiple catalogues of heliospheric events, we propose that the likely source of our observed structure is a stream-interaction region originating from a low-latitude coronal hole. This work demonstrates the power of widefield interplanetary scintillation observations to capture detailed features in the heliosphere which are otherwise unresolvable and go undetected.

Atsuki Ishibashi, Hiroshi Hidaka, W. M. C. Sameera, Yasuhiro Oba, Naoki Watanabe

Accepted in ApJ; 30 pages, 10 figures, 1 Table

The CH$_3$O and CH$_2$OH radicals can be important precursors of complex organic molecules (COMs) in interstellar dust. The COMs presumably originating from these radicals were abundantly found in various astronomical objects. Because each radical leads to different types of COMs, determining the abundance ratio of CH$_3$O to CH$_2$OH is crucial for a better understanding of the chemical evolution to various COMs. Recent work suggested that the reaction between CH$_3$OH and OH on ice dust plays an important role in forming CH$_3$O and CH$_2$OH radicals. However, quantitative details on the abundance of these radicals have not been presented to date. Herein, we experimentally determined the branching ratio (CH$_3$O/CH$_2$OH) resulting from the CH$_3$OH + OH reaction on the water ice surface at 10 K to be 4.3 $\pm$ 0.6. Furthermore, the CH$_3$O product in the reaction would participate in subsequent diffusive reactions even at a temperature as low as 10 K. This fact should provide critical information for COMs formation models in cold molecular clouds.

Raffaella Morganti, Suma Murthy, Tom Oosterloo, Jay Blanchard, Claire Cook, Zsolt Paragi, Monica Orienti, Hiroshi Nagai, Robert Schulz

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

We present new JVLA and VLBA observations tracing the HI in the central region of 3C84 (Perseus A). This radio source is hosted by the bright cluster galaxy NGC 1275 in the centre of the iconic Perseus cluster. With the JVLA, we detected broad (FWHM~500 km/s) HI absorption at arcsecond resolution (~300 pc) centred at the systemic velocity of NGC 1275 against the bright radio continuum, suggesting that the detected gas is very close to the supermassive black hole (SMBH). However, we did not detect any absorption in the higher-resolution VLBA data against the parsec-scale radio core and jet. Based on a comparison of the properties of the HI absorption with those of the molecular circum-nuclear disc (CND) known to be present in NGC 1275, we argue that the HI seen in absorption arises from HI in this fast-rotating CND, and that neutral atomic hydrogen is present as close as ~20 pc from the SMBH. The radio continuum providing the background for absorption arises from non-thermal synchrotron emission from the star formation activity in the CND, whose presence has been reported by earlier VLBA studies. We did not detect any signature that the HI gas is affected by an interaction with the radio jet. Thus, at this stage of the evolution of the source, the impact of the radio jet on the gas in NGC 1275 mainly creates cavities on much larger galaxy scales. Overall, the properties of the CND in Perseus A present strong similarities with Mrk 231, suggesting that, unlike often assumed, HI absorption can arise against the radio emission from star formation in a CND. With the JVLA, we serendipitously detected a new, faint absorbing system that is redshifted by ~2660 km/s, in addition to the already known high-velocity absorption system that is redshifted 2850 km/s with respect to NGC 1275. We identify this new system as gas that is stripped from a foreground galaxy falling into the Perseus cluster.

We use a large set N-body/hydrodynamical simulations to study the physical properties of the merging cluster El Gordo. We found that the observed X-ray structures, along with other data, can be fairly matched by simulations with collision velocities 2,000 kms <= V <= 2,500 kms and impact parameters 600 kpc <= P <= 800 kpc. The mass of the primary is constrained to be between 10^{15} M_sun and ~ 1.6 10^{15} M_sun, in accord with recent lensing-based mass measurements. Moreover, a returning, post-apocenter, scenario is not supported by our head-on simulations. We also consider merger models that incorporate dark matter self-interactions. The simulation results show that the observed spatial offsets between the different mass components are well reproduced in self-interacting dark matter models with an elastic cross-section in the range \sigma_DM/m_X ~ 4 -5 cm^2/gr. In addition, the mean relative line-of-sight radial velocity between the two brightest cluster galaxies is found to be of the order of several hundreds of km/s. We argue that these findings provide an unambiguous signature of a dark matter behavior that exhibits collisional properties in a very energetic high-redshift cluster collision. The range of allowed values we found for sigma_DM/m_X is however inconsistent with present upper limits. To resolve this tension we suggest the possibility that the self-interacting dark matter model used here should be considered as only a low order approximation, and that the underlying physical processes that describe the interaction of dark matter in major cluster mergers are more complex than can be adequately represented by the commonly assumed approach based on scattering of dark matter particles.

Variations of the azimuthal magnetic fields of the Sun in the 23-25 activity cycles of the activity cycles are considered. To identify azimuthal magnetic fields, the analysis of daily observations of LOS magnetic fields from the regions near the solar limb was performed. It is shown that with a sufficiently large averaging of the data, large-scale structures are distinguished that can be interpreted by horizontal magnetic fields directed along the East-West line. Azimuthal magnetic fields are visible both in the low-latitude zone and at high latitudes. Azimuthal fields at the same latitudes have opposite directions in the northern and southern hemispheres, and also change sign in even and odd cycles of activity. The mechanism of formation of global azimuthal magnetic fields and their role in the cycle of solar activity is discussed. The near-surface azimuthal magnetic field is closely related to the activity cycle. Apparently, the azimuthal field is formed from U-shaped flux tubes of active regions (AR). Due to the presence of the tilt angle AR during differential rotation, the subsurface magnetic fields are pulled in the azimuthal direction. The role of azimuthal magnetic fields in solar activity cycles is considered. A scheme for generating a magnetic field according to a scheme different from Babcock-Layton dynamo models is proposed.

Kyu-Hyun Chae

14 pages, 12 figures, submitted to ApJ (this new work complements the paper ApJ, 952, 128 [ arXiv:2305.04613 ] in an important way)

It is found that Gaia DR3 binary stars selected with stringent requirements on astrometric measurements and radial velocities naturally satisfy Newtonian dynamics without hidden close companions when projected separation $s > 2$ kau, showing that pure binaries can be selected. It is then found that pure binaries selected with the same criteria show a systematic deviation from the Newtonian expectation when $s < 2$ kau. When both proper motions and parallaxes are required to have precision better than 0.003 and radial velocities better than 0.2, I obtain 1558 statistically pure binaries within a 'clean' $G$-band absolute magnitude range. From this sample, I obtain an observed to Newtonian predicted kinematic acceleration ratio of $\gamma_g=g_{\rm{obs}}/g_{\rm{pred}}=1.43^{+0.23}_{-0.19}$ for acceleration $< 10^{-10}$ m s$^{-2}$, in excellent agreement with a recent finding $1.43\pm 0.06$ for a much larger general sample with the amount of hidden close companions self-calibrated. I also investigate the radial profile of stacked sky-projected relative velocities without a deprojection to the 3D space. The observed profile matches the Newtonian predicted profile for $s < 2$ kau without any free parameters but shows a clear deviation at a larger separation with a significance of $4.6\sigma$. The projected velocity boost factor for $s>8$ kau is measured to be $\gamma_{v_p} = 1.18\pm 0.06$ matching $\sqrt{\gamma_g}$. Finally, for a small sample of 23 binaries with exceptionally precise radial velocities (precision $<0.0043$) the directly measured relative velocities in the 3D space also show a boost at larger separations. These results robustly confirm the recently reported gravitational anomaly at low acceleration for a general sample.

C. Gieser, H. Beuther, E. F. van Dishoeck, L. Francis, M. L. van Gelder, L. Tychoniec, P. J. Kavanagh, G. Perotti, A. Caratti o Garatti, T. P. Ray, P. Klaassen, K. Justtanont, H. Linnartz, W. R. M. Rocha, K. Slavicinska, L. Colina, M. Güdel, Th. Henning, P.-O. Lagage, G. Östlin, B. Vandenbussche, C. Waelkens, G. Wright

15 pages, 7 figures, accepted for publication in A&A

(abridged) We study and compare the warm (>100 K) and cold (<100 K) material toward the high-mass star-forming region IRAS 23385+6053 (IRAS 23385 hereafter) combining high angular resolution observations in the mid-infrared (MIR) with the JWST Observations of Young protoStars (JOYS) project and with the NOEMA at mm wavelengths at angular resolutions of 0.2"-1". The spatial morphology of atomic and molecular species is investigated by line integrated intensity maps. The temperature and column density of different gas components is estimated using H2 transitions (warm and hot component) and a series of CH3CN transitions as well as 3 mm continuum emission (cold component). Toward the central dense core in IRAS 23385 the material consists of relatively cold gas and dust (~50 K), while multiple outflows create heated and/or shocked H2 and show enhanced temperatures (~400 K) along the outflow structures. An energetic outflow with enhanced emission knots of [Fe II] and [Ni II] hints at J-type shocks, while two other outflows have enhanced emission of only H2 and [S I] caused by C-type shocks. The latter two outflows are also more prominent in molecular line emission at mm wavelengths (e.g., SiO, SO, H2CO, and CH3OH). Even higher angular resolution data are needed to unambiguously identify the outflow driving sources given the clustered nature of IRAS 23385. While most of the forbidden fine structure transitions are blueshifted, [Ne II] and [Ne III] peak at the source velocity toward the MIR source A/mmA2 suggesting that the emission is originating from closer to the protostar.

Takayuki. Saito, K. Hashiyama, H. Iwasaki, H. Kubo, M. Mizote, A. Okumura, H. Tajima, T. Yamamoto

Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:submit/ 2309.08219 )

The recent improvements in the performance of the silicon photomultipliers (SiPMs) made them attractive options as photo sensors of imaging atmospheric Cherenkov telescopes (IACTs). In fact, they are already adopted in some IACTs such as FACT and the Small-Sized Telescopes of the Cherenkov Telescope Array (CTA). However, the application to the Large-Sized Telescopes (LSTs) of CTA requires additional studies. As the pixel size of LSTs is larger than the nominal size of SiPMs, the signal from multiple sensors must be summed up. Also, the high detection efficiency of the night sky background (NSB) photons may degrade the telescope performance. To overcome this, the pulse width must be as small as 3 ns and the detection efficiency for NSB photons must be suppressed as much as possible. Heat generation and gain stabilization are also issues. We studied different types of SiPMs from Hamamatsu photonics and characterized them for the LST application, addressing the previous points. Also, to prove the SiPM performance in LST, we are developing a SiPM module which can be installed in the exisiting LST camera. Here we present the results of this evaluation and the status of the test bench module development.

Mamadou N'Diaye, Arthur Vigan, Byron Engler, Markus Kasper, Serban Leveratto, Johan Floriot, Michel Marcos, Christophe Bailet, Kjetil Dohlen

17 pages, 10 figures, pre-print of the proceeding of the AO4ELT7 conference held in June 2023 in Avignon, France

We propose to explore a cascade extreme Adaptive optics (ExAO) approach with a second stage based on a Zernike wavefront sensor (ZWFS) for exoplanet imaging and spectroscopy. Most exoplanet imagers currently use a single-stage ExAO to correct for the effects of atmospheric turbulence and produce high-Strehl images of observed stars in the near-infrared. While such systems enable the observation of warm gaseous companions around nearby stars, adding a second-stage AO enables to push the wavefront correction further and possibly observe colder or smaller planets. This approach is currently investigated in different exoplanet imagers (VLT/SPHERE, Mag-AOX, Subaru/SCExAO) by considering a Pyramid wavefront sensor (PWFS) in the second arm to measure the residual atmospheric turbulence left from the first stage. Since these aberrations are expected to be very small (a few tens of nm in the near-infrared domain), we propose to investigate an alternative approach based on the ZWFS. This sensor is a promising concept with a small capture range to estimate residual wavefront errors thanks to its large sensitivity, simple phase reconstruction and easiness of implementation. In this contribution, we perform preliminary tests on the GHOST testbed at ESO to validate this approach experimentally. Additional experiments with petalling effects are also showed, giving promising wavefront correction results. Finally, we briefly discuss a first comparison between PWFS-based and ZWFS-based second-stage AO to draw preliminary conclusions on the interests of both schemes for exoplanet imaging and spectroscopy with the upgrade of the current exoplanet imagers and the envisioned ExAO instruments for ELTs.

I. Contopoulos, D. Ntotsikas, K. N. Gourgouliatos

The pulsar magnetosphere is divided into a corotating region of closed field lines surrounded by open field lines that emanate from the two poles of the star, extend to infinity and are separated by an equatorial current sheet. The three regions meet at a magnetospheric Y-point. In steady-state solutions of the ideal force-free magnetosphere, the Y-point may lie at any distance inside the light cylinder. Time-dependent force-free simulations, however, develop closed-line regions that extend all the way to the light cylinder. On the other hand, particle (PIC) solutions consistently develop smaller closed-line regions. In order to understand this effect, we solve the pulsar equation with an improved numerical method. We show that the total electromagnetic energy stored in the ideal force-free magnetosphere manifests a subtle minimum when the closed-line region extends to only 90% of the light cylinder, and thus argue that the system will spontaneously choose this particular configuration. Furthermore, we argue that the intersection of the corotating region with the equatorial current sheet is at right angles, literally leading to a T-point.

We present a comparative study of the effect of low-temperature opacities on stellar models up to the Red Giant branch (RGB), computed with the GARching STellar Evolution Code. We have used two sets of low-temperature opacities; {\AE}SOPUS ({\AE}) from the University of Padova and those from the Wichita State University group (F05). In the relevant range of temperatures for this study, log \k{appa}{\AE} < log \k{appa}F 05. Therefore, to compare stellar evolutionary tracks, we performed a solar calibration of the {\alpha}mlt, for each set of low-temperature opacities. After carrying such a calibration, we find that stellar evolutionary tracks are almost unaffected by the choice of low-temperature opacities, with largest variations of 25-30 K at the latest evolutionary stages of the RGB phase.

Paul C. W. Lai, Matteo Agostini, Foteini Oikonomou, Beatrice Crudele, Ellis R. Owen, Kinwah Wu

8 pages, 2 figures, Presented at the 38th International Cosmic Ray Conference (ICRC2023); PoS (ICRC2023), 1069

In this work, we consider the possible presence of a large population of millisecond pulsars in the Galactic Centre. Their direct detection would be challenging due to severe pulse broadening caused by scattering of radiation. We propose a new method to constrain their population with neutrino imaging of the Galactic Centre. Millisecond pulsars are proposed cosmic-ray accelerators. The high-energy protons they produce will collide with the baryonic matter in the central molecular zone to create charged and neutral pions that decay into neutrinos and $\gamma$-rays, respectively. The specific neutrino and $\gamma$-ray fluxes must be below their corresponding observed values, allowing us to put a conservative upper limit on the millisecond pulsar population of N_MSP < 10,000 within a galacto-centric radius of 20 pc. This upper limit is sensitive to the proton acceleration efficiency of the pulsars, but is less dependent on the particle injection spectral index and the choice of mass tracers. The population will be better constrained when high resolution neutrino observations of the Galactic Centre become available. The presence of these millisecond pulsars can account for the $\gamma$-ray excess in the Galactic Centre.

M. Tavani, G. Piano, A. Bulgarelli, L. Foffano, A. Ursi, F. Verrecchia, C. Pittori, C. Casentini, A. Giuliani, F. Longo, G. Panebianco, A. Di Piano, L. Baroncelli, V. Fioretti, N. Parmiggiani, A. Argan, A. Trois, S. Vercellone, M. Cardillo, L. A. Antonelli, G. Barbiellini, P. Caraveo, P. W. Cattaneo, A. W. Chen, E. Costa, E. Del Monte, G. Di Cocco, I. Donnarumma, Y. Evangelista, M. Feroci, F. Gianotti, C. Labanti, F. Lazzarotto, P. Lipari, F. Lucarelli, M. Marisaldi, S. Mereghetti, A. Morselli, L. Pacciani, A. Pellizzoni, F. Perotti, P. Picozza, M. Pilia, M. Rapisarda, A. Rappoldi, A. Rubini, P. Soffitta, M. Trifoglio, V. Vittorini, F. D'Amico

Accepted for publication in ApJL on September 19, 2023

Gamma-ray emission in the MeV-GeV range from explosive cosmic events is of invaluable relevance to understanding physical processes related to the formation of neutron stars and black holes. Here we report on the detection by the AGILE satellite in the MeV-GeV energy range of the remarkable long-duration gamma-ray burst GRB 221009A. The AGILE onboard detectors have good exposure to GRB 221009A during its initial crucial phases. Hard X-ray/MeV emission in the prompt phase lasted hundreds of seconds, with the brightest radiation being emitted between 200 and 300 seconds after the initial trigger. Very intense GeV gamma-ray emission is detected by AGILE in the prompt and early afterglow phase up to 10,000 seconds. Time-resolved spectral analysis shows time-variable MeV-peaked emission simultaneous with intense power-law GeV radiation that persists in the afterglow phase. The coexistence during the prompt phase of very intense MeV emission together with highly non-thermal and hardening GeV radiation is a remarkable feature of GRB 221009A. During the prompt phase, the event shows spectrally different MeV and GeV emissions that are most likely generated by physical mechanisms occurring in different locations. AGILE observations provide crucial flux and spectral gamma-ray information regarding the early phases of GRB 221009A during which emission in the TeV range was reported.

Hassan Firouzjahi, Antonio Riotto

6 pages, no figure

The abundance of primordial black holes changes in the presence of local non-Gaussianity. A positive non-linear parameter $f_{NL}$ increases the abundance while a negative one reduces it. We show that in non-attractor single-field models of inflation which enhance the curvature power spectrum and may give rise to primordial black holes, $f_{NL}$ is always positive, when computed in correspondence of the peak of the curvature power spectrum where the primordial black hole abundance has its maximum. This implies that the interpretation of the recent pulsar timing arrays data from scalar-induced gravitational waves generated at primordial black hole formation may not be supported by invoking non-Gaussianity within non-attractor single-field models.

Bikash R. Dinda, Narayan Banerjee

19 pages, 3 figures, 1 table. Comments are most welcome

We consider a particular k-essence scalar field model for the late-time cosmic acceleration in which the sound speed, parametrized as $c_s$ is constant. We compute the relevant background and perturbation quantities corresponding to the observables like cosmic microwave background, type Ia supernova, cosmic chronometers, baryon acoustic oscillations, and the $f\sigma_8$. We put constraints on the $c_s^2$ parameter from these observations along with other parameters. We find lower values of $c_s^2$ which are close to zero are tightly constrained. Particularly, we find mean value of $\log_{\rm 10} (c_s^2)$ to be $-0.61$ and $c_s^2 \leq 10^{-3}$ is more than 3$\sigma$ away from this mean value. This means these observations favor a homogeneous dark energy component compared to the clustering one.

Manjeet Kaur (DU), Debottam Nandi (DU), Sharath Raghavan B

14 pages, 6 figures, 1 table

The conventional background solution for the evolution of a single canonical inflaton field performs admirably in extreme scenarios such as the slow-roll phase (where the slow-roll parameter is much less than one) and deep reheating era (where the Hubble parameter is much smaller than the effective mass of the potential and the field oscillates around the minimum of the potential), but fails to accurately depict the dynamics of the Universe near the end of inflation and the initial oscillatory phases. This article proposes a single, unified, model-independent analytical solution for such a model that bridges the gap between these extremes, providing a comprehensive description of the evolution of the Universe. This novel strategy has the potential to substantially enhance both quantitative and qualitative cosmological observational predictions.

K. Hervella Seoane, C. Ramos Almeida, J.A. Acosta Pulido, G. Speranza, C.N. Tadhunter, P.S. Bessiere

26 pages, 9 figures. Accepted for publication in A&A

We present a detailed study of the kinematics of 19 QSO2s in the range 0.3<z<0.41 and [OIII] luminosities $L_{[OIII]} > 10^{8.5}$L$_{\odot}$. We aim at advancing our understanding of the AGN feedback phenomenon by correlating outflow properties with the presence of young stellar populations (YSPs) with ages <100 Myr, the optical morphology and the environment of the galaxies, and the radio luminosity. We characterize the ionized gas kinematics using the [OIII]$\lambda$5007$\r{A}$ profiles, through three different outflow detection methods: multi-component parametric and flux-weighted and peak-weighted non-parametric. We detect ionized outflows in 18 QSO2s using the parametric analysis, and in all of them using the non-parametric methods. We find higher outflow masses using the parametric analysis (log M$_{OF}$(M$_{\odot}$)=6.47$\pm$0.50), and larger mass rates and kinetic powers with the flux-weighted non-parametric method (\.M$_{OF}$=4.0$\pm$4.4 M$_{\odot}$ yr$^{-1}$ and log(\.E$_{kin}$)=41.9$\pm$0.6 erg~s$^{-1}$). However, it is when we use the parametric method and the maximum outflow velocities that we measure the highest outflow mass rates and kinetic energies (23$\pm$35 M$_{\odot}$ yr$^{-1}$ and 42.9$\pm$0.6 erg s$^{-1}$). We do not find any significant correlation between the outflow properties and the previously mentioned galaxy properties. 4 out of 5 QSO2s without a YS<100 Myr show highly disturbed kinematics, whereas only 5 out of the 14 QSO2s with YSPs show similarly asymmetric [OIII] profiles. This might be indicative of negative feedback. The lack of correlation between the outflow properties and the presence of mergers in different interaction stages might be due to their different dynamical timescales. Lastly, the small radio luminosity range covered by our sample may be impeding the detection of any correlation between radio emission and outflow properties.

F. Regnault, N. Al-Haddad, N. Lugaz, C. J. Farrugia, W. Yu, E. E. Davies, A. B. Galvin, B. Zhuang

Accepted for publication at ApJ. 18 pages, 6 figures, 2 tables

In situ measurements from spacecraft typically provide a time series at a single location through coronal mass ejections (CMEs) and they have been one of the main methods to investigate CMEs. CME properties derived from these in situ measurements are affected by temporal changes that occur as the CME passes over the spacecraft, such as radial expansion and ageing, as well as spatial variations within a CME. This study uses multi-spacecraft measurements of the same CME at close separations to investigate both the spatial variability (how different a CME profile is when probed by two spacecraft close to each other) and the so-called ageing effect (the effect of the time evolution on in situ properties). We compile a database of 19 events from the past four decades measured by two spacecraft with a radial separation <0.2 au and an angular separation <10{\deg}. We find that the average magnetic field strength measured by the two spacecraft differs by 18% of the typical average value, which highlights non-negligible spatial or temporal variations. For one particular event, measurements taken by the two spacecraft allow us to quantify and significantly reduce the ageing effect to estimate the asymmetry of the magnetic field strength profile. This study reveals that single-spacecraft time series near 1 au can be strongly affected by ageing and that correcting for self-similar expansion does not capture the whole ageing effect.

Jean Hayoz, Gabriele Cugno, Sascha P. Quanz, Polychronis Patapis, Eleonora Alei, Markus J. Bonse, Felix A. Dannert, Emily O. Garvin, Timothy D. Gebhard, Björn S. Konrad, Lia F. Sartori

The investigation of the atmospheres of closely separated, directly imaged gas giant exoplanets is challenging due to the presence of stellar speckles that pollute their spectrum. To remedy this, the analysis of medium- to high-resolution spectroscopic data via cross-correlation with spectral templates (cross-correlation spectroscopy) is emerging as a leading technique. We aim to define a robust Bayesian framework combining, for the first time, three widespread direct-imaging techniques, namely photometry, low-resolution spectroscopy, and medium-resolution cross-correlation spectroscopy in order to derive the atmospheric properties of close-in directly imaged exoplanets. Our framework CROCODILE (cross-correlation retrievals of directly imaged self-luminous exoplanets) naturally combines the three techniques by adopting adequate likelihood functions. To validate our routine, we simulated observations of gas giants similar to the well-studied $\beta$~Pictoris~b planet and we explored the parameter space of their atmospheres to search for potential biases. We obtain more accurate measurements of atmospheric properties when combining photometry, low- and medium-resolution spectroscopy into atmospheric retrievals than when using the techniques separately as is usually done in the literature. We find that medium-resolution ($R \approx 4000$) K-band cross-correlation spectroscopy alone is not suitable to constrain the atmospheric properties of our synthetic datasets; however, this problem disappears when simultaneously fitting photometry and low-resolution ($R \approx 60$) spectroscopy between the Y and M bands. Our framework allows the atmospheric characterisation of directly imaged exoplanets using the high-quality spectral data that will be provided by the new generation of instruments such as VLT/ERIS, JWST/MIRI, and ELT/METIS.

A. Armeni, B. Stelzer, R. A. B. Claes, C. F. Manara, A. Frasca, J. M. Alcalá, F. M. Walter, Á. Kóspál, J. Campbell-White, M. Gangi, K. Mauco, L. Tychoniec

15 pages, 14 figures. Accepted for publication in A&A

HM Lup is a young M-type star that accretes material from a circumstellar disk through a magnetosphere. Our aim is to study the inner disk structure of HM Lup and to characterize its variability. We used spectroscopic data from HST/STIS, X-Shooter, and ESPRESSO taken in the framework of the ULLYSES and PENELLOPE programs, together with photometric data from TESS and AAVSO. The 2021 TESS light curve shows variability typical for young stellar objects of the "accretion burster" type. The spectra cover the temporal evolution of the main burst in the 2021 TESS light curve. We compared the strength and morphology of emission lines from different species and ionization stages. We determined the mass accretion rate from selected emission lines and from the UV continuum excess emission at different epochs, and we examined its relation to the photometric light curves. The emission lines in the optical spectrum of HM Lup delineate a temperature stratification along the accretion flow. While the wings of the H I and He I lines originate near the star, the lines of species such as Na I, Mg I, Ca I, Ca II, Fe I, and Fe II are formed in an outer and colder region. The shape and periodicity of the 2019 and 2021 TESS light curves, when qualitatively compared to predictions from magnetohydrodynamic models, suggest that HM Lup was in a regime of unstable ordered accretion during the 2021 TESS observation due to an increase in the accretion rate. Although HM Lup is not an extreme accretor, it shows enhanced emission in the metallic species during this high accretion state that is produced by a density enhancement in the outer part of the accretion flow.

R Ignace, J E Bjorkman, A-N Chene, C Erba, L Fabiani, A F J Moffat, R Sincennes, N St-Louis

This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Massive star winds are structured both stochastically ("clumps") and often coherently (Co-rotation Interaction Regions, or CIRs). Evidence for CIRs threading the winds of Wolf-Rayet (WR) stars arises from multiple diagnostics including linear polarimetry. Some observations indicate changes in polarization position angle across optical recombination emission lines from a WR star wind but limited to blueshifted Doppler velocities. We explore a model involving a spherical wind with a single conical CIR stemming from a rotating star as qualitative proof-of-concept. To obtain a realistic distribution of limb polarization and limb darkening across the pseudo-photosphere formed in the optically thick wind of a WR star, we used Monte Carlo radiative transfer (MCRT). Results are shown for a parameter study. For line properties similar to WR 6 (EZ CMa; HD 50896), the combination of the MCRT results, a simple model for the CIR, and the Sobolev approximation for the line formation, we were able to reproduce variations in both polarization amplitude and position angle commensurate with observations. Characterizing CIRs in WR~winds has added importance for providing stellar rotation periods since the v sin i values are unobtainable because the pseudo-photosphere forms in the wind itself.

Giulia Ciacci, Andrea Barucci, Sara Di Ruzza, Elisa Maria Alessi

In this work, we explore how to classify asteroids in co-orbital motion with a given planet using Machine Learning. We consider four different kinds of motion in mean motion resonance with the planet, nominally Tadpole, Horseshoe and Quasi-satellite, building 3 datasets defined as Real (taking the ephemerides of real asteroids from the JPL Horizons system), Ideal and Perturbed (both simulated, obtained by propagating initial conditions considering two different dynamical systems) for training and testing the Machine Learning algorithms in different conditions. The time series of the variable theta (angle related to the resonance) are studied with a data analysis pipeline defined ad hoc for the problem and composed by: data creation and annotation, time series features extraction thanks to the tsfresh package (potentially followed by selection and standardization) and the application of Machine Learning algorithms for Dimensionality Reduction and Classification. Such approach, based on features extracted from the time series, allows to work with a smaller number of data with respect to Deep Learning algorithms, also allowing to define a ranking of the importance of the features. Physical Interpretability of the features is another key point of this approach. In addition, we introduce the SHapley Additive exPlanations for Explainability technique. Different training and test sets are used, in order to understand the power and the limits of our approach. The results show how the algorithms are able to identify and classify correctly the time series, with a high degree of performance.

Xinwu Cao (1), Bei You (2), Xing Wei (3) (1.Zhejiang Univ., 2. Wuhan Univ., 3. Beijing Norm. Univ.)

11 pages, accepted by MNRAS

1ES 1927+654 was known as a type 2 Seyfert galaxy, which exhibited drastic variability recently in ultraviolet (UV)/optical and X-ray bands. An UV/optical outburst was observed in the end of 2017, and it reached the peak luminosity $\sim 50$ days later. The high-cadence observations showed a rapid X-ray flux decline with complete disappearance of the power-law hard X-ray component when the soft X-ray thermal emission reached its lowest level about $150$ days after the UV/optical peak. The power law X-ray component reappeared with thermal X-ray emission brightening from its lowest flux within next $\sim$ 100~days. We assume an episodic accretion event taking place in the outer region of the disc surrounding a central black hole (BH), which is probably due to a red giant star tidally disrupted by the BH. The inner thin disc with corona is completely swept by the accretion event when the gas reaches the innermost circular stable orbit. The field threading the disrupted star is dragged inwards by the disc formed after the tidal disruption event, which accelerates outflows from the disc. The disc dimmed since a large fraction of the energy released in the disc is tapped into the outflows. The accretion rate of the episodic accretion event declines, and ultimately it turns out to be a thin disc, which is inefficient for field advection, and the outflows are switched off. A thin disc with corona reappears later after the outburst.

Winds of massive stars have density inhomogeneities (clumping) that may affect the formation of spectral lines in different ways, depending on their formation region. Most of previous and current spectroscopic analyses have been performed in the optical or ultraviolet domain. However, massive stars are often hidden behind dense clouds rendering near-infrared observations necessary. Our objective is to investigate whether a spectroscopic analysis using either optical or infrared observations results in the same stellar parameters with comparable accuracy, and whether clumping affects them in different ways. We analyzed optical and near-infrared observations of a set of massive O stars with spectral types O4-O9.5 and all luminosity classes. We obtain similar stellar parameters in the optical and the infrared, although with larger uncertainties in the near-infrared, both with and without clumping, albeit with some individual deviating cases. We find that the inclusion of clumping improves the fit to H$_\alpha$ or HeII 4686 in the optical for supergiants, as well as that of Br$_\gamma$ in the near-infrared, but it sometimes worsens the fit to HeII 2.18$\mu$m. Globally, there are no significant differences when using the clumping laws tested in this work. The infrared can be used for spectroscopic analyses, giving similar parameters as from the optical, though with larger uncertainties. The best fits to different lines are obtained with different (linear) clumping laws, indicating that the wind structure may be more complex than adopted in the present work. No clumping law results in a better global fit, or improves the consistency between optical and infrared stellar parameters. Our work shows that the optical and infrared lines are not sufficient to break the dichotomy between the mass-loss rate and clumping factor.

Laura Marcelli (for the JEM-EUSO Collaboration)

18 pages, 10 figures, ICRC2023

Mini-EUSO is a telescope observing the Earth in the ultraviolet band (290-430 nm) since 2019, through a nadir-facing UV-transparent window in the Russian Zvezda module of the International Space Station. The main camera has an optical system composed of two 25 cm diameter Fresnel lenses and a focal surface consisting of 36 multi-anode photomultiplier tubes, 64 pixels each, for a total of 2304 channels. The instrument has a square field of view with a side of 44 degrees, a spatial resolution of about 6.3 km on the Earth surface and a sampling time of 2.5 microseconds. Mini- EUSO has also two cameras in the near infrared and visible ranges and silicon photomultiplier sensors to complement the UV observations. Mini-EUSO has been designed as a small-size version of the original JEM-EUSO space telescope to demonstrate its observational principle. Mini-EUSO is in fact potentially capable of observing extensive air showers generated by ultra-high-energy cosmic rays with an energy above 10 21 eV and of detecting artificial showers generated with lasers from the ground. Other main scientific objectives of the mission are the study of atmospheric phenomena (transient luminous events such as ELVES and sprites), the observation of meteors and among them the search for interstellar meteors and nuclearites such as strange quark matter. Moreover, Mini-EUSO can map night-time UV Earth emissions, both anthropogenic and natural. In this work, we will discuss results and performance of the telescope during its first four years of activity.

F. Sainsbury-Martinez, P. Tremblin, M. Mancip, S. Donfack, E. Honore, M. Bourenane

Accepted for publication in ApJ

In-order to understand the results of recent observations of exoplanets, models have become increasingly complex. Unfortunately this increases both the computational cost and output size of said models. We intend to explore if AI-image-recognition can alleviate this burden. We used DYNAMICO to run a series of HD209458-like models with different orbital-radii. Training data for a number of features of interest was selected from the initial outputs of these models. This was used to train a pair of multi-categorisation convolutional-neural-networks (CNN), which we applied to our outer-atmosphere-equilibrated models. The features detected by our CNNs revealed that our models fall into two regimes: models with a shorter orbital-radii exhibit significant global mixing which shapes the entire atmospheres dynamics. Whereas, models with longer orbital-radii exhibit negligible mixing except at mid-pressures. Here, the initial non-detection of any trained features revealed a surprise: a night-side hot-spot. Analysis suggests that this occurs when rotational influence is sufficiently weak that divergent flows from the day-side to the night-side dominate over rotational-driven transport, such as the equatorial jet. We suggest that image-classification may play an important role in future, computational, atmospheric studies. However special care must be paid to the data feed into the model, from the colourmap, to training the CNN on features with enough breadth and complexity that the CNN can learn to detect them. However, by using preliminary-studies and prior-models, this should be more than achievable for future exascale calculations, allowing for a significant reduction in future workloads and computational resources.

Björn J. R. Davidsson

20 pages, 11 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Planetesimals in the primordial disc may have experienced a collisional cascade. If so, the comet nuclei later placed in the Kuiper belt, scattered disc, and Oort Cloud would primarily be fragments and collisional rubble piles from that cascade. However, the heating associated with the collisions cannot have been strong enough to remove the hypervolatiles that are trapped within more durable ices, because comet nuclei are rich in hypervolatiles. This places constraints on the diameter of the largest bodies allowed to participate in collisional cascades, and limits the primordial disc lifetime or population size. In this paper, the thermophysical code NIMBUS is used to study the thermal evolution of planetesimals before, during, and after catastrophic collisions. The loss of CO during segregation of $\mathrm{CO_2:CO}$ mixtures and during crystallisation of amorphous $\mathrm{H_2O}$ is calculated, as well as mobilisation and internal relocation of $\mathrm{CO_2}$. If an amorphous $\mathrm{H_2O}$ host existed, and was protected by a $\mathrm{CO_2:CO}$ heat sink, only diameter $D<20\,\mathrm{km}$ (inner disc) and $D<64\,\mathrm{km}$ (outer disc) bodies could have been involved in a collisional cascade. If $\mathrm{CO_2}$ was the only CO host, the critical diameters drop to $D<20$-$32\mathrm{km}$. Avoiding disruption of larger bodies requires a primordial disc lifetime of $<9\,\mathrm{Myr}$ at $15\,\mathrm{au}$ and $<50$-$70\,\mathrm{Myr}$ at $30\,\mathrm{au}$. Alternatively, if a $450\,\mathrm{Myr}$ disc lifetime is required to associate the primordial disc disruption with the Late Heavy Bombardment, the disc population size must have been 6-60 times below current estimates.

Klaus Werner, Nicole Reindl, Roberto Raddi, Massimo Griggio, Luigi R. Bedin, María E. Camisassa, Alberto Rebassa-Mansergas, Santiago Torres, Peter Goodhew

Accepted for publication in A&A

Planetary nebulae in Galactic open star clusters are rare objects; only three are known to date. They are of particular interest because their distance can be determined with high accuracy, allowing one to characterize the physical properties of the planetary nebula and its ionizing central star with high confidence. Here we present the first quantitative spectroscopic analysis of a central star in an open cluster, namely the faint nucleus of IPHASX J055226.2$+$323724 in M37. This cluster contains 14 confirmed white dwarf members, which were previously used to study the initial-to-final-mass relation of white dwarfs, and six additional white dwarf candidates. We performed an atmosphere modeling of spectra taken with the 10m Gran Telescopio Canarias. The central star is a hot hydrogen-deficient white dwarf with an effective temperature of 90,000 K and spectral type PG1159 (helium- and carbon-rich). We know it is about to transform into a helium-rich DO white dwarf because the relatively low atmospheric carbon abundance indicates ongoing gravitational settling of heavy elements. The star belongs to a group of hot white dwarfs that exhibit ultrahigh-excitation spectral lines possibly emerging from shock-heated material in a magnetosphere. We find a relatively high stellar mass of $M= 0.85^{+0.13}_{-0.14}$ M$_\odot$. This young white dwarf is important for the semi-empirical initial-final mass relation because any uncertainty related to white-dwarf cooling theory is insignificant with respect to the pre-white-dwarf timescale. Its post-asymptotic-giant-branch age of $170,000-480,000$ yr suggests that the extended planetary nebula is extraordinarily old. We also performed a spectroscopic analysis of the six other white dwarf candidates of M37, confirming one as a cluster member.

Arshia Maria Jacob

Invited contribution, 44 pages, 16 figures, 3 tables, accepted for publication in Astrophysics and Space Sciences

Formed at an early stage of gas-phase ion-molecule chemistry, hydrides -- molecules containing a heavy element covalently bonded to one or more hydrogen atoms -- play an important role in interstellar chemistry as they are the progenitors of larger and more complex species in the interstellar medium. In recent years, the careful analysis of the spectral signatures of hydrides have led to their use as tracers of different constituents, and phases of the interstellar medium and in particular the more diffuse environments. Diffuse clouds form an essential link in the stellar gas life-cycle as they connect both the late and early stages of stellar evolution. As a result, diffuse clouds are continuously replenished by material which makes them reservoirs for heavy elements and hence ideal laboratories for the study of astrochemistry. This review will journey through a renaissance of hydride observations detailing puzzling hydride discoveries and chemical mysteries with special focus carbon-bearing hydrides to demonstrate the big impact of these small molecules and ending with remarks on the future of their studies.

Carla Arce-Tord, Simon Casassus, William R. F. Dent, Sebastián Pérez, Miguel Cárcamo, Philipp Weber, Natalia Engler, Lucas A. Cieza, Antonio Hales, Alice Zurlo, Sebastian Marino

accepted in MNRAS

Warps have often been used to explain disc properties, but well characterised examples are important due to their role in disc evolution. Scattered light images of discs with central gaps have revealed sharp warps, such that the outer rings are shadowed by tilted inner discs. The near-IR intensity drops along the ring around TTauri star DoAr44 have been interpreted in terms of a central warp. We report new ALMA observations of DoAr44 in the continuum at 230 GHz and 350 GHz (at ~10 au), along with a new epoch of SPHERE/IRDIS differential polarised imaging taken during excellent weather conditions. The ALMA observations resolve the ring and confirm the decrements proposed from deconvolution of coarse 336 GHz data. The scattered light image constrains the dips, which correspond to a misaligned inner disc with a relative inclination $\xi$ = 21.4 $^{+6.7}_{-8.3}$ deg. The SPHERE intensity profile shows a morphological change compared to a previous epoch that may be interpreted as a variable orientation of the inner disc, from $\xi$ ~30 deg to $\xi$ ~20 deg. The intensity dips probably correspond to temperature decrements, as their mm-spectral index, $\alpha^{230 GHz}_{350 GHz}$ ~2.0 $\pm$ 0.1, is indicative of optically thick emission. The azimuth of the two temperature decrements are leading clockwise relative to the IR-dips, by $\eta$ = 14.95 deg and $\eta$ = 7.92 deg. For a retrograde disc, such shifts are expected from a thermal lag and imply gas surface densities of $\Sigma_g$ = 117 $\pm$ 10 g/cm$^2$ and $\Sigma_g$ = 48 $\pm$ 10 g/cm$^2$. A lopsided disc, with contrast ratio $f_r$=2.4 $\pm$ 0.5, is also consistent with the large continuum crescent.

Eric C. Bellm, Yuankun Wang, Jan van Roestel, Rebecca A. Phillipson, Michael W. Coughlin, John A. Tomsick, Steven L. Groom, Brian Healy, Josiah Purdum, Ben Rusholme, Jesper Sollerman, Peter Bealo, Stefano Lora, Eddy Muyllaert, Ivo Peretto, Erik J. Schwendeman

12 pages, 6 figures, accepted for publication in ApJ

Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolometric luminosity was only $2\times 10^{35}$ erg s$^{-1}$, implying an Eddington fraction of about $3\times10^{-4}$. The source remained in the hard spectral state throughout the outburst. From optical spectroscopy measurements we estimate an outer disk radius of 10$^{11}$ cm. The low observed X-ray luminosity is not sufficient to irradiate the entire disk, but we observe a surprising exponential decline in the X-ray lightcurve. These observations highlight the potential of optical and infrared (O/IR) synoptic surveys to discover low-luminosity activity from X-ray transients.

M. Albekioni, T.V. Zaqarashvili, V. Kukhianidze, E. Gurgenashvili, P. Bourdin

Observations by recent space missions reported the detection of Rossby waves (r-modes) in light curves of many stars (mostly A, B, and F spectral types) with outer radiative envelope. This paper aims to study the theoretical dynamics of Rossby-type waves in such stars. Hydrodynamic equations in a rotating frame were split into horizontal and vertical parts connected by a separation constant (or an equivalent depth). Vertical equations were solved analytically for a linear temperature profile and the equivalent depth was derived through free surface boundary condition. It is found that the vertical modes are concentrated in the near-surface layer with a thickness of several tens of surface density scale height. Then with the equivalent width, horizontal structure equations were solved, and the corresponding dispersion relation for Rossby, Rossby-gravity, and inertia-gravity waves was obtained. The solutions were found to be confined around the equator leading to the equatorially trapped waves. It was shown that the wave frequency depends on the vertical temperature gradient as well as on stellar rotation. Therefore, observations of wave frequency in light curves of stars with known parameters (radius, surface gravity, rotation period) could be used to estimate the temperature gradient in stellar outer layers. Consequently, the Rossby mode may be considered as an additional tool in asteroseismology apart from acoustic and gravity modes.

Relative astrometric shifts between multiply lensed images provide a valuable tool to investigate haloes in the intergalactic space. In strong lens systems in which a single lens plays the primary role in producing multiple images, the gravitational force exerted by line-of-sight (LOS) haloes can slightly change the relative positions of multiply lensed images produced by the dominant lens. In such cases, a LOS halo positioned sufficiently far from the dominant lens along the LOS can create a pattern in the reduced deflection angle that corresponds to the B-mode (magnetic or divergence-free mode). By measuring both the B-mode and E-mode (electric or rotation-free mode), we can determine the LOS distance ratios, as well as the 'bare' convergence and shear perturbations in the absence of the dominant lens. However, scale variations in the distance ratio lead to mass-sheet transformations in the background lens plane, introducing some uncertainty in the distance ratio estimation. This uncertainty can be significantly reduced by measuring the time delays between the lensed images. Additionally, if we obtain the redshift values of both the dominant and perturbing haloes, along with the time delays between the multiply lensed images that are affected by the haloes, the B-mode can break the degeneracy related to mass-sheet transformations in both the foreground and background lens planes. Therefore, measuring the astrometric lensing B-mode has the potential to substantially decrease the uncertainty in determining the Hubble constant.

Numerous simulations indicate that a large number of subhalos should be hosted by the Milky Way. The potential existence of a nearby subhalo could have important implications for our understanding of dark matter (DM) annihilation. In this study, we investigate the hypothetical presence of a nearby subhalo and set the upper limits on the DM annihilation cross section by analyzing the cosmic-ray antiproton spectrum. By presenting the ratios of annihilation cross section limits for scenarios with and without a nearby subhalo, we can quantitatively evaluate the potential impact of the nearby subhalo on the limits of the DM annihilation cross section. The impacts of the concentration model and the subhalo probability distribution have been considered. We explore the antiproton contribution of the potential nearby DM subhalo accounting for the DAMPE $e^\pm$ spectrum at $\sim 1.4$ TeV and find that the current AMS-02 antiproton results do not place the constraint on this contribution.

Isabella A. Gerrard, Christoph Federrath, Nickolas M. Pingel, Naomi M. McClure-Griffiths, Antoine Marchal, Gilles Joncas, Susan E. Clark, Snežana Stanimirović, Min-Young Lee, Jacco Th. van Loon, John Dickey, Helga Dénes, Yik Ki Ma, James Dempsey, Callum Lynn

20 pages, 16 figures, accepted to MNRAS

Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($\sigma_{\rho/\rho_0}$) to the turbulent sonic Mach number ($\mathcal{M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterised by the turbulence driving parameter: $b=\sigma_{\rho/\rho_0}/\mathcal{M}$. The turbulence driving parameter ranges from $b=1/3$ (purely solenoidal) to $b=1$ (purely compressive), with $b=0.38$ characterising the natural mixture (1/3~compressive, 2/3~solenoidal) of the two driving modes. Here we present a new method for recovering $\sigma_{\rho/\rho_0}$, $\mathcal{M}$, and $b$, from observations on galactic scales, using a roving kernel to produce maps of these quantities from column density and centroid velocity maps. We apply our method to high-resolution HI emission observations of the Small Magellanic Cloud (SMC) from the GASKAP-HI survey. We find that the turbulence driving parameter varies between $b\sim 0.3$ and $b\sim 1.0$ within the main body of the SMC, but the median value converges to $b\sim0.51$, suggesting that the turbulence is overall driven more compressively ($b>0.38$). We observe no correlation between the $b$ parameter and HI or H$\alpha$ intensity, indicating that compressive driving of HI turbulence cannot be determined solely by observing HI or H$\alpha$ emission density, and that velocity information must also be considered. Further investigation is required to link our findings to potential driving mechanisms such as star-formation feedback, gravitational collapse, or cloud-cloud collisions.

Stream interaction regions (SIRs) are often thought to be responsible for the generation of suprathermal population in the interplanetary medium. Despite the source being same, wide variations in spectral indices of suprathermal populations are observed at 1 au during SIRs. This poses significant uncertainty in understanding the generation of suprathermal ion populations by SIRs and indicates interplay of multiple source mechanisms. In the present work, by analyzing variations in suprathermal 4He, O, and Fe for 20 SIR events recorded by STEREO-A during 2007 - 2014, we find that the spectral indices of these elements vary in the range of 2.06-4.08, 1.85-4.56, and 2.11-4.04 respectively for 19 events. However, in one special case, all the three suprathermal elements show nearly identical (1.5) spectral indices. We offer possible mechanisms, which can cause significant variations in the spectral indices of suprathermal particles. More importantly, we show the possible role of merging and/or contraction of small-scale magnetic islands near 1 au in producing nearly identical spectral indices for three different elements with different first ionization potential and mass-to-charge ratio. The occurrence of these magnetic islands near 1 au also supports the minimal modulation in spectral indices of these particles. The role of a possible solar flare in generating these magnetic islands near the heliospheric current sheet is also suggested.

A Thorne-\.{Z}ytkow object can be formed when a neutron star is absorbed by the envelope of its giant companion star, spirals toward the center of the giant star due to the drag from the surrounding envelope, and merges with the star's core. During in-spiral, dynamical friction is the main drag force acting on a neutron star inside the star envelope. However, dynamical friction is caused by the neutron star's gravitational interaction with its own gravitationally induced wake. Therefore, exactly the same gravitational drag force will act on a mirror neutron star, and if a significant part of the dark matter is in the form of mirror matter, then we expect the formation of Thorne-\.{Z}ytkow objects with a mirror neutron star inside.

Doyeon Avery Kim, Susan E Clark, Mary E Putman, Larry Li

15 pages, 16 figures

We characterize the kinematic and magnetic properties of HI filaments located in a high Galactic latitude region ($165^\circ < \alpha < 195^\circ$ and $12^\circ < \delta < 24^\circ$). We extract three-dimensional filamentary structures using \texttt{fil3d} from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) survey 21-cm emission data. Our algorithm identifies coherent emission structures in neighboring velocity channels. Based on the mean velocity, we identify a population of local and intermediate velocity cloud (IVC) filaments. We find the orientations of the local (but not the IVC) HI filaments are aligned with the magnetic field orientations inferred from Planck 353 GHz polarized dust emission. We analyze position-velocity diagrams of the velocity-coherent filaments, and find that only 15 percent of filaments demonstrate significant major-axis velocity gradients with a median magnitude of 0.5 km s$^{-1}$ pc$^{-1}$, assuming a fiducial filament distance of 100 pc. We conclude that the typical diffuse HI filament does not exhibit a simple velocity gradient. The reported filament properties constrain future theoretical models of filament formation.

V. Loaiza-Tacuri, K. Cunha, D. Souto, V. V. Smith, R. Guerço, C. Chiappini, J. V. Sales Silva, D. Horta, C. Allende Prieto, R. Beaton, D. Bizyaev, S. Daflon, P. Frinchaboy, S. Hasselquist, C. R. Hayes, J. A. Holtzman, H. Jönsson, S. R. Majewski, S. Mészáros, D. L. Nidever, M. Pinsonneault, G. Zasowski

17 pages, 10 figures, 6 tables

Previous results in the literature have found the young inner-disk open cluster NGC 6705 to be mildly $\alpha$-enhanced. We examined this possibility via an independent chemical abundance analysis for 11 red-giant members of NGC 6705. The analysis is based on near-infrared APOGEE spectra and relies on LTE calculations using spherical model atmospheres and radiative transfer. We find a mean cluster metallicity of $\rm [Fe/H] = +0.13 \pm 0.04$, indicating that NGC 6705 is metal-rich, as may be expected for a young inner-disk cluster. The mean $\alpha$-element abundance relative to iron is $\rm \langle [\alpha/Fe]\rangle =-0.03 \pm 0.05$, which is not at odds with expectations from general Galactic abundance trends. NGC 6705 also provides important probes for studying stellar mixing, given its turn-off mass of M$\sim$3.3 M$_\odot$. Its red giants have low $^{12}$C abundances ([$^{12}$C/Fe]=$-$0.16) and enhanced $^{14}$N abundances ([$^{14}$N/Fe]=+0.51), which are key signatures of the first dredge-up on the red giant branch. An additional signature of dredge-up was found in the Na abundances, which are enhanced by [Na/Fe]=+0.29, with a very small non-LTE correction. The $^{16}$O and Al abundances are found to be near-solar. All of the derived mixing-sensitive abundances are in agreement with stellar models of approximately 3.3 M$_{\odot}$ evolving along the red giant branch and onto the red clump. As found in young open clusters with similar metallicities, NGC 6705 exhibits a mild excess in the s-process element cerium, with $\rm [Ce/Fe] = +0.13\pm0.07$.

K. Zioutas, V. Anastassopoulos, A. Argiriou, G. Cantatore, S. Cetin, H. Fischer, A. Gardikiotis, H. Haralambous, D. H. H. Hoffmann, S. Hofmann, M. Karuza, A. Kryemadhi, M. Maroudas, A. Mastronikolis, C. Oikonomou, K. Ozbozduman, Y. K. Semertzidis

8 pages, 7 figures, Published in the proceedings of the "15th International Conference on Meteorology, Climatology and Atmospheric Physics (COMECAP 2021)" see this https URL

We are investigating the possible origin of small-scale anomalies, like the annual stratospheric temperature anomalies. Unexpectedly within known physics, their observed planetary "dependency", does not match concurrent solar activity, whose impact on the atmosphere is unequivocal; this points at an additional energy source of exo-solar origin. A viable concept behind such observations is based on possible gravitational focusing by the Sun and its planets towards the Earth of low-speed invisible streaming matter; its influx towards the Earth gets temporally enhanced. Only a somehow "strongly" interacting invisible streaming matter with the small upper atmospheric screening can be behind the observed temperature excursions. Ordinary dark matter (DM) candidates like axions or WIMPs, cannot have any noticeable impact. The associated energy deposition is $\mathcal{O}(\sim 1000\, \mathrm{GeV}/{{\mathrm{cm}}^2}/\mathrm{sec})$. The atmosphere has been uninterruptedly monitored for decades. Therefore, the upper atmosphere can serve as a novel (low-threshold) detector for the dark Universe, with built-in spatiotemporal resolution while the solar system gravity acts temporally as a signal amplifier. Interestingly, the anomalous ionosphere shows a relationship with the inner earth activity like earthquakes. Similarly investigating the transient sudden stratospheric warmings within the same reasoning, the nature of the assumed "invisible streams" could be deciphered.

Amit Mondal, Biswajit Pandey

15 pages, 9 figures, comments welcome

We analyze the spatial anisotropy and the velocity anisotropy in a set of mock stellar halos from the Aquarius simulations. The spatial anisotropy in each mock stellar halo rises progressively with the increasing distance from the halo centre, eventually reaching a maximum near the periphery. Excluding the bound satellites leads to a significant reduction of the spatial anisotropy in each halo. We compare the measured anisotropy in the mock stellar halos with that from their sphericalized versions where all the shape and substructure induced anisotropies are erased. The growth of spatial anisotropy persists throughout the entire halo when the bound satellites are present but remains limited within the inner halo ($<60 \, h^{-1}\, {\rm kpc}$) after their exclusion. This indicates that the spatial anisotropy in the inner halo is induced by the diffuse substructures and the halo shape whereas the outer halo anisotropy is dominated by the bound satellites. We find that the outer parts of the stellar halo are kinematically colder than the inner regions. The stellar orbits are predominantly radial but they become rotationally dominated at certain radii that are marked by the prominent $\beta$ dips. Most of the $\beta$ dips disappear after the removal of the satellites. A few shallow and broad $\beta$ dips arise occasionally due to the presence of diffuse streams and clouds. Our analysis suggests that a combined study of the spatial and velocity anisotropies can reveal the structure and the assembly history of the stellar halos.

Stellar rotation is crucial for studying stellar evolution since it provides information about age, angular momentum transfer, and magnetic fields of stars. In the case of the Sun, due to its proximity, detailed observation of sunspots at various latitudes and longitudes allows the precise estimate of the solar rotation period and its differential rotation. Here, we present for the first time an analysis of stellar differential rotation using starspot transit mapping as a means of detecting differential shear in solar-type and M stars. The aim of this study is to investigate the relationship between rotational shear, $\Delta\Omega$, with both the star's effective temperature ($T_{\text{eff}}$) and average rotation period ($P_{\text{r}}$). We present differential rotation profiles derived from previously collected spot transit mapping data for 13 slowly rotating stars ($P_{\text{rot}} \geq 4.5$ days), with spectral types ranging from M to F, which were observed by the Kepler and CoRoT satellites. Our findings reveal a significant negative correlation between rotational shear and the mean period of stellar rotation (correlation coefficient of -0.77), which may be an indicator of stellar age. On the other hand, a weak correlation was observed between differential rotation and the effective temperature of the stars. Overall, the study provides valuable insights into the complex relationship between stellar parameters and differential rotation, which may enhance our understanding of stellar evolution and magnetic dynamos.

Flux transfer events (FTEs) are a type of magnetospheric phenomena that exhibit distinctive observational signatures from the in-situ spacecraft measurements across the Earth's magnetopause. They are generally believed to possess a magnetic field configuration of a magnetic flux rope and formed through magnetic reconnection at the dayside magnetopause, sometimes accompanied with enhanced plasma convection in the ionosphere. We examine two FTE events under the condition of southward interplanetary magnetic field (IMF) with a dawn-dusk component at the magnetopause by applying the Grad-Shafranov (GS) reconstruction method to the in-situ measurements by the Magnetospheric Multiscale (MMS) spacecraft to derive the magnetic flux contents associated with the FTE flux ropes. In particular, given a cylindrical magnetic flux rope configuration derived from the GS reconstruction, the magnetic flux content can be characterized by both the toroidal (axial) and poloidal fluxes. We then estimate the amount of magnetic flux (i.e., the reconnection flux) encompassed by the area ``opened" in the ionosphere, based on the ground-based Super Dual Auroral Radar Network (SuperDARN) observations. We find that for event 1, the FTE flux rope is oriented in the approximate dawn-dusk direction, and the amount of its poloidal magnetic flux agrees with the corresponding reconnection flux. For event 2, the agreement among the estimates of the magnetic fluxes is uncertain. We provide a detailed description about our interpretation for the topological features of the FTE flux ropes, based on a formation scenario of sequential magnetic field reconnection between adjacent field lines, consistent with our results.

Maxim Pospelov, Anupam Ray

13 pages, 4 figures. Comments are welcome

A sub-component of dark matter with a short collision length compared to a planetary size leads to efficient accumulation of dark matter in astrophysical bodies. We analyze possible neutrino signal from the annihilation of such dark matter and conclude that in the optically thick regime for dark matter capture, the Earth provides the largest neutrino flux. Using the results of the existing searches, we consider two scenarios for neutrino flux, from stopped mesons and prompt higher-energy neutrinos. In both cases we exclude some previously unexplored parts of the parameter space (dark matter mass, its abundance, and the scattering cross section on nuclei) by recasting the existing neutrino searches.

Jordan Koechler

Contribution to the 34th Rencontres de Blois on Particle Physics and Cosmology (Blois 2023), the XVIII International Conference on Topics in Astroparticle and Underground Physics (TAUP 2023) and TeV Particle Astrophysics (TeVPA) 2023

We present updated constraints on 'light' Dark Matter (DM) particles with masses between 1 MeV and 5 GeV. In this range, we can expect DM-produced $e^\pm$ pairs to up-scatter low-energy ambient photons in the Milky Way via the Inverse Compton process, and produce a flux of X-rays that can be probed by a range of space observatories. Using diffuse X-ray data from XMM-Newton, INTEGRAL, NuSTAR and Suzaku, we compute the strongest constraints to date on annihilating DM for 200 MeV $< m_{\rm DM} <$ 5 GeV and decaying DM for 100 MeV $< m_{\rm DM} <$ 5 GeV.

Andrew Ross Wilson, Massimiliano Vasile

18 pages, 3 figures

The recent growth of the space sector, spurred by a surge in private actors, has led to a sharp increase in our ability to address societal challenges through space data. However, this has exacerbated an already critical situation in space: the proliferation of space debris and a critical expansion of space traffic which is leading to high levels of orbital congestion. In parallel, increased levels of spacecraft production and orbital launches are also heightening the environmental footprint of the sector. This might lead to a paradoxical situation whereby the use of space to support the Sustainable Development Goals (SDGs) becomes unsustainable from the perspective of both the Earth and space environment. This situation can be described as the space sustainability paradox. This paper presents this concept for the first time and argues that existing policies and remediation actions are not a long-term sustainable solution for tackling this issue and may actually intensify the problem. This places an added importance upon addressing space sustainability in a more coherent, strategic and responsible manner, potentially based on the doughnut economic model of social and planetary boundaries. Doing so may prevent the sector from falling victim to a tragedy of the commons type of scenario and avoid negative trends from becoming the norm. As a result, this would ensure that outer space can continue to be used by future generations to address societal challenges, without priming severe and enduring damage to the Earth and space environment.

Bradley W. Hindman, Keith Julien

18 pages, 2 figures, submitted to the Astrophysical Journal

Starting from the fully compressible fluid equations in a plane-parallel atmosphere, we demonstrate that linear internal gravity waves are naturally pseudo-incompressible in the limit that the wave frequency $\omega$ is much less than that of surface gravity waves, i.e., $\omega \ll \sqrt{g k_h}$ where $g$ is the gravitational acceleration and $k_h$ is the horizontal wavenumber. We accomplish this by performing a formal expansion of the wave functions and the local dispersion relation in terms of a dimensionless frequency $\varepsilon = \omega / \sqrt{g k_h}$. Further, we show that in this same low-frequency limit, several forms of the anelastic approximation, including the Lantz-Braginsky-Roberts (LBR) formulation, poorly reproduce the correct behavior of internal gravity waves. The pseudo-incompressible approximation is achieved by assuming that Eulerian fluctuations of the pressure are small in the continuity equation. Whereas, in the anelastic approximation Eulerian density fluctuations are ignored. In an adiabatic stratification, such as occurs in a convection zone, the two approximations become identical. But, in a stable stratification, the differences between the two approximations are stark and only the pseudo-incompressible approximation remains valid.

The analyses for the NICER data imply $R_{2.0M_\odot}=12.41^{+1.00}_{-1.10}$ km and $R_{1.4M_\odot}=12.56^{+1.00}_{-1.07}$ km, indicating the lack of significant variation of the radii from $1.4 M_\odot$ to $2.0 M_\odot$. This feature cannot be reproduced by the hadronic matter due to the softening of equation of state (EoS) by hyperon mixing, indicating the possible existence of quark phases in neutron-star interiors. % Two models are used for quark phases: In the quark-hadron transition (QHT) model, quark deconfinement phase transitions from a hadronic-matter EoS are taken into account so as to give reasonable mass-radius ($MR$) curves by adjusting the quark-quark repulsions and the density dependence of effective quark mass. % In the quarkyonic model, the degrees of freedom inside the Fermi sea are treated as quarks and neutrons exist at the surface of the Fermi sea, where $MR$ curves are controlled mainly by the thickness of neutron Fermi layer. % The QHT and quarkyonic EoSs can be adjusted so as to reproduce radii, tidal deformabilities, pressure and central densities inferred from the NICER analysis better than the nucleonic matter EoS, demonstrating the clear impacts of quark phases. Then, the maximum mass for the quakyonic-matter EoS is considerably larger than that for the QHT-matter EoS.

Research on the optical appearance of black holes, both in general relativity and modified gravity, has been in full swing since the Event Horizon Telescope Collaboration announced photos of M87$^{*}$ and Sagittarius A$^{*}$. Nevertheless, limited attention has been given to the impact of tilted accretion disks on black hole images. This paper investigates the $230$ GHz images of non-rotating hairy black holes illuminated by tilted, thin accretion disks in Horndeski gravity with the aid of a ray tracing method. The results indicate that reducing the scalar hair parameter effectively diminishes image luminosity and extends both the critical curve and the inner shadow. This trend facilitates the differentiation between hairy black holes and Schwarzschild black holes. Furthermore, we observe that the inclination of the tilted accretion disk can mimic the observation angle, consequently affecting image brightness and the morphology of the inner shadow. In specific parameter spaces, the disk inclination has the ability to shift the position of the light spot in the images of hairy black holes. This finding may provide potential theoretical evidence for the formation of three flares at different positions in the Sagittarius A$^{*}$ image. Additionally, by examining the images of hairy black holes surrounded by two thin accretion disks, we report the obscuring effect of the accretion environment on the inner shadow of the black hole.

Atefeh Khoshkhahtinat, Ali Zafari, Piyush M. Mehta, Nasser M. Nasrabadi, Barbara J. Thompson, Michael S. F. Kirk, Daniel da Silva

Accepted to IEEE 22$^{nd}$ International Conference on Machine Learning and Applications 2023 (ICMLA) - Selected for Oral Presentation

NASA's Solar Dynamics Observatory (SDO) mission collects large data volumes of the Sun's daily activity. Data compression is crucial for space missions to reduce data storage and video bandwidth requirements by eliminating redundancies in the data. In this paper, we present a novel neural Transformer-based video compression approach specifically designed for the SDO images. Our primary objective is to efficiently exploit the temporal and spatial redundancies inherent in solar images to obtain a high compression ratio. Our proposed architecture benefits from a novel Transformer block called Fused Local-aware Window (FLaWin), which incorporates window-based self-attention modules and an efficient fused local-aware feed-forward (FLaFF) network. This architectural design allows us to simultaneously capture short-range and long-range information while facilitating the extraction of rich and diverse contextual representations. Moreover, this design choice results in reduced computational complexity. Experimental results demonstrate the significant contribution of the FLaWin Transformer block to the compression performance, outperforming conventional hand-engineered video codecs such as H.264 and H.265 in terms of rate-distortion trade-off.