Papers for Friday, May 02 2025

Small satellite technologies have enhanced the potential and feasibility of geodesic missions, through simplification of design and decreased costs allowing for more frequent launches. On-satellite data acquisition systems can benefit from the implementation of machine learning (ML), for better performance and greater efficiency on tasks such as image processing or feature extraction. This work presents convolutional autoencoders for implementation on the payload of small satellites, designed to achieve dual functionality of data compression for more efficient off-satellite transmission, and at-source anomaly detection to inform satellite data-taking. This capability is demonstrated for a use case of disaster monitoring using aerial image datasets of the African continent, offering avenues for both novel ML-based approaches in small satellite applications along with the expansion of space technology and artificial intelligence in Africa.

In this paper we consider the extent to which a lack of observations from SETI may be used to quantify the Fermi paradox. Building on previous research, we construct a geometrical model to compute the probability of at least one detection of an extraterrestrial electromagnetic (EM) signal of galactic origin, as a function of the number $N$ of communicative civilizations. We show how this is derivable from the probability of detecting a single signal; the latter is $\approx 0.6 \delta/R$, where $\delta$ is the distance between the initial and final EM signals and $R$ is the radius of the Milky Way, for $\delta/R \ll 1$. We show how to combine this analysis with the Drake equation $N = \mathscr{N} \delta /c$, where $c$ is the speed of light; this implies, applying a simplified toy model as an example, that the probability of detecting at least one signal is $>99 \%$ for $\delta / c \gtrsim 10^{2.8}$ years, given that $\mathscr{N} = 1$. Lastly, we list this toy model's significant limitations, and suggest ways to ameliorate them in more realistic future models.

Ideal spatial demultiplexing (SPADE) is proven to be a quantum-optimal tool for exoplanet detection, i.e., asymmetric source discrimination. However, recent investigations into the related problems of separation estimation and symmetric source discrimination showed its efficiency to be limited in the presence of noise. In this work, we use analytical tools to scrutinize the practical applicability of SPADE and derive the associated optimal decision strategy for exoplanet detection in the presence of experimental imperfections. On the one hand, we find that the probability of detection of noisy SPADE has the same scaling with planet-star separation and relative brightness as conventional techniques, such as direct imaging and coronagraphs. On the other hand, we prove that, due to a superior scaling coefficient under realistic noise conditions, SPADE remains the most efficient method for practical exoplanet detection in the sub-Rayleigh regime.

Recent observations estimate that 30% of early B and O-type stars are found in triple systems. So far, the evolution of triple star systems has mainly been modeled using fast stellar codes. Their accuracy decreases with increasing mass, limiting their reliability for predicting the evolutionary pathways of massive triple systems. We coupled Tres, which by default uses Seba (fast stellar code) to Mesa to perform the first simulations of triple systems that combine a triple secular evolutionary code with a detailed, on-the-fly stellar code. After examining the differences between the stellar evolution predicted by the two codes, we simulate the evolution of a set of triple systems and compare their predicted evolutionary pathways. The predicted stellar tracks become increasingly divergent with increasing mass and wind mass loss efficiency. The maximal radial extent, crucial for determining whether the components of the triple systems interact, differ by up to two orders of magnitude between the two stellar codes in the considered mass range. This leads to divergences in the triples evolutionary pathways predicted by mesa and seba. Using mesa instead of seba, the minimum period for avoiding inner mass transfer is reduced by three orders of magnitude. This has important consequences for the formation of GW sources through the triple compact object channel. Our simulations offer new insights into the physics of triple systems, as key processes (mass loss, radial expansion, precession) are treated self-consistently. They indicate that the results of triple systems population synthesis studies must be interpreted cautiously, in particular when the considered masses are outside the range of the grid the fast codes are based on and when significant stellar winds are considered.

The existence of billion-solar-mass black holes hosted in luminous quasars within the first gigayear of cosmic history poses a challenge to our understanding of supermassive black hole (SMBH) growth. The problem is further exacerbated by the very short quasar lifetimes of $t_{\rm Q}\lesssim 10^6$ years, as derived from the extent of their proximity zone (PZ) sizes observed in the quasars' rest-UV spectra. However, the quasar lifetime estimates based on the extents of the proximity zones may be underestimated, as time-variable obscuration effects might have limited the quasars' emission along our sightline in the past. In this work, we present independent quasar lifetime measurements for six quasars at $z \sim 6$ leveraging the extended nebular emission perpendicular to our line-of-sight. We use observations from the Very Large Telescope/Multi-Unit Spectroscopic Explorer (MUSE) to search for extended Ly$\alpha$ emission in the circumgalactic medium around quasars with small proximity zones and estimate their lifetimes as the light travel time between the SMBH and the outer edge of the nebula. We find agreement between the independent lifetime estimates. For one object we find a proximate absorption system prematurely truncating the extent of the quasar's proximity zone, which thus results in an expected discrepancy between the lifetime estimates. Our results provide further evidence that the quasars' current accretion episode has only recently begun, challenging our models of SMBH growth.

AT 2019aalc is a peculiar sequence of highly variable emission events observed towards the nucleus of the broad-line AGN SDSS J152416.66+045119.0. The system exhibited two distinct UV/optical flares (the first detected in 2019, the second one in 2023). Spectra obtained following the detection of the second flare revealed prominent Bowen fluorescence and high-ionization coronal emission lines, which were much weaker, if at all detectable, in a spectrum taken following the first flare. We present and analyze a large set of multi-wavelength, multi-epoch data for this source, with particular emphasis on optical spectroscopic monitoring conducted with the Las Cumbres Observatory network. During the relatively slow dimming that followed the second optical flare, the UV/optical light-curve shows a sequence of minor rebrightening events, while the Bowen fluorescence and the coronal lines vary (roughly) in tandem with these "bumps" in the broad-band light-curve. Most of the observed behavior of AT 2019aalc links it to the growing class of Bowen fluorescence flares (BFFs) while setting it apart from canonical tidal disruption events. However, AT 2019aalc has some outstanding peculiarities, including two short flares seen in its soft X-ray light-curve during the dimming phase of the second optical flare, and which do not seem to be linked to the emission line variations. We discuss the optical and X-ray properties of the source and possible scenarios of the origin of the flare, in particular radiation pressure instabilities in the (pre-existing) AGN accretion disk.

