Papers for Thursday, Jun 08 2023

Capturing planets in the act of losing their atmospheres provides rare opportunities to probe their evolution history. Such analysis has been enabled by observations of the helium triplet at 10833 \AA, but past studies have focused on the narrow time window right around the planet's optical transit. We monitored the hot Jupiter HAT-P-32 b using high-resolution spectroscopy from the Hobby-Eberly Telescope covering the planet's full orbit. We detected helium escaping HAT-P-32 b at a $14\sigma$ significance, with extended leading and trailing tails spanning a projected length over 53 times the planet's radius. These tails are among the largest known structures associated with an exoplanet. We interpret our observations using three-dimensional hydrodynamic simulations, which predict Roche Lobe overflow with extended tails along the planet's orbital path.

We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain.

We derive optimal estimators for the binned two-, three-, and four-point correlators of statistically isotropic tensor fields defined on the sphere, in the presence of arbitrary beams, inpainting, and masking. This is a conceptually straightforward extension of the associated scalar field Philcox (2023), but upgraded to include spin-$2$ fields such as Cosmic Microwave Background polarization and galaxy shear, and parity-violating physics in all correlators. All estimators can be realized using spin-weighted spherical harmonic transforms and Monte Carlo summation and are are implemented in the public code PolyBin, with computation scaling, at most, with the total number of bins. We perform a suite of validation tests verifying that the estimators are unbiased and, in limiting regimes, minimum variance. These facilitate general binned analyses of higher-point functions, and allow constraints to be placed on various pheomena, such as non-separable inflationary physics (novelly including polarized trispectra), non-linear evolution in the late Universe, and cosmic parity-violation.

We present a study of the connection between the escape fraction of Lyman-alpha (Ly$\alpha$) and Lyman-continuum (LyC) photons within a sample of N=152 star-forming galaxies selected from the VANDELS survey at $3.85<z_{spec}<4.95$. By combining measurements of H$\alpha$ equivalent width $(W_\lambda(\rm{H\alpha}))$ derived from broad-band photometry with Ly$\alpha$ equivalent width $(W_\lambda(Ly\alpha))$ measurements from VANDELS spectra, we individually estimate $f_{\rm{esc}}^{Ly\alpha}$ for our full sample. In agreement with previous studies, we find a positive correlation between $W_\lambda(Ly\alpha)$ and $f_{\rm{esc}}^{Ly\alpha}$, increasing from $f_{\rm{esc}}^{Ly\alpha}\simeq0.04$ at $W_\lambda(Ly\alpha)=10$\r{A} to $f_{\rm{esc}}^{Ly\alpha}\simeq0.1$ at $W_\lambda(Ly\alpha)=25$\r{A}. For the first time at $z\sim4-5$, we investigate the relationship between $f_{\rm{esc}}^{Ly\alpha}$ and $f_{\rm{esc}}^{\rm{LyC}}$ using $f_{\rm{esc}}^{\rm{LyC}}$ estimates derived using the equivalent widths of low-ionization, FUV absorption lines in composite VANDELS spectra. Our results indicate that $f_{\rm{esc}}^{\rm{LyC}}$ rises monotonically with $f_{\rm{esc}}^{Ly\alpha}$, following a relation of the form $f_{\rm{esc}}^{\rm{LyC}}\simeq 0.15 f_{\rm{esc}}^{Ly\alpha}$. Based on composite spectra of sub-samples with roughly constant $W_\lambda(Ly\alpha)$, but very different $f_{\rm{esc}}^{Ly\alpha}$, we show that the $f_{\rm{esc}}^{\rm{LyC}}-f_{\rm{esc}}^{Ly\alpha}$ correlation is not driven by a secondary correlation between $f_{\rm{esc}}^{Ly\alpha}$and $W_\lambda(Ly\alpha)$. The $f_{\rm{esc}}^{\rm{LyC}}-f_{\rm{esc}}^{Ly\alpha}$ correlation is in good qualitative agreement with theoretical predictions and provides further evidence that estimates of $f_{\rm{esc}}^{\rm{LyC}}$ within the Epoch of Reionization should be based on proxies sensitive to neutral gas density/geometry and dust attenuation.

We present a polarization variability analysis of a sample of 26 $\gamma$-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio quasars (FSRQs). Moreover, BL Lacs display a distribution of their polarization angle typically oriented in a preferential direction, contrary to the rather random distribution of FSRQs. For the latter blazar type, as well as those sources showing a bright stellar emission, we take into account the depolarizing effect introduced by the broad line region and the host galaxy on the measured polarization degree. In this sample we also observe that BL Lacs present an uncorrelated evolution of the flux and the polarization. Contrary, FSRQs show a correlation before the depolarization correction, that is lost however after considering this effect. In addition, we study the behaviour of the polarization angle, searching for angle rotations in its long-term evolution. We derive that the FSRQs studied here show rotations more frequently than BL Lac objects by a factor $\sim$1.5. During these periods we also observe a systematic decrease of the polarization fraction, as well as a marginal flux increase, not significant however to connect rotations with optical flares. We interpret these results within the extended shock-in-jet scenario, able to explain the overall features observed here for the polarization of the blazar sample.

Despite being thought to play an important role in galaxy evolution, the true impact of outflows driven by active galactic nuclei (AGN) on their host galaxies is unclear. In part, this may be because electron densities of outflowing gas are often underestimated: recent studies that use alternative diagnostics have measured much higher densities than those from commonly used techniques, and consequently find modest outflow masses and kinetic powers. Furthermore, outflow ionisation mechanisms - which are often used to probe acceleration mechanisms - are also uncertain. To address these issues, we have analysed archival HST/STIS spectra of the inner regions (r<160pc) of the nearby prototypical Seyfert galaxies NGC 1068 and NGC 4151, which show evidence of warm-ionised outflows driven by the central AGN. We derive high electron densities ($10^{3.6}$<$n_e$<$10^{4.8}$cm$^{-3}$) using the transauroral [OII] and [SII] emission lines ratios for the first time with spatially-resolved observations. Moreover, we find evidence that the gas along the radio axis in NGC 1068 has a significant AGN-photoionised matter-bounded component, and there is evidence for shock-ionisation and/or radiation-bounded AGN-photoionisation along the radio axis in NGC 4151. We also note that the outflow extents are similar to those of the radio structures, consistent with acceleration by jet-induced shocks. Taken together, our investigation demonstrates the diversity of physical and ionisation conditions in the narrow line regions of Seyfert galaxies, and hence reinforces the need for robust diagnostics of outflowing gas densities and ionisation mechanisms.

In this study we present a novel Monte-Carlo code, referred to as GNC, which enables the investigation of dynamical relaxation in clusters comprising multiple mass components in the vicinity of supermassive black holes at the centers of galaxies. Our method is based on two-dimensional Fokker-Planck equations in the energy and angular momentum space, and allows the evolution of multiple mass components, including stars and compact objects. The code demonstrates remarkable flexibility to incorporate additional complex dynamics, such as resonant relaxations and gravitational wave orbital decay. By employing a weighting method, we effectively enhance the statistical accuracy of rare particle results. In this initial publication, we present the fundamental version of our method, focusing on two-body relaxations and loss cone effects. Through comparisons with previous studies, we establish consistent outcomes in terms of relaxation processes, energy and angular momentum distributions, density profiles, and loss cone consumption rates. We consistently observe the development of tangential anisotropy within the cluster, while the outer regions tend to retain near-isotropic characteristics. Moving forward, GNC holds great promise for exploring a wide range of intriguing phenomena within galactic nuclei, in particular relativistic stellar dynamics, providing detailed and insightful outcomes.

A future Galactic Supernova (SN) explosion can lead to a gamma-ray signal induced by ultralight Axion-Like Particles (ALPs) thermally produced in the SN core and converted into high-energy photons in the Galactic magnetic field. The detection of such a signal is in the reach of the Large Area Telescope aboard the \emph{Fermi} Gamma-Ray Space Telescope. The observation of gamma-ray emission from a future SN has a sensitivity to $g_{a\gamma}\gtrsim 4\times 10^{-13}$ GeV$^{-1}$ for a SN at fiducial distance of $10$ kpc and would allow us to reconstruct the ALP-photon coupling within a factor of $\sim2$, mainly due to the uncertainties on the modeling of the Galactic magnetic field.

