Papers for Monday, Aug 07 2023

Classical Cepheids (CCs) are solid distance indicators and tracers of young stellar populations. Our aim is to provide iron, oxygen, and sulfur abundances for the largest and most homogeneous sample of Galactic CCs ever analyzed. The current sample covers a wide range in Galactocentric distances (RG), pulsation modes and periods. High-resolution and high S/N spectra collected with different spectrographs were adopted to estimate the atmospheric parameters. Individual distances are based on Gaia trigonometric parallaxes or on near-infrared Period-Luminosity relations. We found that Fe and alpha-element radial gradients based on CCs display a well-defined change in the slope for RG larger than 12 kpc. Radial gradients based on open clusters, covering a wide range in age, display similar trends, meaning that the flattening in the outer disk is an intrinsic feature of the radial gradients since it is independent of age. Empirical evidence indicates that the radial gradient for S is steeper than for Fe. The difference in the slope is a factor of two in the linear fit. We also found that S is, on average, under-abundant compared with O. We performed a detailed comparison with Galactic chemical evolution models and we found that a constant Star Formation Efficiency for RG larger than 12 kpc takes account for the flattening in both Fe and alpha-elements. To further constrain the impact that predicted S yields for massive stars have on radial gradients, we adopted a "toy model" and we found that the flattening in the outermost regions requires a decrease of a factor of four in the current S predictions. Sulfur photospheric abundances, compared with other alpha-elements, have the key advantage of being a volatile element. Therefore, stellar S abundances can be directly compared with nebular S abundances in external galaxies.

Meteorites provide an important source of information about the formation and composition of asteroids, because the level of accuracy of studies and analyses performed in a laboratory cannot be achieved by any ground or space based observation. To better understand what asteroid types a meteorite represents, it is crucial to identify the body they originated from. In this paper, we aim to determine possible parent bodies for the known meteorite falls among the known population of near-Earth asteroids (NEAs). By using the similarity criterion $D_N$, based on geocentric quantities, we found 20 possible NEA-meteorite pairs. By performing additional numerical simulations of the backward dynamics, we found that 12 of these pairs may be associated with a possible separation event from the progenitor NEA, while the remaining 8 pairs appear to be ambiguous or random associations. The most interesting are the Pribram and Neuschwanstein meteorites, which are dynamically associated with (482488) 2012 SW20 with a common separation age dating back to about 20$-$30 kyr ago, and the Motopi Pan meteorite, that has three candidate parent bodies: (454100) 2013 BO73, 2017 MC3, and 2009 FZ4. The average time of separation between our meteorite list and the progenitor body appears to be about 10 kyr, a time consistent with what is expected from the collision frequency of small NEAs. Based on our results, we suggest that about 25 per cent of meteorites do not originate in the main belt, but mainly from little collision events happening between NEAs in the inner Solar System.

We present AutoPhot, a fast, powerful, and user-friendly Python based astronomical image photometry solver. AutoPhot incorporates automatic differentiation and GPU (or parallel CPU) acceleration, powered by the machine learning library PyTorch. Everything: AutoPhot can fit models for sky, stars, galaxies, PSFs, and more in a principled Chi^2 forward optimization, recovering Bayesian posterior information and covariance of all parameters. Everywhere: AutoPhot can optimize forward models on CPU or GPU; across images that are large, multi-band, multi-epoch, rotated, dithered, and more. All at once: The models are optimized together, thus handling overlapping objects and including the covariance between parameters (including PSF and galaxy parameters). A number of optimization algorithms are available including Levenberg-Marquardt, Gradient descent, and No-U-Turn MCMC sampling. With an object-oriented user interface, AutoPhot makes it easy to quickly extract detailed information from complex astronomical data for individual images or large survey programs. This paper outlines novel features of the AutoPhot code and compares it to other popular astronomical image modeling software. AutoPhot is open-source, fully Python based, and freely accessible here: https://github.com/Autostronomy/AutoPhot

Using $N$-body simulations, we demonstrate that satellite dwarf galaxy pairs which undergo significant mixing ($\sim 6$ Gyr) can have their respective most bound particles separated great distances upon subsequently merging with a more massive host. This may provide an explanation as to the origin of the complex globular cluster NGC 2149, which is found within the tail of the Sagittarius dwarf spheroidal galaxy, yet separated from its central remnant by over 100 kpc. Dynamical investigations could support the chemical evidence which already points to the NGC 2419 being a nuclear star cluster. Motivated by the distinct nature of NGC 2419, we run a suite of simulations whereby an initial pre-infall merger of two satellites is followed by a post-infall merger of the remnant into a MW-like host potential. We present a striking example from our suite in this work, in which this separation is reproduced by the most bound particles of the two pre-infall satellites. Additionally, this double merger scenario can induce unusual on-sky features in the tidal debris of the post-infall merger, such as clouds, overdensities, and potentially new arms that could resemble the bifurcation observed in Sagittarius.