Pebble accretion is the leading theory for the formation of exoplanets more massive than the Earth. Many parameters influence planet growth in the pebble accretion models. In this paper, we study the influence of pebble fragmentation velocity, turbulence strength, stellar metallicity, stellar mass, and planet location on the growth of planets located within 1 au from their parent stars. Analysing the close-in planets from NASA's Exoplanet Archive, we find that the turbulence strength influences planet growth more than the pebble fragmentation velocity does. Planets orbiting stars with higher metallicity have an overall higher probability of reaching their pebble isolation mass than those orbiting lower metallicity stars, but the impact of metallicity is not as high as that of stellar mass, orbital separation, and most importantly disc turbulence.

Free-floating planets comprise one of the most enigmatic populations of exoplanets in the Galaxy. Though ground-based observations point to a large abundance of these worlds, little is known about their origins and demographics. In the coming years, the Nancy Grace Roman Space Telescope's Galactic Bulge Time Domain Survey is expected to detect several hundred free-floating planets, providing the first opportunity to characterize these worlds at the population level. We present a first study of Roman's prospects for reconstructing the mass distribution of free-floating planets through population-level statistical inference. We find that depending on the true underlying mass distribution of free-floating planets in the Galaxy, Roman will be able to improve upon existing estimates of the abundance by orders of magnitudes in the largely unexplored mass range below that of Earth. When applied to Roman's observations, the methodology we present herein will be capable of discriminating between different hypothesized mass distributions at high statistical significance, opening a new window into the origins of these rogue worlds.

