Papers for Wednesday, Jan 07 2026

RS Oph is a recurrent nova, a kind of cataclismic variable that shows bursts in a period approximately shorter than a century. Persistent homology, a technique from topological data analysis, studies the evolution of topological features of a simplicial complex composed of the data points or an embedding of them, as some distance parameter is varied. For this work I trained a supervised learning model based on several featurizations, namely persistence landscapes, Carlsson coordinates, persistent images, and template functions, of the persistence diagrams of sections of the lightcurve of RS Oph. A tenfold cross validation of the model based on one of the featurizations, persistence landscapes, consistently shows high recalls and accuracies. This method serves the purpose of predicting whether RS Oph is bursting within a year.

A model one-dimensional self consistent steady state collisionless self-gravitating system in which all the particles have the same energy is presented. This has the remarkable property that the position and velocity of the particles orbiting in their own self consistent potential are given exactly, in terms of time, by the truncations of sine and cosine functions to the first two terms in their respective Taylor series. The potential and density also have simple analytic expressions in terms of time as parameter. It is not being claimed that this system has any direct astronomical application. However, it does motivate a conjecture about the behaviour of the density, potential, and orbits near caustics in simulations of cold collisionless dark matter. It is a rather surprising result which might interest practitioners of stellar dynamics and serve as an elementary example in teaching the subject.

The majority of stars more massive than the Sun is found in binary or multiple star systems and many of them will interact during their evolution. Specific interactions, where progenitors and post-mass transfer (MT) systems are clearly linked, can provide yet missing observational constraints. Volume-complete samples of progenitor and post-MT systems are well suited to study those processes. To compile them, we need to determine the parameters of thousands of binary systems with periods spanning several orders of magnitude. The bottleneck are the orbital parameters, because accurate determinations require a good coverage of the orbital phases. The next generation of time-resolved spectroscopic surveys should be optimized to follow-up and solve whole populations of binary systems in an efficient way. To achieve this, a scheduler predicting the best times of the next observation for any given target in real time should be combined with a flexibly schedulable multi-object spectrograph or ideally a network of independent telescopes.

We validate Leonid Space's satellite lifetime prediction pipeline through comprehensive backtesting against 934 non-maneuvering satellites that deorbited from LEO between 1961 and 2024. This represents the first large-scale validation of lifetime prediction tooling using forecasted space weather conditions rather than historical hindsight. Our toolchain combines ballistic coefficient estimation from on-orbit data with probabilistic orbit propagation under varying environmental conditions. Using TLE data and space weather records spanning six solar cycles, our three-stage validation approach progressively removes hindsight bias to arrive at fully predictive operational conditions. We achieve 1-year prediction accuracy (median continuously ranked probability score) of 6.0 days (1.6%) under perfect knowledge conditions, 18.6 days (5.1%) with estimated ballistic coefficients and known space weather, and 45.5 days (12.4%) under fully predictive conditions. Comparison against ESA's standard DRAMA & DISCOS toolchain demonstrates a 4x improvement in state-of-the-art accuracy for well-characterized satellites. A custom semianalytic propagator provides a >3500x speedup over Orekit and 4.5x speedup over DRAMA, enabling rapid Monte Carlo analysis across large satellite populations. Our analysis reveals that solar cycle forecasting dominates error budgets after ballistic coefficient estimation, with higher-fidelity propagators and atmosphere models providing marginal benefit. These results establish a validated performance baseline for operational lifetime prediction services supporting LEO mission planning and regulatory compliance.

We study a minimal type-I seesaw framework in which a first-order phase transition (FOPT), driven by a singlet scalar, produces right-handed neutrinos (RHNs) through bubble collisions, realizing a cosmic-scale collider that probes ultra-high energy scales. The resulting RHN distribution sources novel low-frequency gravitational-waves (GWs) in addition to the standard bubble-collision contribution. A stable lightest RHN can account for the observed dark matter (DM) relic abundance for masses as low as $M_{1} \equiv m_{\rm DM} \gtrsim 10^{6}\,\mathrm{GeV}$, with the associated novel GW signal accessible in LISA, ET and upcoming LVK detectors. If the RHNs are unstable, their CP-violating decays generate the observed baryon asymmetry via leptogenesis for $M_{1} \gtrsim 10^{11}\,\mathrm{GeV}$ and phase transition temperatures $T_* \gtrsim 10^{6}\,\mathrm{GeV}$, for which the novel GW spectrum is detectable in ET, BBO and upcoming LVK. If RHN decays also populate a dark-sector fermion with mass $m_{\chi} \in [10^{-4},10^{4}],\mathrm{GeV}$, successful co-genesis of baryons and asymmetric dark matter occurs for $T_* \gtrsim 10^{7}\,\mathrm{GeV}$ and $M_{1} \gtrsim 10^{9}\,\mathrm{GeV}$, naturally explaining $\Omega_{\rm DM} \simeq 5\Omega_{\rm B}$. The corresponding GW signals are testable with LISA, ET, and BBO. Finally, we analyze a UV-complete multi-Majoron model, based on a global $U(1)_N \times U(1)_{\rm B-L}$ extension, motivated from the hierarchy of lepton masses, which we dub as Mojaron collider. The corresponding FOPT in this model leaves a distinctive GW signature arising from RHN production during $U(1)_N$ symmetry breaking detectable by BBO, ET and upcoming LVK. Successful leptogenesis is realized for heaviest RHN mass $M_3 \sim 10^{10}\,\mathrm{GeV}$ and a $U(1)_N$ breaking vev $v_2 \sim \mathcal{O}(\mathrm{TeV})$, which sets the seesaw scale.

Relativistic magnetic reconnection is thought to power various multi-wavelength emission signatures from neutron stars and black holes. Relativistic resistive magnetohydrodynamics (RRMHD) offers the simplest model of reconnection. However, a small uniform resistivity underestimates the reconnection rate compared to first-principles kinetic models. By employing an effective resistivity based on kinetic models - which connects the reconnection electric field to the charge-starved current density - we show that RRMHD can reproduce the increased reconnection rate of kinetic models, both in local current sheets and in global black hole magnetospheres.

Determining where, when, and how neutrino flavor oscillations must be included in large-scale simulations of hot and dense astrophysical environments is an enduring challenge that must be tackled to obtain accurate predictions. Using an angular moment-based linear stability analysis framework, we examine the different kinds of flavor instabilities that can take place in the context of the post-processing of a neutron star merger simulation, with a particular focus on the collisional flavor instability and a careful assessment of several commonly used approximations. First, neglecting anisotropies of the neutrino field, we investigate the extent to which commonly used monoenergetic growth rates reproduce the results obtained from a full multi-energy treatment. Contrary to the large discrepancies found in core-collapse supernova environments, we propose a simple combination of energy-averaged estimates that reproduces the multi-energy growth rates in our representative simulation snapshot. We then quantify the impact of additional physical effects, including nuclear many-body corrections, scattering opacities, and the inclusion of the vacuum term in the neutrino Hamiltonian. Finally, we include the neutrino distribution anisotropies, which allows us to explore, for the first time in a multi-energy setting, the interplay between collisional, fast, and slow modes in a moment-based neutron star merger simulation. We find that despite a dominance of the fast instability in most of the simulation volume, certain regions only exhibit a collisional instability, while others, especially at large distances, exhibit a slow instability that is largely underestimated if anisotropic effects are neglected.

Physically motivated Gaussian process (GP) kernels for stellar variability, like the commonly used damped, driven simple harmonic oscillators that model stellar granulation and p-mode oscillations, quantify the instantaneous covariance between any two points. For kernels whose timescales are significantly longer than the typical exposure times, such GP kernels are sufficient. For time series where the exposure time is comparable to the kernel timescale, the observed signal represents an exposure-averaged version of the true underlying signal. This distinction is important in the context of recent data streams from Extreme Precision Radial Velocity (EPRV) spectrographs like fast readout stellar data of asteroseismology targets and solar data to monitor the Sun's variability during daytime observations. Current solar EPRV facilities have significantly different exposure times per-site, owing to the different design choices made. Consequently, each instrument traces different binned versions of the same "latent" signal. Here we present a GP framework that accounts for exposure times by computing integrated forms of the instantaneous kernels typically used. These functions allow one to predict the true latent oscillation signals and the exposure-binned version expected by each instrument. We extend the framework to work for instruments with significant time overlap (i.e., similar longitude) by including relative instrumental drift components that can be predicted and separated from the stellar variability components. We use Sun-as-a-star EPRV datasets as our primary example, but present these approaches in a generalized way for application to any dataset where exposure times are a relevant factor or combining instruments with significant overlap.

We present the largest direct-method abundance catalogue of galaxies to date, containing measurements of 49$\,$959 star-forming galaxies at $z < 0.96$ from DESI Data Release 2. By directly measuring electron temperatures across multiple ionization zones, we provide constraints on a number of electron temperature relations finding good consistency with previous literature relations. Using these temperature measurements, we derive reliable abundances for N, O, Ne, S and Ar and measure the evolution of abundances and abundance ratios of as a function of metallicity and other galaxy properties. Our measurements include direct oxygen abundances for 49$\,$766 galaxies, leading to the discovery of the two most metal-poor galaxies in the nearby Universe, with oxygen abundances of $\rm 12+\log(O/H) = 6.77_{-0.03}^{+0.03}~\rm dex$ (1.2\% $\rm Z_{\odot}$) and $\rm 12+\log(O/H) = 6.81_{-0.04}^{+0.04}~\rm dex$ (1.3\% $\rm Z_{\odot}$). We identify a rare outlier population of 139 galaxies with high N/O ratios at low metallicity, reminiscent of galaxy abundances observed in the early Universe. We find these high N/O galaxies are more massive than typical galaxies at the same metallicity. We find the Ne/O ratio is constant at low metallicity but increases significantly at $\rm 12+log(O/H) > 8.105\pm0.004$ dex. We show that the S/O and Ar/O abundance ratios are strongly correlated, consistent with the expected additional Type Ia enrichment channel for S and Ar. In this work we present an initial survey of the key properties of the sample, with this dataset serving as a foundation for extensive future work on galaxy abundances at low redshift.

Measuring the density profile and mass concentration of dark-matter haloes is a key test of the standard cold dark matter paradigm. Such objects are dark and thus challenging to characterise, but they can be studied via gravitational lensing. Recently, a million-solar-mass object was discovered superposed on an extended and extremely thin gravitational arc. Here we report on extensive tests of various assumptions for the mass density profile and redshift of this object. We find models that best describe the data have two components: an unresolved point-mass of radius $\leq10$ pc centred on an extended mass distribution with an almost constant surface density out to a truncation radius of 139 pc. These properties do not resemble any known astronomical object. However, if the object is dark matter-dominated, its structure is incompatible with cold dark matter models, but may be compatible with a self-interacting dark matter halo where the central region has collapsed to form a black hole. This detection could thus carry substantial implications for our current understanding of dark matter.