We propose that the polarization of the emission from the X-ray binary Cygnus X-1, measured using the Imaging X-ray Polarimetry Explorer (IXPE), is imprinted by bulk Comptonization of coronal emission in a mildly relativistic wind or jet with a hollow-cone geometry. Models based on scattering in a static corona overlying a thin accretion disc have difficulty reproducing the relatively high polarization degree (PD ~4%) concurrently with the low inclination (~30 deg) of the binary orbit. We show that bulk outflow with a Lorentz factor > 1.5 is adequate to reproduce the observed PD, with position angle parallel to the large-scale jet, provided that the scattering occurs in a conical sheath offset from the jet axis and our line of sight aligns roughly with the opening angle of the cone. Physically, this flow geometry could represent the entrainment of dense material near the base of an accelerating jet as it passes through the disc corona, or a slow (but still relativistic) sheath around a fast jet. If similar outflows are present in other X-ray binaries at higher inclination, we might expect to see still higher degrees of linear polarization ~10% with an orientation perpendicular to the jet direction.

The Tidal Disruption Event (TDE) AT 2018hyz exhibited a delayed radio flare almost three years after the stellar disruption. Here we report new radio observations of the TDE AT 2018hyz with the AMI-LA and ATCA spanning from a month to more than four years after the optical discovery and 200 days since the last reported radio observation. We detected no radio detection from 30-220 days after the optical discovery in our observations at 15.5 GHz down to a $3\sigma$ level of < 0.14 mJy. The fast-rising, delayed, radio flare is observed in our radio data set and continues to rise almost ~1580 days after the optical discovery. We find that the delayed radio emission, first detected $972$ days after optical discovery, evolves as $t^{4.2 \pm 0.9}$, at 15.5 GHz. Here, we present an off-axis jet model that can explain the full set of radio observations. In the context of this model, we require a powerful narrow jet with an isotropic equivalent kinetic energy $E_{\rm k,iso} \sim 10^{55}$ erg, an opening angle of $ \rm \sim 7^{\circ}$, and a relatively large viewing angle of $ \rm \sim 42^{\circ}$, launched at the time of the stellar disruption. Within our framework, we find that the minimal collimated energy possible for an off-axis jet from AT 2018hyz is $E_k \geq 3 \times 10^{52}$ erg. Finally, we provide predictions based on our model for the light curve turnover time, and for the proper motion of the radio emitting source.

Circumbinary discs (CBDs) arise in many astrophysical settings, including young stellar binaries and supermassive black hole binaries. Their structure is mediated by gravitational torques exerted on the disc by the central binary. The spatial distribution of the binary torque density (so-called excitation torque density) in CBDs is known to feature global large-amplitude, quasi-periodic oscillations, which are often interpreted in terms of the local resonant Lindblad torques. Here we investigate the nature of these torque oscillations using 2D, inviscid hydrodynamic simulations and theoretical calculations. We show that torque oscillations arise due to the gravitational coupling of the binary potential to the density waves launched near the inner cavity and freely propagating out in the disc. We provide analytical predictions for the radial periodicity of the torque density oscillations and verify them with simulations, showing that disc sound speed and the multiplicity of the density wave spiral arms are the key factors setting the radial structure of the oscillations. Resonant Lindblad torques play no direct role in determining the radial structure and periodicity of the torque oscillations and manifest themselves only by driving the density waves in the disc. We also find that vortices forming at the inner edge of the disc can provide a non-trivial contribution to the angular momentum transport in the CBD. Our results can be applied to understanding torque behaviour in other settings, e.g. discs in cataclysmic variables and X-ray binaries.

Embedded, Class 0/I protostellar disks represent the initial condition for planet formation. This calls for better understandings of their bulk properties and the dust grains within them. We model multi-wavelength dust continuum observations of the disk surrounding the Class I protostar TMC1A to provide insight on these properties. The observations can be well fit by a gravitationally self-regulated (i.e., marginally gravitationally unstable and internally heated) disk model, with surface density $\Sigma \sim 1720 (R/10au)^{-1.96} g/cm^2$ and midplane temperature $T_{mid} \sim 185 (R/10au)^{-1.27} K$. The observed disk contains a $m=1$ spiral substructure; we use our model to predict the spiral's pitch angle and the prediction is consistent with the observations. This agreement serves as both a test of our model and strong evidence of the gravitational nature of the spiral. Our model estimates a maximum grain size $a_{max}\sim 196(R/10au)^{-2.45} \mu m$, which is consistent with grain growth being capped by a fragmentation barrier with threshold velocity $\sim 1 m/s$. We further demonstrate that observational properties of TMC1A are typical among the observed population of Class 0/I disks, which hints that traditional methods of disk data analyses based on Gaussian fitting and the assumption of the optically thin dust emission could have systematically underestimated disk size and mass and overestimated grain size.

