Papers for Monday, Jun 26 2023

I report on Chandra X-ray observations of SN 2018cqj, a low luminosity Type Ia supernova that showed an H$\alpha$ line in its optical spectrum. No X-ray emission was detected at the location of the SN, with an upper limit to the X-ray luminosity of 2 $\times 10^{39}$ erg s$^{-1}$.

We explore reprocessing models for a sample of 17 hypervariable quasars, taken from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, which all show coordinated optical luminosity hypervariability with amplitudes of factors $\gtrsim 2$ between 2014 and 2020. We develop and apply reprocessing models for quasar light curves in simple geometries that are likely to be representative of quasar inner environments. In addition to the commonly investigated thin-disk model, we include the thick-disk and hemisphere geometries. The thick-disk geometry could, for instance, represent a magnetically-elevated disk, whereas the hemisphere model can be interpreted as a first-order approximation for any optically-thick out-of-plane material caused by outflows/winds, warped/tilted disks, etc. Of the 17 quasars in our sample, eleven are best-fit by a hemisphere geometry, five are classified as thick disks, and both models fail for just one object. We highlight the successes and shortcomings of our thermal reprocessing models in case studies of four quasars that are representative of the sample. While reprocessing is unlikely to explain all of the variability we observe in quasars, we present our classification scheme as a starting point for revealing the likely geometries of reprocessing for quasars in our sample and hypervariable quasars in general.

The recent discovery of two detached black hole-star (BH-star) binaries from Gaia's third data release has sparkled interest in understanding the formation mechanisms of these systems. We investigate the formation of these systems by dynamical processes in young open star clusters (SCs) and via isolated binary (IB) evolution, using a combination of direct $N$-body models and population synthesis simulations. By comparing dynamical and isolated systems created using the same model of binary stellar evolution, we find that dynamical formation in SCs is nearly 40 times as efficient per unit of star formation at producing BH-star binaries compared to IB evolution. We expand this analysis to the full Milky Way (MW) using a FIRE-2 hydrodynamical simulation of a MW-mass galaxy. Even assuming that only $10\%$ of star formation produces SCs with masses $> 1000\,\mathrm{M_{\odot}}$, we find that the MW contains $\sim 2 \times 10^5$ BH-star systems, with approximately 4 out of every 5 systems being formed dynamically. Many of these dynamically-formed systems have larger orbital periods, eccentricities, and black hole masses than their isolated counterparts. For binaries older than 100 Myr, we show that any detectable system with $e\gtrsim0.5$ or $M_{\rm BH}\gtrsim 10\,\mathrm{M_{\odot}}$ can only be formed through dynamical processes. Our MW model predicts between 61 and 210 such detections from the complete DR4 Gaia catalog, with the majority of systems being dynamically formed in massive and metal-rich SCs. Finally, we compare our populations to the recently discovered Gaia BH1 and Gaia BH2, and conclude that the dynamical scenario is the most favorable formation pathway for both systems.

We study signatures of macroscopic dark matter (DM) in current and future gravitational wave (GW) experiments. Transiting DM with a mass of $\sim10^5-10^{15}$ kg that saturates the local DM density can be potentially detectable by GW detectors, depending on the baseline of the detector and the strength of the force mediating the interaction. In the context of laser interferometers, we derive the gauge invariant observable due to a transiting DM, including the Shapiro effect, and adequately account for the finite photon travel time within an interferometer arm. In particular, we find that the Shapiro effect can be dominant for short-baseline interferometers such as Holometer and GQuEST. We also find that proposed experiments such as Cosmic Explorer and Einstein Telescope can constrain a fifth force between DM and baryons, at the level of strength $\sim 10^3$ times stronger than gravity for, e.g., kg mass DM with a fifth-force range of $10^6$ m.

Using a representative sample of Milky Way (MW)-like galaxies from the TNG50 cosmological-volume simulation, we investigate physical processes driving the formation of galactic disks. A disk forms as a result of the interplay between inflow and outflow carrying angular momentum in and out of the galaxy. Interestingly, the inflow and outflow have remarkably similar distributions of angular momentum, suggesting an exchange of angular momentum and/or outflow recycling, leading to continuous feeding of pre-aligned material from the co-rotating circumgalactic medium. We show that disk formation in TNG50 is correlated with stellar bulge formation, in qualitative agreement with a recent theoretical model of disk formation facilitated by steep gravitational potentials. Disk formation is also correlated with the formation of a hot circumgalactic halo with a significant fraction of the inflow occurring at sub- and transonic velocities. In the context of recent theoretical works connecting disk settling and hot halo formation, our results imply that the subsonic part of the inflow may settle into a disk while the remaining supersonic inflow will perturb this disk via the chaotic cold accretion. We find that disks tend to form when the host halos become more massive than $\sim (1-2) \times 10^{11} M_\odot$, consistent with previous theoretical findings and observational estimates of the pre-disk protogalaxy remnant in the MW. Our results do not prove that either co-rotating outflow recycling, gravitational potential steepening, or hot halo formation cause disk formation but they show that all these processes occur concurrently and may play an important role in disk growth.

Curvatons are light (compared to the Hubble scale during inflation) spectator fields during inflation that potentially contribute to adiabatic curvature perturbations post-inflation. They can alter CMB observables such as the spectral index $n_s$, the tensor-to-scalar ratio $r$, and the local non-Gaussianity $\;f_{\rm NL}^{\rm (loc)}$. We systematically explore the observable space of a curvaton with a quadratic potential. We find that when the underlying inflation model does not satisfy the $n_s$ and $r$ observational constraints but can be made viable with a significant contribution from what we call a savior curvaton, a large $\;f_{\rm NL}^{\rm (loc)}$ is inevitable. On the other hand, when the underlying inflation model already satisfies the $n_s$ and $r$ observational constraints, so significant curvaton contribution is forbidden, a large $\;f_{\rm NL}^{\rm (loc)}$ is possible in the exceptional case when the isocurvature fluctuation in the curvaton fluid is much greater than the global curvature fluctuation.