Weak gravitational lensing (WL) is a key cosmological probe that requires precise measurement of galaxy images to infer shape distortions, or shear, and constrain cosmology. Accurate estimation of the Point Spread Function (PSF) is crucial for shear measurement, but the wavelength dependence of the PSF introduces chromatic biases that can systematically impact shear inference. We focus on biases arising from spectral energy distribution (SED) differences between stars, used for PSF modeling, and galaxies, used for shear measurement. We investigate these effects in $\textit{Roman's}$ four design reference mission WL bands (Y106, J129, H158, F184) and wide filter (W146). Using $\textit{Roman}$-like image simulations, we quantify the induced shear biases and compare them to requirements on those biases. Multiplicative biases over all galaxies hover around $\sim$0.2% in the WL bands and 2% in the wide filter, exceeding the mission requirement of $|m| < 0.032\%$ and relaxed requirement of $|m| < 0.1\%$. In individual redshift bins, biases can reach 0.4$\unicode{x2013}$0.9% for the WL bands and 3$\unicode{x2013}$6% for the wide filter. Additive biases remain acceptable in the WL bands but exceed systematic limits in the wide filter. We develop and test PSF-level corrections, showing that a first-order correction reduces biases within survey requirements for the WL bands; however, higher-order terms are necessary for the wide filter. Our results highlight the necessity of chromatic corrections for precision WL with $\textit{Roman}$ and provide a framework for mitigating these biases. Finally, we compare analytical color-based corrections to self-organizing maps (SOMs) and find that both methods effectively reduce biases.

Metallicity ($Z$) estimates based on ultraviolet (UV) emission lines from the narrow-line regions (NLRs) of active galactic nuclei (AGNs) have been found to differ from those derived from optical lines. However, the origin of this discrepancy ($ZR$) remains poorly understood. To investigate the source of $ZR$, we compiled from the literature the fluxes of narrow near-UV ($1000 < \lambda(\angstrom) < 2000)$ and optical ($3000 < \lambda(\angstrom) < 7000)$ emission line measurements for a sample of 11 AGNs (9 at $z<0.4$ and 2 at $z\sim2.4$). Metallicity values for our sample were derived using a semi-empirical calibration based on the $C43$=log[(\ion{C{iv}$\lambda$1549+\ion{C{iii}]$\lambda$1909)/\ion{He}{ii}$\lambda$1640] emission-line ratio and compared with those obtained via direct measurement of the electron temperature ($T_{\rm e}$-method) and via calibrations based on optical emission-lines. The source of the discrepancy was investigated in terms of the ionization parameter ($U$), electron density ($N_{\rm e}$), and carbon abundance (C/H). We found a weak correlation between $ZR$, $U$ and $N_{\rm e}$. However, a moderate correlation was observed between $ZR$ and direct estimates of C/H, suggesting that the previously assumed (C/O)-$Z$ relations in photoionization models used to derive UV carbon-line calibrations may not be valid for AGNs. By combining a large set of abundance estimates for local star-forming regions with those of our AGN sample, we derived a new (C/O)-$Z$ relation. Comparisons between the results of photoionization models that assume this new abundance relation and the UV observational data of our sample produce $Z$ values derived from the $C43$ index that are consistent with those obtained using the $T_{\rm e}$-method.

We present the results of our identification of 11 X-ray sources detected on the half of the sky $0^\circ<l<180^\circ$ in the 4-12 keV energy band on the combined map of the first five all-sky surveys with the Mikhail Pavlinsky ART-XC telescope onboard the SRG observatory. All these sources were also detected by the SRG/eROSITA telescope in the 0.2-8 keV energy band, whose data have allowed us to improve their positions and to investigate their X-ray spectra. Five of them have been detected in X-rays for the first time, while the remaining ones have already been known previously, but their nature has remained unknown. We have taken optical spectra for nine sources with the 1.6-m AZT-33IK telescope at the Sayan Observatory (the Institute of Solar-Terrestrial Physics, the Siberian Branch of the Russian Academy of Sciences); for two more objects we have analyzed the archival spectra from SDSS and the 6dF survey. The objects are classified as Seyfert galaxies (seven Sy1, three Sy1.9, and one Sy2) at redshifts $z$=0.029-0.258. Our analysis of the X-ray spectra has revealed a noticeable intrinsic absorption ($N_{\rm H} \sim 10^{22}$ cm$^{-2}$) in two of the four Seyfert 2 galaxies (Sy1.9-2). The spectrum of one more of them (SRGA J000132.9+240237) cannot be described within the model of an absorbed Comptonization continuum, which may point to a strong absorption and a significant contribution of the reflected radiation. However, the available SRG all-sky survey data are not enough to obtain reliable constraints on the absorption column density in this object, which is also interesting in that it is radio loud. Longer X-ray observations are required to refine the physical properties of this active galactic nucleus.

We performed astrophysics model calculations with updated nuclear data to identify a possible bypass of the $^{64}{\rm Ge}$ waiting-point, a defining feature of the rapid-proton capture (rp-) process that powers type-I x-ray bursts on accreting neutron stars. We find that the rp-process flow through the $^{64}{\rm Ge}$ bypass could be up to 36\% for astrophysically relevant conditions. Our results call for new studies of $^{65}{\rm Se}$, including the nuclear mass, $\beta$-delayed proton emission branching, and nuclear structure as it pertains to the $^{64}{\rm As}(p,\gamma)$ reaction rate at x-ray burst temperatures.

WR 31a (Hen 3-519) is likely a post-luminous blue variable (LBV) star that is evolving to become a classical Wolf-Rayet star. Multicolor (UBVR) photopolarimetric observations of WR 31a were obtained over nine nights in early 2007. The linear polarization data of WR 31a trace a "loop" structure in a Stokes Q-U diagram, which is similar in all four passbands. After mean subtraction, the four loops align to form a single overall pattern. Such loops can be expected to arise from binary systems. We test the binary hypothesis with two models. The data are fit for a strictly circular orbit to derive an orbital period of 16.7 d, requiring a high inclination perspective of $i\sim 80^\circ$. We also consider an elliptical orbit under simplifying assumptions, yielding a match for $i\sim 75^\circ$ with eccentricity $e\sim 0.5$ and a longer orbital period of about 70 d. The prevalence of binarity among massive stars is well-known; the prospect of detecting a binary companion during the post-LBV stage of WR 31a would add to an emerging narrative of diverse interactions between massive multiple components as a function of evolutionary stage. However, if the loop originates because of a co-rotating interaction region (CIR), then the rotation period could be 8.5 d or 17 d. This would give an estimated equatorial rotation speed of 95 or 190 km/s. Either of these is a significant fraction of the estimated critical speed of rotational break-up at 320 km/s (for an Eddington factor of $\Gamma=0$).

State-of-the-art long-term solar system integrations include several second order effects such as the Sun's quadrupole moment J2 and a contribution from asteroids (plus the Moon and general relativity). We recently showed that including 10 asteroids and a reduced J2 in our astronomical solutions provides the best match with geologic data to -58 Myr. However, the rationale for the reduced J2 remained ambiguous and may suggest that parameters for long-term integrations compatible with geologic observations are not fully compatible with our knowledge of the current solar system (specifically J2). Here we show that a reduced J2 compensates for a diminished asteroid population in long-term solar system integrations, which may appear surprising. We present an analysis and offer a mechanism for the long-term compensating effects of J2 and asteroid mass in the solar system (not planetary systems in general). Our analysis suggests that "differential effects" on specific secular frequencies involved in resonant terms (i.e., (g4-g3) and (s4-s3)), are critical in the long term, rather than short-term effects on the orbital elements of individual planetary orbits across the board. Also, our results indicate that if long-term intergrations including the full asteroid population were computationally feasible, a J2 value (within errors) compatible with our current knowledge of the solar system could be used. Attempts to improve the long-term accuracy of astronomical solutions by, e.g., tinkering with initial conditions using current/future astronomical observations are futile unless asteroid deficiencies in the solar system model are addressed.

Future high-precision X-ray and gravitational wave observations of neutron stars (NSs) are expected to measure NS radii to better than $\sigma=0.1$ km accuracy, providing unprecedented opportunities to extract novel information about the nature and equation of state (EOS) of supradense matter in NS cores. Within a Bayesian framework using a meta-model for NS EOS encapsulating a first-order hadron-quark phase transition and satisfying all known constraints from both nuclear physics and astrophysics, we investigate how NS radius data with higher precision may better inform us about (1) the NS crust-core transition density $\rho_{cc}$, (2) the hadron-quark transition density $\rho_t$, quark matter fraction $F_{\rm{QM}}$ and its radius $R_{\rm{QM}}$, and (3) high-density NS EOS parameters. Using fiducial NS radius data with mocked precisions varying from $\sigma=1.0$ km to 0.1 km, we found (a) the most probable crust-core transition density $\rho_{cc}$ and its 68\% confidence boundaries are essentially unaffected by $\sigma$ especially for massive NSs; (b) our answers to the questions (2) and (3) listed above depend sensitively on the prior range of hadron-quark transition density $\rho_t$ assumed. Using its fiducial range of $(1.0-6.0)\rho_0$, the posterior PDF($\rho_t$) has a major peak around $(1.7-2.0)\rho_0$ that is sufficient but unnecessary in describing all existing NS radius data, and a minor peak around $(3.0-5.0)\rho_0$ consistent with the indication about $\rho_t$ of recent Beam Energy Scan Experiments at RHIC. Narrowing down the prior range of $\rho_t$ to $(3.0-6.0)\rho_0$, NS radius data with smaller $\sigma$ can constrain more stringently the posterior PDF($\rho_t$), $F_{\rm{QM}}$, $R_{\rm{QM}}$ and several high-density hadronic EOS parameters.

The first simultaneous observations of the Fe XVIII 974.86 Å and Fe XX 721.56 Å forbidden lines from the Spectral Imaging of the Coronal Environment (SPICE) spectrograph on Solar Orbiter are presented. The lines were observed from the post-flare loops of an M2.5 class solar flare that peaked at 23:49 UT on 2024 March 23. The Fe XX/Fe XVIII ratio is used to derive a temperature space-time map for the flaring period, with values ranging from 8 to 20 MK. The map reveals repeated episodes of heating at the SPICE slit location over a 30 min period. For one location with the brightest emission, the plasma cooled from 10 MK to 8 MK in 5 min which is longer than the expected conductive cooling time of 170 s, suggesting continued background heating during the cooling period. Doppler shifts of between 0 and +10 km/s were obtained with precisions of 1-4 km/s, but the accuracies are lower due to uncertainties over the absolute wavelength calibration hence we can not conclude there are plasma flows in the flare loops. The widths of the two lines were found to be close to the instrumental widths with no evidence of non-thermal broadening, although this result is limited by the instrument resolution. The Fe XVIII and Fe XX lines have high signal-to-noise with only a 5 s exposure time, demonstrating that the lines will be valuable for high-cadence flare studies with SPICE.

There are several astrophysical configurations where one is interested only in the long-term dynamical evolution. Although the first-order version of this approximation is usually sufficient in applications, second-order corrections may be relevant, too. Here we use the Hamiltonian formalism to show how such higher-order terms lead to the long-term evolution of the semi-major axis.

In this letter, we investigate cosmology within the framework of modified $f(Q, L_m)$ gravity using the non-linear model $f(Q, L_m) = -Q + \alpha L_m^n + \beta$, where $\alpha$, $\beta$, and $n$ are free parameters. The modified Friedmann equations are derived for a matter-dominated universe, and an analytical solution is obtained. Using Hubble, Pantheon+, and joint datasets, we constrain the model parameters and examine their implications. The results indicate that the model accommodates varying $H_0$ values, contributing to the Hubble tension, while $n \neq 1$ suggests deviations from general relativity (GR). The deceleration parameter confirms a transition to acceleration at $z_t \approx 0.6 - 0.8$, with present values supporting cosmic acceleration. This model offers a viable alternative to GR-based cosmology.

When a neutron star (NS) intercepts gas from a non-degenerate star, e.g., in a tidal disruption event, a common-envelope phase, or the collapsing core of a massive star, photons become trapped in the hot flow around the NS. This gas forms a radiatively inefficient accretion flow (RIAF) until the density and temperature close to the NS surface grow large enough for binding energy to be converted to neutrinos. Here we present three-dimensional, general-relativistic, magnetohydrodynamic simulations of accretion onto a non-rotating, unmagnetized NS. These connect, for the first time, an extended accretion disk with a self-consistent hydrostatic atmosphere around the star. The impact of different seed magnetic fields and accretion rates is studied by approximating the radiation-pressure dominated flow as an ideal gas with an adiabatic index of $4/3$, coupled to a variable neutrino emissivity. At low accretion rates, the hydrostatic atmosphere shows slow rotation and weak magnetization, transitioning to an outer RIAF structure. A toroidal magnetic field mediates the inward flow of energy and angular momentum through the atmosphere, which reaches a steady state when neutrino emission balances the accretion power. We develop a one-dimensional analytical model connecting these results with more general initial conditions and describing the main features of the flow. Our results have implications for the spin and mass evolution of hypercritically accreting NSs.

The prestellar core Barnard 68 (B68) is a prototypical source to study the initial conditions and chemical processes of star formation. A previous numerical simulation suggested the southeastern bullet is impacting on the main body of B68. In order to obtain more observational evidence, mapping observations of the ground state SO ($1_0-0_1$) emission line at 30 GHz were made with the Effelsberg 100 m telescope. Based on the velocity field and channel maps derived from SO, three velocity components were clearly detected. The velocity field of the main body indicates rotation and is well fitted by a solid-body rotation model. The measured radial velocity difference between the bullet and the main core is about 0.4 km s$^{-1}$, which is almost equal to the velocity obtained by the previous numerical simulation. Therefore, the bullet is most likely impacting onto the rotating main body of B68. A 1D spherical non-LTE Monte-Carlo radiation transfer RATRAN code is performed to derive the radial abundance profile of SO by analyzing the observed velocity-integrated intensity. SO is depleted inside a 60$^{\prime\prime}$ (0.02 pc) radius from the core. The abundance stays constant at 2.0$\times$10$^{-9}$ for radii larger than 60$^{\prime\prime}$ from the center of the main core. The abundance is enhanced at the interface of the bullet and the main core indicating that shock waves were produced by the collision between the bullet and the main core. In conclusion, based on the kinematical and chemical analysis, our observational results support the previously proposed core-core collision scenario in B68.

In $f(R)$ gravity, the scalaron$\unicode{x2014}$a scalar degree of freedom arising from modification of General Relativity$\unicode{x2014}$could account for all dark matter if its mass lies in the meV$\unicode{x2013}$MeV range. In this work, we revisit the scalaron's interactions with Standard Model particles, assuming their minimal coupling to gravity. In particular, we refine the scalaron decay rate into photons. Assuming the scalaron constitutes all of dark matter, we calculate the average cosmological background radiation produced by these decays. We also estimate the contribution of primordial thermal scalarons to the present dark matter density and find it to be negligible. This supports the original scenario in which the scalaron dark matter behaves as a coherently oscillating field.

Huairou Solar Observing Station of the National Astronomical Observatories of the Chinese Academy of Sciences has been in operation since 1987. Successful observations of the solar vector magnetic field have been conducted during its operation. Based on the achievements at Huairou, we analyze the methods of observing the solar magnetic field, including discussions of the approximation of the transfer theory of the solar magnetic field in the atmosphere, wide field of view polarized observation, and some questions on the inversion of solar magnetic field data. We also present relevant challenges for further research.

(abridged) Weak gravitational lensing (WL) is the unique and powerful probe into the large-scale structures of the Universe. Removing the shape noise from the observed WL field, i.e., denoising, enhances the potential of WL by accessing information at small scales where the shape noise dominates without denoising. We utilise two machine learning (ML) models for denosing: generative adversarial network (GAN) and diffusion model (DM). We evaluate the performance of denosing with GAN and DM utilising the large suite of mock WL observations, which serve as the training and test data sets. We apply denoising to 1,000 maps with GAN and DM models trained with 39,000 mock observations. Both models can fairly well reproduce the true convergence map on large scales. Then, we measure cosmological statistics: power spectrum, bispectrum, one-point probability distribution function, peak and minima counts, and scattering transform coefficients. We find that DM outperforms GAN in almost all statistics and recovers the correct statistics down to small scales within roughly $0.3 \sigma$ level at the scales accessible from current and future WL surveys. We also conduct the stress tests on the trained model; denoising the maps with different characteristics, e.g., different source redshifts, from the data used in training. The performance degrades at small scales, but the statistics can still be recovered at large scales. Though the training of DM is more computationally demanding compared with GAN, there are several advantages: numerically stable training, higher performance in the reconstruction of cosmological statistics, and sampling multiple realisations once the model is trained. It has been known that DM can generate higher-quality images in real-world problems than GAN, the superiority has been confirmed as well in the WL denoising problem.

Recently, deep space exploration, especially focusing on halo orbits, the periodic orbits of the Moon, has been widely studied. The spacecraft in halo orbits performs periodic orbital motion, which affects the attitude motion by periodic disturbances. The conventional attitude control method, PD control, is widely used, but its application to periodic disturbance attenuation is inefficient. To address these challenges, this study proposes a predictive Repetitive Control (RC) approach that addresses periodic disturbances, particularly GG torque, by exploiting the periodic nature of the system dynamics. The proposed method is also applied to the case of using a Reaction Wheel (RW) as an attitude control actuator. Despite the inherent challenges posed by RW limitations, including saturation torque and transmission delay, our predictive RC approach effectively mitigates these effects. Numerical simulations demonstrate the robust performance of the proposed method in maintaining attitude control for spacecraft traversing halo orbits near the Earth-Moon $L_2$ point, validating its potential for future deep space exploration missions.

Planet formation via core accretion involves the growth of solids that can accumulate to form planetary cores. There are a number of barriers to the collisional growth of solids in protostellar discs, one of which is the drift, or metre, barrier. Solid particles experience a drag force that will tend to cause them to drift towards the central star in smooth, laminar discs, potentially removing particles before they grow large enough to decouple from the disc gas. Here we present 3-dimensional, shearing box simulations that explore the dynamical evolution of solids in a protostellar disc that is massive enough for the gravitational instability to manifest as spiral density waves. We expand on earlier work by considering a range of particle sizes and find that the spirals can still enhance the local solid density by more than an order of magnitude, potentially aiding grain growth. Furthermore, if solid particles have enough mass, and the particle size distribution extends to sufficiently large particle sizes, the solid component of the disc can undergo direct gravitational collapse to form bound clumps with masses typically between $1$ and $10$ M$_\oplus$. Thus, the concentration of dust in a self-gravitating disc could bypass the size barrier for collisional growth and directly form planetary cores early in the lifetime of the disc.

The power spectra of the redshifted 21-cm signal from the Epoch of Reionization (EoR) contain information about the ionization and thermal states of the intergalactic medium (IGM), and depend on the properties of the EoR sources. Recently, Mertens et al 2025 has analysed 10 nights of LOFAR high-band data and estimated upper limits on the 21-cm power spectrum at redshifts 8.3, 9.1 and 10.1. Here we use these upper limit results to constrain the properties of the IGM at those redshifts. We focus on the properties of the ionized and heated regions where the temperature is larger than that of the CMB. We model the 21-cm power spectrum with the code GRIZZLY, and use a Bayesian inference framework to explore the source parameters for uniform priors on their ranges. The framework also provides information about the IGM properties in the form of derived parameters. In a model which includes a radio background in excess of the CMB, the 95 (68) per cent credible intervals of disfavoured models at redshift 9.1 for the chosen priors correspond to IGM states with averaged ionization and heated fraction below 0.46 ($\lesssim 0.05$), an average gas temperature below 44 K (4 K), and a characteristic size of the heated region $\lesssim 14 ~h^{-1} ~\mathrm{Mpc}$ ($\lesssim 3 ~h^{-1} ~\mathrm{Mpc}$). The 68 per cent credible interval suggests an excess radio background which is more than 100 per cent of the CMB at 1.42 GHz, while the 95 per cent credible interval of the radio background efficiency parameter spans the entire prior range. The behaviour of the credible intervals is similar at all redshifts. The models disfavoured by the LOFAR upper limits are extreme ones, as they are mainly driven by rare and large ionized or heated regions.

Comet--asteroid transition (CAT) objects are small solar system bodies in the process of evolving from cometary nuclei into asteroids, as they gradually lose volatile substances due to solar heating. The volatile material is mainly water ice, and the time required for its complete depletion is called the desiccation time. Estimating the desiccation time is important for examining the formation and evolution of small solar system bodies. Here, we propose a new theoretical model for evaluating the desiccation time as a function of orbital elements, considering the contraction of the entire cometary nucleus due to ice sublimation. First, we performed numerical calculations of the thermal evolution of a cometary nucleus in an eccentric orbit, considering the seasonal variation in the solar heating rate. Next, we derived the desiccation time analytically as a function of orbital elements based on a steady-state model considering the solar heating rate averaged over the seasons. We compared the numerical solutions for the desiccation time with the analytical solutions and clarified the conditions under which the analytical model can be applied. Additionally, based on the analytical model, we derived formulae for estimating the emission rates of water vapor and dust on the surface of the cometary nucleus, the maximum size of the emitted dust, and the dust emission velocity, by assuming the amount of ice remaining inside the nucleus. Using these analytical solutions, we considered the internal structure and evolution process of typical CAT objects. Our analytical model was generally consistent with that of the results of earlier observations of these objects. Our model provides a theoretical guideline for discussing the evolution of cometary nuclei and the possibility of retaining internal ice in asteroids.

The radio emission of the quiet Sun in the metric and decametric bands has not been well studied historically due to limitations of existing instruments. It is nominally dominated by thermal brehmsstrahlung of the solar corona, but may also include significant gyrosynchrotron emission, usually assumed to be weak under quiet conditions. In this work, we investigate the expected gyrosynchrotron contribution to solar radio emission in the lowest radio frequencies observable by ground instruments, for different regions of the low and middle corona. We approximate the coronal conditions by a synoptic magnetohydrodynamic (MHD) model. The thermal emission is estimated from a forward model based on the simulated corona. We calculate the expected gyrosynchrotron emission with the Fast Gyrosynchrotron Codes framework by Fleishman & Kuznetsov (2010). The model emissions of different coronal regions are compared with quiet-time observations between 20-90 MHz by the LOw Frequency ARray (LOFAR) radio telescope. The contribution of gyrosynchrotron radiation to low frequency solar radio emission may shed light on effects such as the hitherto unexplained brightness variation observed in decametric coronal hole emission, and help constrain measurements of the coronal magnetic fields. It can also improve our understanding of electron populations in the middle corona and their relation to the formation of the solar wind.

The accuracy of a mass model in the strong lensing analysis is crucial for unbiased predictions of physical quantities such as magnifications and time delays. While the mass model is optimized by changing parameters of the mass model to match predicted positions of multiple images with observations, positional uncertainties of multiple images often need to be boosted to take account of the complex structure of dark matter in lens objects, making the interpretation of the chi-square value difficult. We introduce the Jackknife method as a new method to validate strong lens mass models, specifically focusing on cluster-scale mass modeling. In this approach, we remove multiple images of a source from the fitting and optimize the mass model using multiple images of the remaining sources. We then calculate the multiple images of the removed source and quantitatively evaluate how well they match the observed positions. We find that the Jackknife method performs effectively in simulations using a simple model. We also demonstrate our method with mass modeling of the galaxy cluster MACS J0647.7+7015. We discuss the potential of using the Jackknife method to validate the error estimation of the physical quantities by the Markov Chain Monte Carlo.

Pairwise velocities of the large-scale structure encode valuable information about the growth of structure. They can be observed indirectly through redshift-space distortions and the kinetic Sunyaev-Zeldovich effect. Whether it is Gaussian or non-Gaussian, pairwise velocity has a broad distribution, but the cosmologically useful information lies primarily in the mean - the streaming velocities; the dispersion around the mean is often treated as a nuisance and marginalized over. This conventional approach reduces the constraining power of our observations. Here, we show that this does not have to be the case, provided the physics behind the dispersion is understood. We demonstrate that by splitting the halo/galaxy samples according to their density environments and measuring the streaming velocities separately, the total signal-to-noise is several times greater than in conventional global measurements of the pairwise velocity distribution (PVD). This improvement arises because the global PVD is a composite of a series of near-Gaussian distributions with different means and dispersions, each determined by its local density environment. Around underdense and overdense regions, the mean streaming velocities are positive and negative, respectively. By splitting the data, we avoid cancellation between these opposing velocities, effectively turning what would be considered dispersion in the global PVD into a signal. Our findings indicate substantial potential for improving the analysis of PVD observations using the kinetic Sunyaev-Zeldovich effect and redshift-space distortions.

Simulation-based inference provides a powerful framework for extracting rich information from nonlinear scales in current and upcoming cosmological surveys, and ensuring its robustness requires stringent validation of forward models. In this work, we recast forward model validation as an out-of-distribution (OoD) detection problem within the framework of machine learning (ML)-based simulation-based inference (SBI). We employ probability density as the metric for OoD detection, and compare various density estimation techniques, demonstrating that field-level probability density estimation via continuous time flow models (CTFM) significantly outperforms feature-level approaches that combine scattering transform (ST) or convolutional neural networks (CNN) with normalizing flows (NFs), as well as NF-based field-level estimators, as quantified by the area under the receiver operating characteristic curve (AUROC). Our analysis shows that CTFM not only excels in detecting OoD samples but also provides a robust metric for model selection. Additionally, we verified CTFM maintains consistent efficacy across different cosmologies while mitigating the inductive biases inherent in NF architectures. Although our proof-of-concept study employs simplified forward modeling and noise settings, our framework establishes a promising pathway for identifying unknown systematics in the cosmology datasets.

Kilohertz quasi-periodic oscillations (kHz QPOs) are believed to originate from the orbital timescales of the inner accretion flow, reflecting the dynamics of the innermost disk regions under strong gravitational forces. Despite numerous radiative and geometric models proposed so far, a comprehensive explanation of the observed properties of these variability components remains elusive. This study systematically examines kHz QPOs, their variability, and their connection to spectral properties in $4U 1636-536$ using AstroSat data. Our analysis tracks the source transition from hard to soft states in the hardness-intensity diagram. Broad spectral analysis (0.7-25 keV) using SXT and LAXPC data indicates a spectrum shaped by reflection from a thermal corona, with contributions from boundary layer emission and a soft disk component. We find significant changes in optical depth, blackbody temperature, and inner disk temperature that likely drive state transitions. Power density spectra reveal three variability types: a low frequency QPO (LFQPOs) (~30 Hz), and two simultaneous kHz QPOs. The LFQPOs and the upper kHz QPOs appear more prominently in soft spectral states. The presence of LFQPOs and twin kHz QPOs in soft spectral states enable us to estimate the neutron star mass at (2.37 $\pm$ 0.02) $M_\odot$ using the relativistic precession model (RPM). Additionally, time-lag and root mean square (rms) analysis provide insights into the size of the corona and the radiative origin of these variability components.

We present JWST observations of the near-Earth asteroid (3200) Phaethon using the Near-Infrared Camera (NIRCam), Near-Infrared Spectrograph (NIRSpec), and Mid-Infrared Instrument (MIRI) to further investigate the composition of Phaethon's surface. Our NIRSpec data confirms that Phaethon's surface is dehydrated, showing no evidence of hydrated minerals in the 3-$\mu$m region. We estimate an upper limit on the hydrogen content in phyllosilicates of 0.06 wt%. Comparisons with laboratory spectra of carbonaceous chondrites suggest that Phaethon's surface composition is best matched by thermally metamorphosed samples of the CM chondrite Murchison (heated to 1000$^{\circ}$C), rather than CY meteorites as previous work suggested. We find no evidence of ongoing surface evolution due to recent perihelion passages. A comparison of the mid-infrared spectra of Phaethon and Bennu shows distinct spectral differences that are consistent with their different thermal histories. Our findings further refine our understanding of Phaethon's current surface composition and evolution and provide additional insights for the upcoming DESTINY+ mission.

Gamma-ray bursts (GRBs) are powerful probes of the high-redshift universe. However, the proportion of collapsar GRBs among long GRBs and their event rate relative to the star formation rate (SFR) remain contentious issues. We assume that long GRBs with $z\geq 2$ are all collapsar GRBs and construct the luminosity function using a high-redshift sample from the Swift satellite spanning 2004 to 2019. We model the luminosity function with a broken power-law form and consider three scenarios: no evolution, luminosity evolution, and density evolution. Our results are as follows: 1) The no-evolution model can be ruled out. 2) The fitting results indicate that to adequately explain the observations, a significant redshift evolution in either luminosity (evolution index $\delta = 1.54^{+0.21}_{-0.22}$) or density ($\delta = 2.09^{+0.29}_{-0.26}$) is required. This excludes the possibility that the evolution of long GRBs with redshift is due to contamination from non-collapsar GRBs. 3) The luminosity evolution model predicts that the number of collapsar GRBs with $z<2$ and $P \geq 1 \, \text{ph} \, \text{cm}^{-2} \, \text{s}^{-1}$ is 138.6, accounting for 82.5% of the observed long GRBs with $z<2$ and $P \geq 1 \, \text{ph} \, \text{cm}^{-2} \, \text{s}^{-1}$. The density evolution model predicts that the number of collapsar GRBs with $z<2$ and $P \geq 1 \, \text{ph} \, \text{cm}^{-2} \, \text{s}^{-1}$ is 80.2, accounting for 47.7% of the observation. Regardless of the model, a substantial portion of the long GRBs are not collapsar GRBs.