Above $\sim$3 keV, the X-ray spectrum of the active galactic nuclei (AGN) is characterized by the intrinsic continuum and compton reflection features. For type-1 AGN, several regions could contribute to the reflection. In order to investigate the nature of the reflecting medium, we perform a systematic analysis of the reflector using XMM-Newton and NuSTAR observations of a sample of 22 type-1 AGN. We create a baseline model which includes Galactic absorption and an intrinsically absorbed power-law plus a reflection model. We test a set of nine reflection models in a sub-sample of five objects. Based on these results, we select three models to be tested on the entire sample, accounting for distinct physical scenarios: neutral/distant reflection, ionized/relativistic reflection, and neutral/distant+ionized/relativistic reflection, namely hybrid model. We find that 18 sources require the reflection component to fit their spectra. Among them, 67$\%$ prefer the hybrid model. Neutral and ionized models are equally preferred by three sources. We conclude that both the neutral/distant reflector most probably associated with the inner edges of the torus and the ionized/relativistic reflector associated with the accretion disk are required to describe the reflection in type-1 AGN.

We announce the detection of escaping helium from TOI 2134b, a mini-Neptune a few Gyr old. The average in-transit absorption spectrum shows a peak of 0.37 +- 0.05% and an equivalent width of $W_{\rm avg}=3.3 \pm 0.3$ m$\r{A}$. Among all planets with helium detections, TOI 2134b is the only mature mini-Neptune, has the smallest helium signal, and experiences the lowest XUV flux. Putting TOI 2134b in the context of all other helium detections, we report the detection of a strong (p=3.0e-5) and theoretically expected correlation between $F_{\rm XUV}/\rho_{\rm XUV}$ (proportional to the energy-limited mass loss rate) and $R_* W_{\rm avg}$ (roughly proportional to the observationally inferred mass loss rate). Here, $W_{\rm avg}$ is the equivalent width of the helium absorption and $\rho_{\rm XUV}$ is the density of the planet within the XUV photosphere, but the correlation is similarly strong if we use the optical photosphere. TOI 2134b anchors the relation, having the lowest value on both axes. We encourage further observations to fill in missing regions of this parameter space and improve estimates of $F_{\rm XUV}$.

Large High Altitude Air Shower Observatory~(LHAASO) opens the window of ultra-high-energy~(UHE) photon detection, broadens the path of testing basic physical concept such as Lorentz symmetry, and brings possibility of potential high-energy physical phenomenon research such as photon decay and electron decay. Currently, the UHE photons from LHAASO observation set strict constraints on photon and electron Lorentz symmetry violation~(LV) effects. To obtain a global impression of the photon-electron LV parameter plane, we make a detailed analysis for photon decay and electron decay. Our discussion gives the corresponding decay thresholds and energy-momentum distributions in different LV parameter configurations. We get corresponding constraints on photon LV parameter, electron LV parameter and the photon-electron LV parameter plane from LHAASO observation. For the space allowed for LV effect, that is beyond relativity, we also provide corresponding boundaries from LHAASO observation.

The Nancy Grace Roman Space Telescope will survey a large area of the sky at near-infrared wavelengths with its Wide Field Instrument (WFI). The performance of the 18 WFI H4RG-10 detectors will need to be well-characterized and regularly monitored in order for Roman to meet its science objectives. Weak lensing science goals are particularly sensitive to instrumental distortions and patterns that might masquerade as astronomical signals. We apply the wavelet scattering transform in order to analyze localized signals in Roman WFI images that have been taken as part of a dark image test suite. The scattering transform quantifies shapes and clustering information by reducing images into non-linear combinations of wavelet modes on multiple size scales. We show that these interpretable scattering statistics can separate rare correlated patterns from typical noise signals, and we discuss the results in context of power spectrum analyses and other computer vision methods.

Context: Potentially hazardous asteroids (PHA) in Earth-crossing orbits pose a constant threat to life on Earth. Several mitigation methods have been proposed, and the most feasible technique appears to be the disintegration of the impactor and the generation of a fragment cloud by explosive penetrators at interception. Mitigation analyses, however neglect the effect of orbital dynamics on fragments trajectory. Aims: We aim at studying the effect of orbital dynamics of the impactor's cloud on the number of fragments that hit the Earth assuming different interception dates. The effect of self-gravitational cohesion and the axial rotation of the impactor are also investigated. Methods: The orbits of 10^5 fragments are computed with a high-precision direct N-body integrator of the 8th order, running on GPUs. We consider orbital perturbations from all large bodies in the Solar System and the self-gravity of the cloud fragments. Results: Using a series of numerical experiments, we show that orbital shear causes the fragment cloud to adopt the shape of a triaxial ellipsoid. The shape and alignment of the triaxial ellipsoid are strongly modulated by the cloud's orbital trajectory, and hence the impact cross-section of the cloud with respect to the Earth. Therefore, the number of fragments hitting the Earth is strongly influenced by the orbit of the impactor and the time of interception. A minimum number of impacts occurs for a well-defined orientation of the impactor rotational axis, depending on the date of interception. Conclusions: To minimise the lethal consequences of an PHA's impact, a well-constrained interception timing is necessary. Too early interception may not be ideal for PHAs in the Apollo or Aten groups. The best time to intercept PHA is when it is at the pericentre of its orbit.