The Gaia mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple stellar systems are common in our galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orbit about the inner binary. In such a system, the two orbits torque each other via the eccentric Kozai-Lidov mechanism (EKL), which can alter the orbital configuration of the inner binary. We simulate thousands of triple stellar systems for over 10 Gyr, tracking gravitational interactions, tides, general relativity, and stellar evolution up to their white dwarf fate. As demonstrated here, three-body dynamics coupled with stellar evolution is a critical channel to form tight WD binaries or lead to a merger of a WD and a companion. Amongst these triples, we explore their manifestations as cataclysmic variables, Type Ia supernovae, and gravitational-wave events. The simulated systems are then compared to a sample of WD triples selected from the Gaia catalog. We find that including the effect of mass loss-induced kicks is crucial for producing a distribution of the inner binary-tertiary separations that is consistent with Gaia observations. Lastly, leveraging this consistency, we estimate that ~30% of solar-type stars in the local 200 parsec were born in triples.

The Milky Way has distinct structural stellar components linked to its formation and subsequent evolution, but disentangling them is nontrivial. With the recent availability of high-quality data for a large numbers of stars in the Milky Way, it is a natural next step for research in the evolution of the Galaxy to perform automated explorations with unsupervised methods of the structures hidden in the combination of large-scale spectroscopic, astrometric, and asteroseismic data sets. We determine precise stellar properties for 21,076 red giants, mainly spanning 2-15 kpc in Galactocentric radii, making it the largest sample of red giants with measured asteroseismic ages available to date. We explore the nature of different stellar structures in the Galactic disc by using Gaussian mixture models as an unsupervised clustering method to find substructure in the combined chemical, kinematic, and age subspace. The best-fit mixture model yields four distinct physical Galactic components in the stellar disc: the thin disc, the kinematically heated thin disc, the thick disc, and the stellar halo. We find hints of an age asymmetry between the Northern and Southern hemisphere and we measure the vertical and radial age gradient of the Galactic disc using the asteroseismic ages extended to further distances than previous studies.

AM CVn systems are ultra-compact binaries where a white dwarf accretes from a helium-rich degenerate or semi-degenerate donor. Some AM CVn systems will be among the loudest sources of gravitational waves for the upcoming Laser Interferometer Space Antenna (LISA), yet the formation channel of AM CVns remains uncertain. We report the study and characterisation of a new eclipsing AM CVn, SRGeJ045359.9+622444 (hereafter SRGeJ0453), discovered from a joint SRG/eROSITA and ZTF program to identify cataclysmic variables (CVs). We obtained optical photometry to confirm the eclipse of SRGeJ0453 and determine the orbital period to be $P_\textrm{orb} = 55.0802 \pm 0.0003$ min. We constrain the binary parameters by modeling the high-speed photometry and radial velocity curves and find $M_\textrm{donor} = 0.044 \pm0.024 M_{\odot}$ and $R_\textrm{donor}=0.078 \pm 0.012 R_{\odot}$. The X-ray spectrum is approximated by a power-law model with an unusually flat photon index of $\Gamma\sim 1$ previously seen in magnetic CVs with SRG/eROSITA, but verifying the magnetic nature of SRGeJ0453 requires further investigation. Optical spectroscopy suggests that the donor star of SRGeJ0453 could have initially been a He star or a He white dwarf. SRGeJ0453 is the ninth eclipsing AM CVn system published to date, and its lack of optical outbursts have made it elusive in previous surveys. The discovery of SRGeJ0453 using joint X-ray and optical surveys highlights the potential for discovering similar systems in the near future.

Binary pulsars can be used to probe Galactic potential gradients through calculating their line-of-sight accelerations. We present the first data release of direct line-of-sight acceleration measurements for 26 binary pulsars. We validate these data with a local acceleration model, and compare our results to those extracted from indirect methods. We find evidence for an acceleration gradient in agreement with these values, with our results indicating a local disk density of $\rho_\mathrm{d} = 0.064_{-0.033}^{+0.025} \ \mathrm{M_\odot}\mathrm{pc}^{-3}$. We also find evidence for unmodeled noise of unknown origin in our data set.

The rapidly increasing number of delta Scuti stars with regular patterns among their pulsation frequencies necessitates modelling tools to better understand the observations. Further, with a dozen identified modes per star, there is potential to make meaningful inferences on stellar structure using these young $\delta$ Sct stars. We compute and describe a grid of $>$800,000 stellar models from the early pre-main-sequence to roughly one third of the main-sequence lifetime, and calculate their pulsation frequencies. From these, we also calculate asteroseismic parameters and explore how those parameters change with mass, age, and metal mass fraction. We show that the large frequency separation, $\Delta\nu$, is insensitive to mass at the zero-age main sequence. In the frequency regime observed, the $\Delta\nu$ we measure (from modes with $n\sim5$--9) differs from the solar scaling relation by $\sim$13%. We find that the lowest radial order is often poorly modelled, perhaps indicating that the lower-order pressure modes contain further untapped potential for constraining the stellar structure. We also show that different nuclear reaction networks available in MESA can affect the pulsation frequencies of young $\delta$ Sct stars by as much as 5%. We apply the grid to five newly modelled stars, including two pre-main-sequence stars each with 15+ modes identified, and we make the grid available as a community resource.

Binary neutron star mergers probe the dense-matter equation of state (EoS) across a wide range of densities and temperatures, from the cold conditions of the inspiral to the high-temperature matter of the massive neutron star remnant. In this paper, we explore the sensitivity of neutron star mergers to uncertainties in the finite-temperature part of the EoS with a series of merger simulations performed in full general relativity. We expand on our previous work to explore the interplay between the thermal prescription and the stiffness of the zero-temperature EoS, which determines the compactness of the initial neutron stars. Using a phenomenological model of the particle effective mass, $M^*$, to calculate the finite-temperature part of the EoS, we perform merger simulations for a range of thermal prescriptions, together with two cold EoSs that predict either compact or large-radius initial neutron stars. We report on how the choice of $M^*$-parameters influences the thermal properties of the post-merger remnant, and how this varies for stars with different initial stellar compactness. We characterize the post-merger gravitational wave signals, and find differences in the peak frequencies of up to 190 Hz depending on the choice of $M^*$-parameters. Finally, we find that the total dynamical ejecta is in general only weakly sensitive to the thermal prescription, but that a particular combination of $M^*$-parameters, together with a soft cold EoS, can lead to significant enhancements in the ejecta.