In gravitational-wave astronomy, as in other scientific disciplines, ``exceptional'' sources attract considerable interest because they challenge our current understanding of the underlying (astro)physical processes. Crucially, ``exceptionality'' is defined only relative to the rest of the detected population. For instance, among all gravitational-wave events detected so far, GW231123 is the binary black hole with the largest total mass, while GW241110 is the binary black hole with the most strongly misaligned spin relative to the orbital angular momentum. Mandel [this http URL. 996 (2026) 1, L4] argued that apparent ``exceptionality'' may reflect measurement error rather than an extreme true value, and suggested that the total mass of GW231123 may be significantly overestimated. Here we present a quantitative analysis that supports this conceptual point. We find that claims of ``exceptionality'' obtained under agnostic priors should be critically questioned whenever measurement uncertainties are comparable to the width of the underlying population. Specifically, we find that the total mass of GW231123 is unlikely to be meaningfully affected by this effect while the spin of GW241110 is far less likely to be anti-aligned than initially claimed: about 70% of realizations that appear to yield an ``exceptionally anti-aligned'' spin are in fact consistent with either nonspinning or aligned configurations.

The James Webb Space Telescope (JWST) has spectroscopically confirmed galaxies up to $z\sim14$, 300 Myr after the Big Bang, and several candidates have been discovered at $z\sim15-25$, with one candidate as high as $z\sim30$, only 100 Myr after the Big Bang. Such objects are unexpected, since theoretical studies have not predicted the existence of detectable galaxies at $z\sim30$. While any $z\sim30$ candidates may be contaminants at lower redshifts, we explore whether such extreme redshift sources could be consistent with hyper-energetic transient events linked to the formation of the first, metal-free, stars. Specifically, we consider pair-instability supernovae (PISNe), a predicted class of extreme thermonuclear explosions that leave no remnant behind. Using cosmological simulations, we investigate an overdense cosmic region, where star formation and subsequent PISNe occur at $z\sim30-40$, even within standard cosmology. Assessing the likelihood of such a region, the corresponding number of PISNe at $z\gtrsim20$, and their observed flux, we find that JWST has a non-negligible chance to detect a PISN event at extremely high redshifts. If a transient event were confirmed at $z\sim30$, this would provide a direct glimpse into the epoch of first star formation, dramatically extending the empirical reach of astronomy.

A star entering the tidal sphere of a supermassive black hole (SMBH) can be partially stripped of mass, resulting in a partial tidal disruption event (TDE). Here we develop an analytical model for properties of these events, including the peak fallback rate, $\dot{M}_{\rm peak}$, the time at which the peak occurs, $t_{\rm peak}$, and the amount of mass removed from the star, $\Delta M$, for any star and any pericenter distance associated with the stellar orbit about the black hole. We compare the model predictions to 1276 hydrodynamical simulations of partial TDEs of main-sequence stars by a $10^6 M_\odot$ SMBH. The model yields $t_{\rm peak}$ predictions that are in good agreement (to within tens of percent) with the numerical simulations for any stellar mass and age. The agreement for $\dot{M}_{\rm peak}$ is weaker due to the influence of self-gravity on the debris stream dynamics, which remains dynamically important for partial TDEs; the agreement for $\dot{M}_{\rm peak}$ is, however, to within a factor of $\sim 2-3$ in the majority of cases considered, with larger differences for low-mass stars ($M_\star \lesssim 0.5 M_\odot$) on grazing orbits with small mass loss. We show that partial TDE lightcurves for disruptions caused by $\sim 10^6M_\odot$ SMBHs can span $\sim 20-100$ day peak timescales, whereas grazing encounters of high-mass stars with high-mass SMBHs can yield longer peak timescales ($t\gtrsim 1000$ days), associated with some observed transients. Our model provides a significant step toward an analytical prescription for TDE lightcurves and luminosity functions.

Planet formation is inherently linked to protoplanetary disc evolution, which recent developments suggest is driven by magnetised winds rather than turbulent viscosity. We study planet formation in magnetohydrodynamic (MHD) wind-driven discs, assuming accretion occurs in a laminar surface layer above a weakly turbulent midplane. Our goal is to assess the global consequences of recent hydrodynamical results, including inefficient midplane heating and the existence of two Type II migration regimes: slow viscosity-dominated and fast wind-driven migration. We perform single-embryo planetary population syntheses with varying initial disc conditions (i.e. disc mass, size and angular momentum transport), and embryo starting locations, testing different prescriptions for the accretion layer thickness $\Sigma_\text{active}$. Thin ($\lesssim0.01\mathrm{g\,cm^{-2}}$) or fast ($\gtrsim12\%$ sonic velocity) accretion layers result in slow, viscosity-dominated regime which strongly limits the extent of Type II migration. For thick ($\gtrsim1\mathrm{g\,cm^{-2}}$) or slow ($\lesssim3\%$ sonic velocity) accretion layers, fast wind-driven Type II migration occurs frequently, leading to long-range inward migration that sets in once planets reach masses sufficient to block the accreting layer. Disk-limited gas accretion is also strongly affected by deep and early gap opening, limiting maximum giant planet masses. These effects strongly influence the final mass-distance distribution. For thin layers, giant planets form nearly in situ once they have entered Type II migration, which happens already at a few Earth masses, while thick layers lead to numerous migrated Hot Jupiters. Overall, we find that while the global properties of the emerging planet population are strongly modified relative to classical viscous discs, key properties of the observed population can be reproduced within this new paradigm.

We present the largest catalogue to date of Galaxies On Top Of Quasars (GOTOQs), systems where the sightline to a background quasar passes directly through or very close to a foreground galaxy. Using $\approx$1.1 million quasar spectra from SDSS DR16 and the DESI Early Data Release, we identify 1345 unique GOTOQs over the redshift range $0.0045 \leq z \leq 1.09$, more than quadrupling the number previously known. The catalogue combines both absorption agnostic and absorption selected searches, enabling a statistically robust characterization of gas in galaxies at projected separations $D \lesssim 20$ kpc. The GOTOQ emission line ratios indicate that their host galaxies are predominantly normal, star forming disks with typical dust extinctions of $A_V \approx 0.5$ mag. We measure the MgII covering fraction at $D \lesssim 20$ kpc and find it to be remarkably high, $f_c = 98.2^{+1.3}_{-4.5}$ per cent, indicating that these systems trace gas rich, metal enriched regions at the disk halo interface. The median colour excess towards the quasar line of sight, $E(B-V) = 0.087$, is significantly higher than that of typical MgII absorbers, underscoring the dust rich nature of GOTOQ sightlines. From a high signal to noise composite spectrum, we report the first statistical detection of the diffuse interstellar band at $\lambda 4428$ ($W_r = 0.055 \pm 0.011$ angstrom), revealing the presence of complex organic molecules at the disk halo interface. This catalogue provides a powerful reference sample for future multi wavelength studies of gas flows, cold gas, and dust evolution in the inner circumgalactic medium.

Fe K$\alpha$ line emission from Galactic center molecular clouds can be produced either via fluorescence after illumination by an X-ray source or by cosmic ray ionization. Unparalleled high-resolution X-ray spectroscopy obtained by XRISM-Resolve for the galactic center molecular cloud G0.11-0.11 resolves its Fe K$\alpha$ line complex for the first time, and points to a new method for discrimination between the X-ray reflection and cosmic ray ionization models. The Fe K$\alpha$ line complex is resolved into Fe K$\alpha_1$ at $E_{1} = 6.4040 \: \rm{keV}$ and Fe K$\alpha_2$ at $E_{2}= 6.3910 \:\rm{keV}$. Both lines have non-instrumental FWHM of $\approx 3 \:\rm{eV}$, close to the predicted quantum mechanical width of the lines, suggesting scant other sources of line broadening other than instrumental and quantum effects. We measure a radial velocity of $v_{\rm{LSR}} = 50 \pm 12_{fit} \pm 14_{scale} \:\rm{km/s}$ for G0.11-0.11, achieving the same precision reached by radio observations of such clouds. The high-resolution spectrum tests for the presence of secondary Fe K$\alpha$ lines, expected as a signature of cosmic ray proton/ion ionization. The absence of the secondary lines argues against the cosmic ray ionization model for G0.11-0.11. In the preferred X-ray reflection model, if the illuminating source is Sgr A$^{\star}$, the required luminosity for an X-ray outburst about 200 years ago is $L_8 \approx 10^{38} \:\rm{erg/s}$ in an $8\:\rm{keV}$-wide band at $8\:\rm{keV}$.

Context. Cosmic rays and magnetic fields play an important role for the formation and evolution of galaxies. Radio continuum observations allows us to study them in the haloes of edge-on galaxies. Aims. We explore the frequency dependence of the radio scale height which depends on cosmic ray transport and electron cooling. We test the influence of fundamental galaxy properties, such as star-formation rate (SFR), mass and size. Methods. We used radio continuum data of 16 edge-on galaxies from the Continuum HAloes in Nearby Galaxies - an EVLA Survey (CHANG-ES). We included maps from the LOw Frequency ARray at 144 MHz and from the Jansky Very Large Array at 3 GHz with 7" angular resolution. We extracted vertical intensity profiles within the effective radio radius and fitted beam-convolved double- exponential models to separate thin and thick discs. For the thick radio discs, we computed mean spectral indices and scale-height ratios between 144 MHz and 3 GHz. Results. We find a mean scale-height ratio of 1.26 \pm 0.16. This is much lower than what we would expect for either cosmic ray diffusion or advection if synchrotron and inverse Compton losses dominate for the electrons. There is a moderate positive correlation between the ratio and spectral index of the thick disc: galaxies with high ratios have flat radio spectra. The ratio does not depend on any other galaxy parameter. The radio spectrum of the thick disc, as indicated by the radio spectral index, steepens with total mass (strong correlation) and flattens with SFR-to-mass surface density (moderate correlation). Conclusions. Galaxies with galactic winds have flat radio continuum spectra and large scale heights at low frequencies. This shows effective transport of cosmic rays in such systems.

This manuscript investigates the impact of key dust evolution parameters on dust retention and trapping in protoplanetary discs. Using models with and without pressure bumps, combined with radiative transfer simulations, images of the dust continuum emission at (sub-)millimeter wavelengths, their fluxes and observed disc sizes are presented. For discs without pressure bumps (smooth discs), significant dust mass can only be retained over Myr timescales when dust fragmentation velocities are low (1m/s) and with viscosity values of $\alpha=10^{-3}$. For such a combination of fragmentation velocity and viscosity, the synthetic images show a bright inner emission follow by a shallow emission with potential gaps if they are present in the gas profile as well. At higher fragmentation velocities (5-10m/s), most dust is lost due to radial drift at million-year timescales unless pressure traps are present, in which case dust masses can increase by orders of magnitude and structures are observed in synthetic images. The viscosity parameter strongly shapes observable features, with low $\alpha$ producing sharper, potentially asymmetric inner wall cavities in inclined discs due to optically thick emission. High $\alpha$ favors the appearance of shoulders around the predominant rings that dust trapping produces. However, distinguishing between different fragmentation velocities observationally remains challenging. The inferred dust disc sizes from synthetic observations do not always correspond directly to dust model sizes or to the location of pressure bumps. Finally, we discuss implications for pebble fluxes and the delivery of volatiles to the inner disc. These results emphasize the strong degeneracies among dust evolution parameters and highlight the need for multi-wavelength, high-resolution observations to disentangle the processes shaping the formation of planets in protoplanetary discs.