We present the results of our analysis of multiwavelength observations for the long gamma-ray burst GRB 200829A. The burst redshift $z \approx 1.29 \pm 0.04$ has been determined photometrically at the afterglow phase. In gamma rays the event is one of the brightest (in isotropic equivalent), $E_{iso} \gtrsim 10^{54}$ erg. The multicolor light curve of the GRB 200829A afterglow is characterized by chromatic behavior and the presence of a plateau gradually transitioning into a power-law decay that can also be interpreted as a quasi-synchronous inhomogeneity (flare). We assume that the presence of a chromatic inhomogeneity in the early afterglow is consistent with the model of a structured jet.

We report direct observational evidence for a latitudinal dependence of dust cloud opacity in ultracool dwarfs, indicating that equatorial latitudes are cloudier than polar latitudes. These results are based on a strong positive correlation between the viewing geometry and the mid-infrared silicate absorption strength in mid-L dwarfs using mid-infrared spectra from the Spitzer Space Telescope and spin axis inclination measurements from available information in the literature. We confirmed that the infrared color anomalies of L dwarfs positively correlate with dust cloud opacity and viewing geometry, where redder objects are inclined equator-on and exhibit more opaque dust clouds while dwarfs viewed at higher latitudes and with more transparent clouds are bluer. These results show the relevance of viewing geometry to explain the appearance of brown dwarfs and provide insight into the spectral diversity observed in substellar and planetary atmospheres. We also find a hint that dust clouds at similar latitudes may have higher opacity in low-surface gravity dwarfs than in higher-gravity objects.

We review the iid2022 workshop on statistical methods for X-ray and $\gamma$-ray astronomy and high--energy astrophysics event data in astronomy, held in Guntersville, AL, on Nov. 15-18 2022. New methods for faint source detection, spatial point processes, variability and spectral analysis, and machine learning are discussed. Ideas for future developments of advanced methodology are shared.

As Physics-Informed Neural Networks and other methods for full-vector-field construction or analysis become more prominent, a need has developed for a large set of simulated active regions for training, validation and testing purposes. We use a state-of-the-art magnetohydrostatic extrapolation method to develop a public dataset of over five thousand data cubes based on the Spaceweather HMI Active Region Patch (SHARP) library of active region magnetogram images. Each cube resolves the magnetic field vector and plasma forcing at approximately 100,000 scattered points that are adaptively clustered near the high-flux regions of the domain. This paper describes the methodology of construction of the Plasma-prescribed Active Region Static Extrapolation (PARSE) dataset, as well as its structure and how to access it.

Interstellar dust extinction law is essential for interpreting observations. In this work, we investigate the ultraviolet (UV)--mid-infrared (IR) extinction law of the Taurus molecular cloud and its possible variations. We select 504,988 dwarf stars (4200 K < Teff < 8000 K) and 4,757 giant stars (4200 K < Teff < 5200 K) based on the stellar parameters of Gaia DR3 as tracers. We establish the Teff--intrinsic color relations and determine the intrinsic color indices and color excesses for different types of stars. In the determination of color excess ratios (CERs), we analyze and correct the curvature of CERs and derive the UV--mid-IR CERs of 16 bands. We consider different effective wavelengths for different types of stars when converting CERs to relative extinction, and obtain the extinction law with a better wavelength resolution. In addition, we analyze the possible regional variation of extinction law and derive the average extinction law of Rv=3.13+-0.32 for the Taurus molecular cloud. Only 0.9% of subregions have deviations >3sigma, indicating limited regional variation in the extinction law. We also discuss the effect of Gaia Teff overestimation on the determination of the Taurus extinction law and find that the effect is negligible.

RR Lyrae (RR Lyr) stars are a well-known and useful distance indicator for old stellar populations such as globular clusters and dwarf galaxies. Fundamental-mode RR Lyr (RRab) stars are commonly used to measure distances, and the accuracy of the determined distance is strongly constrained by metallicity. Here, we investigate the metallicity dependence in the period-luminosity (PL) relation of double-mode RR Lyr (RRd) stars. We find and establish a linear relation between metallicity and period or period ratio for RRd stars. This relation can predict the metallicity as accurately as the low-resolution spectra. Based on this relation, we establish a metallicity-independent PL relation for RRd stars. Combining the distance of the Large Magellanic Cloud and Gaia parallaxes, we calibrate the zero point of the derived PL relation to an error of 0.022 mag. Using RRd stars, we measure the distances of globular clusters and dwarf galaxies with an accuracy of 2-3% and 1-2%, respectively. In the future, RRd stars could anchor galaxy distances to an accuracy of 1.0% and become an independent distance ladder in the Local Group.