SDO/HMI observations reveal a class of solar flares with substantial energy and momentum impacts in the photosphere, concurrent with white-light emission and helioseismic responses, known as sunquakes. Previous radiative hydrodynamic modeling has demonstrated the challenges of explaining sunquakes in the framework of the standard flare model of `electron beam' heating. One of the possibilities to explain the sunquakes and other signatures of the photospheric impact is to consider additional heating mechanisms involved in solar flares, for example, via flare-accelerated protons. In this work, we analyze a set of single-loop RADYN+FP simulations where the atmosphere is heated by non-thermal power-law-distributed proton beams which can penetrate deeper than the electron beams into the low atmospheric layers. Using the output of the RADYN models, we calculate synthetic Fe I 6173 A line Stokes profiles and from those the line-of-sight (LOS) observables of the SDO/HMI instrument, as well as the 3D helioseismic response and compare them with the corresponding observational characteristics. These initial results show that the models with proton beam heating can produce the enhancement of the HMI continuum observable and explain qualitatively generation of sunquakes. The continuum observable enhancement is evident in all models but is more prominent in ones with $E_{c}\geq$500 keV. In contrast, the models with $E_{c}\leq$100 keV provide a stronger sunquake-like helioseismic impact according to the 3D acoustic modeling, suggesting that low-energy (deka- and hecto-keV) protons have an important role in the generation of sunquakes.

Upcoming studies at the Large Hadron Collider (LHC) aim to extend the rapidity coverage in measurements of the production cross section of antinuclei ${\rm \bar d}$ and $\overline{^3\rm He}$. We illustrate the impact of such studies on cosmic ray (CR) flux predictions, important, in turn, for the interpretation of results from CR experiments. We show that, in terms of the rapidity effect, covering the range $|y|<1.5$ at the LHC should be sufficient for the astrophysical CR calculation. Important extrapolation remains in other aspects of the problem, notably $\sqrt{s}$.

We present results exploring the role that probabilistic deep learning models play in cosmology from large-scale astronomical surveys through photometric redshift (photo-z) estimation. Photo-z uncertainty estimates are critical for the science goals of upcoming large-scale surveys such as LSST, however common machine learning methods typically provide only point estimates and lack uncertainties on predictions. We turn to Bayesian neural networks (BNNs) as a promising way to provide accurate predictions of redshift values with uncertainty estimates. We have compiled a new galaxy training data set from the Hyper Suprime-Cam Survey with grizy photometry, which is designed to be a smaller scale version of large surveys like LSST. We use this data set to investigate the performance of a neural network (NN) and a probabilistic BNN for photo-z estimation and evaluate their performance with respect to LSST photo-z science requirements. We also examine the utility of photo-z uncertainties as a means to reduce catastrophic outlier estimates. The BNN model outputs the estimate in the form of a Gaussian probability distribution. We use the mean and standard deviation as the redshift estimate and uncertainty, respectively. We find that the BNN can produce accurate uncertainties. Using a coverage test, we find excellent agreement with expectation -- 67.2% of galaxies between 0 < 2.5 have 1-$\sigma$ uncertainties that cover the spectroscopic value. We find the BNN meets 2/3 of the LSST photo-z science requirements in the range 0 < z < 2.5 and generally outperforms the alternative photo-z methods considered here on the same data.

Galaxy cluster masses derived from observations of weak lensing suffer from a number of biases affecting the accuracy of mass-observable relations calibrated from such observations. In particular, the choice of the cluster center plays a prominent role in biasing inferred masses. In the past, empirical miscentring distributions have been used to address this issue. Using hydro-dynamical simulations, we aim to test the accuracy of weak lensing mass bias predictions based on such miscentring distributions by comparing the results to mass biases computed directly using intra-cluster medium (ICM)-based centers from the same simulation. We construct models for fitting masses to both centered and miscentered Navarro-Frenk-White profiles of reduced shear, and model the resulting distributions of mass bias with normal and log-normal distributions. We find that the standard approach of using miscentring distributions leads to an over-estimation of cluster masses at levels of between 2\% and 6\% when compared to the analysis in which actual simulated ICM centers are used, even when the underlying miscentring distributions match in terms of the miscentring amplitude. We find that neither log-normal nor normal distributions are generally reliable for approximating the shapes of the mass bias distributions, regardless of whether a centered or miscentered radial model is used.

Binary and dual active galactic nuclei (AGN) are an important observational tool for studying the formation and dynamical evolution of galaxies and supermassive black holes (SMBHs). An entirely new method for identifying possible AGN pairs makes use of the exquisite positional accuracy of Gaia to detect astrometrically-variable quasars, in tandem with the high spatial resolution of the Karl G. Jansky Very Large Array (VLA). We present a new pilot study of radio observations of 18 quasars (0.8 < z < 2.9), selected from the SDSS DR16Q and matched with the Gaia DR3. All 18 targets are identified by their excess astrometric noise in Gaia. We targeted these 18 quasars with the VLA at 2-4 GHz (S-band) and 8-12 GHz (X-band), providing resolutions of 0.65" and 0.2", respectively, in order to constrain the origin of this variability. We combine these data with ancillary radio survey data and perform radio spectral modeling. The new observations are used to constrain the driver of the excess astrometric noise. We find that ~39% of the target sample is likely to be either candidate dual AGN or gravitationally lensed quasars. Ultimately, we use this new strategy to help identify and understand this sample of astrometrically-variable quasars, demonstrate the potential of this method for systematically identifying kpc-scale dual quasars.

The evolution of exoplanetary systems with a close-in planet is ruled by the tides mutually raised on the two bodies and by the magnetic braking of the host star. This paper deals with consequences of this evolution and some features that can be observed in the distribution of the systems two main periods: the orbital periods and the stars rotational periods. The results of the simulations are compared to plots showing both periods as determined from the light curves of a large number of Kepler objects of interest. These plots show important irregularities as a dearth of systems in some regions and accumulations of hot Jupiters in others. It is shown that the accumulation of short-period hot Jupiters around stars with rotation periods close to 25 days results from the evolution of the systems under the joint action of tides and braking, and requires a relaxation factor for solar-type stars of around $10 \, s^{-1}$.

Using CORSIKA simulations, we investigate the mass sensitivity of cosmic-ray air-shower observables for sites at the South Pole and Malarg\"ue, Argentina, the respective locations of the IceCube Neutrino Observatory and the Pierre Auger Observatory. Exact knowledge of observables from air-shower simulations was used to study the event-by-event mass separation between proton, helium, oxygen, and iron primary cosmic rays with a Fisher linear discriminant analysis. Dependencies on the observation site as well as the energy and zenith angle of the primary particle were studied in the ranges from $10^{16.0}-10^{18.5}\,$eV and $0^\circ$ to $60^\circ$: they are mostly weak and do not change the qualitative results. Promising proton-iron mass separation is achieved using combined knowledge of all studied observables, also when typical reconstruction uncertainties are accounted for. However, even with exact measurements, event-by-event separation of intermediate-mass nuclei is challenging and better methods than the Fisher discriminant and/or the inclusion of additional observables will be needed. As an individual observable, high-energy muons ($> 500\,$GeV) provide the best event-by-event mass discrimination, but the combination of muons of any energy and $X_{\text{max}}$ provides already a high event-by-event separation between proton-iron primaries at both sites. We also confirm that the asymmetry and width parameters of the air-shower longitudinal profile, $R$ and $L$, are mass sensitive. Only $R$ seems to be suitable for event-by-event mass separation, but $L$ can potentially be used to statistically determine the proton-helium ratio. Overall, our results motivate the coincident measurement of several air-shower observables, including at least $X_{\text{max}}$ and the sizes of the muonic and electromagnetic shower components, for the next generation of air-shower experiments.

Astronomers have studied the Vela pulsar (PSR~J0835$-$4510) for decades. This study analyses almost one hundred hours of single-pulse data collected over three consecutive days from 2016 and 2020. The work investigates the fascinating phenomena of the earlier arrival of brighter pulses with their increase in peak intensity. We found a hyperbolic relation between them by constructing integrated pulse profiles using flux density intervals and examining the relationship between pulse arrival time and intensity. We identified a phase limit of $-0.85~\pm 0.0109$~ms for the earliest arrival of the brightest pulses. This study offers exciting prospects for further exploring the emission regions responsible for the Vela pulsar's regular and giant micro-pulses.

Proximity observations of (162173) Ryugu by the telescopic Optical Navigation Camera onboard Hayabusa2 and (101955) Bennu by MapCam onboard Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer found opposite spectral trends of space weathering on these carbonaceous asteroids. Whether the space weathering trends on these asteroids evolved from the same starting spectra would place an important constraint for understanding their relation. However, systematic error between data obtained by the two imagers needed to be reduced for accurate comparison. To resolve this problem, we cross calibrated albedo and color data using the Moon as the common standard. We show that the cross-calibrated reflectance can be obtained by upscaling the pre-cross-calibrated reflectance of Bennu by 12 +/- 2% at v-band, reducing the systematic errors down to 2%. The cross-calibrated data show that Bennu is brighter by 16 +/- 2% at v-band and bluer in spectral slope by 0.19 +/- 0.05 (/um) than Ryugu. The spectra of fresh craters on Ryugu and Bennu before cross calibration appeared to follow two parallel trend lines with offset, but they converged to a single trend after cross calibration. Such a post-cross-calibration perspective raise the possibility that Ryugu and Bennu evolved from materials with similar visible spectra but evolved in diverging directions by space weathering. The divergent evolution can be caused by the difference in space weathering dose/process and/or composition of the starting material. Thus, comparing the composition of samples returned from Ryugu and Bennu may change the way we interpret the spectral variation of C-complex asteroids.

We describe a 2.5D MHD simulation describing the evolution of cool jets triggered by initial vertical velocity perturbations in the solar chromosphere. We implement random velocity pulses of amplitude 20-50 km/s between 1 Mm and 1.5 Mm, along with various switch-off periods between 50 s and 300 s. The applied vertical velocity pulses create a series of magnetoacoustic shocks steepening above TR. These shocks interact with each other in the inner corona, leading to complex localized velocity fields. The upward propagation of such perturbations creates low-pressure regions behind them, which propel a variety of cool jets and plasma flows. We study the transverse oscillations of a representative cool jet J1 , which moves up to the height of 6.2 Mm above the TR from its origin point. During its evolution, the plasma flows make the spine of jet J1 radially inhomogeneous, which is visible in the density and Alfv\'en speed smoothly varying across the jet. The highly dense J1 supports the propagating transverse wave of period of approximately 195 s with a phase speed of about 125 km/s. In the distance-time map of density, it is manifested as a transverse kink wave. However, the careful investigation of the distance-time maps of the x- and z-components of velocity reveals that these transverse waves are actually the mixed Alfv\'enic modes. The transverse wave shows evidence of damping in the jet. We conclude that the cross-field structuring of the density and characteristic Alfv\'en speed within J1 causes the onset of the resonant conversion and leakage of the wave energy outward to dissipate these transverse oscillations via resonant absorption. The wave energy flux is estimated as approximately of 1.0 x 10^6 ergs cm^{-2} s^{-1}. This energy, if it dissipates through the resonant absorption into the corona where the jet is propagated, is sufficient energy for the localized coronal heating.

Richardson-Lucy (RL) deconvolution is one of the classical methods widely used in X-ray astronomy and other areas. Amid recent progress in image processing, RL deconvolution still leaves much room for improvement under a realistic situations. One direction is to include the positional dependence of a point-spread function (PSF), so-called RL deconvolution with a spatially variant PSF (RL$_{\rm{sv}}$). Another is the method of estimating a reliable number of iterations and their associated uncertainties. We developed a practical method that incorporates the RL$_{\rm{sv}}$ algorithm and the estimation of uncertainties. As a typical example of bright and high-resolution images, the Chandra X-ray image of the supernova remnant Cassiopeia~A was used in this paper. RL$_{\rm{sv}}$ deconvolution enables us to uncover the smeared features in the forward/backward shocks and jet-like structures. We constructed a method to predict the appropriate number of iterations by using statistical fluctuation of the observed images. Furthermore, the uncertainties were estimated by error propagation from the last iteration, which was phenomenologically tested with the observed data. Thus, our method is a practically efficient framework to evaluate the time evolution of the remnants and their fine structures embedded in high-resolution X-ray images.

Context. Dynamically linking a meteor shower with its parent body is challenging, and chaos in the dynamics of meteoroid streams may be one of the reasons. For a robust identification of parent bodies, it is therefore necessary to quantify the amount of chaos involved in the evolution of meteoroid streams. Aims. The characterisation of chaos in meteoroid streams thanks to chaos maps is still a new field of study. We aim to study two very different meteoroid streams, the Draconids and the Leonids, in order to obtain a general view of this topic. Methods. We use the method developed in a previous paper dedicated to Geminids, drawing chaos maps with the orthogonal fast Lyapunov indicator. We choose four particle size ranges to investigate the effect of non-gravitational forces. As the dynamics is structured by mean-motion resonances with planets, we compute the locations and widths of the resonances at play. We use semi-analytical formulas valid for any eccentricity and inclination and an arbitrary number of planets. Results. We pinpoint which mean-motion resonances with Jupiter play a major role in the dynamics of each meteoroid stream. We show how those resonances tend to trap mostly large particles, preventing them from meeting with Jupiter. We also study particles managing to escape those resonances, thanks to the gravitational perturbation of Saturn for example. Finally, we explain why non-gravitational forces do not disturb the dynamics much, contrary to what is observed for the Geminids.

Upcoming wide-field spectroscopic surveys will observe galaxies in a range of cosmic web environments in and around galaxy clusters. In this paper, we test and quantify how successfully we will be able to identify the environment of individual galaxies in the vicinity of massive galaxy clusters, reaching out to $\sim5R_{200}$ into the clusters' infall region. We focus on the WEAVE Wide Field Cluster Survey (WWFCS), but the methods we develop can be easily generalised to any similar spectroscopic survey. Using numerical simulations of a large sample of massive galaxy clusters from \textsc{TheThreeHundred} project, we produce mock observations that take into account the selection effects and observational constraints imposed by the WWFCS. We then compare the `true' environment of each galaxy derived from the simulations (cluster core, filament, and neither core nor filament, {``NCF''}) with the one derived from the observational data, where only galaxy sky positions and spectroscopic redshifts will be available. We find that, while cluster core galaxy samples can be built with a high level of completeness and moderate contamination, the filament and NCF galaxy samples will be significantly contaminated and incomplete due to projection effects exacerbated by the galaxies' peculiar velocities. We conclude that, in the infall regions surrounding massive galaxy clusters, associating galaxies with the correct cosmic web environment is highly uncertain. However, with large enough spectroscopic samples like the ones the WWFCS will provide (thousands of galaxies per cluster, {out to $5R_{200}$}), and the correct statistical treatment that takes into account the probabilities we provide here, we expect we will be able to extract robust and well-quantified conclusions on the relationship between galaxy properties and their environment.