The widely used Milky Way dust reddening map, the Schlegel, Finkbeiner, & Davis (1998, SFD) map, was found to contain extragalactic large-scale structure (LSS) imprints (Chiang & M\'enard 2019). Such contamination is inherent in maps based on infrared emission, which pick up not only Galactic dust but also the cosmic infrared background (CIB). When SFD is used for extinction correction, over-correction occurs in a spatially correlated and redshift-dependent manner, which could impact precision cosmology using galaxy clustering, lensing, and supernova Ia distances. Similarly, LSS imprints in other Galactic templates can affect intensity mapping and cosmic microwave background experiments. This paper presents a generic way to remove LSS traces in Galactic maps and applies it to SFD. First, we measure descriptive summary statistics of the CIB in SFD by cross-correlating the map with spectroscopic galaxies and quasars in SDSS tomographically as functions of redshift and angular scale. To reconstruct the LSS on the map level, however, additional information on the phases is needed. We build a large set of 180 overcomplete, full-sky basis template maps from the density fields of over 450 million galaxies in WISE and find a linear combination that reproduces all the high-dimensional tomographic two-point statistics of the CIB in SFD. After subtracting this reconstructed LSS/CIB field, the end product is a full-sky Galactic dust reddening map that supersedes SFD, carrying all Galactic features therein, with maximally suppressed CIB. We release this new dust map dubbed CSFD, the Corrected SFD, to the community at https://idv.sinica.edu.tw/ykchiang/CSFD.html.

We present Balmer decrement H$\alpha$/ H$\beta$ measurements for a sample of 51 galaxies at redshifts z = 4-7 observed with the JWST/NIRSpec MSA, as part of the JADES survey. Leveraging 28-hour long exposures and the efficiency of the prism/clear configuration (but also using information from the medium-resolution gratings), we are able to probe directly the low-mass end of the galaxy population, reaching stellar masses Mstar as low as 10^7 Msun . We find that the correlation between Balmer decrement and Mstar is already established at these high redshifts, indicating a rapid build up of dust in moderately massive galaxies at such early epochs. The lowest-mass galaxies in our sample (Mstar = 1-3 x 10^7 Msun ) display a remarkably low Balmer decrement of 2.88 $\pm$ 0.08, consistent with case B, suggesting very little dust content. However, we warn that such a low observed Balmer decrement may also partly be a consequence of an intrinsically lower H$\alpha$/ H$\beta$, resulting from the extreme conditions of the ionized gas in these primeval and unevolved systems. We further compare the Balmer decrement to continuum-derived star-formation rates (SFR), finding tentative evidence of a correlation, which likely traces the underlying connection between SFR and mass of cold gas. However, we note that larger samples are required to distinguish between direct and primary correlations from indirect and secondary dependencies at such high redshifts.

The Magellanic Stream is unique to sample the MW potential from ~50 kpc to 300 kpc, and is also unique in constraining the LMC mass, an increasingly important question for the Local Group/Milky Way modeling. Here we compare strengths and weaknesses of the two types of models (tidal and ram-pressure) of the Magellanic Stream formation. I will present our modeling for the formation of the Magellanic System, including those of the most recent discoveries in the Stream, in the Bridge and at the outskirts of Magellanic Clouds. This model has been successful in predicting most recent observations in both properties of stellar and gas phase. It appears that it is an over-constrained model and provides a good path to investigate the Stream properties. In particular, this model requires a LMC mass significantly smaller than 110^{11} Msun

We analyze the internal cross-section structure and dynamics of inter-halo cosmic-web filaments that connect to massive high-$z$ galaxies. Our analysis is based on a high-resolution \texttt{AREPO} cosmological simulation zooming in on a volume encompassing three filaments of $\sim 1 \, \Mpc$ length feeding three massive haloes of $\sim 10^{12}\,\msun$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a strong cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \kpc$, within which the gas is in virial equilibrium in the potential well set by an isothermal dark-matter filament. The cross-section of the filament can be divided into three radial zones. In the outer ``thermal'' (\textbf{T}) zone, $r \geq 0.65 \, r_{\rm shock}$, the inward gravity and ram-pressure forces from the inflowing gas are over-balanced by outward thermal pressure forces, causing the inflowing gas to decelerate and the shock to expand outwards. In the intermediate ``vortex'' (\textbf{V}) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outward force is provided mainly by centrifugal forces associated with these vortices and, to a lesser extent, by global rotation and thermal pressure, while shear and turbulent forces that are also associated with the vortex structure add to the force inward. The inner ``stream'' (\textbf{S}) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core defining the cold streams that feed the galaxies, with $T\sim 3 \times 10^4 \, \Kdegree$ and $n_{\rm H}\sim 0.01 \cmc$. The core is formed by an isobaric cooling flow that results from post-shock cooling on a free-fall time, and is associated with a decrease in the outward forces, though the gas within the core exhibits both inflows and outflows. [abridged]

Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto-neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process ($r$-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical (GRMHD) simulations of PNS winds which include the effects of a dynamically strong ($B \gtrsim 10^{15}$ G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a ``closed'' zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of $\sim$60 ms, consistent with the predictions of \citet{Thompson03}. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient to enable a successful 2nd-peak $r$-process in at least a modest $\gtrsim 1\%$ of the equatorial wind ejecta.

Early-type B stars with strong magnetic fields and rapid rotation form centrifugal magnetospheres (CMs), as the relatively weak stellar wind becomes magnetically confined and centrifugally supported above the Kepler co-rotation radius. CM plasma is concentrated at and above the Kepler co-rotation radius at the intersection between the rotation and magnetic field axis. Stellar rotation can cause these clouds of material to intersect the viewer's line-of-sight, leading to photometric eclipses. However, for stars with strong ($\sim 10\,{\rm kG}$) magnetic fields and rapid rotation, CMs can become optically thick enough for emission to occur via electron scattering. Using high-precision space photometry from a sample of stars with strong H$\alpha$ emission, we apply simulated light curves from the Rigidly Rotating Magnetosphere model to directly infer magnetic and rotational properties of these stars. By comparing the values inferred from photometric modelling to those independently determined by spectropolarimetry, we find that magnetic obliquity angle $\beta$, viewer inclination $i$ and critical rotation fraction $W$ can be approximately recovered for 3 of the 4 stars studied here. However, there are large discrepancies between the optical depth at the Kepler radius $\tau_{\rm K}$ expected from magnetometry, and the values required to match the observations. We show that $\tau_{\rm K}$ of order unity is needed to reasonably match the light curve morphology of our sample stars.

Magnetic fields on small scales are ubiquitous in the universe. Though they can often be observed in detail, their generation mechanisms are not fully understood. One possibility is the so-called small-scale dynamo (SSD). Prevailing numerical evidence, however, appears to indicate that an SSD is unlikely to exist at very low magnetic Prandtl numbers ($Pr_M$) such as are present in the Sun and other cool stars. We have performed high-resolution simulations of isothermal forced turbulence employing the lowest $Pr_M$ values so far achieved. Contrary to earlier findings, the SSD turns out to be not only possible for $Pr_M$ down to 0.0031, but even becomes increasingly easier to excite for $Pr_M$ below $\simeq\,$0.05. We relate this behaviour to the known hydrodynamic phenomenon referred to as the bottleneck effect. Extrapolating our results to solar values of $Pr_M$ indicates that an SSD would be possible under such conditions.

We introduce the Virgo Consortium's FLAMINGO suite of hydrodynamical simulations for cosmology and galaxy cluster physics. To ensure the simulations are sufficiently realistic for studies of large-scale structure, the subgrid prescriptions for stellar and AGN feedback are calibrated to the observed low-redshift galaxy stellar mass function and cluster gas fractions. The calibration is performed using machine learning, separately for three resolutions. This approach enables specification of the model by the observables to which they are calibrated. The calibration accounts for a number of potential observational biases and for random errors in the observed stellar masses. The two most demanding simulations have box sizes of 1.0 and 2.8 Gpc and baryonic particle masses of $1\times10^8$ and $1\times10^9 \text{M}_\odot$, respectively. For the latter resolution the suite includes 12 model variations in a 1 Gpc box. There are 8 variations at fixed cosmology, including shifts in the stellar mass function and/or the cluster gas fractions to which we calibrate, and two alternative implementations of AGN feedback (thermal or jets). The remaining 4 variations use the unmodified calibration data but different cosmologies, including different neutrino masses. The 2.8 Gpc simulation follows $3\times10^{11}$ particles, making it the largest ever hydrodynamical simulation run to $z=0$. Lightcone output is produced on-the-fly for up to 8 different observers. We investigate numerical convergence, show that the simulations reproduce the calibration data, and compare with a number of galaxy, cluster, and large-scale structure observations, finding very good agreement with the data for converged predictions. Finally, by comparing hydrodynamical and `dark-matter-only' simulations, we confirm that baryonic effects can suppress the halo mass function and the matter power spectrum by up to $\approx20$ per cent.

We investigate the diversity in the sizes and average surface densities of the neutral atomic hydrogen (HI) gas discs in ~280 nearby galaxies detected by the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). We combine the uniformly observed, interferometric HI data from pilot observations of the Hydra cluster and NGC 4636 group fields with photometry measured from ultraviolet, optical and near-infrared imaging surveys to investigate the interplay between stellar structure, star formation and HI structural parameters. We quantify the HI structure by the size of the HI relative to the optical disc and the average HI surface density measured using effective and isodensity radii. For galaxies resolved by >1.3 beams, we find that galaxies with higher stellar masses and stellar surface densities tend to have less extended HI discs and lower HI surface densities: the isodensity HI structural parameters show a weak negative dependence on stellar mass and stellar mass surface density. These trends strengthen when we limit our sample to galaxies resolved by >2 beams. We find that galaxies with higher HI surface densities and more extended HI discs tend to be more star forming: the isodensity HI structural parameters have stronger correlations with star formation. Normalising the HI disc size by the optical effective radius (instead of the isophotal radius) produces positive correlations with stellar masses and stellar surface densities and removes the correlations with star formation. This is due to the effective and isodensity HI radii increasing with mass at similar rates while, in the optical, the effective radius increases slower than the isophotal radius. Our results demonstrate that with WALLABY we can begin to bridge the gap between small galaxy samples with high spatial resolution HI data and large, statistical studies using spatially unresolved, single-dish data.

Dual quasars at small physical separations are an important precursor phase of galaxy mergers, ultimately leading to the coalescence of the two supermassive black holes. Starting from a sample of dual/lensed quasar candidates discovered using astrometric jitter in Gaia data, we present a pilot case study of one of the most promising yet puzzling candidate dual quasars at cosmic noon (z$\sim$1.8). Using multi-wavelength imaging and spectroscopy from X-ray to radio, we test whether the J0823+2418 system is two individual quasars in a bound pair at separation$\sim$0.64'', or instead a single quasar being gravitationally lensed by a foreground galaxy. We find consistent flux ratios ($\sim$1.25-1.45) between the two sources in optical, NIR, UV, and radio, and thus similar spectral energy distributions, suggesting a strong lensing scenario. However, differences in the radio spectral index, as well as changing X-ray fluxes, hint at either a dual quasar with otherwise nearly identical properties, or perhaps lensing-based time lag of $\sim$3 days paired with intrinsic variability. We find with lens mass modeling that the relative NIR positions and magnitudes of the two quasars and a marginally detected central galaxy are consistent with strong lensing. Archival SDSS spectra likewise suggest a foreground absorber via Mg II absorption lines. We conclude that J0823+2418 is likely a lensed quasar, and therefore that the VODKA sample contains a population of these lensed systems (perhaps as high as 50%) as well as dual quasars.

A fundamental prediction of relativistic cosmologies is that, due to the expansion of space, observations of the distant cosmos should be time dilated and appear to run slower than events in the local universe. Whilst observations of cosmological supernovae unambiguously display the expected redshift-dependent time dilation, this has not been the case for other distant sources. Here we present the identification of cosmic time dilation in a sample of 190 quasars monitored for over two decades in multiple wavebands by assessing various hypotheses through Bayesian analysis. This detection counters previous claims that observed quasar variability lacked the expected redshift-dependent time dilation. Hence, as well as demonstrating the claim that the lack of the redshift dependence of quasar variability represents a significant challenge to the standard cosmological model, this analysis further indicates that the properties of quasars are consistent with them being truly cosmologically distant sources.

Astronomy is presently experiencing profound growth in the deployment of machine learning to explore large datasets. However, transient quasi-periodic oscillations (QPOs) which appear in power density spectra of many X-ray binary system observations are an intriguing phenomena heretofore not explored with machine learning. In light of this, we propose and experiment with novel methodologies for predicting the presence and properties of QPOs to make the first ever detections and characterizations of QPOs with machine learning models. We base our findings on raw energy spectra and processed features derived from energy spectra using an abundance of data from the NICER and RXTE space telescope archives for two black hole low mass X-ray binary sources, GRS 1915+105 and MAXI J1535-571. We advance these non-traditional methods as a foundation for using machine learning to discover global inter-object generalizations between - and provide unique insights about - energy and timing phenomena to assist with the ongoing challenge of unambiguously understanding the nature and origin of QPOs. Additionally, we have developed a publicly available Python machine learning library, QPOML, to enable further Machine Learning aided investigations into QPOs.

Cygnus X-1, the well-known accreting black hole system, exhibits several observational features hinting at an intricate interplay between the accretion disk, its atmosphere known as the corona and the putative relativistic jet. It has been extensively studied using all available observational methods, including using the newly available technique of sensitive X-ray polarimetry. X-ray polarization characteristics are distinct for coronal and jet emissions. The low X-ray polarization measured below $\sim$100 keV is understood as arising from the corona. In contrast, the high polarization measurements reported above $\sim$400 keV required a separate jet-dominated spectral component, which spectroscopy does not demonstrate conclusively. Here we report precise polarization measurements in the 100-380 keV region made during three different sub-classes of spectral states of the source using the CZTI instrument onboard {\em AstroSat}. A high polarization (23$\pm$4 \%) is found mainly in the Intermediate Hard State of the source, and the energy-resolved measurements smoothly connect the coronal and the jet regimes. When high polarization is observed, the simultaneous spectral data hints at a separate power law component above 100 keV. We examine the possible sources of this energy-dependent high polarization in Cygnus X-1.

Observations of ionized AGN outflows have provided compelling evidence that the radiation field transfers both momentum and energy to the plasma. At parsec scale distances in AGN, energy transfer can dominate, in which case the only force needed to launch an outflow is that from gas pressure. Much closer to the black hole, gravity dominates thermal energy due to insufficient heating by the radiation and the gas is in the so-called `cold' regime. Only magnetic or radiation forces can then lead to outflow, but it is unclear at what temperature and ionization state the radiation force weakens, as these properties depend on the spectral energy distribution (SED). In this work, we survey the parameter space of radiation forces due to spectral lines resulting from blackbody SEDs with varying temperatures in the range $\sim 10^4 - 10^6$~K to identify the radiation temperature at which line-driving begins to lose efficiency. We find that the temperature $\lesssim4\times10^5$~K marks the transition to inefficient line driving. We also self-consistently compute the heating and cooling balance to estimate the gas temperature, so that our parameter survey covers the transition where thermal driving goes from negligible to comparable to line driving. We summarize a large set of hydrodynamical simulations of radial flows to illustrate how the wind properties change during the transition and the dependence of these properties on the assumed SED and governing flow parameters.

Understanding the collisional behavior of dust aggregates is essential in the context of planet formation. It is known that low-velocity collisions of dust aggregates result in bouncing rather than sticking when the filling factor of colliding dust aggregates is higher than a threshold value. However, a large discrepancy between numerical and experimental results on the threshold filling factor was reported so far. In this study, we perform numerical simulations using soft-sphere discrete element methods and demonstrate that the sticking probability decreases with increasing aggregates radius. Our results suggest that the large discrepancy in the threshold filling factor may reflect the difference in the size of dust aggregates in earlier numerical simulations and laboratory experiments.

The eccentricity of a sub-stellar companion is an important tracer of its formation history. Directly imaged companions often present poorly constrained eccentricities. A recently developed prior framework for orbit fitting called ''observable-based priors'' has the advantage of improving biases in derived orbit parameters for objects with minimal phase coverage, which is the case for the majority of directly imaged companions. We use observable-based priors to fit the orbits of 21 exoplanets and brown dwarfs in an effort to obtain the eccentricity distributions with minimized biases. We present the objects' individual posteriors compared to their previously derived distributions, showing in many cases a shift toward lower eccentricities. We analyze the companions' eccentricity distribution at a population level, and compare this to the distributions obtained with the traditional uniform priors. We fit a Beta distribution to our posteriors using observable-based priors, obtaining shape parameters $\alpha = 1.09^{+0.30}_{-0.22}$ and $\beta = 1.42^{+0.33}_{-0.25}$. This represents an approximately flat distribution of eccentricities. The derived $\alpha$ and $\beta$ parameters are consistent with the values obtained using uniform priors, though uniform priors lead to a tail at high eccentricities. We find that separating the population into high and low mass companions yields different distributions depending on the classification of intermediate mass objects. We also determine via simulation that the minimal orbit coverage needed to give meaningful posteriors under the assumptions made for directly imaged planets is $\approx$ 15% of the inferred period of the orbit.

We present a study of the structural properties and metallicity distribution of the nearby peculiar galaxy NGC3077. Using data from our survey of the M81 Group with the Hyper Suprime-Cam on the Subaru Telescope, we construct deep color-magnitude diagrams that are used to probe the old red giant branch population of NGC3077. We map these stars out to and beyond the nominal tidal radius, which allows us to derive the structural properties and stellar content of the peripheral regions. We show that NGC3077 has an extended stellar halo and pronounced ``S-shaped" tidal tails that diverge from the radial profile of the inner region. The average metallicity of the old population in NGC3077 is estimated from individual RGBs to be $\rm{[M/H]}=-0.98 \pm 0.26$, which decreases with the distance from the galaxy center as $\rm{[M/H]}=-0.17$ dex $\rm{R_{h}}^{-1}$. The metallicity of the S-shaped structure is similar to that of the regions lying at $r\sim4\times\rm{R_{h}} (\sim 30$~kpc), indicating that the stellar constituents of the tidal tails have come from the outer envelope of NGC3077. These results suggest that this peculiar galaxy was probably a rather normal dwarf elliptical galaxy before the tidal interaction with M81 and M82. We also examine the evidence in our dataset for the six recently-reported ultra-faint dwarf candidates around NGC3077. We recover a spatial overdensity of sources coinciding with only one of these.

We investigate the magnetic flux budget in a previously developed dynamo model of Solar Cycles 23 and 24. The mean-field 3D MHD model simulates the global dynamo process and the magnetic buoyancy instability that leads to the formation of bipolar magnetic regions (BMR) on the solar surface. The initial perturbations of the instability correspond to the distribution of active regions observed during the solar cycles. The toroidal and poloidal flux budgets are calculated by applying the Stokes theorem. The results are compared with a baseline 2D dynamo model without bipolar magnetic regions and with the synoptic observations from Kitt Peak Observatory and SoHO (Solar and Heliospheric Observatory) and SDO (Solar Dynamics Observatory) space missions. We find that while the regions of the high radial rotational shear at the boundaries of the convection zone are important for maintaining the dynamo process, the toroidal magnetic flux that results in the formation of BMR is generated by the latitudinal differential rotation. The toroidal flux generation reaches the maximum at the poles during the solar minima, which explains the correlation between the polar field strength and the subsequent sunspot maximum. However, we find that the generation rate of this flux strongly depends on the radial magnetic field distribution near the solar poles. Our results suggest that while the surface magnetic activity contributes to the poloidal magnetic flux budget, a significant part of the poloidal flux is generated in the deep convection zone, contrary to the Babcock-Leighton solar cycle scenario.

The data on spectroscopic galaxy clustering collected by the Dark Energy Spectroscopic Instrument (DESI) will allow the significant detection of subtle features in the galaxy two-point correlation in redshift space, beyond the "standard" redshift-space distortions. Here we present an independent assessment of the detectability of the relativistic dipole in the cross-correlation of two populations of galaxies if they would be selected from the Bright Galaxy Survey (BGS) of DESI. We build synthetic galaxy catalogues with the characteristics of the BGS using the light cone of a relativistic $N$-body simulation. Exploring different ways of splitting the populations of galaxies we find that with an unequal split with more bright galaxies than faint galaxies the detectability is significantly boosted, reaching 19 $\sigma$ in the redshift bin $0.2 \lesssim z \lesssim 0.3$ and expected to be even higher at lower redshift. Moreover, we find that the measured dipole agrees very well with the prediction of relativistic effects from linear theory down to separations of $\sim$ 30 Mpc/$h$.

We study the solar flare index (SFI) for the solar cycles 18\,--\,24. We find that SFI has deeper Gnevyshev gap (GG) in its first principal component than other atmospheric parameters. The GG is extremely clear especially in the even cycles. The GG of the SFI appears about a half year later as a drop in the interplanetary magnetic field near the Earth and in the geomagnetic Ap-index. The instantaneous response of the magnetic field to solar flares, however, shows about two to three days after the eruption as a high, sharp peak in the cross-correlation of the SFI and Ap-index and as a lower peak in SFI vs. IMF B cross-correlation. We confirm these rapid responses using superposed-epoch analysis. The most active flare cycles during 1944-2020 are the Cycles 19 and 21. The Cycle 18 has very strong SFI days as much as Cycle 22, but it has least nonzero SFI days in the whole interval. Interestingly Cycle 20 can be compared to the Cycles 23 and 24 in its low flare activity, although it locates between the most active SFI cycles.

Context. The structural scaling relations (SSRs) of galaxies, i.e. the observed correlations between effective radius, effective surface intensity and velocity dispersion, are important tools for understanding how evolution proceeds. Aims. In this paper we aim to demonstrate that the evolution of the SSRs back in time is governed by the combination of the virial theorem (VT) and the relation $L=L'_0 \sigma^{\beta(t)}$, where the parameters $\beta$ and $L'_0$ vary with time and from galaxy to galaxy. Methods. Using the WINGS database for the galaxies at redshift $z=0$ and the Illustris-1 and Illustris-TNG databases of artificial galaxies, for the galaxies up to redshift $z=4$, we analyse the SSRs back in time and, by means of simple algebraic expressions for $L'_0$ and $\beta$ (functions of time and other physical quantities), we derive the expected paths followed by galaxies in the various SSRs toward the distributions observed at $z=0$. Results. The distribution of galaxies in the SSRs is ultimately related to the evolution in luminosity and velocity dispersion that are empirically mirrored by the $L=L'_0 \sigma^{\beta(t)}$ law. Furthermore, the $\beta$ parameter works as a thermometer of the virialization of a galaxy. This parameter can assume either positive or negative values, and its absolute value attains high values when the galaxy is close to the virial condition, while it tends to zero when the galaxy is far from it. Conclusions. As the SSRs change with time, the method we are proposing allows us to decipher the temporal evolution of galaxies.

We present results from a radio survey for variable and transient sources on 15-min timescales, using the Australian SKA Pathfinder (ASKAP) pilot surveys. The pilot surveys consist of 505 h of observations conducted at around 1 GHz observing frequency, with a total sky coverage of 1476 deg$^2$. Each observation was tracked for approximately 8-10h, with a typical rms sensitivity of $\sim$30 $\mu$jy/beam and an angular resolution of $\sim$12 arcsec. The variability search was conducted within each 8-10h observation on a 15-min timescale. We detected 38 variable and transient sources. Seven of them are known pulsars, including an eclipsing millisecond pulsar, PSR J2039$-$5617. Another eight sources are stars, only one of which has been previously identified as a radio star. For the remaining 23 objects, 22 are associated with active galactic nuclei or galaxies (including the five intra-hour variables that have been reported previously), and their variations are caused by discrete, local plasma screens. The remaining source has no multi-wavelength counterparts and is therefore yet to be identified. This is the first large-scale radio survey for variables and transient sources on minute timescales at a sub-mJy sensitivity level. We expect to discover $\sim$1 highly variable source per day using the same technique on the full ASKAP surveys.

The experimental determination of the reaction rate coefficients for production and destruction of $\text{HCN}^+$ and $\text{HNC}^+$ in collisions with $\text{H}_2$ is presented. A variable temperature 22 pole radio frequency ion trap was used to study the reactions in the temperature range of $17 - 250\;\text{K}$. The obtained rate coefficients for the reaction of $\text{CN}^+$ and of $\text{HCN}^+$ with $\text{H}_2$ are close to the collisional (Langevin) value, whereas that for the reaction of $\text{HNC}^+$ with $\text{H}_2$ is quickly decreasing with increasing temperature. The product branching ratios for the reaction of $\text{CN}^+$ with $\text{H}_2$ are also reported and show a notable decrease of $\text{HNC}^+$ product with respect to $\text{HCN}^+$ product with increasing temperature. These measurements have consequences for current astrochemical models of cyanide chemistry, in particular for the $\text{HCNH}^+$ cation.

A large sample of long-period giant planets has been discovered thanks to long-term radial velocity surveys, but only a few dozen of these planets have a precise radius measurement. Transiting gas giants are crucial targets for the study of atmospheric composition across a wide range of equilibrium temperatures and for shedding light on the formation and evolution of planetary systems. Indeed, compared to hot Jupiters, the atmospheric properties and orbital parameters of cooler gas giants are unaltered by intense stellar irradiation and tidal effects. We identify long-period planets in the Transiting Exoplanet Survey Satellite (TESS) data as duo-transit events. To solve the orbital periods of TESS duo-transit candidates, we use the CHaracterising ExOPlanet Satellite (CHEOPS) to observe the highest-probability period aliases in order to discard or confirm a transit event at a given period. We also collect spectroscopic observations with CORALIE and HARPS in order to confirm the planetary nature and measure the mass of the candidates. We report the discovery of a warm transiting Neptune-mass planet orbiting TOI-5678. After four non-detections corresponding to possible periods, CHEOPS detected a transit event matching a unique period alias. Joint modeling reveals that TOI-5678 hosts a 47.73 day period planet. TOI-5678 b has a mass of 20 (+-4) Me and a radius of 4.91 (+-0.08 Re) . Using interior structure modeling, we find that TOI-5678 b is composed of a low-mass core surrounded by a large H/He layer with a mass of 3.2 (+1.7, -1.3) Me. TOI-5678 b is part of a growing sample of well-characterized transiting gas giants receiving moderate amounts of stellar insolation (11 Se). Precise density measurement gives us insight into their interior composition, and the objects orbiting bright stars are suitable targets to study the atmospheric composition of cooler gas giants.

We investigate the impact of binary orbital motions on the dynamical modeling of dwarf galaxies with intrinsic line-of-sight velocity dispersions ($\sigma_{v_r}$) of 1 to 9 km/s. Using dwarf galaxies from the Auriga level-2 and level-3 simulations, we apply the Jeans Anisotropic Multi-Gaussian Expansion modelling to tracer stars before and after including binaries to recover the dynamical masses. The recovered total masses within the half-mass radius of tracers, $M(<r_\mathrm{half})$, are always inflated due to binary motions, with greater inflations occurring for smaller $\sigma_{v_r}$. However, many dwarf galaxies experience central density deflated due to binary motions, with little dependences on $\sigma_{v_r}$. This is due to the negative radial gradients in the velocity dispersion profiles, with the fractional inflation in $\sigma_{v_r}$ due to binaries more significant in outskirts. An extreme binary fraction of 70% can lead to central density deflation of up to 10-20% at 3 km/s$<\sigma_{v_r}<$8 km/s, with $M(<r_\mathrm{half})$ inflated by 4% at 9 km/s and up to 15% at 3 km/s. A lower binary fraction of 36% leads to similar deflations, with the inflations decreasing to approximately 10% at 3 km/s and becoming statistically insignificant. The choice of binary orbit distribution models does not result in significant differences, and observational errors tend to slightly weaken the deflations in the recovered central density. Two observations separated by one year to exclude binaries lead to almost zero inflations/deflations for a binary fraction of 36% over 3 km/s$<\sigma_{v_r}<$9 km/s. For $\sigma_{v_r}\sim$1 km/s to 3 km/s, a binary fraction of 70% (36%) still results in 60% (30%) to 10% (1%) of inflations in $M(<r_\mathrm{half})$, even with two-epoch observation.

One key question in tidal disruption events theory is that how much of the fallback debris can be accreted to the black hole. Based on radiative hydrodynamic simulations, we study this issue for efficiently `circularized' debris accretion flow. We find that for a black hole disrupting a solar type star, $15\%$ of the debris can be accreted for a $10^7$ solar mass ($M_\odot$) black hole. While for a $10^6M_\odot$ black hole, the value is $43\%$. We find that wind can be launched in the super-Eddington accretion phase regardless of the black hole mass. The maximum velocity of wind can reach $0.7c$ (with $c$ being speed of light). The kinetic power of wind is well above $10^{44} {\rm erg \ s^{-1}}$. The results can be used to study the interaction of wind and the circumnuclear medium around quiescent super-massive black holes.

We construct an efficient parameterization of the pure neutron-matter equation of state (EoS) that incorporates the uncertainties from both chiral effective field theory ($\chi$EFT) and phenomenological potential calculations. This parameterization yields a family of EoSs including and extending the forms based purely on these two calculations. In combination with an agnostic inner core EoS, this parameterization is used in a Bayesian inference pipeline to obtain constraints on the e os parameters using multi-messenger observations of neutron stars. We specifically considered observations of the massive pulsar J0740+6620, the binary neutron star coalescence GW170817, and the NICER pulsar J0030+0451. Constraints on neutron star mass-radius relations are obtained and compared. The Bayes factors for the different EoS models are also computed. While current constraints do not reveal any significant preference among these models, the framework developed here may enable future observations with more sensitive detectors to discriminate them.

It is well known that metal-poor red giant branch (RGB) stars show variations in some elemental abundances, including carbon, due to the internal mixing accompanied by their own in situ CN cycle in the hydrogen burning shell. With our new photometric carbon abundance measurements of RGB stars in M22 and other globular clusters (GCs) in our previous studies, M5, M3, and M92, we derive the carbon depletion rates against the $V$ magnitude, $d\mathrm{[C/Fe]}/M_V$, for individual populations in each GC. We find the metallicity dependence of the carbon depletion rates, $d\mathrm{[C/Fe]}/M_V$ $\propto$ $-$0.25[Fe/H]. Our results also suggest that the carbon depletion rates of the second generation (SG) of stars are larger than those of the first generation (FG) of stars in our sample GCs, most likely due to different internal temperature profiles with different initial helium abundances between the FG and SG. Our results can provide critical constraints both on understanding the mixing efficiency in the theoretical models, which is largely unknown, and on interpretation of the observational carbon abundance evolution of the bright halo RGB stars.

We present the results of a high-cadence and high-precision radial velocity (RV) monitoring of 3 late-type dwarf stars hosting long-period giants with well-measured orbits, in order to search for short-period sub-Neptunes (SN, $M \sin i < 30$ M$_\oplus$). Building on the results and expertise of our previous studies, we carry out combined fits of our HARPS-N data with literature RVs, using MCMC analyses and Gaussian Process regression. We then use the results of our survey to estimate the frequency of sub-Neptunes in systems hosting cold-Jupiters, $f(SN|CJ)$, and compare it with the frequency around field M-dwarfs, $f(SN)$. We identify a new short-period low-mass planet orbiting GJ 328, GJ 328\,c, with $P_c = 241.8^{+1.3}_{-1.7}$ d and $M_c \sin i = 21.4^{+ 3.4}_{- 3.2}$ M$_\oplus$. We moreover identify and model the chromospheric activity signals and rotation periods of GJ 649 and GJ 849, around which no additional planet is found. Then, taking into account also planetary system around the previosuly-analyzed low-mass star BD-11 4672, we derive an estimate of the frequencies of inner planets in such systems. In particular $f(SN|CJ) = 0.25^{+0.58}_{-0.07}$ for mini-Neptunes ($10$ M$_\oplus < M \sin i < 30$ M$_\oplus$, $P < 150$ d), marginally larger than $f(SN)$. For lower-mass planets ($M \sin i < 10$ M$_\oplus$) instead $f(SN|CJ) <0.69$, compatible with $f(SN)$. In light of the newly detected mini-Neptune, we find tentative evidence of a positive correlation between the presence of those planets and that of inner low-mass planets, $f(SN|CJ) > f(SN)$. This might indicate that cold Jupiters have an opposite influence in the formation of inner sub-Neptunes around late-type dwarfs as opposed to their solar-type counterparts, boosting the formation of mini-Neptunes instead of impeding it.

Our aim is to detect and characterise long-period companions around main sequence stars (spectral types late F to early M). We use the RV method to search for exoplanets around stars. The RV variations are measured with HARPS at the ESO 3.6 metre telescope. The true mass and inclination of our heavier companions are provided by astrometry, for which we use proper motions from Hipparcos and Gaia. Five Jupiter-mass exoplanets are reported to orbit HIP54597, BD-210397 (x2), HD74698, and HD94771 with 8.9 yr, 5.2 yr, 17.4 yr, 9.4 yr, and 5.9 yr orbits, and to have minimum masses of $2.01 \pm 0.03$, $0.7 \pm 0.1$, $2.4^{+1.5}_{-0.2}$, $0.40 \pm 0.06$, and $0.53 \pm 0.03 M_J$ respectively. HD74698 also hosts a highly irradiated Neptune in a 15 day orbit with a minimum mass of $0.07 \pm 0.01 M_J$. The mass of HIP54597 b can maximally increase by 10% - 30%, the minimum mass of HD74698 c is likely equal to its true mass, and BD-210397 c has a mass of $2.66^{+0.63}_{-0.32} M_J$. HD62364 hosts a brown dwarf with a true mass of $18.77^{+0.66}_{-0.63} M_J$ in an orbit of 14 yr. HD56380B, HD221638B, and HD33473C have minimum masses within the brown dwarf limits, in orbits of 8.9 yr, 16.6 yr, and 50 yr respectively; however, astrometric measurements reveal them to be stellar binaries, with masses of $375.3^{+8.6}_{-8.4}$, $110.0^{+3.9}_{-3.7}$, and $271.0^{+3.9}_{-3.8} M_J$. The orbits of the stellar binaries HD11938 and HD61383 are incomplete. The preliminary result for HD61383 is a 0.190 $M_{\odot}$ binary in a 39 yr orbit. The secondary of the binary system HD11938 has a mass of 0.33 $M_{\odot}$ - which is confirmed by a secondary peak in the CCF - and a preliminary period of 35 yr. The origin of the 3.0 yr RV signal of HD3964 is uncertain as it shows entanglement with the magnetic cycle of the star. We finally report one more star, HD11608, with a magnetic cycle that mimics a planetary signal.

HIP 9618 (HD 12572, TOI-1471, TIC 306263608) is a bright ($G=9.0$ mag) solar analogue. TESS photometry revealed the star to have two candidate planets with radii of $3.9 \pm 0.044$ $R_\oplus$ (HIP 9618 b) and $3.343 \pm 0.039$ $R_\oplus$ (HIP 9618 c). While the 20.77291 day period of HIP 9618 b was measured unambiguously, HIP 9618 c showed only two transits separated by a 680-day gap in the time series, leaving many possibilities for the period. To solve this issue, CHEOPS performed targeted photometry of period aliases to attempt to recover the true period of planet c, and successfully determined the true period to be 52.56349 d. High-resolution spectroscopy with HARPS-N, SOPHIE and CAFE revealed a mass of $10.0 \pm 3.1 M_\oplus$ for HIP 9618 b, which, according to our interior structure models, corresponds to a $6.8\pm1.4\%$ gas fraction. HIP 9618 c appears to have a lower mass than HIP 9618 b, with a 3-sigma upper limit of $< 18M_\oplus$. Follow-up and archival RV measurements also reveal a clear long-term trend which, when combined with imaging and astrometric information, reveal a low-mass companion ($0.08^{+0.12}_{-0.05} M_\odot$) orbiting at $26^{+19}_{-11}$ au. This detection makes HIP 9618 one of only five bright ($K<8$ mag) transiting multi-planet systems known to host a planet with $P>50$ d, opening the door for the atmospheric characterisation of warm ($T_{\rm eq}<750$ K) sub-Neptunes.

Multi-planet systems are important sources of information regarding the evolution of planets. However, the long-period planets in these systems often escape detection. HD 22946 is a bright star around which 3 transiting planets were identified via TESS photometry, but the true orbital period of the outermost planet d was unknown until now. We aim to use CHEOPS to uncover the true orbital period of HD 22946d and to refine the orbital and planetary properties of the system, especially the radii of the planets. We used the available TESS photometry of HD 22946 and observed several transits of the planets b, c, and d using CHEOPS. We identified 2 transits of planet d in the TESS photometry, calculated the most probable period aliases based on these data, and then scheduled CHEOPS observations. The photometric data were supplemented with ESPRESSO radial velocity data. Finally, a combined model was fitted to the entire dataset. We successfully determined the true orbital period of the planet d to be 47.42489 $\pm$ 0.00011 d, and derived precise radii of the planets in the system, namely 1.362 $\pm$ 0.040 R$_\oplus$, 2.328 $\pm$ 0.039 R$_\oplus$, and 2.607 $\pm$ 0.060 R$_\oplus$ for planets b, c, and d, respectively. Due to the low number of radial velocities, we were only able to determine 3$\sigma$ upper limits for these respective planet masses, which are 13.71 M$_\oplus$, 9.72 M$_\oplus$, and 26.57 M$_\oplus$. We estimated that another 48 ESPRESSO radial velocities are needed to measure the predicted masses of all planets in HD 22946. Planet c appears to be a promising target for future atmospheric characterisation. We can also conclude that planet d, as a warm sub-Neptune, is very interesting because there are only a few similar confirmed exoplanets to date. Such objects are worth investigating in the near future, for example in terms of their composition and internal structure.

The ExoMol database currently provides comprehensive line lists for modelling the spectroscopic properties of molecules in hot atmospheres. Extending the spectral range of the data provided to ultraviolet (UV) wavelengths brings into play three processes not currently accounted for in the ExoMol data structure, namely photodissociation, which is an important chemical process in its own right,the opacity contribution due to continuum absorption and predissociation which can lead to significant and observable line broadening effects. Data structures are proposed which will allow these processes to be correctly captured and the (strong) temperature-dependent effects predicted for UV molecular photoabsorption in general and photodissociation in particular to be represented.

Asteroseismic modelling will be part of the pipeline of the PLATO mission and will play a key role in the mission precision requirements on stellar mass, radius and age. It is therefore crucial to compare how current modelling strategies perform, and discuss the limitations and remaining challenges for PLATO, such as the so-called surface effects, the choice of physical ingredients, and stellar activity. In this context, we carried out a systematic study of the impact of surface effects on the estimation of stellar parameters. In this work, we demonstrated how combining a mean density inversion with a fit of frequencies separation ratios can efficiently damp the surface effects and achieve precise and accurate stellar parameters for ten Kepler LEGACY targets, well within the PLATO mission requirements. We applied and compared two modelling approaches, directly fitting the individual frequencies, or coupling a mean density inversion with a fit of the ratios, to six synthetic targets with a patched 3D atmosphere from Sonoi et al. (2015) and ten actual targets from the LEGACY sample. The fit of the individual frequencies is unsurprisingly very sensitive to surface effects and the stellar parameters tend to be biased, which constitutes a fundamental limit to both accuracy and precision. In contrast, coupling a mean density inversion and a fit of the ratios efficiently damps the surface effects, and allows us to get both precise and accurate stellar parameters. The average statistical precision of our selection of LEGACY targets with this second strategy is 1.9% for the mass, 0.7% for the radius, and 4.1% for the age, well within the PLATO requirements. Using the mean density in the constraints significantly improves the precision of the mass, radius and age determinations, on average by 20%, 33%, and 16%, respectively.

We report the discovery of two warm sub-Neptunes transiting the bright (G = 9.5 mag) K-dwarf HD 15906 (TOI 461, TIC 4646810). This star was observed by the Transiting Exoplanet Survey Satellite (TESS) in sectors 4 and 31, revealing two small transiting planets. The inner planet, HD 15906 b, was detected with an unambiguous period but the outer planet, HD 15906 c, showed only two transits separated by $\sim$ 734 days, leading to 36 possible values of its period. We performed follow-up observations with the CHaracterising ExOPlanet Satellite (CHEOPS) to confirm the true period of HD 15906 c and improve the radius precision of the two planets. From TESS, CHEOPS and additional ground-based photometry, we find that HD 15906 b has a radius of 2.24 $\pm$ 0.08 R$_\oplus$ and a period of 10.924709 $\pm$ 0.000032 days, whilst HD 15906 c has a radius of 2.93$^{+0.07}_{-0.06}$ R$_\oplus$ and a period of 21.583298$^{+0.000052}_{-0.000055}$ days. Assuming zero bond albedo and full day-night heat redistribution, the inner and outer planet have equilibrium temperatures of 668 $\pm$ 13 K and 532 $\pm$ 10 K, respectively. The HD 15906 system has become one of only six multiplanet systems with two warm ($\lesssim$ 700 K) sub-Neptune sized planets transiting a bright star (G $\leq$ 10 mag). It is an excellent target for detailed characterisation studies to constrain the composition of sub-Neptune planets and test theories of planet formation and evolution.

The Virgo cluster is the nearest massive galaxy cluster and thus a prime target to study astrophysical processes in dense large-scale environments. In the radio band, we can probe the non-thermal components of the inter-stellar medium (ISM), intracluster medium (ICM) and of active galactic nuclei (AGN). With the ViCTORIA (Virgo Cluster multi-Telescope Observations in Radio of Interacting galaxies and AGN) project, we are carrying out multiple wide-field surveys of the Virgo cluster at different frequencies. We aim to investigate the impact of the environment on the evolution of galaxies and the contribution of AGN to the ICM-heating, from the inner cluster regions out to beyond the virial radius. We present a survey of the cluster at 120-168 MHz using LOFAR. We image a 132 deg$^2$ region of the cluster, reaching an order of magnitude greater sensitivity than existing wide-field radio surveys of this field at three times higher spatial resolution compared to other low-frequency observations. We developed a tailored data processing strategy to subtract the bright central radio galaxy M87 from the data. This allowed us to correct for the systematic effects due to ionospheric variation as a function of time and direction. In the final mosaic with a resolution of 9"x5", we reach a median noise level of 140 ${\mu}$Jy/beam inside the virial radius and 280 ${\mu}$Jy/beam for the full area. We detect 112 Virgo member galaxies and 114 background galaxies. In at least 18 cases, the radio morphology of the cluster member galaxies shows clear signs of ram-pressure stripping. This includes three previously unreported candidates. In addition, we reveal for the first time 150 kpc long tails from a previous epoch of AGN activity for NGC 4472 (M 49). While no cluster-scale diffuse radio sources are discovered, we find the presence of an extended radio signature of the W$'$-group.

The technique of frequency phase transfer (FPT), enabled by multiband receivers with shared optical path (SOP), is set to become a true backbone of VLBI operations at frequencies above 22 GHz. The FPT has been successfully implemented at the Korean VLBI Network (KVN), while gaining ever more prominent attention worldwide. Over the next few years, FPT VLBI at 22/43/86 GHz will become feasible at more than ten telescopes in Eurasia and Australia. This development would bring order of magnitude improvements of sensitivity and dynamic range of VLBI imaging at 86 GHz and deliver astrometric measurements with an accuracy of one microsecond of arc. The resulting exceptional discovery potential would strongly impact a number of scientific fields ranging from fundamental cosmology and black hole physics to stellar astrophysics and studies of transient phenomena. It is now the right moment for establishing a Science Working Group and a Technical Working Group for FPT VLBI in order to actively focus and coordinate the relevant activities at all stakeholder institutes and ultimately to realize the first global FPT VLBI instrument operating at 22/43/86 GHz.

New MiZo line lists are presented for the D$_2$H$^+$ and D$_3^+$ isotopologues of H$_3^+$. These line lists plus the existing H$_3^+$ MiZATeP and the Sochi H$_2$D$^+$ line lists are updated using empirical energy levels generated using the MARVEL procedure for H$_3^+$, H$_2$D$^+$ and D$_2$H$^+$, and effective Hamiltonian energies for D$_3^+$ for which there is significantly less laboratory data available. These updates allow accurate frequencies for far infrared lines for these species to be predicted. Assignments of the energy levels of H$_3^+$ and D$_3^+$ are extended using a combination of high accuracy variational calculations and analysis of transition intensities. All line lists are made available via www.exomol.com.

We present the properties of 16 faint Lyman-$\alpha$ emitting galaxies (LAEs) at $z>5.8$ from the JWST Advanced Deep Extragalactic Survey (JADES) spectroscopic data in the Hubble Ultra Deep Field/GOODS-S. These LAEs span a redshift range $z\approx5.8-8.0$ and UV magnitude range $M_{\textrm{UV}} \approx -17$ to $-20.6$, with Ly$\alpha$ equivalent width (EW) in the range $\approx 25-350$ \AA. The detection of other rest-optical emission lines in the spectra of these LAEs enables the determination of accurate systemic redshifts and Ly$\alpha$ velocity offsets, as well as the physical and chemical composition of their stars and interstellar media. These faint LAEs are consistent with metal-poor systems with high ionization parameters, similar to the general galaxy population at $z>6$. We measure an average ionizing photon production efficiency, log($\xi_{\textrm{ion}}$/erg$^{-1}$ Hz) $\approx 25.56$ across our LAEs, which does not evolve strongly with redshift. We report an anti-correlation between Ly$\alpha$ escape fraction and velocity offset from systemic, consistent with model expectations. We further find that the strength and velocity offset of Ly$\alpha$ are not correlated with galaxy spectroscopic properties nor with $\xi_{\textrm{ion}}$. We find a decrease in Ly$\alpha$ escape fractions with redshift, indicative of decreasing sizes of ionized bubbles around LAEs at high redshifts. We use a range of galaxy properties to predict Lyman continuum escape fractions for our LAEs, finding that the ionizing photon output into the intergalactic medium from our LAEs remains roughly constant across the observed UV magnitude and Ly$\alpha$ equivalent width, showing a mild increase with redshift. We derive correlations between the ionizing photon output from LAEs and UV magnitude Ly$\alpha$ strengths and redshift, which can be used to build realistic reionization models.

The recent observations revealed that the cosmic-ray (CR) proton spectrum showed a complex structure: the hardening at $\rm \sim 200\,GeV$ and softening at $\rm \sim 10\,TeV$. However, so far the physical origins of this spectral feature remain strongly debated. In this work, we simulate the acceleration of cosmic-ray protons in a nearby Supernova remnant (SNR) by solving numerically the hydrodynamic equations and the equation for the quasi-isotropic CR momentum distribution in the spherically symmetrical case to derive the spectrum of protons injected into the interstellar medium (ISM), and then simulate the propagation process of those accelerated CR particles to calculate the proton fluxes reaching the Earth. Besides, we use the DRAGON numerical code to calculate the large-scale cosmic-ray proton spectrum. Our simulated results are in good agreement with the observed data (including the observed data of proton fluxes and dipole anisotropy). We conclude that the spectral feature of cosmic-ray protons in this energy band may originate from the superposition of the distribution from the nearby SNR and background diffusive cosmic-ray component. We find that the release of particles from this nearby SNR has a time delay. Besides, it can be found that the nonlinear response of energetic particles, release time of CR protons, and age of the local SNR can leave strong signatures in the spectrum of the resulting CR proton fluxes.

The Near Infrared Imager and Slitless Spectrograph instrument (NIRISS) is the Canadian Space Agency (CSA) contribution to the suite of four science instruments of JWST. As one of the three NIRISS observing modes, the Single Object Slitless Spectroscopy (SOSS) mode is tailor-made to undertake time-series observations of exoplanets to perform transit spectroscopy. The SOSS permits observing point sources between 0.6 and 2.8 um at a resolving power of 650 at 1.25 um using a slit-less cross-dispersing grism while its defocussing cylindrical lens enables observing targets as bright as J=6.7 by spreading light across 23 pixels along the cross-dispersion axis. This paper officially presents the design of the SOSS mode, its operation, characterization, and its performance, from ground-based testing and flight-based Commissioning. On-sky measurements demonstrate a peak photon conversion efficiency of 55% at 1.2 um. The first time-series on the A-type star BD+60o1753 achieves a flux stability close to the photon-noise limit, so far tested to a level of 20 parts per million on 40-minute time-scales after simply subtracting a long-term trend. Uncorrected 1/f noise residuals underneath the spectral traces add an extra source of noise equivalent to doubling the readout noise. Preliminary analysis of a HAT-P-14b transit time-series indicates that it is difficult to remove all the noise in pixels with partially saturated ramps. Overall, the SOSS delivers performance at the level required to tackle key exoplanet science programs such as detecting secondary atmospheres on terrestrial planets and measuring abundances of several chemical species in gas giants.

The recent discovery of objects with redshift $z>10$ with the help of James Webb Space Telescope (JWST) poses serious challenges to the $\Lambda$CDM cosmological model, which has been in vogue for some time now. The new data indicate that galaxy formation must have taken place much earlier than expected in this model. Another viable class of cosmological models is that of the so-called coasting models, in which the scale factor of the universe varies proportionately with time. In these models, the universe at redshift $z=12$ has ample time ($\sim 1070$ Myrs) for galaxy formation. The earliest such model is the one proposed by E.A. Milne, based on his `kinematic relativity', but it is considered unrealistic for not treating gravity as relevant at cosmological scales. A closed version of an eternal coasting FLRW model was proposed by the present authors even before SNe Ia data began to pour in. Subsequently we developed a more general model of the same class, which is valid for all the three possible geometries, with open, closed or flat spatial sections. In the nonrelativistic era, this model makes the falsifiable prediction that the ratio of matter density to dark energy density is 2. This avoids the cosmic coincidence problem. Moreover, this eternal coasting model allows room for creation of matter from dark energy, that may speed up galaxy and structure formation at the early epochs, as implied by the JWST data. The paper also attempts to review some similar coasting models, but emphasizes the eternal coasting cosmology as the most probable candidate model capable of explaining the presence of high redshift galaxies discovered by JWST.

Gravitational waves (GW) emanating from unstable quasi-normal modes in Neutron Stars (NS) could be accessible with the improved sensitivity of the present gravitational wave (GW) detectors or with the next-generation GW detectors and therefore employed to study the NS interior. By taking into account potential GW candidates detectable by A+ and Einstein Telescope (ET) originating from f-modes excited by glitches in isolated pulsars, we demonstrate the inverse problem of NS asteroseismology in a Bayesian formalism to constrain the nuclear parameters within a relativistic mean field (RMF) description of NS interior. We find that for a single detected GW event from the Vela pulsar in A+ and ET, with the considered RMF model, the nucleon effective mass ($m^*$) can be restricted (within $90\%$ credible interval) within $10\%$ and $5\%$, respectively. With the considered RMF model, the incompressibility ($K$) and the slope of the symmetry energy ($L$) are only loosely constrained. With a single observed event in A+ and ET, the f-mode frequency of a $1.4M_{\odot}$ ($f_{1.4M_{\odot}}$) inside a 90\% symmetric credible interval (SCI) can be confined to 100 Hz and 50 Hz, respectively. Additionally, we consider multiple GW candidates in our analysis. For detecting multiple (ten) events with A+ and ET, $m^*$ can be constrained to $3\%$ and $2\%$, respectively. All the other nuclear saturation parameters get well constrained. In particular, $K$ and $L$ can be constrained within $10\%$ and $20\%$ (< $90\%$ SCI), respectively. Within the 90\% SCI, $f_{1.4M_{\odot}}$ can be estimated within 50 Hz and 20 Hz in A+ and ET, respectively. Uncertainty of other NS properties such as radius of a $1.4M_{\odot}$ ($R_{1.4M_{\odot}}$), f-mode damping time of a $1.4M_{\odot}$ ($\tau_{1.4M_{\odot}}$) and few equations of state (EOS) properties including squared speed of sound ($c_s^2$) are also estimated.

We present a study of the environments of 16 Lyman-$\alpha$ emitting galaxies (LAEs) in the reionisation era ($5.8 < z < 8$) identified by $\mathit{JWST}$/NIRSpec as part of the $\mathit{JWST}$ Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Lyman-$\alpha$ emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Lyman-$\alpha$ velocity offsets ($\Delta v_\mathrm{Ly\alpha} < 300 \, \mathrm{km \, s^{-1}}$) with moderately high Lyman-$\alpha$ escape fractions ($f_\mathrm{esc, \, Ly\alpha} > 5\%$) observed in our sample of LAEs, indicating the presence of ionised ``bubbles'' with physical sizes of the order of $0.1 \, \mathrm{pMpc} \lesssim R_\text{ion} \lesssim 1 \, \mathrm{pMpc}$ in a patchy reionisation scenario where the bubbles are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at $z \sim 5.8$-$5.9$ and $z \sim 7.3$, suggesting Lyman-$\alpha$ transmission is strongly enhanced in such overdense regions, and underlining the importance of LAEs as tracers of the first large-scale ionised bubbles. Considering only spectroscopically confirmed galaxies, we find our sample of UV-faint LAEs ($M_\text{UV} \gtrsim -20 \, \mathrm{mag}$) and their direct neighbours are generally not able to produce the required ionised regions based on the Lyman-$\alpha$ transmission properties, suggesting lower-luminosity sources likely play an important role in carving out these bubbles. These observations demonstrate the combined power of $\mathit{JWST}$ multi-object and slitless spectroscopy in acquiring a unique view of the early stages of Cosmic Reionisation via the most distant LAEs.

Large-scale asymmetries (i.e. lopsidedness) are a common feature in the stellar density distribution of nearby disk galaxies both in low- and high-density environments. In this work, we characterize the present-day lopsidedness in a sample of 1435 disk-like galaxies selected from the TNG50 simulation. We find that the percentage of lopsided galaxies (10%-30%) is in good agreement with observations if we use similar radial ranges to the observations. However, the percentage (58%) significantly increases if we extend our measurement to larger radii. We find a mild or lack of correlation between lopsidedness amplitude and environment at z=0 and a strong correlation between lopsidedness and galaxy morphology regardless of the environment. Present-day galaxies with more extended disks, flatter inner galactic regions and lower central stellar mass density (i.e. late-type disk galaxies) are typically more lopsided than galaxies with smaller disks, rounder inner galactic regions and higher central stellar mass density (i.e. early-type disk galaxies). Interestingly, we find that lopsided galaxies have, on average, a very distinct star formation history within the last 10 Gyr, with respect to their symmetric counterparts. Symmetric galaxies have typically assembled at early times (~8-6 Gyr ago) with relatively short and intense bursts of central star formation, while lopsided galaxies have assembled on longer timescales and with milder initial bursts of star formation, continuing building up their mass until z=0. Overall, these results indicate that lopsidedness in present-day disk galaxies is connected to the specific evolutionary histories of the galaxies that shaped their distinct internal properties.

Detecting and coherently characterizing thousands of gravitational-wave signals is a core data-analysis challenge for the Laser Interferometer Space Antenna (LISA). Transient artifacts, or "glitches", with disparate morphologies are expected to be present in the data, potentially affecting the scientific return of the mission. We present the first joint reconstruction of short-lived astrophysical signals and noise artifacts. Our analysis is inspired by glitches observed by the LISA Pathfinder mission, including both acceleration and fast displacement transients. We perform full Bayesian inference using LISA time-delay interferometric data and gravitational waveforms describing mergers of massive black holes. We focus on a representative binary with a detector-frame total mass of $6 \times 10^7 M_\odot$ at redshift $7$, yielding a signal lasting $\sim 30~\mathrm{h}$ in the LISA sensitivity band. We explore two glitch models of different flexibility, namely a fixed parametric family and a shapelet decomposition. In the most challenging scenario, we report a complete loss of the gravitational-wave signal if the glitch is ignored; more modest glitches induce biases on the black-hole parameters. On the other hand, a joint inference approach fully sanitizes the reconstruction of both the astrophysical and the glitch signal. We also inject a variety of glitch morphologies in isolation, without a superimposed gravitational signal, and show we can identify the correct transient model. Our analysis is an important stepping stone toward a realistic treatment of LISA data in the context of the highly sought-after "global fit".

We investigate the inflationary cosmology involving an SU(5) GUT (grand unified theory) singlet scalar with non-minimal coupling to the Ricci scalar. In this scenario the scale of grand unification is set by the inflaton vev when the inflaton rolls down its potential towards its minimum $v$, thereby relating inflationary dynamics to GUT symmetry breaking with a prediction $r \simeq 0.025$ for the tensor-to-scalar ratio to be tested by the next generation CMB experiments. We show in this inflationary framework involving inflection-point how a suitable choice of parameters in $SU(5)$ leads to a bump in the scalar power spectrum with production of Primordial Blackholes (PBH) of masses $10^{17}-10^{18}$g ($ 10 - 100 M_\odot$). We derive the constraints on the self quartic and mixed quartic couplings of the inflaton in SU(5) that are consistent with the inflationary analysis. Moreover, we also show that this scenario leads to large amplitude induced second-order tensor perturbations propagating as Gravitational Waves (GW) with amplitude $\Omega_{\rm GW}h^2 \sim 10^{-17}$ and peak frequency $f_{\rm peak} \sim$ (0.1 - 1) Hz, which can be detected in the next generation GW observatories like LISA, BBO, ET, etc. Thus, we unify the $SU(5)$ framework with PBH via inflection-point inflation showing how the upcoming measurements of PBH and GW will enable us to probe the scale of $SU(5)$ symmetry breaking, and thereby complementing the laboratory based experiments. We also discuss scenarios involving the Pati-Salam and Trinification gauge groups and its impact on quartic and mixed-quartic couplings that may lead to PBH and detectable GW signals.

We demonstrate that a model with extra dimensions formulated in Phys. Rev. D, 62, 045015 , which fatefully reproduces Friedmann-Robertson-Walker (FRW) equations on the brane, allows for an apparent superluminal propagation of massless signals. Namely, a massive brane curves the spacetime and affects the trajectory of a signal in a way that allows a signal sent from the brane through the bulk to arrive (upon returning) to a distant point on the brane faster than the light can propagate along the brane. In particular, the signal sent along the brane suffers a greater gravitational time delay than the bulk signal due to the presence of matter on the brane. While the bulk signal never moves with the speed greater than the speed of light in its own locality, this effect still enables one to send signals faster than light from the brane observer's perspective. For example, this effect might be used to resolve the cosmological horizon problem. In addition, one of the striking observational signatures would be arrival of the same gravitational wave signal at two different times, where the first signals arrives before its electromagnetic counterpart. We used GW170104 gravitational wave event to impose a strong limit on the model with extra dimensions in question.

We present a model of Quarkyonic matter that is explicitly dual between quarks and baryons. The baryon number and energy densities are expressed as functionals of either the baryon momentum distribution, $f_{\rm B}$, or the quark distribution, $f_{\rm Q}$, which are subject to the constraints on fermions, $0 \le f_{\rm B,Q} \le 1$. The theory is ideal in the sense that the confinement of quarks into baryons is reflected in the duality relation between $f_{\rm Q}$ and $f_{\rm B}$, while other possible interactions among quarks and baryons are all neglected. The variational problem with the duality constraints is formulated and we explicitly construct analytic solutions, finding two distinct regimes: A nuclear matter regime at low density and a Quarkyonic regime at high density. In the Quarkyonic regime, baryons underoccupy states at low momenta but form a momentum shell with $f_{\rm B}=1$ on top of a quark Fermi sea. Such a theory describes a rapid transition from a soft nuclear equation of state to a stiff Quarkyonic equation of state. At this transition, there is a rapid increase in the sound speed.

GW170817 - GRB 170817A provided the first observation of gravitational waves from a neutron star merger with associated transient counterparts across the entire electromagnetic spectrum. This discovery demonstrated the long-hypothesized association between short gamma-ray bursts and neutron star mergers. More joint detections are needed to explore the relation between the parameters inferred from the gravitational wave and the properties of the gamma-ray burst signal. We developed a joint multi-messenger analysis of LIGO, Virgo, and Fermi/GBM data designed for detecting weak gravitational-wave transients associated with weak gamma-ray bursts. As such, it does not start from confident (GWTC-1) events only. Instead, we take the full list of existing compact binary coalescence triggers generated with the PyCBC pipeline from the second Gravitational-Wave Observing Run (O2), and reanalyze the entire set of public Fermi/GBM data covering this observing run to generate a corresponding set of gamma-ray burst candidate triggers. We then search for coincidences between the gravitational-wave and gamma-ray burst triggers without requiring a confident detection in any channel. The candidate coincidences are ranked according to a statistic combining each candidate's strength in gravitational-wave and gamma-ray data, their time proximity, and the overlap of their sky localization. The ranking is then converted to a false-alarm rate using time shifts between the gravitational-wave and gamma-ray burst triggers. We present the results using O2 triggers which allowed us to check the validity of our method against GW170817 - GRB 170817A. We also discuss the different configurations tested to maximize the significance of the joint detection.

Gravitational waves (GWs) from stellar-mass compact binary coalescences (CBCs) are expected to be strongly lensed when encountering large agglomerations of matter, such as galaxies or clusters. Searches for strongly lensed GWs have been conducted using data from the first three observing runs of the LIGO-Virgo GW detector network. Although no confirmed detections have been reported, interesting candidate lensed pairs have been identified. In this work, we delineate a preliminary analysis that rapidly identifies pairs to be further analyzed by more sophisticated Bayesian parameter estimation (PE) methods. The analysis relies on the Gaussian/Fisher approximation to the likelihood and compares the corresponding approximate posteriors on the chirp masses of the candidate pair. It additionally cross-correlates the rapidly produced localization sky areas (constructed by Bayestar sky-localization software). The analysis was used to identify pairs involving counterparts from targeted sub-threshold searches to confidently detected super-threshold CBC events. The most significant candidate ``super-sub'' pair deemed by this analysis was subsequently found, by more sophisticated and detailed joint-PE analyses, to be among the more significant candidate pairs, but not sufficiently significant to suggest the observation of a lensed event [1].

Palatini $f(R)$ gravity is probably the simplest extension of general relativity (GR) and the simplest realization of a metric-affine theory. It has the same number of degrees of freedom as GR and, in vacuum, it is straightforwardly mapped into GR with a cosmological constant. The mapping between GR and Palatini $f(R)$ inside matter is possible but at the expense of reinterpreting the meaning of the matter fields. The physical meaning and consequences of such mapping will depend on the physical context. Here we consider three such cases within the weak field limit: Solar System dynamics, planetary internal dynamics (seismology), and galaxies. After revising our previous results on the Solar System and Earth's seismology, we consider here the possibility of $f(R)$ Palatini as a dark matter candidate. For any $f(R)$ that admits a polynomial approximation in the weak field limit, we show here, using SPARC data and a recent method that we proposed, that the theory cannot be used to replace dark matter in galaxies. We also show that the same result applies to the Eddington-inspired Born-Infeld gravity. Differently from the metric $f(R)$ case, the rotation curve data are sufficient for this conclusion. This result does not exclude a combination of modified gravity and dark matter.

Ultra-high frequency gravitational waves in the MHz to THz regime promise a unique possibility to probe the very early universe, particle physics at very high energies and exotic astrophysical objects - but achieving the sensitivity required for detection is an immense challenge. This is a brief summary of recent progress in electromagnetic high-frequency gravitational wave searches, which are based on classical electromagnetism in a space-time perturbed by gravitational waves. A particular focus is given to synergies with axion searches and atomic precision measurements. This article was prepared as proceedings for Moriond EW 2023.