Many interstellar complex organic molecules (COMs) are believed to be produced on the surfaces of icy grains at low temperatures. Atomic carbon is considered responsible for the skeletal evolution processes, such as C-C bond formation, via insertion or addition reactions. Before reactions, C atoms must diffuse on the surface to encounter reaction partners; therefore, information on their diffusion process is critically important for evaluating the role of C atoms in the formation of COMs. In situ detection of C atoms on ice was achieved by a combination of photostimulated desorption and resonance enhanced multiphoton ionization methods. We found that C atoms weakly bound to the ice surface diffused approximately above 30 K and produced C2 molecules. The activation energy for C-atom surface diffusion was experimentally determined to be 88 meV (1,020 K), indicating that the diffusive reaction of C atoms is activated at approximately 22 K on interstellar ice. The facile diffusion of C at T > 22 K atoms on interstellar ice opens a previously overlooked chemical regime where the increase in complexity of COMs as driven by C atoms. Carbon addition chemistry can be an alternative source of chemical complexity in translucent clouds and protoplanetary disks with crucial implications in our current understanding on the origin and evolution of organic chemistry in our Universe.

Magnetohydrodynamic kink waves naturally form as a consequence of perturbations to a structured medium, for example transverse oscillations of coronal loops. Linear theory has provided many insights in the evolution of linear oscillations, and results from these models are often applied to infer information about the solar corona from observed wave periods and damping times. However, simulations show that nonlinear kink waves can host the Kelvin-Helmholtz instability (KHi) which subsequently creates turbulence in the loop, dynamics which are beyond linear models. In this paper we investigate the evolution of KHi-induced turbulence on the surface of a flux tube where a non-linear fundamental kink-mode has been excited. We control our numerical experiment so that we induce the KHi without exciting resonant absorption. We find two stages in the KHi turbulence dynamics. In the first stage, we show that the classic model of a KHi turbulent layer growing $\propto t$is applicable. We adapt this model to make accurate predictions for damping of the oscillation and turbulent heating as a consequence of the KHi dynamics. In the second stage, the now dominant turbulent motions are undergoing decay. We find that the classic model of energy decay proportional to $t^{-2}$ approximately holds and provides an accurate prediction of the heating in this phase. Our results show that we can develop simple models for the turbulent evolution of a non-linear kink wave, but the damping profiles produced are distinct from those of linear theory that are commonly used to confront theory and observations.

An open supercluster (OSC) is defined as a cluster of at least six open clusters (OCs) born from the same giant molecular cloud (GMC). We survey the recent catalogs of OCs based on Gaia data and relevant literature to find 17 OSCs of the third galactic quadrant, along with 190 likely members of them. OSCs are frequent enough to be considered an extra class of objects in the hierarchy of star formation. Some of these supersystems are new and most of them contain more members than previously thought. The detailed study of some OSCs has leaded to the discovery of four new young clusters that are members of them, named Casado-Hendy 2 to 5. In certain instances, subgroups with distinct PMs or 3D positions have been found within an OSC, suggesting the presence of multiple generations of stars formed from several bursts of star formation within the same GMC. OSCs are typically unbound and tend to disintegrate on timescales of 0.1 Gyr. The present results support the Primordial Group hypothesis (Casado 2022) and suggest that globular clusters are not formed from the accretion of open superclusters, at least in the local Universe at late times.

In the absence of numerous gravitational-wave detections with confirmed electromagnetic counterparts, the "dark siren" method has emerged as a leading technique of gravitational-wave cosmology. The method allows redshift information of such events to be inferred statistically from a catalogue of potential host galaxies. Due to selection effects, dark siren analyses necessarily depend on the mass distribution of compact objects and the evolution of their merger rate with redshift. Informative priors on these quantities will impact the inferred posterior constraints on the Hubble constant ($H_0$). It is thus crucial to vary these unknown distributions during an $H_0$ inference. This was not possible in earlier analyses due to the high computational cost, restricting them to either excluding galaxy catalogue information, or fixing the gravitational-wave population mass distribution and risking introducing bias to the $H_0$ measurement. This paper introduces a significantly enhanced version of the Python package GWCOSMO, which allows joint estimation of cosmological and compact binary population parameters. This thereby ensures the analysis is now robust to a major source of potential bias. The gravitational-wave events from the Third Gravitational-Wave Transient Catalogue are reanalysed with the GLADE+ galaxy catalogue, and an updated, more reliable measurement of $H_0=69^{+12}_{-7}$ km s$^{-1}$ Mpc$^{-1}$ is found (maximum a posteriori probability and 68% highest density interval). This improved method will enable cosmological analyses with future gravitational-wave detections to make full use of the information available (both from galaxy catalogues and the compact binary population itself), leading to promising new independent bounds on the Hubble constant.

Chariklo has two narrow and dense rings, C1R and C2R, located at 391 km and 405 km, respectively. In the light of new stellar occultation data, we study the stability around Chariklo. We also analyse three confinement mechanisms, to prevent the spreading of the rings, based on shepherd satellites in resonance with the edges of the rings. This study is made through a set of numerical simulations and the Poincar\'e surface of section technique. From the numerical simulation results we verify that, from the current parameters referring to the shape of Chariklo, the inner edge of the stable region is much closer to Chariklo than the rings. The Poincar\'e surface of sections allow us to identify the first kind periodic and quasi-periodic orbits, and also the resonant islands corresponding to the 1:2, 2:5, and 1:3 resonances. We construct a map of a versus e space which gives the location and width of the stable region and the 1:2, 2:5, and 1:3 resonances. We found that the first kind periodic orbits family can be responsible for a stable region whose location and size meet that of C1R, for specific values of the ring particles' eccentricities. However, C2R is located in an unstable region if the width of the ring is assumed to be about 120 m. After analysing different systems we propose that the best confinement mechanism is composed of three satellites, two of them shepherding the inner edge of C1R and the outer edge of C2R, while the third satellite would be trapped in the 1:3 resonance.