The High-Energy Ray Observatory (HERO) is a space experiment based on a heavy ionization calorimeter for direct study of cosmic rays. The effective geometrical factor of the apparatus varies from 12 to 60 m$^2$sr for protons depending on the weight of the calorimeter from 10 to 70 tons. During the exposure for $\sim$5 years this mission will make it possible to measure energy spectra of all abundant cosmic ray nuclei in the knee region ($\sim$3 PeV) with individual resolution of charges with energy resolution better than 30\% and provide useful information to solve the puzzle of the cosmic ray knee origin. HERO mission will make it also possible to measure energy spectra of cosmic rays nuclei for energies 1-1000 TeV with very high precision and energy resolution (up to 3\% for calorimeter 70 tons) and study the fine structure of the spectra. The planned experiment launch is no earlier than 2029.

The fine structure of solar prominences is made by thin threads that outline the magnetic field lines. Observations show that transverse waves of Alfv\'enic nature are ubiquitous in prominence threads. These waves are driven at the photosphere and propagate to prominences suspended in the corona. Heating due to Alfv\'en wave dissipation could be a relevant mechanism in the cool and partially ionized prominence plasma. We explore the construction of 1D equilibrium models of prominence thin threads that satisfy energy balance between radiative losses, thermal conduction, and Alfv\'en wave heating. We assume the presence of a broadband driver at the photosphere that launches Alfv\'en waves towards the prominence. An iterative method is implemented, in which the energy balance equation and the Alfv\'en wave equation are consecutively solved. From the energy balance equation and considering no wave heating initially, we compute the equilibrium profiles along the thread of the temperature, density, ionisation fraction. We use the Alfv\'en wave equation to compute the wave heating rate, which is then put back in the energy balance equation to obtain new equilibrium profiles. The process is repeated until convergence to a self-consistent thread model heated by Alfv\'en waves is achieved. We have obtained equilibrium models composed of a cold and dense thread, a extremely thin PCTR, and an extended coronal region. The length of the cold thread decreases with the temperature at the prominence core and increases with the Alfv\'en wave energy flux. Equilibrium models are not possible for sufficiently large wave energy fluxes when the wave heating rate inside the cold thread becomes larger than radiative losses. The maximum value of the wave energy flux that allows an equilibrium depends on the prominence core temperature. This constrains the existence of equilibria in realistic conditions.