Using a model-independent analysis method, we analyze $H(z)$ parameter data with some restrictive conditions. We find that: (a) the Universe may experience an accelerated expansion with a confidence level greater than 5 $\sigma$ at redshift $z_{101}\in (0, 0.36)$ and greater than 1.9 $\sigma$ at redshifts $z_{3835}\in (1.3, 1.53)$ and $z_{3836}\in (1.43, 1.53)$; (b) the Universe may experience a decelerated expansion with a confidence level greater than 1.5 $\sigma$ at redshift $z_{2012}\in (0.40, 0.52)$; (c) $w_{\rm{x}}\leq w_{\rm{t}}<-1$ with confidence level great than 1.6 $\sigma$ at redshift $z_{3836}\in (1.43, 1.53)$.

The dispersion in the speed of gravitational waves is a novel way to test the general theory of relativity and understand whether the origin of cosmic acceleration is due to any alternative theory of gravity. Several alternative theories of gravity predict dispersion in the gravitational wave signal in a frequency-dependent deviation from the speed of light at lower frequencies than accessible from current ground-based detectors. We show how a multi-band observation of gravitational wave signal combining deci-Hertz gravitational wave signal from LGWA (Lunar Gravitational Wave Antenna) with ground-based detectors such as Cosmic Explorer or Einstein Telescope, and also including LISA (Laser Interferometer Space Antenna), we can probe the energy scale associated with effective theory of gravity with a precision of approximately $8.6\%$ by combining only $\mathcal{O}(10)$ high SNR multi-band gravitational wave events. This precision will further improve with the inclusion of more events as $\sqrt{N}$. In the future, this measurement will shed light on an unexplored domain of fundamental physics and will bring deeper insights into the phenomenon of cosmic acceleration. The operation of the gravitational wave detector in the deci-Hertz frequency band is key to exploring this frontier of fundamental physics.