Multiple models have been suggested over the years to explain the structure and support of accretion disks around supermassive black holes, from the standard thin thermal-pressure-dominated $\alpha$-disk model to more recent models that describe geometrically thicker radiation or magnetic or turbulence-dominated disks. In any case, objects embedded in the disk (e.g. compact objects, stars, gas, dust) can undergo gravitational and hydrodynamic interactions with each other leading to interesting processes such as binary interaction/capture, gravitational wave merger events, dynamical friction, accretion, gap opening, etc. It has long been argued that disks of active galactic nuclei (AGN) can enhance the rates for many of these events; however, almost all of that analysis has assumed specific thin-disk models (with aspect ratios $H/R \lesssim 0.01$). We show here that the rates for processes such as these that are mediated by gravitational cross-sections has a very strong inverse dependence on the thickness $H/R$ (scaling as steeply as $(H/R)^{-8}$), and $H/R$ can vary in the outer disk (where these processes are often invoked) by factors $\gtrsim 1000$ depending on the assumed source of pressure support in the disk. This predicts rates that can be lower by tens of orders-of-magnitude in some models, demonstrating that it is critical to account for disk parameters such as aspect ratio and different sources of disk pressure when computing any meaningful predictions for these rates. For instance, if magnetic pressure is important in the outer disk, as suggested in recent work, capture rates would be suppressed by factors $\sim 10^{10}-10^{20}$ compared to previous studies where magnetic pressure was ignored.

The intermediate neutron-capture process (i-process) can occur during proton ingestion events (PIEs), which may take place in the early evolutionary phases of asymptotic giant branch (AGB) stars. We investigate the impact of rotational and magnetic mixing on i-process nucleosynthesis in low-metallicity, low-mass AGB stars. We computed AGB models with [Fe/H] = $-2.5$ and $-1.7$ and initial masses of 1 and 1.5 $M_{\odot}$ using the STAREVOL code, including a network of 1160 nuclei coupled to transport equations. Rotating models incorporate a calibrated Tayler-Spruit (TS) dynamo to account for core rotation rates inferred from asteroseismic observations of solar-metallicity sub-giants and giants. Initial rotation velocities of 0, 30, and 90 km s$^{-1}$ were considered, along with varying assumptions for magnetic mixing. We find that rotation without magnetic fields strongly suppresses the i-process due to the production of primary $^{14}$N, which is subsequently converted into $^{22}$Ne $-$ a potent neutron poison during the PIE. Including magnetic fields via the TS dynamo restores the models close to their non-rotating counterparts: strong core-envelope coupling suppresses shear mixing and prevents primary $^{14}$N synthesis, yielding i-process nucleosynthesis similar to non-rotating models. We also find that rotational mixing during the AGB phase is insufficient to affect the occurrence of PIEs. Proton ingestion event-driven nucleosynthesis proceeds similarly in asteroseismic-calibrated magnetic rotating AGB stars and non-rotating stars, producing identical abundance patterns.

We constrain the rotation and turbulent support of the intracluster medium (ICM) in Abell 2029 (A2029), using dynamical equilibrium models and a combination of state-of-the-art X-ray datasets. We reduce and conduct the spectral analysis of the XRISM/Resolve data. The rotating, turbulent ICM in the model has a composite polytropic distribution in equilibrium in a spherically-symmetric, cosmologically motivated dark halo. The profile of rotation velocity and the distribution of turbulent velocity dispersion are described with flexible functional forms, consistent with the properties of synthetic clusters formed in cosmological simulations. Adopting realistic profiles for the metallicity distribution of the ICM and for the point spread function of XRISM and XMM-Newton, we tune via a Markov chain Monte Carlo algorithm the observables of the intrinsic quantities of the plasma in our model to reproduce the radial profiles of the thermodynamic quantities as derived from the spectral analysis of the XMM-Newton and Planck maps and the measurements of the line-of-sight (LOS) non-thermal velocity dispersion and redshift (probing the LOS velocity) in the XRISM pointings. Our model accurately reproduces the measurements of redshift and LOS non-thermal velocity dispersion, as further demonstrated by simulating and analyzing synthetic counterparts of the XRISM spectra, in accordance with the posterior distribution of our model. We find turbulence-to-total pressure ratio $\approx$ 2% across the (0 - 650) kpc radial range, and a rotation-to-dispersion velocity ratio peaking at 0.15 between 200 - 600 kpc. The hydrostatic-to-total mass ratio is $\approx$ 0.97 at r2500, the radius enclosing an overdensity of 2500 times the average value.

Many Wolf-Rayet (WR) stars have optically thick winds that cloak the hydrostatic layers of the underlying star. In these cases, traditional spectral analysis methods are plagued by degeneracies that make it difficult to constrain parameters such as the stellar radius and the deeper density and velocity structure of the atmosphere. Focussing on the regime of nitrogen-rich WN4-stars with strong emission lines, we employ hydrodynamically-consistent modelling using the PoWR-HD code branch to perform a next generation spectral analysis. The inherent coupling of the stellar and wind parameters enables us to break parameter degeneracies, constrain the wind structure, and get a mass estimate. With this information, we can draw evolutionary implications and test current mass-loss descriptions for WR stars. We selected a sample of six Galactic WN4b stars. Applying updated parallaxes from Gaia DR3 and calculating PoWR-HD models that sufficiently resemble most of their spectral appearance, we obtain new values for the stellar and wind parameters of the WN4b sample. We compare our results to previous studies employing grid models with a beta = 1 velocity structure and cross-check our derived parameters with stellar structure predictions from GENEC and FRANEC evolution tracks. For all six targets, we obtain a narrow range of stellar temperatures T~140 kK, in contrast to previous grid-model analyses. We confirm the existence of WRs with luminosities as low as log L/Lsol = 5.0 and M~5 Msol. All derived velocity fields include a plateau feature at ~85% of the terminal velocity. Both the distance updates and the switch to dynamically-consistent atmospheres lead to substantial parameter adjustments compared to earlier grid-based studies. A comparison of the derived mass-loss rates favours a different description for the WN4b sample than for WN2 stars analysed with the same methodology.

This work investigates variations in the star formation rate during galaxy collisions when the initial conditions of velocity and gas mass are altered. For this purpose, hydrodynamic simulations were performed using the GADGET-4 code, with initial conditions generated by the Galstep and SnapshotJoiner programs. Systems of two galaxies on a head-on collision course were modeled with relative initial velocities ranging from 100 km/s to 1000 km/s, considering two scenarios: the first with identical galaxies, and the second with galaxies of different sizes. In simulations of systems with higher initial relative velocities, both found more intense peaks in the star formation rate, triggered by the first contact of the collision, followed by a strong decline caused by gas dispersion. In contrast, for systems with lower initial velocities, mergers between galaxies were observed, leading to multiple peaks in the star formation rate. A greater initial distance between galaxies has also been linked to whether or not the galaxy system merges, since it implies longer timescales for gravitational action, which leads to higher relative velocities at the moment of collision. Furthermore, the star formation rate in galaxies was found to have a clear dependence on initial gas content. Furthermore, the initial gas content in galaxies was found to have a clear dependence on star formation rates. Overall, our results show that the relative impact velocity, the initial distance between the galaxies, and the gas content are important parameters for analyzing the star formation rate in colliding galaxies.

One of the Milky Way's characteristic features is a strongly bimodal distribution of $\alpha$-process elements, such as Mg, at fixed [Fe/H] in stellar abundances. We examine patterns in [Mg/Fe] versus [Fe/H] in FIRE-2 simulations of Milky Way-mass galaxies. Out of 16 galaxies, 4 are capable of producing a strongly bimodal distribution. In all four galaxies, the high-$\alpha$ population corresponds to an older, radially-compact, thick disk, and the low-$\alpha$ population corresponds to a younger, radially-extended, thin disk, similar to the this http URL transition from high- to low-$\alpha$ took $0.3-1.2\Gyr$ and began $5.5-6.5\Gyr$ ago. [Mg/Fe] decreased at relatively fixed [Fe/H], both in the galaxy overall and at fixed radii: Fe enrichment nearly balanced gas accretion (and therefore dilution), but Mg enrichment was weaker. Importantly, this transition occurred during a period of relatively low gas fraction ($5-15\%$), immediately after a rapid decline in star formation (halving within a few hundred Myr), which caused an increase in Fe-producing white-dwarf supernovae relative to Mg-producing core-collapse supernovae. Only one case coincided with a major merger coalescence. We find similar trends in measuring stars by their current radius and by their birth radius, therefore, radial redistribution did not play a dominant role in the formation of a bimodality or its spatial dependence today. Overall, in FIRE-2, strong $\alpha$-element bimodalities are relatively uncommon ($\sim25\%$), often not associated with a major merger, and arise primarily from a rapid decline in star formation during relatively low gas fraction.

Astronomical measurements are often integrated over finite exposures, which can obscure latent variability on comparable timescales. Correctly accounting for exposure integration with Gaussian Processes (GPs) in such scenarios is essential but computationally challenging: once exposure times vary or overlap across measurements, the covariance matrix forfeits any quasiseparability, forcing O($N^2$) memory and O($N^3$) runtime costs. Linear Gaussian state space models (SSMs) are equivalent to GPs and have well-known O($N$) solutions via the Kalman filter and RTS smoother. In this work, we extend the GP-SSM equivalence to handle integrated measurements while maintaining scalability by augmenting the SSM with an integral state that resets at exposure start times and is observed at exposure end times. This construction yields exactly the same posterior as a fully integrated GP but in O($N$) time on a CPU, and is parallelizable down to O($N/T + \log T$) on a GPU with $T$ parallel workers. We present smolgp (State space Model for O(Linear/log) GPs), an open-source Python/JAX package offering drop-in compatibiltiy with tinygp while supporting both standard and exposure-aware GP modeling. As SSMs provide a framework for representing general GP kernels via their series expansion, smolgp also brings scalable performance to many commonly used covariance kernels in astronomy that lack quasiseparability, such as the quasiperiodic kernel. The substantial performance boosts at large $N$ will enable massive multi-instrument cross-comparisons where exposure overlap is ubiquitous, and unlocks the potential for analyses with more complex models and/or higher dimensional datasets.

The third interstellar object to be discovered, 3I/ATLAS, has a unique and continually unfolding story to tell about its nature and origin as it is monitored by telescopes on Earth, orbiting Earth and around the Solar System. Previous research into missions using chemical propulsion have only really addressed the direct case, where the opportunity to launch already expired before 3I/ATLAS's discovery. In contrast, investigations herein exploit 'Optimum Interplanetary Trajectory Software' to simulate an alternative indirect option for chemical propulsion, namely the Solar Oberth Manoeuvre (SOM). For a SOM, a low perihelion burn provides maximum benefit from the Oberth Effect, and accelerates the spacecraft rapidly towards the receding 3I/ATLAS. Though in principle feasible, results indicate this option presents significant challenges. For possible launch years between 2031 and 2037 inclusive, a 2035 launch permits the most efficient transfer to 3I/ATLAS. The reference mission requires a SOM at 3.2 Solar Radii from the Sun's centre, with an intercept after 35-50 years. It is found the SOM can leverage spacecraft masses up to $\sim{500}$ kg. Two or three solid propellant boosters could deliver the required SOM $\Delta$V, and furthermore a refuelled Starship Block 3 in LEO has sufficient performance for such a mission. As inevitable with a SOM, some of the payload mass would be needed for a heat shield to protect against the high solar flux at low perihelion.