The leading spectrographs used for exoplanets' sceince offer online data reduction softwares (DRS) that yield as an ancillary result the full-width at half-maximum (FWHM) of the cross-correlation function (CCF) that is used to estimate the radial velocity of the host star. The FWHM also contains information on the stellar projected rotational velocity vsini We wanted to establish a simple relationship to derive the vsini directly from the FWHM computed by the HARPS-N DRS in the case of slow-rotating solar-like stars. This may also help to recover the stellar inclination i, which in turn affects the exoplanets' parameters. We selected stars with an inclination of the spin axis compatible with 90 deg by looking at exoplanetary transiting systems with known small sky-projected obliquity: for these stars, we can presume that vsini is equal to stellar equatorial velocity veq. We derived their rotational periods from photometric time-series and their radii from SED fitting. This allowed us to recover their veq, which we could compare to the FWHM values of the CCFs obtained both with G2 and K5 spectral type masks. We obtained an empirical relation for each mask, useful for slow rotators (FWHM < 20 km/s). We applied them to 273 exoplanet-host stars observed with HARPS-N, obtaining homogeneous vsini measurements. We compared our results with the literature ones to confirm the reliability of our work, and we found a good agreement with the values found with more sophisticated methods for stars with log g > 3.5. We also tried our relations on HARPS and SOPHIE data, and we conclude that they can be used also on FWHM derived by HARPS DRS with G2 and K5 mask, and they may be adapted to the SOPHIE data as long as the spectra are taken in the high-resolution mode. We were also able to recover or constrain i for 12 objects with no prior vsini estimation.