Modern gravitational-wave science demands increasingly accurate and computationally intensive numerical relativity (NR) simulations. The Python-based, open-source NRPy framework generates optimized C/C++ code for NR, including the complete NR code BlackHoles@Home (BH@H), which leverages curvilinear coordinates well-suited to many astrophysical scenarios. Historically, BH@H was limited to single-node OpenMP CPU parallelism. To address this, we introduce superB, an open-source extension to NRPy that enables automatic generation of scalable, task-based, distributed-memory Charm++ code from existing BH@H modules. The generated code partitions the structured grids used by NRPy/BH@H, managing communication between them. Its correctness is validated through bit-identical results with the standard OpenMP version on a single node and via a head-on binary black hole simulation in cylindrical-like coordinates, accurately reproducing quasi-normal modes (up to $\ell=8$). The superB/NRPy-generated code demonstrates excellent strong scaling, achieving an $\approx 45$x speedup on 64 nodes (7168 cores) compared to the original single-node OpenMP code for a large 3D vacuum test. This scalable infrastructure benefits demanding simulations and lays the groundwork for future multi-patch grid support, targeting long inspirals, extreme parameter studies, and rapid follow-ups. This infrastructure readily integrates with other NRPy/BH@H-based projects, enabling performant scaling for the general relativistic hydrodynamics code GRoovy, and facilitating future coupling with GPU acceleration via the NRPy-CUDA project.