The IceCube Observatory provides our highest-statistics picture of the cosmic-ray arrival directions in the Southern Hemisphere, with over 700 billion cosmic-ray-induced muon events collected between May 2011 and May 2022. Using the larger data volume, we find an improved significance of the PeV cosmic ray anisotropy down to scales of $6^\circ$. In addition, we observe a variation in the angular power spectrum as a function of energy, hinting at a relative decrease in large-scale features above 100 TeV. The data-taking period covers a complete solar cycle, providing new insight into the time variability of the signal. We present preliminary results using this up-to-date event sample.

We present the results of N-body simulations meant to reproduce the long-term effects of mutually inclined exoplanets on debris disks, using the HD 202628 system as a proxy. HD 202628 is a Gyr-old solar-type star that possesses a directly observable, narrow debris ring with a clearly defined inner edge and non-zero eccentricity, hinting at the existence of a sculpting exoplanet. The eccentric nature of the disk leads us to examine the effect on it over Gyr timescales from an eccentric and inclined planet, placed on its orbit through scattering processes. We find that, in systems with dynamical timescales akin to that of HD 202628, a planetary companion is capable of completely tilting the debris disk. This tilt is preserved over the Gyr age of the system. Simulated observations of our models show that an exoplanet around HD 202628 with an inclination misalignment $\gtrsim\,10$ degrees would cause the disk to be observably diffuse and broad, which is inconsistent with ALMA observations. With these observations, we conclude that if there is an exoplanet shaping this disk, it likely had a mutual inclination of less than 5 degrees with the primordial disk. Conclusions of this work can be either applied to debris disks appearing as narrow rings (e.g., Fomalhaut, HR 4796), or to disks that are vertically thick at ALMA wavelengths (e.g., HD 110058).

We investigate the disk formation process in the TNG50 simulation, examining the profiles of SFR surface density ($\Sigma_{\rm SFR}$), gas inflow and outflow, and the evolution of the angular momentum of inflowing gas particles. The TNG50 galaxies tend to have larger star-forming disks, and also show larger deviations from exponential profiles in $\Sigma_{\rm SFR}$ when compared to real galaxies in the MaNGA (Mapping Nearby Galaxies at APO) survey. The stellar surface density of TNG50 galaxies show good exponential profiles, which is found to be the result of strong radial migration of stars over time. However, this strong radial migration of stars in the simulation produces flatter age profiles in TNG50 disks compared to observed galaxies. The star formation in the simulated galaxies is sustained by a net gas inflow and this gas inflow is the primary driver for the cosmic evolution of star formation, as expected from simple gas-regulator models of galaxies. There is no evidence for any significant loss of angular momentum for the gas particles after they are accreted on to the galaxy, which may account for the large disk sizes in the TNG50 simulation. Adding viscous processes to the disks, such as the magnetic stresses from magneto-rotational instability proposed by Wang & Lilly 2022, will likely reduce the sizes of the simulated disks and the tension with the sizes of real galaxies, and may produce more realistic exponential profiles.

The origin of rapid rotation in massive stars remains debated, although binary interactions are now often advocated as a cause. However, the broad and shallow lines in the spectra of fast rotators make direct detection of binarity difficult. In this paper, we report on the discovery and analysis of multiplicity for three fast-rotating massive stars: HD25631 (B3V), HD191495 (B0V), and HD46485 (O7V). They display strikingly similar TESS light curves, with two narrow eclipses superimposed on a sinusoidal variation due to reflection effects. We complement these photometric data by spectroscopy from various instruments (X-Shooter, Espadons, FUSE...), to further constrain the nature of these systems. The detailed analyses of these data demonstrates that the companions of the massive OB stars have low masses (~1Msol) with rather large radii (2-4 Rsol) and low temperatures (<15 kK). These companions display no UV signature, which would exclude a hot subdwarf nature, but disentangling of the large set of X-Shooter spectra of HD25631 revealed the typical signature of chromospheric activity in the companion's spectrum. In addition, despite the short orbital periods (P=3-7d), the fast-rotating OB-stars still display non-synchronized rotation and all systems appear young (<20Myr). This suggests that, as in a few other cases, these massive stars are paired in those systems with non-degenerate, low-mass PMS companions, implying that fast rotation would not be a consequence of a past binary interactions in their case.