The thermal Sunyaev-Zeldovich (SZ) effect serves as a direct potential probe of the energetic outflows from quasars that are responsible for heating the intergalactic medium. In this work, we use the GIZMO meshless finite mass hydrodynamic cosmological simulation SIMBA (Dave et al. 2019), which includes different prescriptions for quasar feedback, to compute the SZ effect arising from different feedback modes. From these theoretical simulations, we perform mock observations of the Atacama Large Millimeter Array (ALMA) in four bands (320 GHz, 135 GHZ, 100 GHz and 42 GHz) to characterize the feasibility of direct detection of the quasar SZ signal. Our results show that for all the systems we get an enhancement of the SZ signal, when there is radiative feedback, while the signal gets suppressed when the jet mode of feedback is introduced in the simulations. Our mock ALMA maps reveal that, with the current prescription of jet feedback, the signal goes below the detection threshold of ALMA. We also find that the signal is higher for high redshift systems, making it possible for ALMA and cross SZ-X-ray studies to disentangle the varying modes of quasar feedback and their relative importance in the cosmological context.

Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum. The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas. With nearly $60,000$ antennas envisioned to be used for air shower studies at the Square Kilometre Array (SKA), simulating all of these would come at unreasonable costs. We present an interpolation algorithm to reconstruct the full pulse time series at any position in the radio footprint, from a set of antennas simulated on a polar grid. Relying on Fourier series representations and cubic splines, it significantly improves on existing linear methods. We show that simulating about 200 antennas is sufficient for high-precision analysis in the SKA era, including e.g. interferometry which relies on accurate pulse shapes and timings. We therefore propose the interpolation algorithm and its implementation as a useful extension of radio simulation codes, to limit computational effort while retaining accuracy.

The design of a downlink communication system for returning scientific data from an interstellar flyby probe is reviewed in this tutorial white paper. It its assumed that the probe is ballistic, and data is downloaded during a period following encounter with the target star and its exoplanet(s). Performance indices of interest to scientific investigators include the total launch-to-completion data latency and the total volume of data reliably recovered. Issues considered include the interaction between the speed and mass of the probe and the duration of downlink transmission. Optical communication using pulse-position modulation (PPM) with error-correction coding (ECC) is assumed. A very large receiver collection area on or near Earth is composed of individual incoherently-combined diffraction-limited apertures. Other important issues in the design including transmit and receive pointing accuracy and beam size and receiver field of view are reviewed. Numerical examples assume a mission to Proxima Centauri (the nearest star to our Sun) initially launched by directed-energy propulsion from the vicinity of Earth.

The G31.41+0.31 Unbiased ALMA sPectral Observational Survey (GUAPOS) project targets the hot molecular core (HMC) G31.41+0.31 (G31), to unveil the complex chemistry of one of the most chemically rich high-mass star-forming regions outside the Galactic Center (GC). In the third paper of the project, we present a study of nine O-bearing (CH$_3$OH, $^{13}$CH$_3$OH, CH$_3^{18}$OH, CH$_3$CHO, CH$_3$OCH$_3$, CH$_3$COCH$_3$ , C$_2$H$_5$OH, aGg'-(CH$_2$OH)$_2$, and gGg'-(CH$_2$OH)$_2$) and six N-bearing (CH$_3$CN, $^{13}$CH$_3$CN, CH$_3^{13}$CN, C$_2$H$_3$CN, C$_2$H$_5$CN, and C$_2$H$_5^{13}$CN) complex organic molecules toward G31. The aim of this work is to characterize the abundances in one of the most chemically-rich hot molecular cores outside the GC and to search for a possible chemical segregation between O-bearing and N-bearing species in G31, which hosts four compact sources as seen with higher angular resolution data. Observations were carried out with the interferometer ALMA and covered the entire Band 3 from 84 to 116 GHz ($\sim 32$ GHz bandwidth) with an angular resolution of $1.2''$ ($\sim4400\,\mathrm{au}$). The spectrum has been analyzed with the tool SLIM of MADCUBA to determine the physical parameters of the emitting gas. Moreover, we have analyzed the morphology of the emission of the molecular species. We have compared the abundances w.r.t methanol of COMs in G31 with other twenty-seven sources, including other hot molecular cores inside and outside the Galactic Center, hot corinos, shocked regions, envelopes around young stellar objects, and quiescent molecular clouds, and with chemical models. Different species peak at slightly different positions, and this, together with the different central velocities of the lines obtained from the spectral fitting, point to chemical differentiation of selected O-bearing species.

The asteroseismic analysis of stellar power density spectra is often computationally expensive. The models used in the analysis may use several dozen parameters to accurately describe features in the spectra caused by oscillation modes and surface granulation. Many parameters are often highly correlated, making the parameter space difficult to quickly and accurately sample. They are, however, all dependent on a smaller set of parameters, namely the fundamental stellar properties. We aim to leverage this to simplify the process of sampling the model parameter space for the asteroseismic analysis of solar-like oscillators, with an emphasis on mode identification. Using a large set of previous observations, we applied principal component analysis to the sample covariance matrix to select a new basis on which to sample the model parameters. Selecting the subset of basis vectors that explains the majority of the sample variance, we redefine the model parameter prior probability density distributions in terms of a smaller set of latent parameters. We are able to reduce the dimensionality of the sampled parameter space by a factor of two to three. The number of latent parameters needed to accurately model the stellar oscillation spectra cannot be determined exactly but is likely only between four and six. Using two latent parameters, the method is able to describe the bulk features of the oscillation spectrum, while including more latent parameters allows for a frequency precision better than $\approx10\%$ of the small frequency separation for a given target. We find that sampling a lower-rank latent parameter space still allows for accurate mode identification and parameter estimation on solar-like oscillators over a wide range of evolutionary stages. This allows for the potential to increase the complexity of spectrum models without a corresponding increase in computational expense.

In this paper, we presented the 23.3 years of pulsar timing results of PSR J1456-6413 based on the observation of Parkes 64m radio telescope. We detected two new glitches at MJD 57093(3) and 59060(12) and confirmed its first glitch at MJD 54554(10). Using the "Cholesky" timing analysis method, we have determined its position, proper motion, and two-dimensional transverse velocities from the data segments before and after the second glitch, respectively. Furthermore, we detected exponential recovery behavior after the first glitch, with a recovery time scale of approximately 200 days and a corresponding exponential recovery factor Q of approximately 0.15(2), while no exponential recovery was detected for the other two glitches. More interestingly, we found that the leading component of the integral pulse profile after the second glitch became stronger, while the main component became weaker. Our results will expand the sample of pulsars with magnetosphere fluctuation triggered by the glitch event.