The equilibrium configurations of slowly rotating anisotropic self-gravitating fluids are computed using the extended Hartle structure equations, including anisotropic effects, derived in our previous paper. We focus on the so-called $\mathcal{C}$-star, whose anisotropic pressure follows a fully covariant equation of state (EoS), while a standard polytrope describes the radial pressure. We determine surface and integral properties, such as the moment of inertia, mass change, mass quadrupole moment, and ellipticity. Notably, for certain values of the compactness parameter, highly anisotropic $\mathcal{C}$-stars exhibit a prolate shape rather than the typical oblate form, an intriguing behavior also observed in other anisotropic systems like Bowers-Liang spheres and stars governed by a quasi-local EoS. Although the $\mathcal{C}$-stars considered in this study are limited by stability criteria and cannot sustain compactness beyond $M/R\approx0.38$, we found indications that certain rotational perturbations exhibit similarities to those observed in other ultracompact systems approaching the black hole limit.

Periodic sediment patterns have been observed on Earth in riverbeds and sand and snow deserts, but also in other planetary environments. One of the most ubiquitous patterns, familiar wind or 'impact' ripples, adorns sand beaches and arid regions on Earth. The observation of aeolian impact ripples on Mars the same size as their terrestrial counterparts despite a thinner atmosphere raises questions about their formation. Here we show in a numerical simulation that the emergent wavelength of impact ripples is controlled by the mechanics of grain-bed impacts and not the characteristic trajectories of grains above the bed. We find that the distribution of grain trajectories in transport is essentially scale-free, invoking the proximity of a critical point and precluding a transport-related length scale that selects ripple wavelengths. By contrast, when a grain strikes the bed, the process leading to grain ejections introduces a collective granular length scale that determines the scale of the ripples. We propose a theoretical model that predicts a relatively constant ripple size for most planetary conditions. In addition, our model predicts that for high-density atmospheres, such as on Venus, or for sufficiently large sand grains on Earth, impact ripples propagate upwind. Although wind-tunnel and field experiments are needed to confirm the existence of such 'antiripples', we suggest that our quantitative model of wind-blown sediment transport may be used to deduce geological and environmental conditions on other planets from the sizes and propagation speeds of impact ripples.