We present a study of the tilt-to-length coupling noise during the LISA Pathfinder mission and how it depended on the system's alignment. Tilt-to-length coupling noise is the unwanted coupling of angular and lateral spacecraft or test mass motion into the primary interferometric displacement readout. It was one of the major noise sources in the LISA Pathfinder mission and is likewise expected to be a primary noise source in LISA. We demonstrate here that a recently derived and published analytical model describes the dependency of the LISA Pathfinder tilt-to-length coupling noise on the alignment of the two freely falling test masses. This was verified with the data taken before and after the realignments performed in March (engineering days) and June 2016, and during a two-day experiment in February 2017 (long cross-talk experiment). The latter was performed with the explicit goal of testing the tilt-to-length coupling noise dependency on the test mass alignment. Using the analytical model, we show that all realignments performed during the mission were only partially successful and explain the reasons why. In addition to the analytical model, we computed another physical tilt-to-length coupling model via a minimising routine making use of the long cross-talk experiment data. A similar approach could prove useful for the LISA mission.

Here we show how $H_0$ can be derived purely from the gravitational waves (GW) of neutron star-black hole (NSBH) mergers. This new method provides an estimate of $H_0$ spanning the redshift range, $z<0.25$ with current GW sensitivity and without the need for any afterglow detection. We utilise the inherently tight neutron star mass function together with the NSBH waveform amplitude and frequency to estimate distance and redshift respectively, thereby obtaining $H_0$ statistically. Our first estimate is $H_0 = 86^{+55}_{-46}$ km s$^{-1}$ Mpc$^{-1}$ for the secure NSBH events GW190426 and GW200115. We forecast that soon, with 10 more such NSBH events we can reach competitive precision of $\delta H_0/H_0 \lesssim 20\%$.

Cepheids have been the cornerstone of the extragalactic distance scale for a century. With high-quality data, these luminous supergiants exhibit a small dispersion in their Leavitt (period-luminosity) relation, particularly at longer wavelengths, and few methods rival the precision possible with Cepheid distances. In these proceedings, we present an overview of major observational programs pertaining to the Cepheid extragalactic distance scale, its progress and remaining challenges. In addition, we present preliminary new results on Cepheids from the James Webb Space Telescope (JWST). The launch of JWST has opened a new chapter in the measurement of extragalactic distances and the Hubble constant. JWST offers a resolution three times that of the Hubble Space Telescope (HST) with nearly 10 times the sensitivity. It has been suggested that the discrepancy in the value of the Hubble constant based on Cepheids compared to that inferred from measurements of the cosmic microwave background requires new and additional physics beyond the standard cosmological model. JWST observations will be critical in reducing remaining systematics in the Cepheid measurements and for confirming if new physics is indeed required. Early JWST data for the galaxy, NGC 7250 show a decrease in scatter in the Cepheid Leavitt law by a factor of two relative to existing HST data and demonstrate that crowding/blending effects are a significant issue in a galaxy as close as 20 Mpc.

High-contrast imaging of debris disk systems permits us to assess the composition and size distribution of circumstellar dust, to probe recent dynamical histories, and to directly detect and characterize embedded exoplanets. Observations of these systems in the infrared beyond 2--3 $\mu$m promise access to both extremely favorable planet contrasts and numerous scattered-light spectral features -- but have typically been inhibited by the brightness of the sky at these wavelengths. We present coronagraphy of the AU Microscopii (AU Mic) system using JWST's Near Infrared Camera (NIRCam) in two filters spanning 3--5 $\mu$m. These data provide the first images of the system's famous debris disk at these wavelengths and permit additional constraints on its properties and morphology. Conducting a deep search for companions in these data, we do not identify any compelling candidates. However, with sensitivity sufficient to recover planets as small as $\sim 0.1$ Jupiter masses beyond $\sim 2^{\prime\prime}$ ($\sim 20$ au) with $5\sigma$ confidence, these data place significant constraints on any massive companions that might still remain at large separations and provide additional context for the compact, multi-planet system orbiting very close-in. The observations presented here highlight NIRCam's unique capabilities for probing similar disks in this largely unexplored wavelength range, and provide the deepest direct imaging constraints on wide-orbit giant planets in this very well studied benchmark system.

I consider the current sample of galaxy nuclei producing quasiperiodic eruptions (QPEs). If the quasiperiod results from the orbital motion of a star around the central black hole, the dearth of associated black hole masses $\gtrsim 10^6\msun$ places tight constraints on models. It disfavours those assuming wide orbits and small eccentricities, because there is ample volume within pericentre to allow significantly more massive holes in QPE systems than are currently observed. If instead the orbiting star is assumed to pass close to the black hole, the same lack of large black hole masses strongly suggests that the stellar orbits must be significantly eccentric, with $1 - e \lesssim {\rm few}\times 10^{-2}$. This favours a tidal disruption near-miss picture where QPEs result from repeated accretion from an orbiting star (in practice a white dwarf) losing orbital angular momentum to gravitational radiation, even though this is not assumed in deriving the eccentricity constraint. Given the tight constraints resulting from the current small observed sample, attempts to find QPE systems in more massive galaxies are clearly important.