The production mechanism of repeating fast radio bursts (FRBs) is still a mystery, and correlations between burst occurrence times and energies may provide important clues to elucidate it. While time correlation studies of FRBs have been mainly performed using wait time distributions, here we report the results of a correlation function analysis of repeating FRBs in the two-dimensional space of time and energy. We find the universal laws on temporal correlations by analyzing nearly 7,000 bursts reported in the literature for the three most active sources of FRB 20121102A, 20201124A, and 20220912A. A clear power-law signal of the correlation function is seen, extending to the typical burst duration ($\sim$ 10 msec) toward shorter time intervals ($\Delta t)$. The correlation function indicates that every single burst has about a 10--60\% chance of producing an aftershock at a rate decaying by a power-law as $\propto (\Delta t)^{-p}$ with $p =$ 1.5--2.5, like the Omori-Utsu law of earthquakes. The correlated aftershock rate is stable regardless of source activity changes, and there is no correlation between emitted energy and $\Delta t$. We demonstrate that all these properties are quantitatively common to earthquakes, but different from solar flares in many aspects, by applying the same analysis method for the data on these phenomena. These results suggest that repeater FRBs are a phenomenon in which energy stored in rigid neutron star crusts is released by seismic activity. This may provide a new opportunity for future studies to explore the physical properties of the neutron star crust.

For existing and future projects dedicated to measuring precise radial velocities (RVs), we have created an open-source, flexible data reduction software to extract RVs from \'echelle spectra via the iodine (I$_2$) absorption cell method. The software, called $pyodine$, is completely written in Python and has been built in a modular structure to allow for easy adaptation to different instruments. We present the fundamental concepts employed by $pyodine$, which build on existing I$_2$ reduction codes, and give an overview of the software's structure. We adapted $pyodine$ to two instruments, Hertzsprung SONG located at Teide Observatory (SONG hereafter) and the Hamilton spectrograph at Lick Observatory (Lick hereafter), and demonstrate the code's flexibility and its performance on spectra from these facilities. Both for SONG and Lick data, the $pyodine$ results generally match the RV precision achieved by the dedicated instrument pipelines. Notably, our code reaches a precision of roughly $0.69 \,m\,s^{-1}$ on a short-term solar time series of SONG spectra, and confirms the planet-induced RV variations of the star HIP~36616 on spectra from SONG and Lick. Using the solar spectra, we also demonstrate the capabilities of our software in extracting velocity time series from single absorption lines. A probable instrumental effect of SONG is still visible in the $pyodine$ RVs, despite being a bit damped as compared to the original results. With $pyodine$ we prove the feasibility of a highly precise, yet instrument-flexible I$_2$ reduction software, and in the future the code will be part of the dedicated data reduction pipelines for the SONG network and the Waltz telescope project in Heidelberg.

We analyze the X-ray properties for a sample of 23 high probability AGN candidates with ultraviolet variability identified in Wasleske et al. (2022). Using data from the Chandra X-ray Observatory and the XMM-Newton Observatory, we find 11/23 nuclei are X-ray detected. We use SED modeling to compute star formation rates and show that the X-ray luminosities are typically in excess of the X-ray emission expected from star formation by at least an order of magnitude. Interestingly, this sample shows a diversity of optical spectroscopic properties. We explore possible reasons for why some objects lack optical spectroscopic signatures of black hole activity while still being UV variable and X-ray bright. We find that host galaxy stellar emission and obscuration from gas and dust are all potential factors. We study where this sample falls on relationships such as $\alpha_{\rm OX}-L_{2500}$ and $L_{X}-L_{IR}$ and find that some of the sample falls outside the typical scatter for these relations, indicating they differ from the standard quasar population. With the diversity of optical spectroscopic signatures and varying impacts of dust and stellar emissions on our sample, these results emphasizes the strength of variability in selecting the most complete set of AGN, regardless of other host galaxy properties.

The presence of a plethora of light spin 0 and spin 1 fields is motivated in a number of BSM scenarios, such as the axiverse. The study of the interactions of such light bosonic fields with the Standard Model has focused mostly on interactions involving only one such field, such as the axion ($\phi$) coupling to photons, $\phi F \tilde F$, or the kinetic mixing between photon and the dark photon, $ F F_D$. In this work, we continue the exploration of interactions involving two light BSM fields and the standard model, focusing on the mixed axion-photon-dark-photon interaction $\phi F \tilde F_D$. If either the axion or dark photon are dark matter, we show that this interaction leads to conversion of the CMB photons into a dark sector particle, leading to a distortion in the CMB spectrum. We present the details of these unique distortion signatures and the resulting constraints on the $\phi F \tilde F_D$ coupling. In particular, we find that for a wide range of masses, the constraints from these effect are stronger than on the more widely studied axion-photon coupling.

The QCD axion has been postulated to exist because it solves the strong CP problem. Furthermore, if it exists axions should be created in the early Universe and could account for all the observed dark matter. In particular, axion masses of order $10^{-10}$ to $10^{-7}$ eV correspond to axions in the vicinity of the GUT-scale. In this mass range many experiments have been proposed to search for the axion through the standard QED coupling parameter $g_{a\gamma\gamma}$. Recently axion electrodynamics has been expanded to include two more coupling parameters, $g_{aEM}$ and $g_{aMM}$, which could arise if heavy magnetic monopoles exist. In this work we show that both $g_{aMM}$ and $g_{aEM}$ may be searched for using a high voltage capacitor. Since the experiment is not sensitive to $g_{a\gamma\gamma}$, it gives a new way to search for effects of heavy monopoles if the GUT-scale axion is shown to exist, or to simultaneously search for both the axion and the monopole at the same time.

We revisit stellar constraints on dark photons. We undertake dynamical stellar evolution simulations which incorporate the resonant and off-resonant production of transverse and longitudinal dark photons. We compare our results with observables derived from measurements of globular cluster populations, obtaining new constraints based on the luminosity of the tip of the red-giant branch (RGB), the ratio of populations of RGB to horizontal branch (HB) stars (the $R$-parameter), and the ratio of asymptotic giant branch to HB stars (the $R_2$-parameter). We find that previous bounds derived from static stellar models do not capture the effects of the resonant production of light dark photons leading to overly conservative constraints, and that they over-estimate the effects of heavier dark photons on the RGB-tip luminosity. This leads to differences in the constraints of up to an order of magnitude in the kinetic mixing parameter.

Ultralight dark matter exhibits an order-one density fluctuation over the spatial scale of its wavelength. These fluctuations gravitationally interact with gravitational wave interferometers, leading to an additional noise floor or signals. We investigate the ultralight dark matter-induced effects in the gravitational wave interferometers. We perform a systematic computation of the power spectrum of ultralight dark matter in interferometers. We show that the ultralight dark matter-induced effect is most relevant for the interferometers with long baseline and that it only constitutes a sub-leading noise floor compared to the estimated noise level in the case of Laser Interferometer Space Antenna or future interferometers with an arm-length comparable to a few astronomical units. Gravitational wave interferometers can then place upper limits on the ultralight dark matter density in the solar system. We find that, under certain assumptions, future interferometers with AU-scale arm-length might probe the dark matter density a few hundred times the local dark matter density, which is measured over a much larger spatial scale.