Our present knowledge of the nuclear equation of state is briefly reviewed in this article intended for a wider readership. Particular emphasis is given to the asymmetric-matter equation of state required for modeling neutron stars, neutron-star mergers, and r-process nucleosynthesis. Recent analyses based on combining information obtained from nuclear theory, heavy-ion collisions and astrophysical observations confine the obtained radii of the canonical 1.4-solar-mass neutron star to values between 12 km and 13 km. The remaining uncertainty is primarily related to missing information in the density interval between nuclear saturation density and about twice that value which, however, is accessible with laboratory experiments.

We investigate the impact of spacetime non-commutativity on the tidal deformability of compact objects and explore the feasibility of detecting non-commutative (NC) effects through gravitational wave (GW) observations. We considered NC modifications to spacetime geometry based on de Sitter gauge theory of gravity and calculate their impact on tidal deformability. While several types of compact objects have been proposed as candidates for probing spacetime non-commutativity, particularly at the horizon scales, our study showed analytically that, for compact objects with non-singular metric at their surface (such as neutron stars and boson stars), the NC correction to their tidal deformability converge to a finite value at the black-hole-compactness limit, eliminating infinite enhancement at the horizon scales. We then compute the NC corrections for neutron stars and boson stars, considering several different models, and analyze their imprints on the GW signals. By comparing the results, we assess the scale of NC effects across different compactness regimes and discuss the conditions under which these NC effects can be amplified. While our findings suggest that the leading-order NC correction dominates the tidal deformability of a compact object near the black-hole-compactness limit, we demonstrate that neutron stars and boson stars are not viable candidates to constrain spacetime non-commutativity, while relying on the tidal deformability through GW observations.