Using VLT/MUSE data, we study the ionised-gas kinematics in a mass- and volume-limited ($M_* \geq 10^{10} M_\odot$, $D \leq 20$\,Mpc) sample of 21 nearby galaxies to identify signatures of extended shocks within their inner kiloparsec, which appear as coherent velocity jumps in kinematic maps. By removing angular momentum, shocks in gas cause inflows, which can trigger nuclear star formation and fuel AGN activity. To identify the signatures of extended shocks, we examine residual velocity fields after subtracting a modelled rotating disc, and we study velocity difference between various gas tracers. Combining our kinematic analysis with BPT ionisation diagnostic maps and dust morphology, we find that 11 of 21 galaxies ($\sim$52%) show extended shock signatures with velocity jumps consistent with models of bar-driven shocks. This is likely a lower limit, as three additional galaxies ($\sim$15%) exhibit shocks along bars, potentially reaching the nucleus but obscured by AGN outflows. We trace shock signatures inwards close to the resolution limit, which suggests that shocks may be the prevailing mechanism of inflow in the central kpc of galaxies. The only two unbarred galaxies in our sample are also the only systems with unperturbed kinematics and no shocks, strongly linking the perturbed gas dynamics in centres of galaxies to the presence of bars. All galaxies with inner bars show LINER- or Seyfert-like nuclear emission, whereas galaxies without inner bars exhibit all emission types, indicating that regardless of gas supply, inner bars suppress star formation in galactic nuclei.

Pebble drift is an important mechanism for supplying the materials needed to build planets in the inner region of protoplanetary disks. Thus, constraining pebble drift's timescales and mass flux is essential to understanding planet formation history. Current pebble drift models suggest pebble fluxes can be constrained from the enhancement of gaseous volatile abundances when icy pebbles sublimate after drifting across key snowlines. In this work, we present ALMA observations of spatially resolved $^{13}$C$^{18}$O J=2-1 line emission inside the midplane CO snowline of the HD 163296 and MWC 480 protoplanetary disks. We use radiative transfer and thermochemical models to constrain the spatial distribution of CO gas column density. We find that both disks display centrally peaked CO abundance enhancement of up to ten times of ISM abundance levels. For HD 163296 and MWC 480, the inferred enhancements require 250-350 and 480-660 Earth Masses of pebbles to have drifted across their CO snowlines, respectively. These ranges fall within cumulative pebble mass flux ranges to grow gas giants in the interior to the CO snowline. The centrally peaked CO enhancement is unexpected in current pebble drift models, which predict CO enhancement peaks at the CO snowline or is uniform inside the snowline. We propose two hypotheses to explain the centrally-peaked CO enhancement, including a large CO desorption distance and CO trapped in water ice. By testing both hypotheses with the 1D gas and dust evolution code chemcomp, we find that volatile trapping (about 30\%) best reproduces the centrally peaked CO enhancement observed.

The rise of oxygen in the Earth's atmosphere during the Great Oxidation Event (GOE) occurred about 2.3 billion years ago. There is considerably greater uncertainty for the origin of oxygenic photosynthesis, but it likely occurred significantly earlier, perhaps by 700 million years. Assuming this time lag is proportional to the rate of oxygen generation, we can estimate how long it would take for a GOE-like event to occur on a hypothetical Earth-analog planet orbiting the star TRAPPIST-1 (a late M star with Teff 2560 K). Although in the habitable zone, an Earth-analog planet located in TRAPPIST-1e's orbit would receive only 0.9% of the Photosynthetically Active Radiation (PAR) that the Earth gets from the Sun. This is because most of the star's light is emitted at wavelengths longer than the 400-700 nm PAR range. Thus it would take 63 Gyrs for a GOE to occur. But the linear assumption is problematic; as light levels increase, photosynthesis saturates then declines, an effect known as photoinhibition. Photoinhibition varies from species to species and depends on a host of environmental factors. There is also sensitivity to the upper wavelength limit of the PAR: extending just 50 nm increases the number of photons by a factor of 2.5. Including these and other factors greatly reduces the timescale to roughly 1-5 Gyrs for a GOE. However, non-oxygenic photosynthetic bacteria can thrive in low-light environments and can use near-IR light out to 1100 nm, providing 22 times as many photons. With this huge light advantage, and because they evolved earlier, anoxygenic photosynthesizers would likely dominate the ecosystem. On a late M-star Earth-analog planet, oxygen may never reach significant levels in the atmosphere and a GOE may never occur, let alone a Cambrian Explosion. Thus complex animal life is unlikely.