We apply quadratic $f(R)=R+\epsilon R^2$ field equations, where $\epsilon$ has a dimension [L$^2$], to static spherical stellar model. We assume the interior configuration is determined by Krori-Barua ansatz and additionally the fluid is anisotropic. Using the astrophysical measurements of the pulsar PSR J0740+6620 as inferred by NICER and XMM observations, we determine $\epsilon\approx \pm 3$ km$^2$. We show that the model can provide a stable configuration of the pulsar PSR J0740+6620 in both geometrical and physical sectors. We show that the Krori-Barua ansatz within $f(R)$ quadratic gravity provides semi-analytical relations between radial, $p_r$, and tangential, $p_t$, pressures and density $\rho$ which can be expressed as $p_r\approx v_r^2 (\rho-\rho_1)$ and $p_r\approx v_t^2 (\rho-\rho_2)$, where $v_r$ ($v_t$) is the sound speed in radial (tangential) direction, $\rho_1=\rho_s$ (surface density) and $\rho_2$ are completely determined in terms of the model parameters. These relations are in agreement with the best-fit equations of state as obtained in the present study. We further put the upper limit on the compactness, which satisfies the $f(R)$ modified Buchdahl limit. Interestingly, the quadratic $f(R)$ gravity with negative $\epsilon$ naturally restricts the maximum compactness to values lower than Buchdahl limit, unlike the GR or $f(R)$ gravity with positive $\epsilon$ where the compactness can arbitrarily approach the black hole limit $C\to 1$. The model predicts a core density a few times the saturation nuclear density $\rho_{\text{nuc}} = 2.7\times 10^{14}$ g/cm$^3$, and a surface density $\rho_s > \rho_{\text{nuc}}$. We provide the mass-radius diagram corresponding to the obtained boundary density which has been shown to be in agreement with other observations.

The characterization of Earth-like exoplanets and precision tests of cosmological models using next-generation telescopes such as the ELT will demand precise calibration of astrophysical spectrographs in the visible region, where stellar absorption lines are most abundant. Astrocombs--lasers providing a broadband sequence of ultra-narrow, drift-free, regularly spaced optical frequencies on a multi-GHz grid--promise an atomically-traceable, versatile calibration scale, but their realization is challenging because of the need for ultra-broadband frequency conversion of mode-locked infrared lasers into the blue-green region. Here, we introduce a new concept achieving a broad, continuous spectrum by combining second-harmonic generation and sum-frequency-mixing in an aperiodically-poled MgO:PPLN waveguide to generate gap-free 390-520 nm light from a 1 GHz Ti:sapphire laser frequency comb. We lock a low-dispersion Fabry-Perot etalon to extract a sub-comb of bandwidth from 392-472 nm with a spacing of 30 GHz, visualizing the thousands of resulting comb modes on a high resolution cross-dispersion spectrograph. Complementary experimental data and simulations demonstrate the effectiveness of the approach for eliminating the spectral gaps present in second-harmonic-only conversion, in which weaker fundamental frequencies are suppressed by the quadratic \{chi}^((2)) nonlinearity. Requiring only ~100 pJ pulse energies, our concept establishes a practical new route to broadband UV-visible generation at GHz repetition rates.

We study the bulk viscosity of moderately hot and dense, neutrino-transparent relativistic $npe\mu$ matter arising from weak-interaction direct Urca processes. This work parallels our recent study of the bulk viscosity of $npe\mu$ matter with a trapped neutrino component. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations -- DDME2 (which does not allow for the low-temperature direct-Urca process at any density) and NL3 (which allows for low-temperature direct-Urca process above a low-density threshold). We compute the equilibration rates of Urca processes of neutron decay and lepton capture, as well as the rate of the muon decay, and find that the muon decay process is subdominant to the Urca processes at temperatures $T\geq 3$~MeV in the case of DDME2 model and $T\geq 1$~MeV in the case of NL3 model. Thus, the Urca-process-driven bulk viscosity is computed with the assumption that pure leptonic reactions are frozen. As a result the electronic and muonic Urca channels contribute to the bulk viscosity independently and at certain densities the bulk viscosity of $npe\mu$ matter shows a double-peak structure as a function of temperature instead of the standard one-peak (resonant) form. In the final step, we estimate the damping timescales of density oscillations by the bulk viscosity. We find that, \eg, at a typical oscillation frequency $f=1$~kHz, the damping of oscillation is most efficient at temperatures $3\leq T\leq 5$~MeV and densities $n_B\leq 2n_0$ where they can affect the evolution of the post-merger object.

The nanohertz stochastic gravitational wave background (SGWB) is an excellent early universe laboratory for testing the fundamental properties of gravity. In this letter, we elucidate on the full potential of pulsar timing array (PTA) by utilizing cosmic variance-limited, or rather experimental noise-free, correlation measurements to understand the SGWB and by extension gravity. We show that measurements of the angular power spectrum play a pivotal role in the PTA precision era for scientific inferencing. In particular, we illustrate that cosmic variance-limited measurements of the first few power spectrum multipoles enable us to clearly set apart general relativity from alternative theories of gravity.

Astronomical precision spectroscopy underpins searches for life beyond Earth, direct observation of the expanding Universe and constraining the potential variability of physical constants across cosmological scales. Laser frequency combs can provide the critically required accurate and precise calibration to the astronomical spectrographs. For cosmological studies, extending the calibration with such astrocombs to the ultraviolet spectral range is highly desirable, however, strong material dispersion and large spectral separation from the established infrared laser oscillators have made this exceedingly challenging. Here, we demonstrate for the first time astronomical spectrograph calibrations with an astrocomb in the ultraviolet spectral range below 400 nm. This is accomplished via chip-integrated highly nonlinear photonics in periodically-poled, nano-fabricated lithium niobate waveguides in conjunction with a robust infrared electro-optic comb generator, as well as a chip-integrated microresonator comb. These results demonstrate a viable route towards astronomical precision spectroscopy in the ultraviolet and may contribute to unlocking the full potential of next generation ground- and future space-based astronomical instruments.