We construct a set of unified equations of state based on the quark mean field (QMF) model, calibrated to different values of nuclear symmetry energy slope at the saturation density ($L_0$), with the aim of exploring both the static properties and dynamical behavior of neutron stars (NSs), and building a coherent picture of their internal structure. We assess the performance of these QMF models in describing the mass-radius relation, the cooling evolution of isolated NSs and X-ray transients, and the instabilities (e.g., the r-mode). In comparison to relativistic mean field (RMF) models formulated at the hadronic level, the QMF model predicts heavier nuclear clusters and larger Wigner-Seitz cell sizes in the NS crust, while the density of the free neutron gas remains largely similar between the two approaches. For the cooling of isolated NSs, the thermal evolution is found to be insensitive to both the many-body model and the symmetry energy slope in the absence of the direct Urca (dUrca) process. However, when rapid cooling via the dUrca process is allowed, in the case of large $L_0$ values (e.g., $L_0 \gtrsim 80$ MeV) in our study, the QMF model predicts a longer thermal relaxation time. Both the QMF and RMF models can reproduce cooling curves consistent with observations of X-ray transients (e.g., KS 1731--260) during their crustal cooling phase, although stellar parameters show slight variations depending on the model and symmetry energy slope. Within our unified framework, a larger $L_0$ value generally results in a wider instability window, while increasing the stellar mass tends to suppress the instability window. We also provide simple power-law parameterizations that quantify the dependence of bulk and shear viscosities on the symmetry energy slope for nuclear matter at saturation density.

We show that a vector field non-minimally coupled to gravity reproduces the dynamics of an Einstein cluster. Our results suggest a vectorial nature for dark matter.

Gravitational wave data are often contaminated by non-Gaussian noise transients, glitches, which can bias the inference of astrophysical signal parameters. Traditional approaches either subtract glitches in a pre-processing step, or a glitch model can be included from an agnostic wavelet basis (e.g. BayesWave). In this work, we introduce a machine-learning-based approach to build a parameterised model of glitches. We train a normalising flow on known glitches from the Gravity Spy catalogue, constructing an informative prior on the glitch model. By incorporating this model into the Bayesian inference analysis with Bilby, we estimate glitch and signal parameters simultaneously. We demonstrate the performance of our method through bias reduction, glitch identification and Bayesian model selection on real glitches. Our results show that this approach effectively removes glitches from the data, significantly improving source parameter estimation and reducing bias.