We have used 23 years of Hubble Space Telescope ACS/SBC data to study what background levels are encountered in practice and how much they vary. The backgrounds vary considerably, with F115LP, F122M, F125LP, PR110L, and PR130L all showing over an order of magnitude of variation in background between observations, apparently due to changes in airglow. The F150LP and F165LP filters, which are dominated by dark rate, not airglow, exhibit a far smaller variation in backgrounds. For the filters where the background is generally dominated by airglow, the backgrounds measured from the data are significantly lower than what the ETC predicts (as of ETC v33.2). The ETC predictions for `average' airglow are greater than the median of our measured background values by factors of 2.51, 2.64, 105, and 3.64, for F115LP, F122M, F125LP, and F140LP, respectively. A preliminary analysis suggests this could be due to certain OI airglow lines usually being fainter than expected by the ETC. With reduced reduced background levels, the shorter-wavelength SBC filters can conduct background-limited observations much more rapidly than had previously been expected. As of ETC v34.1, a new option will be included for SBC calculations, allowing users to employ empirical background percentiles to estimate required exposure times.

The Lazuli Space Observatory is a 3-meter aperture astronomical facility designed for rapid-response observations and precision astrophysics across visible to near-infrared wavelengths (400-1700 nm bandpass). An off-axis, freeform telescope delivers diffraction-limited image quality (Strehl $>$0.8 at 633 nm) to three instruments across a wide, flat focal plane. The three instruments provide complementary capabilities: a Wide-field Context Camera (WCC) delivers multi-band imaging over a 35' $\times$ 12' footprint with high-cadence photometry; an Integral Field Spectrograph (IFS) provides continuous 400-1700 nm spectroscopy at R $\sim$ 100-500 for stable spectrophotometry; and an ExtraSolar Coronagraph (ESC) enables high-contrast imaging expected to reach raw contrasts of $10^{-8}$ and post-processed contrasts approaching $10^{-9}$. Operating from a 3:1 lunar-resonant orbit, Lazuli will respond to targets of opportunity in under four hours--a programmatic requirement designed to enable routine temporal responsiveness that is unprecedented for a space telescope of this size. Lazuli's technical capabilities are shaped around three broad science areas: (1) time-domain and multi-messenger astronomy, (2) stars and planets, and (3) cosmology. These capabilities enable a potent mix of science spanning gravitational wave counterpart characterization, fast-evolving transients, Type Ia supernova cosmology, high-contrast exoplanet imaging, and spectroscopy of exoplanet atmospheres. While these areas guide the observatory design, Lazuli is conceived as a general-purpose facility capable of supporting a wide range of astrophysical investigations, with open time for the global community. We describe the observatory architecture and capabilities in the preliminary design phase, with science operations anticipated following a rapid development cycle from concept to launch.

We present a proposal for a nanomechanical membrane resonator integrated into a moderate-finesse ($\mathcal{F}\sim 10$) optical cavity as a versatile platform for detecting high-frequency gravitational waves and vector dark matter. Gravitational-wave sensitivity arises from cavity-length modulation, which resonantly drives membrane motion via the radiation-pressure force. This force also enables in situ tuning of the membrane's resonance frequency by nearly a factor of two, allowing a frequency coverage from 0.5 to 40 kHz using six membranes. The detector achieves a peak strain sensitivity of $2\times 10^{-23}/\sqrt{\text{Hz}}$ at 40 kHz. Using a silicon membrane positioned near a gallium-arsenide input mirror additionally provides sensitivity to vector dark matter via differential acceleration from their differing atomic-to-mass number ratios. The projected reach surpasses the existing limits in the range of $2\times 10^{-12}$ to $2\times 10^{-10}$ $\text{eV}/c^2$ for a one-year measurement. Consequently, the proposed detector offers a unified approach to searching for physics beyond the Standard Model, probing both high-frequency gravitational waves and vector dark matter.

Polar-structure galaxies (PSGs) host photometrically and kinematically decoupled components oriented at large angles to one another. These systems, which include polar rings, polar disks, polar halos, polar bulges, polar dust lanes, and polar tidal structures, provide valuable insights into galaxy formation and evolution, although their rarity has limited statistical studies. We aim to construct the largest and most homogeneous catalog of PSGs to date in order to enable robust statistical studies of their properties and occurrence rates in the local Universe. Using DESI Legacy Imaging Surveys (DR10) data, we identified PSG candidates in the Siena Galaxy Atlas (SGA) through visual inspection, convolutional neural network classification, and cross-matching with previously reported systems. Each galaxy was assigned a PSG subtype and host morphology. We analyzed general properties of PSGs and compared them with those for all galaxies in the SGA. Simple image simulations were used to evaluate projection effects. The resulting Catalog of Unusual Galaxies with polar Structures in the DESI Legacy Imaging Surveys (COUGS-DESI) contains 2,989 PSG candidates, including 342 previously known objects. Projection effects from random galaxy overlaps are negligible. The sample spans a wide range of polar morphologies, with 1,113 polar rings, 75 polar bulges, 216 polar halos, 185 polar dust lanes, and 1,315 polar tidal structures. PSGs constitute 2.2% of local non-dwarf galaxies, with polar rings representing 0.7%. Approximately 1% of S0 galaxies in the SGA host polar rings, whereas spirals constitute the most common morphological type among the PSG hosts in our catalog. COUGS--DESI increases the number of known PSG candidates by an order of magnitude and provides a foundation for detailed studies of the formation and evolution of polar structures.

Isotopic abundance ratios in protoplanetary disks are critical for understanding volatile inheritance and chemical evolution in planet-forming environments. We present Atacama Large Millimeter/submillimeter Array observations of the rare isotopologue 13C18O(2-1) at approximately 0.3 arcsec resolution from the disk around the Herbig Ae star HD 163296, combined with archival observations of C17O(2-1), C18O(1-0), and C17O(1-0), to empirically constrain carbon and oxygen isotopic ratios without detailed disk modeling. Both the C17O/13C18O(2-1) and C18O/C17O(1-0) flux ratios rise sharply across the CO snowline and flatten beyond 1.5 arcsec (r >= 150 au), where the emission becomes optically thin. This transition, reflecting a steep drop in CO column density set by the disk's thermal structure, makes HD 163296 an optimal case for isotopic analysis. Using beam-averaged intensities of the four transitions measured in this optically thin region, we derive isotopic ratios of 12C/13C = 75.3 (+14.7/-11.4) and 18O/17O = 3.28 (+0.31/-0.26), both consistent with local interstellar medium values. The 16O/18O ratio remains weakly constrained due to moderate optical depth in the C18O(1-0) line and degeneracy with CO column density. These results demonstrate that rare CO isotopologues can provide robust, empirical constraints on isotopic ratios in disks when sharp structural transitions allow for the identification of optically thin regions, and establish HD 163296 as a benchmark for extending such studies to other systems with resolved snowline structures.

Recent observations of small bodies of the Solar System showed evidence of the presence of refractory (asteroidal) material in the Oort cloud. Different models of the origin of the Solar System predict different numbers of rocky objects in the Oort cloud, meaning that measurement of this population can be used as an observational constraint for cosmogonic models. The aim of our work is to study how the data obtained from meteor observations can be used as a tool for distinguishing among the existing cosmogonic models. We investigated two meteor databases collected by the cameras of the All-Sky Meteor Orbit System (AMOS) located in the Canary Islands and in Chile. We describe methodology and results of the search for unusually strong rocky meteoroids on cometary orbits with the origin in the Oort cloud. These data will be used to calculate the fluxes of meteors of different compositions in order to constrain the ratio of icy and rocky components of the Oort cloud. For the flux determination, we estimate the observational time and effective area of the AMOS system.

Binary neutron star merger (BNSM) ejecta are considered a primary repository of $r$-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. We find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary $r$-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the $r$-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Stellar photospheric heterogeneities (e.g., starspots, faculae) distort the stellar spectrum in transit and imprint wavelength-dependent biases on the planet-to-star radius ratio (Transit Light Source Effect, TLSE). The Rackham-TLSE (R-TLSE) prescription applies a disc-averaged correction based solely on filling factor and spectral contrast, but transmission spectroscopy also depends on limb darkening, active-region distribution, and transit geometry. We include these in a pixel-resolved framework, ECLIPSE-Xlambda, and run idealised noise-free model-model comparisons to R-TLSE. For LHS 1140 b, K2-18 b, and WASP-69 b, disc-averaged corrections differ from the pixel model by up to about 400 ppm in the optical for active hosts and non-equatorial transits, but stay below about 10 ppm in the near-infrared where limb darkening is weak. We then apply both approaches to the JWST/NIRISS SOSS spectrum of LHS 1140 b. With limb darkening set to zero, ECLIPSE-Xlambda recovers stellar-contamination parameters matching the reference R-TLSE solution, confirming consistency in the disc-averaged limit. With wavelength-dependent limb darkening, reproducing the short-wavelength slope via stellar contamination alone requires hot faculae (delta Tfac about 600 K; ffac about 0.35), equivalent to a circular facular region of radius about 0.6 Rstar (about 60% of the stellar radius) on the disc; such an extended unocculted region is physically unlikely even for an active M dwarf. Purely stellar contamination would therefore require extreme faculae, whereas a genuine atmospheric contribution complementing a more modest facular signal is more plausible. These results delineate the validity regime of R-TLSE and underscore the need for geometry-aware stellar-heterogeneity models including limb darkening in high-precision transmission spectroscopy.

We investigate the evolution of the ionization state of metal-enriched gas in and around galaxies near the epoch of reionization using a sample of 488 metal absorption systems at 4.3<z<6.3 from the E-XQR-30 survey. We classify the absorption systems based on whether they display only low-ionization absorption (CII, SiII, MgII), only high-ionization absorption (CIV, SiIV), or both. The percentage of low-ionization-only systems decreases from 24% at $z\sim$6 to 2% at $z\sim$4.3, whilst the fraction of high-ionization-only systems increases from 52% to 82%. For mixed absorbers (with both low and high ionization absorption), we use the column density ratios log(N_CII/N_CIV) and log(N_SiII/N_SiIV) to quantify the average ionization as a function of redshift. The log(N_SiII/N_SiIV) ratio does not change significantly over 5$\lesssim z \lesssim$6.3. We combine the E-XQR-30 log(N_CII/N_CIV) measurements with literature measurements at $z\sim$2-4 and find that the log(N_CII/N_CIV) ratio declines by a factor of $\sim$20 between $z\sim$6 and $z\sim$2. To explore possible drivers of this evolution, we run photoionization models of gas slabs illuminated by a uniform UV background at fixed density, metallicity and HI column density. We find that the increase in the ionization state of metal absorbers towards lower redshifts can likely be explained by some combination of 1) an increase in the metallicity of CIV-absorbing gas and 2) a decrease in the typical HI column densities of the absorbing gas, driven by the declining cosmic mean density and a rapid rise in the strength of the UV background during the final stages of reionization.

Though missions devoted to time series photometry focus primarily on targets far beyond the solar system, their observations can be contaminated by foreground minor planets, especially near the ecliptic plane where solar system objects are most prevalent. Crucially, depending on one's choice of data reduction/background estimation algorithm, these objects can induce both apparent brightening and/or dimming events in processed light curves. To quantify the impact of these objects on archived TESS light curves, we used N-body integrations of all currently known minor planets to postdict all 600,000+ of their interactions with stars selected for high-cadence observations during TESS ecliptic plane sectors. We then created mock images of these moving sources and performed simple aperture photometry using the same target and background apertures used in SPOC processing. Our resulting 10,000+ target-specific light curves, which faithfully model the time-dependent positions and magnitudes of the actual solar system objects that approached each target, reveal that $>95\%$ of high-cadence ecliptic plane targets experience a minor planet crossing within 1 TESS pixel of the source. Additionally, 50% of all $T>13$ mag targets experience at least one instantaneous moment where the contaminating flux from minor planets exceeds 1% of the target flux. We discuss these population-level results and others, and highlight several case studies of bright flybys.

We present photometric and spectroscopic observations of SN 2024abfl, a low-luminosity Type IIP supernova (LLSN) discovered shortly after explosion. The transient reached a peak absolute magnitude of $M_V = -14.9$ and exhibited an extended, flat plateau lasting $\sim$125 days. From the late-time bolometric light curve, we estimate a $^{56}$Ni mass of $\sim0.01~M_\odot$, consistent with other LLSNe. Analytical shock-cooling models fail to reproduce the rapid early rise, indicating that circumstellar material (CSM) interaction contributed to the initial emission. The spectroscopic evolution is typical of LLSNe, with relatively narrow metal lines and low expansion velocities ($\lesssim 3000$ km s$^{-1}$) that decline slowly over time. We detect a broad ``ledge'' feature around 4600 Åwithin three days of explosion, which we interpret as a blend of high-ionization shock-accelerated CSM lines. Multi-peaked H$\alpha$ profiles develop during the plateau phase, consistent with complex ejecta-CSM interactions. As one of the best-observed examples of LLSNe, SN 2024abfl exhibits a weak explosion and signatures of nearby CSM, offering new insights into progenitor properties, pre-explosion mass loss, and the diversity of LLSNe.

Subhalo spin is essential for modeling galaxy formation and controlling systematic uncertainties in intrinsic alignment (IA) studies. However, the physical mechanisms governing subhalo spin acquisition within the tidal environments of host halos remain poorly understood. In this work, we investigate the alignment between subhalo and host halo spins using the high-resolution cosmological $N$-body simulation, Shin-Uchuu. We find that the spin alignment between subhalos and host halos becomes increasingly pronounced toward the central regions. Our analysis reveals that subhalos typically acquire spin in the same direction as their orbital angular momentum. Since the orbital angular momentum of most subhalos is aligned with the host halo spin, an overall alignment between subhalo and host spins emerges. When classified by orbital orientation, however, subhalo spins in the inner regions are found to be oriented perpendicularly or anti-parallel to the host spin for polar and retrograde orbits, respectively. These results provide strong evidence that subhalo spins are acquired through torques exerted by the tidal field of the host halo. Furthermore, we demonstrate that the mass ratio and the radial distance from the host center are the primary parameters governing subhalo spin alignment, while the dependence on the accretion redshift is found to be negligible.

Experiments designed to detect ultra-high energy (UHE) neutrinos using radio techniques are also capable of detecting the radio signals from cosmic-ray (CR) induced air showers. These CR signals are important both as a background and as a tool for calibrating the detector. The Askaryan Radio Array (ARA), a radio detector array, is designed to detect UHE neutrinos. The array currently comprises five independent stations, each instrumented with antennas deployed at depths of up to 200 meters within the ice at the South Pole. In this study, we focus on a candidate event recorded by ARA Station 2 (A2) that shows features consistent with a downward-going CR-induced air shower. This includes distinctive double-pulse signals in multiple channels, interpreted as geomagnetic and Askaryan radio emissions arriving at the antennas in sequence. To investigate this event, we use detailed simulations that combine a modern ice-impacting CR shower simulation framework, FAERIE, with a realistic detector simulation package, AraSim. We will present results for an optimization of the event topology, identified through simulated CR showers, comparing the vertex reconstruction of both the geomagnetic and Askaryan signals of the event, as well as the observed time delays between the two signals in each antenna.

Cosmic dust plays a vital role in stellar and galactic formation and evolution, but its three-dimensional structure in the Milky Way has remained unclear due to insufficient precise reddening and distance measurements. Although early studies typically adopted a single-disk model, we detect two distinct components at Galactocentric distances of 5-14 kpc, enabled by photometric, spectroscopic, and astrometric measurements of over 5 million stars. The thin dust disk's scale height increases radially from 60 to 200 pc, while the thick disk grows from 300 to 800 pc. For the first time, we find the thin and thick dust disk correlates spatially with molecular and atomic hydrogen disk, respectively. The thin, thick, and combined disks have scale lengths of 9.6+1.2-1.1 kpc, 4.2+0.4-0.3 kpc, and 6.6+0.3-0.3 kpc, respectively. The gas-to-dust ratio shows an exponential radial gradient, increasing from around 60 at 5 kpc to around 470 at 14 kpc. These findings provide new insights into dust morphology in the Galaxy and raise fundamental questions that require further investigation.

Recent observations of FRB 20220529 reveal significant variation and a partial reversal in its rotation measure (RM), suggesting the presence of a dynamically evolving magnetized environment, which could be caused by the orbital motion of the magnetar within the binary system. Here we develop the binary model by suggesting that the spin and magnetic axis of the companion star could undergo precession around the orbital axis. It is then investigated how the precession period and the inclination of the magnetic axis, as well as a possible disc wind, can influence the evolution behaviors of the RM and dispersion measure (DM) of FRB emission. As the foremost consequence, the RM variation can be significantly altered on timescales longer than the orbital period, producing super-orbital evolution and complex patterns. Applying this model to FRB 20220529, we find that its RM evolution could be reproduced with a precession period of 182 days and an inclination angle of approximately $19^{\circ}$, while the other binary parameters are fixed at their typical values. Meanwhile, the absence of significant variation of the DM argues against the presence of a dense equatorial disc around the companion star, which would be constrained by future long-term observations.

We present high-resolution near-infrared spectroscopic observations of the newly identified supernova remnant (SNR) candidate G25.8+0.2 obtained with the Immersion Grating Infrared Spectrograph (IGRINS) on the Harlan J. Smith Telescope. The source was originally discovered in the UKIRT Wide-field Infrared Survey for Fe+ (UWIFE; Lee et al. 2014). Our spectra reveal multiple kinematic components in the [Fe II] 1.644 um emission. The high-velocity components exhibit elevated [Fe II]/Br gamma ratios characteristic of shock excitation, while the low-velocity components are dominated by hydrogen and helium recombination lines and are consistent with photoionized gas, indicating an H II-region origin. G25.8+0.2 lies within the G26 complex, a large (~15'x 30', corresponding to ~28 pc x 57 pc at a distance of 6.5 kpc) star-forming region in the inner Galaxy. The shock-excited [Fe II] filaments closely trace the morphology of the bright radio shell that partially encloses centrally filled soft X-ray emission, strongly suggesting recent supernova activity in this region. We discuss the physical nature of G25.8+0.2 and its relationship to the surrounding G26 star-forming complex. In addition, we derive the extinction toward the source using Brackett-line ratios and constrain the spectral types of the dominant ionizing stars from the He I 2.058 um/Br gamma line ratios.

Three-dimensional (3D) magnetic field in the solar atmosphere provides crucial information to understand the explosive phenomenon such as solar flares and coronal mass ejections. Since it is still hard that we determine the 3D magnetic field from direct observation, a nonlinear force-free field (NLFFF) extrapolation is one of the best modeling methods that provides 3D magnetic field. However, the method is based on zero-beta assumption, i.e., the model ignores the gas pressure gradient and gravitational force. The magnetic field based on an NLFFF is not well reconstructed in high-beta region, such as in chromospheric or lower height layer and in weak field region. To overcome this problem, we need to consider the magnetohydrostatic equilibrium. In this study, we developed a finite plasma beta magnetic field extrapolation method based on magnetohydrodynamic relaxation method. In our method, we consider a force balance of the Lorentz force and the gas pressure. We tested three different schemes and extrapolated 3D magnetic field using an observational photospheric vector magnetic field of one solar active region, which is a quadrupole complex sunspot group and well studied with an NLFFF. The verification of three schemes is performed by comparing the residual force, and we concluded that our methods reduce ~4% of residual force of the previous NLFFF. We also examined the plasma beta profile along the height and found that, in the core of the active region, plasma beta reaches a local minimum of ~0.01 in the lower corona with beta ~1 at the photosphere.

We explore the transition from the $\Lambda$CDM to an interacting dark sector by introducing a model with a redshift threshold that controls the onset of the energy transfer between the dark energy and the dark matter. Below the transition redshift, the interaction between dark matter and dark energy becomes active, while at earlier times the cosmological evolution coincides with that of $\Lambda$CDM. This approach allows us to determine the epoch in the comic history where the interacting effects have an impact in the description of the dark sector. We constrain the free parameters of the model using late-time cosmological observations, namely Cosmic Chronometers, DESI DR2 Baryonic Acoustic Oscillations, and Supernova data from the Pantheon Plus, Union3.0, and DES-Dovekie catalogues. The analysis provides an indication of a strong interacting term that describes energy transfer from dark energy to dark matter, which is activated at low redshifts. The PantheonPlus sample provides a threshold of $z_{T}<0.624$, the Union3.0 sample yields $z_{T}=0.400_{-0.23}^{+0.021}$, and the DES-Dovekie sample gives $z_{T}=0.371_{-0.26}^{+0.028}$. The model fits the data in a similar way to the CPL parametrization, without the dark energy to cross the phantom divide line.

We present a new fiber assignment algorithm for a robotic fiber positioner system in multi-object spectroscopy. Modern fiber positioner systems typically have overlapping patrol regions, resulting in the number of observable targets being highly dependent on the fiber assignment scheme. To maximize observable targets without fiber collisions, the algorithm proceeds in three steps. First, it assigns the maximum number of targets for a given field of view without considering any collisions between fiber positioners. Then, the fibers in collision are grouped, and the algorithm finds the optimal solution resolving the collision problem within each group. We compare the results from this new algorithm with those from a simple algorithm that assigns targets in descending order of their rank by considering collisions. As a result, we could increase the overall completeness of target assignments by 10% with this new algorithm in comparison with the case using the simple algorithm in a field with 150 fibers. Our new algorithm is designed for the All-sky SPECtroscopic survey of nearby galaxies (A-SPEC) based on the K-SPEC spectrograph system, but can also be applied to similar fiber-based systems with heavily overlapping fiber positioners.

The 1.3 mm ground-based very long baseline interferometry (VLBI) array, the Event Horizon Telescope (EHT), is limited by the Earth's diameter and can image the supermassive black hole (SMBH) shadows of only M87* and Sgr A*. Extending the array with an assumed lunar-based telescope could achieve $\sim 0.85\ \mu$as angular resolution at 230 GHz, enabling black hole shadow detection for a larger SMBH sample. The concept is motivated by space VLBI missions and lunar exploration, including the ongoing Lunar Orbit VLBI Experiment (LOVEX) aboard QueQiao-2 (Chang'E-7) and the planned International Lunar Research Station (ILRS). We assess shadow detectability for 31 SMBH with predicted large angular sizes, exploring different telescope location and antenna size. Assuming a telescope at the lunar antipode, we simulate the Moon-Earth (u,v) coverage and show that source geometry relative to the Moon's orbit determines whether the primary indicator of shadow, first visibility null, can be sampled. Using a geometric ring model, we identify six high-priority targets: M104, NGC 524, PGC 049940, NGC 5077, NGC 5252, and NGC 1052. Shadows of M104, NGC 5077, and NGC 1052 are detectable with a 5 m lunar-based telescope; PGC 049940 requires 20 m; NGC 524 and NGC 5252 require 100 m. Photon ring detection for Sgr A*, M87*, NGC 1600, and M31 is possible if space telescopes fill the baseline coverage gaps and sensitivity requirements are met. These results provide a clear scientific and technical motivation for lunar-based telescopes in future black hole shadow studies.

We present a numerical investigation of nonlinear cluster lens reconstruction using weak lensing mass mapping. Recent advances in imaging and shear estimation have pushed reliable reduced shear measurements closer to cluster cores, making mass reconstruction accessible in the nonlinear regime. However, the Kaiser-Squires based algorithm becomes unstable in cluster cores, where convergence $\kappa$ significantly deviates from zero and the linear approximation breaks down. To address this limitation, we develop a reconstruction framework with two key modifications: applying smooth masks to these regions and using a model-derived analytical solution as the initial guess, rather than assuming $\kappa = 0$. We validate our framework using simulated cluster lensing data with known mass distributions, incorporating realistic masks that arise from limitations in reduced shear measurements. We show that in the absence of shape noise, our framework yields high-fidelity mass reconstruction in regions of large reduced shear, with the best-performing method achieving residuals below $0.02 \sigma$ in the unmasked regions. This pushes mass reconstruction to higher accuracy in the nonlinear regime.

Procyon A is an F-type main-sequence star in a binary system. It has been the subject of numerous ground-based and space-based observing campaigns, providing precise classical constraints, including a well-determined mass. It was also among the first stars in which individual frequencies were detected, making it a crucial benchmark for F-type stars. Our goal is to investigate the surface effect, namely the discrepancy between observed and model oscillation frequencies due to inadequate modeling of the surface stellar layers, especially important in F-type stars. Using Procyon A as a case study, we aim to understand how different surface correction prescriptions impact the inference of the fundamental properties of this star, and compare the results with those obtained when the surface corrections are neglected. We inferred the fundamental stellar properties employing a grid of models computed with MESA, including gravitational settling, radiative accelerations, and turbulent mixing. We selected the best-fit models using the AIMS code taking into account different methods to fit the individual frequencies. We find that the use of surface corrections can introduce uncertainties up to 7\% in the inferred stellar mass. We identify that the most reliable stellar mass estimates are obtained when using frequency ratios, the Sonoi et al. (2015) surface correction or directly fitting the individual frequencies. Our results indicate that the surface effects in F-type stars differ from those found in the Sun and in solar-like stars, highlighting the need to be careful when considering the surface corrections for these stars.

We propose a data-driven approach to reconstruct the all-sky distribution of the dispersion measure contribution from the Galactic halo ($\mathrm{DM_{halo}}$) through a spherical harmonic expansion, enabling an investigation of its possible anisotropies. Based on the NE2001 model and using 92 localized and 574 unlocalized non-repeating fast radio bursts (FRBs) at Galactic latitudes $|b|>15^\circ$, we find a significant dipole anisotropy in $\mathrm{DM_{halo}}$, pointing toward $(l=130^\circ,\, b=+5^\circ)$ with a $1\sigma$ uncertainty of approximately $28^\circ$. The $\mathrm{DM_{halo}}$ value in this direction is $63\pm9~\mathrm{pc~cm^{-3}}$, exceeding the all-sky mean by about $2.6\sigma$. This result is not significantly affected by the choice of Galactic ISM models. Furthermore, even when using a refined sample of 62 localized FRBs (excluding CHIME detections, repeaters, and unlocalized events), the dipole anisotropic structure persists, with a direction of $(l=141^\circ,\, b=+51^\circ)$ and a larger 1$\sigma$ uncertainty of $\sim 44^\circ$. Model comparisons using the Akaike Information Criterion and Bayesian evidence yield consistent preferences, and together they suggest that current FRB data slightly favor the existence of a dipole structure in $\mathrm{DM_{halo}}$. If this feature is not a statistical fluctuation or systematic error, its physical origin requires further investigation. Future FRB samples with larger sizes and more complete sky coverage will be essential to confirm or refute this possible anisotropic structure.

We present photometric properties of 161 galaxy candidates at $z=7-14$ selected from the second data release (DR2) of BEACON, a JWST Cycle 2 pure-parallel NIRCam imaging program. Carefully selected from 36 independent pointings (corresponding to $\sim350$\,arcmin$^2$ sky coverage), and hence with reduced cosmic variance, our galaxy candidates provide an unbiased sample for investigating galaxy properties over a wide range of environments. We measure the physical properties, including UV continuum slope ($\beta_{\rm UV}$), stellar mass ($M_*$), star formation rate (SFR), and sizes. Our highest redshift galaxy candidate at $z=13.71\pm0.15$ has a remarkably bright UV luminosity of $M_{\rm UV}=-21.19\pm0.08$, making it the brightest galaxy at $z>12$ if spectroscopically confirmed. With an extremely blue UV slope, compact morphology, and high star formation rate surface density ($\Sigma_{\rm SFR}$), this candidate may have extremely low metallicity, high ionizing photon escape fraction, or contributions from an AGN. Among our multiple independent sightlines, we identify three fields of galaxy number overdensity with $>3\sigma$ significance. The properties of galaxies in various environments do not exhibit significant differences, implying either that accelerated galaxy evolution in overdense regions is not yet widespread at $z>7$, or that the current constraints are limited by sample size. Our simulations indicate that increasing the sample by an order of magnitude would allow such environmental trends to be robustly confirmed or ruled out, underscoring the importance of future pure-parallel observations.

We present JWST/MIRI imaging of eight nearby Active Galactic Nuclei (AGN) from the GATOS survey to investigate the physical conditions of extended dust in their narrow line regions (NLRs). In four galaxies (ESO 428-G14, NGC 4388, NGC 3081, and NGC 5728), we detect spatially resolved dust structures extending ~100-200 pc along the NLR. In these systems, we find a strong link between the morphology of the dust, the radio ejecta, and the coronal [Si VI] emission, implying that dust carries imprints of the processes shaping the NLR. Using spatially resolved spectral energy distributions, we show that dust in the NLR has systematically steeper slopes than star forming clumps. This dust emits at temperatures in the range 150 - 220 K, at a distance of ~150 pc from the nucleus. Using simple models, we show that, even under optimistic assumptions of grain size and AGN luminosity, the excess MIR emission cannot be explained by AGN illumination alone. We interpret this excess heating as in-situ. We show that shocks with velocities of $v_{\rm shock} \sim 200- 400 \, \rm km/s$ in dense gas can close this gap, and in some cases even account for the total observed emission. This, combined with multiple lines of evidence for shocks in these regions, supports a scenario in which shocks not only coexist with dust but may be playing a key role in heating it. Our findings reveal shocks may be an important and previously overlooked driver of extended dust emission in the central hundreds of parsecs in AGN.

We propose an innovative definition of the scale at which the Universe becomes homogeneous based on measurements of velocities rather than densities. When using the matter density field, one has to choose an arbitrary scale (e.g. within 1\% of the average density) to define the transition to homogeneity. Furthermore, the resulting homogeneity scale is strongly degenerate with the galaxy bias. By contrast, peculiar velocities (PV) allow us to define an unambiguous scale of homogeneity, namely the distance at which the velocities between pairs of galaxies change from being on-average correlated to anti-correlated. Physically, this relates to when the motion of pairs of galaxies is influenced by the matter density between them, rather than beyond. The disadvantage is that peculiar velocities are more difficult to measure than positions, resulting in smaller samples with larger uncertainties. Nevertheless, we illustrate the potential of this approach using the peculiar velocity correlation functions obtained from the Sloan Digital Sky Survey PV catalog, finding an homogeneity scale of $R_H\approx 133\substack{+28 \\ -52}\, \rm{Mpc/h}$. Finally, we show that more precise measurements are within reach of upcoming peculiar velocity surveys, and highlight this homogeneity scale's potential use as a standard ruler within the standard cosmological model.

Detecting and characterizing individual ionized bubbles during the Epoch of Reionization (EoR) using the redshifted HI 21-cm signal provides a direct probe of the early ionizing sources and the intergalactic medium. We develop and validate a computationally efficient estimator that operates on gridded visibilities to detect ionized bubbles. This serves as an accurate alternative to the more computationally demanding bare estimator that uses all baselines and frequency channels. Further, we employ a non-parametric foreground-subtraction method based on Gaussian process regression, which minimizes loss of the HI 21-cm signal and yields improved signal-to-noise ratios. Our analysis indicates that ionized bubbles at redshifts $z \sim 7 - 8$ can be detected with SNR $\gtrsim 10$ using $\sim 100$ hours of SKA1-Low AA$^*$ and AA4 observations. We further derive a scaling relation that connects the SNR to the bubble radius, redshift, total observing time, and the mean neutral hydrogen fraction of the surrounding IGM. This helps to quickly predict the observational outcome for any planned observations and is, therefore, useful for devising observational strategies. Finally, we apply a Bayesian likelihood framework with Markov Chain Monte Carlo sampling to the residual visibilities to recover ionized bubble properties, including radius, position, and the mean neutral fraction. The resulting posterior distributions demonstrate accurate recovery of the bubble parameters. This confirms the feasibility of robustly characterizing individual ionized regions with the SKA1-Low.

Accurately measuring stellar ages and internal structures is challenging, but the inclusion of asteroseismic observables can substantially improve precision. However, the curse of dimensionality means this comes at a high computational cost when using standard interpolation methods across grids of stellar models. Furthermore, without a rigorous treatment of random uncertainties in grid-based modelling, it is not possible to address systematic errors in stellar models. We present PITCHFORK -- a multilayer perceptron neural network with a branching architecture capable of rapid emulation of both classical stellar observables and individual asteroseismic oscillation modes of solar-like oscillators. PITCHFORK can predict the classical observables $T_{\text{eff}}$, $L$, and $\left[\mathrm{Fe}/\mathrm{H}\right]$ with precisions of $5.88\,\text{K}$, $0.014\,\text{L}_{\odot}$, and $0.001\,\text{dex}$, respectively, and can predict 35 individual radial mode frequencies with a uniform precision of $0.02$ per cent. PITCHFORK is coupled to a vectorised Bayesian inference pipeline to return well-sampled and fully marginalised posterior distributions. We validate our rigorous treatment of the random uncertainties -- including the asteroseismic surface effect -- in an extensive hare-and-hounds exercise. We also demonstrate our ability to infer the stellar properties of benchmark stars -- namely, the Sun and the binary stars 16 Cygni A and B. This work demonstrates a computationally scalable and statistically robust framework for stellar parameter inference of solar-like oscillators using individual asteroseismic mode frequencies. This provides a foundation for the treatment of systematics in preparation for the imminent abundance of asteroseismic data from future missions.

We report the discovery of LAMOST J0041+3948, the most luminous post-AGB Type II Cepheid (TIIC) known, located in the Andromeda Giant Stellar Stream. Its spectral energy distribution (SED) exhibits a strong near-infrared excess, indicating the presence of a circumbinary dusty disk and hence binarity. SED fitting yields an effective temperature of $T_{\rm eff}=6738_{-262}^{+234}\,$K and a post-AGB luminosity of $\log(L/L_{\odot})=4.32_{-0.08}^{+0.07}$. Comparison with theoretical evolutionary tracks suggests a ~$2.0$-$4.0\,M_{\odot}$ progenitor when accounting for a possible scattered-light contribution. ZTF Light curves reveal a pulsation period of 89d that lies close to the period-luminosity relation for long-period RV Tauri stars. Follow-up spectroscopy reveals clear $s$-process enrichment and signatures consistent with an accretion disk around the companion. The inferred progenitor is significantly younger and more massive than a typical stream member, suggesting that an additional mechanism such as a stellar merger is required. We propose a formation channel in which the present post-AGB binary descends from a hierarchical triple system. In this scenario, the inner binary merged after the system was displaced to its current location by the galaxy merger event, and the resulting massive merger remnant subsequently evolved into the extremely luminous post-AGB star observed today.

In standard cosmology, Cosmic Microwave Background photons near the last scattering surface exhibit only linear polarization due to Compton scattering, leading to the assumption that primordial circular polarization is negligible. However, the physics of Lorentz violation (LV), associated with specific operators, can influence these polarization characteristics. This study employs the Boltzmann equation within the framework of the Standard Model Extension (SME) to explore how the background LV tensor $ K_{AF} $ can induce circular polarization in CMB radiation. By computing the transformation of linear polarization into circular polarization and utilizing the Faraday conversion angle, we derive a bound for $ K_{AF} $ on the order of $ 10^{-41} \, \text{GeV} $, aligning with recent findings. Additionally, we consider the total pure photon terms within the SME, demonstrating that LV in the presence of scalar perturbations can also generate cosmic birefringence (CB) in the CMB radiation. Through analysis of best-fitting CB angles, we establish a more stringent bound of approximately $ 10^{-32} \, \text{GeV} $ for $ K_{F} $.

We present new Rosseland and Planck mean opacity tables relevant to the shallow interiors and atmospheres of giant planets. The tables span metallicities from 0.31 to 50 times solar, temperatures from 100 - 6000 K, and pressures from 1e-6 - 1e5 bar, thereby covering a wider parameter space than previous data sets. Our calculations employ the latest molecular and atomic line lists and pressure-broadening treatments, and include contributions from collision-induced absorption, free electrons, and scattering processes. We further provide cloudy mean opacity tables that account for cloud particle extinction across a range of particle sizes and capture the sequential removal of condensates as the gas cools. We benchmark our cloud-free tables against widely used opacity tables and find significant relative differences, exceeding 100% in Rosseland mean opacities at T \gtrsim 3000 K due to the inclusion of additional short-wavelength absorbers. Differences in Planck mean opacities at high temperatures are even larger, in some cases exceeding two orders of magnitude, which is most likely driven by the inclusion of Ca, Mg, and Fe cross-sections and updated Na D and K I resonance line treatments. Cloud opacities substantially increase Rosseland mean opacities for T \lesssim 2800 K, while their effect on Planck mean opacities is weaker. We also discuss limitations of our mean opacities at high pressures, where non-ideal effects become important. This work provides improved cloud-free mean opacity tables for giant planets, as well as the first publicly available cloudy mean opacity tables, which will enable more realistic modeling of their atmospheres and interiors.

Recently, the Dark Energy Spectroscopic Instrument Data Release 2 (DESI DR2) suggests that the dark energy in our universe might be evolving, favoring the Chevallier-Polarski-Linder (CPL) parameterization and a lower Hubble constant. In our previous work, it has been reported that cosmological model with the non-singular bounce inflation (BI) scenario and $\Lambda$CDM might alleviate the Hubble tension into 3$\sigma$ confidence. In this paper, we study the cosmological model of BI with a dynamical dark energy. We find that individual consideration of the CPL parameterization and the data \texttt{DESI DR2} tend to larger Hubble constants for both BI and power law (PL) case with cosmic microwave background (CMB) data. Employing BI with combined CPL parameterization and \texttt{DESI DR2}, we obtain the Hubble constant $H_ 0 = 65.2^{ + 1.8}_{ - 2.2} \ \mathrm{km} \cdot \mathrm{s}^{ -1 } \cdot \mathrm{Mpc}^{ -1 }$, which is larger than $H_ 0 = 64.0 \pm 2.1 \ \mathrm{km} \cdot \mathrm{s}^{ -1 } \cdot \mathrm{Mpc}^{ -1 }$ for the PL case. After considering nontrivial weak lensing effect and spatial curvature as well as adding \texttt{Pantheon+}, BI fits 3.1$\sigma$ confidence of $\Lambda$CDM with $w_ 0 = -0.919 \pm 0.038$ and $w_{ \mathrm{a}} = -0.37 \pm 0.12$, and it prefers evolving dark energy than the PL case with $w_ 0 = -0.960 \pm 0.074$ and $w_{ \mathrm{a}} = -0.15^{ +0.28}_{ -0.25}$.

OBA hypergiants (OBAHGs) are evolved massive stars with notable wind features in their optical spectrum. Positioned at the cool edge of the line-driven wind regime, many are candidate luminous blue variables (LBVs) likely near the Eddington Limit. Although brief, this evolutionary stage deeply impacts their surroundings and subsequent evolution. We study the mechanisms behind OBAHG winds and spectra, covering the temperature range of non-eruptive LBVs. Using the PoWR atmosphere code, we compute models with an Eddington parameter Gamma_e ~ 0.4 and moderate turbulent pressure, typical for cool hypergiants, varying the effective temperature from ~12.5 to ~38.0 kK at solar metallicity. Our models show a complex temperature-dependent mass-loss pattern, with regions of higher/lower rates linked to two wind solutions: "dense" and "airy." Spectra of known OBAHGs and LBVs match models from all solution regions. We find bi-stability jumps -- with sharp mass-loss increases -- at temperatures where Fe IV recombines to Fe III (and Fe III to Fe II). "Drops" in mass-loss also occur when the leading Fe ion changes at wind onset, signaling a switch to airy solutions under insufficient driving opacity. The resulting velocity fields also reflect these different regimes: airy solutions match the empirical terminal velocity vs temperature relation, while dense ones deviate. Turbulent pressure is crucial for wind acceleration at cooler temperatures, especially in airy regimes. We demonstrate that the bi-stability jumps exist in OBAHGs but are part of a broader complex behavior not replicated by current mass-loss recipes. Combining our and other recent results, we suggest that the switch between airy and dense solutions only occurs within a certain proximity to the Eddington Limit. Testing this requires future models with broader parameters and advanced treatments of radiatively-driven turbulence.

Using the Vera C. Rubin Observatory Data Preview 1 (DP1), we present the detection and photometric analysis of a stellar stream extending along the minor axis of a massive disk-galaxy, LEDA 751050, at a redshift of z = 0.1 in the Extended Chandra Deep Field-South (ECDFS). This feature appears to be an analog of the M31 Giant Stellar Stream, spanning around 130 kpc by 15 kpc, with an estimated stellar mass of M_* = 6x10^8 M_Sun, and a mean g-band surface brightness (SB) of mu_g = 28.6 mag/arcsec^2. This discovery demonstrates Rubin's ability to reveal low-SB substructures and foreshadows the future capabilities when the full Legacy Survey of Space and Time (LSST) data set is actualized.

In this study, we present observational constraints on a class of phenomenological interacting dark energy (IDE) models that admit analytical solutions for the Hubble parameter $H(z)$. We consider a set of five linear and three non-linear IDE scenarios, encompassing both interactions proportional to the dark matter and/or dark energy densities, as well as non-linear combinations of the two. For all eight IDE models, we find a better fit than $\Lambda$CDM from a $\Delta\chi^2$ analysis for both combinations of datasets considered. When using the Akaike Information Criterion ($\Delta$AIC), we find a similarly improved fit in all cases, except for one dataset combination in $Q=3H\delta\rho_{\rm de}$. Our analysis also shows a preference for sign-switching interactions, with energy transfer from dark energy to dark matter at low redshift, reversing direction at higher redshift. These results should be interpreted with caution, as the latter direction of energy transfer is accompanied by negative dark energy densities in the past, which may be unphysical. Models that do not allow sign-changing behaviour instead show a preference for energy flow from dark matter to dark energy, and hence negative dark energy densities. The only exceptions are $Q=3H\delta\rho_{\rm de}$ and $Q=3H\delta\left(\tfrac{\rho_{\rm de}^2}{\rho_{\rm dm}+\rho_{\rm de}}\right)$, which exhibit energy flow in the opposite direction. Furthermore, for all interactions considered, we find a phantom-divide crossing in the effective dark energy equation of state $w^{\rm eff}_{\rm de}$, with the dark energy density decreasing ($w^{\rm eff}_{\rm de}>-1$) at present and at low redshift, while increasing ($w^{\rm eff}_{\rm de}<-1$) in the past at high redshift. These results highlight the promising, but problematic, nature of dark sector interactions, as well as the need to extend the analysis using early-time physics and datasets.

Black holes (BHs) in microquasars can launch powerful relativistic jets that have the capacity to travel up to several parsecs from the compact object and interact with the interstellar medium. Recently, the detection of large-scale very-high-energy (VHE) gamma-ray emission around the black hole transient V4641 Sgr and other BH-jet systems suggested that jets from microquasars may play an important role in the production of galactic cosmic rays. V4641 Sgr is known for its superluminal radio jet discovered in 1999, but no radio counterpart of a large-scale jet has been observed. The goal of this work is to search for a radio counterpart of the extended VHE source. We observed V4641 Sgr with the MeerKAT radio telescope at the L and UHF bands and produced deep maps of the field using high dynamic range techniques. We report the discovery of a large-scale (35 pc), bow-tie-shaped, diffuse, radio structure around V4641 Sgr, with similar angular size to the extended X-ray emission discovered by XRISM. However, it is not spatially coincident with the extended VHE emission. After discussing the association of the structure with V4641 Sgr, we investigate the nature of the emission mechanism. We suggest that the bow-tie structure arose from the long-term action of large-scale jets or disk winds from V4641 Sgr. If the emission mechanism is of synchrotron origin, the radio/X-ray extended structure implies acceleration of electrons up to more than 100 TeV as far as tens of parsecs from the black hole.

The exact location of the solar dynamo remains uncertain--whether it operates primarily in the near-surface shear layer, throughout the entire convection zone, or near the tachocline, a region of sharp transition in the solar rotation, located at the base of the convection zone, approximately 200,000 km beneath the surface. Various studies have supported each of these possibilities. Notably, the solar magnetic "butterfly" diagram and the pattern of zonal flows ("torsional oscillations") exhibit strikingly similar characteristics, suggesting a link between magnetic field evolution and solar flows. Since magnetic fields cannot be measured directly in the deep solar interior, torsional oscillations and rotation gradients are employed as diagnostic proxies. Our analysis reveals that the gradient of rotation displays "butterfly"-like behavior near the tachocline, which is similar to the magnetic butterfly diagram at the surface. This result supports the idea that the solar dynamo has a deep-seated origin, likely operating either near the tachocline or throughout the convection zone, thereby disfavoring the recent scenario of a shallow, near-surface dynamo. This finding may also have important implications for understanding how stellar dynamos operate in general.

With the release of the Gravitational-Wave Transient Catalog GWTC-4.0 by the LIGO-Virgo-KAGRA (LVK) collaboration, 218 candidate detections of gravitational waves (GWs) from compact binary coalescences (CBCs) have been reported. This milestone represents a major advancement for GW cosmology, as many methods, particularly those employing the spectral siren approach, critically depend on the number of available sources. We investigate the impact of a novel parametric model describing the full population mass spectrum of CBCs on the estimation of the Hubble constant. This model is designed to test the impact of heavy black holes in GW cosmology. We perform a joint inference of cosmological and population parameters using 142 CBCs from GWTC-4.0 with a false alarm rate smaller than 0.25 per year, using both spectral and dark siren approaches. With spectral sirens, we estimate the Hubble constant to be $H_0 = 78.8^{+19.0}_{-15.3}\,{\rm km s^{-1} Mpc^{-1}}$ (68% CL), while the dark siren method yields $H_0 = 82.5^{+16.8}_{-14.3}\,{\rm km s^{-1} Mpc^{-1}}$ (68% CL). These results improve the uncertainty by approximately 30.4% and 36.2%, respectively. We show that this improvement is linked to the presence of a new mass scale in the binary black hole mass spectrum at $63.3^{+4.8}_{-4.8}\,M_{\odot}$, which provides additional constraints on the Hubble constant. Besides providing the tightest standard-siren constraints on $H_0$, this highlights the importance of a heavy-mass feature in the black hole spectrum.

We present a complete analytic and semi-analytic study of gravitational collapse and primordial black hole (PBH) formation in the quadratic $f(R)$ model $f(R)=R+\alpha R^2$. We first derive the perturbative expansion around General Relativity (GR), working to first order in the small parameter $\alpha$. For a collapsing flat FLRW dust interior we compute the explicit first-order corrections to the scale factor, the stellar radius, and the horizon formation time. We then use these results to obtain the shift in the PBH formation threshold $\delta_c$. The perturbative effect is small for PBHs forming in the deep radiation era, but becomes important when the background curvature is high. To access this early regime we reformulate the theory in the Einstein frame, where the model becomes GR plus the scalaron field $\phi$ with the Starobinsky potential. We provide the complete ODE system governing both the cosmological background and the evolution of an overdense closed FLRW patch. This system can be numerically integrated to obtain the critical overdensity $\delta_c(k)$ for PBH formation near the end of inflation.

Stellar-mass binaries evolving in the vicinity of supermassive black holes (SMBHs) may be common in the universe, either in active galactic nuclei or in other astrophysical environments. Here, we study in detail the resonant excitation of SMBH modes driven by a nearby stellar-mass binary. The resulting resonant energy fluxes vary with the orbital location and frequency of the binary, exhibiting a rich and complex structure. In particular, we find that the total energy flux radiated to infinity is maximized at a gravitational-wave frequency that is close to, but not exactly equal to, the real part of the corresponding quasinormal-mode frequency. Moreover, as the binary is moved farther away from the SMBH, this offset from the mode frequency becomes increasingly pronounced. In addition, for suitable orientations, the binary can effectively ``feed'' the light ring of the SMBH, selectively exciting particular oscillation modes. For rotating (Kerr) black holes, the mode spectrum is significantly more intricate; however, individual modes are also less strongly damped, leading to an enhanced - but more difficult to interpret - resonant response.

The rapid proliferation of LLM agent frameworks has forced developers to choose between vendor lock-in through provider-specific SDKs and complex multi-package ecosystems that obscure control flow and hinder reproducibility. Integrating tool calling across multiple LLM providers remains a core engineering challenge due to fragmented APIs, incompatible message formats, and inconsistent streaming and tool-calling behavior, making it difficult to build portable, reliable agent systems. We introduce Orchestral, a lightweight Python framework that provides a unified, type-safe interface for building LLM agents across major providers while preserving the simplicity required for scientific computing and production deployment. Orchestral defines a single universal representation for messages, tools, and LLM usage that operates seamlessly across providers, eliminating manual format translation and reducing framework-induced complexity. Automatic tool schema generation from Python type hints removes the need for handwritten descriptors while maintaining type safety across provider boundaries. A synchronous execution model with streaming support enables deterministic behavior, straightforward debugging, and real-time interaction without introducing server dependencies. The framework's modular architecture cleanly separates provider integration, tool execution, conversation orchestration, and user-facing interfaces, enabling extensibility without architectural entanglement. Orchestral supports advanced agent capabilities found in larger frameworks, including rich tool calling, context compaction, workspace sandboxing, user approval workflows, sub-agents, memory management, and MCP integration.

We devise a technique to suppress the effect of noise transients occurring at gravitational-wave detectors based on temporal anticoincidence. Searches for gravitational-wave signals in the detector data are prone to spurious disturbances of terrestrial origin. The technique presented here benefits from the fact that such effects are generally non-coincident in time at geographically separated detectors. Therefore, abnormally loud triggers that are not time-coincident can be vetoed. We implement the veto technique in a matched-filter search for transient signals from binary black holes and observe search backgrounds to be generally close to the Gaussian limit. An improvement in the sensitivity of the search is demonstrated using simulated signals. The technique is expected to especially improve the detection efficiency of the search toward short duration transient signals.

The equation-of-state (EOS) parameter $\phi \equiv P/\varepsilon$, defined as the ratio of pressure to energy density, encapsulates the fundamental response of matter under extreme compression. Its value at the center of the most massive neutron star (NS), $\x \equiv \phi_{\rm c} = P_{\rm c}/\varepsilon_{\rm c}$, sets a universal upper bound on the maximum denseness attainable by any form of visible matter anywhere in the Universe. Remarkably, owing to the intrinsically nonlinear structure of the EOS in General Relativity (GR), this bound is forced to lie far below the naive Special Relativity (SR) limit of unity. In this work, we refine the theoretical upper bound on $\x$ in a self-consistent manner by incorporating, in addition to the causality constraint from SR, the mass-sphere stability condition associated with the mass evolution pattern in the vicinity of the NS center. This condition is formulated within the intrinsic-and-perturbative analysis of the dimensionless Tolman--Oppenheimer--Volkoff equations (IPAD-TOV) framework. The combined constraints yield an improved bound, $\x \lesssim 0.385$, which is slightly above but fully consistent with the previously derived causal-only limit, $\x \lesssim 0.374$. We further derive an improved scaling relation for NS compactness and verify its universality across a broad set of 284 realistic EOSs, including models with first-order phase transitions, exotic degrees of freedom, continuous crossover behavior, and deconfined quark cores. The resulting bound on $\x$ thus provides a new, EOS-independent window into the microphysics of cold superdense matter compressed by strong-field gravity in GR.

In this paper, we study two inflationary models, namely, monomial inflation and the simplest $\alpha$-attractor inflation, within extra-dimensional frameworks. We consider three extra-dimensional setups: Dark Dimension, which embeds one flat extra-dimension to explain the observed smallness of the 4D cosmological constant $\Lambda_4$; and the two Randall-Sundrum scenarios with one warped extra-dimension, namely RS1 with two branes and RS2 with one brane. We derive the corresponding Friedmann equations, compute the slow-roll parameters in each case, and we fit the experimental data for ($n_s - 1$, $\alpha$, $\Delta_s^2$, $r$), using Planck, BICEP, and ACT data. We find that monomial inflation is strongly disfavored in all scenarios, while $\alpha$-attractor inflation provides an excellent fit to current observations, with extra-dimensional setups offering additional flexibility compared to the standard 4D case.

In the Standard Model, the Higgs potential develops an instability at high field values when the quartic self-coupling runs negative. Large quantum fluctuations during cosmic inflation could drive the Higgs field beyond the potential barrier, creating regions that would be catastrophic for our observable universe. We point out that the extreme-value statistics describing the peaks (maxima) of the Higgs values is the correct statistics to infer the condition to avoid vacuum instability. Even if this statistics delivers a bound on the Hubble rate during inflation which is only a factor $\sqrt{2}$ stronger than the one commonly adopted in the literature, it is qualitatively distinct and we believe worthwhile communicating it.