Dark matter (DM) could be a pseudo-Dirac thermal relic with a small mass splitting that is coupled off-diagonally to a kinetically mixed dark photon. This model, particularly in the sub-GeV mass range, is a key benchmark for accelerator searches and direct detection experiments. Typically, the presence of even a tiny fraction of pseudo-Dirac DM in the excited state around the time of recombination would be excluded by DM annihilation bounds from the cosmic microwave background (CMB); thus, viable thermal histories must typically feature an exponential suppression of the excited state. We revisit assumptions about the thermal history in the resonant regime, where the dark photon mass is slightly more than twice the DM mass (to within $\sim10\%$), leading to an $s$-channel resonance in the annihilation cross section. This resonance substantially reduces the couplings required for achieving the observed relic abundance, implying that in much of the parameter space, the DM kinetically decouples from the Standard Model well before the final DM relic abundance is achieved. We find that the excited state is not thermally depopulated in this regime. In spite of this, we find that the presence of the excited state does $\textit{not}$ violate CMB bounds, even for arbitrarily small mass splittings. The present-day abundance of the excited state opens up the possibility of signatures that are usually not relevant for pseudo-Dirac DM, including indirect detection, direct detection, and self-interacting DM signatures.

The phenomenon of the Dark matter baffles the researchers: the underlying dark particle has escaped so far the detection and its astrophysical role appears complex and entangled with that of the standard luminous particles. We propose that, in order to act efficiently, alongside with abandoning the current $\Lambda CDM$ scenario, we need also to shift the Paradigm from which it emerged.

Large-scale compressive slow-mode-like fluctuations can cause variations in the density, temperature, and magnetic-field magnitude in the solar wind. In addition, they also lead to fluctuations in the differential flow $U_{\rm p\alpha}$ between $\alpha$-particles and protons ($p$), which is a common source of free energy for the driving of ion-scale instabilities. If the amplitude of the compressive fluctuations is sufficiently large, the fluctuating $U_{\rm p\alpha}$ intermittently drives the plasma across the instability threshold, leading to the excitation of ion-scale instabilities and thus the growth of corresponding ion-scale waves. The unstable waves scatter particles and reduce the average value of $U_{\rm p\alpha}$. We propose that this "fluctuating-beam effect" maintains the average value of $U_{\rm p\alpha}$ well below the marginal instability threshold. We model the large-scale compressive fluctuations in the solar wind as long-wavelength slow-mode waves using a multi-fluid model. We numerically quantify the fluctuating-beam effect for the Alfv\'en/ion-cyclotron (A/IC) and fast-magnetosonic/whistler (FM/W) instabilities. We show that measurements of the proton-$\alpha$ differential flow and compressive fluctuations from the {\it Wind} spacecraft are consistent with our predictions for the fluctuating-beam effect. This effect creates a new channel for a direct cross-scale energy transfer from large-scale compressions to ion-scale fluctuations.

We examine the Sommerfeld enhancement effect for the puffy self-interacting dark matter. We find out two new parameters to classify the self-scattering cross section into the Born, the resonance and the classical regimes for the puffy dark matter. Then we observe that the resonance peaks for the puffy dark matter self-scattering and for the Sommerfeld enhancement effect have the same locations. Further, we find that for a large ratio between $R_{\chi}$ (radius of a puffy dark matter particle) and $1/m_{\phi}$ (force range), the Sommerfeld enhancement factor approaches to 1 (no enhancement). Finally, for the puffy SIDM scenario to solve the small-scale problems, the values of the Sommerfeld enhancement factor are displayed in the allowed parameter regions.

In recent years, $f(Q)$ gravity has enjoyed considerable attention in the literature and important results have been obtained. However, the question of how many physical degrees of freedom the theory propagates -- and how this number may depend on the form of the function $f$ -- has not been answered satisfactorily. In this article we show that a Hamiltonian analysis based on the Dirac-Bergmann algorithm -- one of the standard methods to address this type of question -- fails. We isolate the source of the failure, show that other commonly considered teleparallel theories of gravity are affected by the same problem, and we point out that the number of degrees of freedom obtained in Phys. Rev. D 106 no. 4, (2022) by K. Hu, T. Katsuragawa, and T. Qui (namely eight), based on the Dirac-Bergmann algorithm, is wrong. Using a different approach, we show that the upper bound on the degrees of freedom is seven. Finally, we propose a more promising strategy for settling this important question.

This work outlines a time-domain numerical integration technique for linear hyperbolic partial differential equations sourced by distributions (Dirac $\delta$-functions and their derivatives). Such problems arise when studying binary black hole systems in the extreme mass ratio limit. We demonstrate that such source terms may be converted to effective domain-wide sources when discretized, and we introduce a class of time-steppers that directly account for these discontinuities in time integration. Moreover, our time-steppers are constructed to respect time reversal symmetry, a property that has been connected to conservation of physical quantities like energy and momentum in numerical simulations. To illustrate the utility of our method, we numerically study a distributionally-sourced wave equation that shares many features with the equations governing linear perturbations to black holes sourced by a point mass.