Papers for Tuesday, Feb 24 2026

This paper provides experimental and numerical evidence supporting the occurrence of liposome congregation at the floors of meteor craters on Early Earth. This work builds on our earlier research, which demonstrated that liposomes submerged in a shallow Archean pond are protected from harmful UV radiation. This protection allows them to survive long enough for autocatalytic replication of amphiphiles and for mutation and selection of assemblies that maximize membrane stability. For liposomes to fuse, grow, exchange contents and membranes, and divide, they need to establish a population, which means forming a dense conglomerate that enables close physical contact. The study demonstrates that such a congregation is feasible in bowl-shaped meteor craters on Early Earth, especially under periodic seismic disturbances.

We analyse data on the post-perihelion morphology, including jet position angles (PAs) and coma dominated photometry of the interstellar object 3I/ATLAS. From Hubble Space Telescope (HST) images processed with a Larson Sekanina rotational gradient filter, we measure the PAs of three main persistent jet like features between November 30 and December 27, 2025 and fit a weighted Fourier model in a period scan. The dominant jet PA wobble yields Pjet = 7.20 +/- 0.05 h. An independent Gr (R band) time-series photometry data set, using two different apertures from MPC station L92, analyzed with nightly offsets and 30 minute binning, gives Pphot = 7.136 +/- 0.001 h (formal 1sigma), with a semiamplitude A about 0.311 mag and scatter sigmajit about 0.089 mag. The close agreement beetwen the periods supports a characteristic post-perihelion period of about 7.1 h. We interpret this period as an attitude precession/nutation (non principal axis rotation) traced by jet orientation and coma flux redistribution. The jet structure precesses about the rotation axis with a characteristic angular excursion of order about 20 degs, and the rotation axis is aligned with the sunward direction to within about 20 degs.

I present the first systematic catalogue of simultaneous solar transits of both Phobos and Deimos as observed from the surface of Mars. Using the JPL mar099 ephemeris (Brozovic et al. 2025) and SPICE toolkit, I searched the millennium 1600-2600 CE for epochs at which both Martian moons project onto the solar disc at the same instant for at least one surface location. I identify 8565 grazing double transits, 49 partial-overlap double transits (both moons simultaneously on the solar disc, with at least one partially cut off by the limb), and 17 full double transits in which both moons lie wholly within the solar disc at the same moment. All events cluster tightly around the Martian equinoxes (Ls ~ 0 deg or 180 deg) and within +/- 9 deg of the equator, reflecting the near-equatorial orbital inclinations of both moons. I derive a hard theoretical latitude limit of +/- 13.1 deg beyond which simultaneous double transits are geometrically impossible. The next observable partial double transit, excluding less prominent grazing events, is predicted for 2034 April 17. The next full double transit, with both moons wholly inside the solar disc, with a gap between each silhouette and the solar limb, occurs on 2118 November 20. The geometries for both these events were confirmed with JPL Horizons. I provide uncertainty estimates based on the Brozovic et al. covariance model, with predicted position errors growing from ~1 km for near-term events to ~600 km at the catalogue boundaries, and note that the JAXA MMX mission (~2031) will dramatically reduce uncertainties for all post-2030 predictions.

Using SDSS-V DR19 Milky Way Mapper APOGEE data, we measure the impact that close binarity has on surface chemistry across the Hertzsprung-Russell diagram in a broad set of abundances by studying single-lined spectroscopic binaries (SB1s) in open clusters. We derive binary membership and orbital parameters for 103 SB1s by analysing APOGEE radial velocities with The Joker and UltraNest. We perform a detailed abundance analysis with BACCHUS to derive abundances in fourteen chemical species: Si, Fe, C, N, O, Na, Mg, Al, Ca, Ti, Cr, Ni, Ce, and Nd. Leveraging the assumptions of chemical homogeneity in open clusters, we compare the surface abundances of SB1s to non-binary stars at similar evolutionary states. We find that a subset of binaries with significant UV excess have a $\Delta$[C/N] that is 0.2--0.5 dex higher than expected, resulting in overestimated [C/N]-based ages for those stars. This points to pollution from an evolved companion and has implications for [C/N]-based age studies of the broader Milky Way. At the population level, we find that SB1s in our sample can be treated as statistically chemically homogeneous with their single-star counterparts, and we find no connection between orbital separation and chemical enrichment or depletion. We show that at separations up to ~5 pc, co-eval stars can be considered chemically homogeneous with one another within current abundance precisions, regardless of multiplicity.

We analyze the public DESI full-shape clustering data using simulation-based priors (SBPs). Our priors are obtained by fitting normalizing flows to the distribution of EFT parameters measured from field-level simulations, themselves generated using tailored halo occupation distribution (HOD) models for each tracer. Incorporating SBPs in a power spectrum analysis significantly enhances $\Lambda$CDM cosmological parameter constraints; in combination with BAO information from DESI DR2 and a BBN prior on the baryon density, we find the matter density parameter $\Omega_m=0.2987\pm 0.0066$, the Hubble constant $H_0=68.80\pm 0.35\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, and the mass fluctuation amplitude $\sigma_8 = 0.766\pm 0.015$ (or the lensing parameter $S_8=0.764\pm 0.018$), which are $1\%$, $40\%$ and $50\%$ stronger than the baseline results, though with a notable downwards shift in $\sigma_8$. The SBPs also have a significant impact in extended models, with the dark energy figure-of-merit improving by $70\%$ ($20\%$) in a $w_0w_a$CDM analysis when combining with the CMB (and supernovae). In the SBP analysis, we do not find statistically significant evidence for dynamical dark energy: the equation of state parameters are consistent with a cosmological constant within $2.2\sigma$ ($1.4\sigma$) in analyses without (with) supernovae. The neutrino mass constraints are also enhanced, with the $95\%$ limits $M_\nu<0.073\,\mathrm{eV}$ and $M_\nu<0.090\,\mathrm{eV}$ in $\Lambda$CDM and $w_0w_a$CDM respectively. The latter is the strongest constraint obtained to date and reinforces the preference for the normal neutrino mass hierarchy, regardless of the background dynamics. While our results are sensitive to HOD modeling assumptions, they clearly demonstrate that the inclusion of small-scale information can significantly sharpen cosmological parameter constraints.

Recent JWST observations reveal massive, UV-bright galaxies at $z > 10$ with little apparent dust attenuation, whereas ALMA detections at $z \simeq 7$ show similarly massive systems that are already dust-rich and IR-luminous. This raises a fundamental question: can a single physical model of star formation and dust production explain both populations across cosmic time? We address this using a minimal framework with only two free parameters--the instantaneous star formation efficiency ($\epsilon_\star$) and the dust yield per Type II supernova ($y_d$)--and predict the rest-frame UV and IR luminosity functions (LFs) from $z \simeq 14$ to 7. For a uniform ISM, we find a UV-IR tension at the bright end of the LFs at $z \ge 7$. The UV LF requires low dust yields ($y_d \lesssim 0.01\,M_\odot$), whereas the $z=7$ IR LF requires higher yields ($y_d \sim 0.1\,M_\odot$) unless the star formation efficiency is boosted above $\epsilon_\star \approx 5$-10%. We show that incorporating a porous, turbulent ISM largely resolves this tension: turbulence opens low-column-density sightlines that enhance the UV escape fraction while leaving the total absorbed energy--and thus the IR luminosity--nearly unchanged once radiative-transfer--induced flattening of the attenuation curve is included. Large-grain dust distributions, while reducing UV opacity, play a secondary role once ISM porosity and radiative transfer are taken into account. At $z > 10$, however, even strong turbulence cannot reproduce the bright end of the UV LF at high dust yield. This could be resolved either by efficient dust removal in early massive systems or by substantial ISM dust growth by $z \simeq 7$. Our results highlight dust physics as a key lever for interpreting the rapidly growing UV and IR constraints within the broader context of early galaxy formation.

Galactic binaries are expected to be the most numerous LISA sources and to produce a stochastic gravitational-wave foreground whose spectral shape encodes information about the underlying population. Extracting this information with standard hierarchical methods is challenging due to the high dimensionality of the problem and the computational cost of global-fit analyses. We present a simulation-based inference framework to measure the population properties of galactic binaries directly from the reconstructed foreground. Adopting an astrophysically agnostic parametrization in the observable space -- defined by signal amplitude, frequency, and frequency derivative -- we generate synthetic catalogs and foreground spectra using a global-fit-inspired subtraction algorithm. We then train a neural posterior estimator to map spectra to population parameters. We validate our method on simulated data and recover population parameters with good accuracy, including the total number of binaries. As a by-product, we present a GPU-accelerated version of the subtraction algorithm, which delivers a ~100X speed-up compared to previous implementations in the literature. Our results demonstrate that LISA's stochastic foreground alone carries significant information about the Galactic binary population and provide a practical step toward joint inference from resolved and unresolved sources.

The rotational properties of small near-Earth asteroids (NEAs) provide crucial insights into their internal structure and collisional history. However, systematic surveys targeting metre- to decametre-sized bodies are rare, thus leaving their spin distribution poorly constrained. Our aim was to quantify the prevalence of fast rotation and characterise the spin-rate distribution of small NEAs to constrain their internal strength and evolution. We conducted a dedicated high-cadence photometric survey of 249 NEAs using the Two-meter Twin Telescope (TTT). Rotation periods and amplitudes were derived from dense time series to classify objects as fast or non-fast rotators. We determined the rotation period of 156 new fast rotators (P < 2.2 h), including 87 that rotate faster than 10 min. The prevalence of fast rotators increases with absolute magnitude: from 60.6-80.3 % for 22 < H < 24 to 77.3-89.4 % for 24 < H < 26 and 94.1-96.1 % for H > 26, indicating that fast rotation dominates in the small NEA population. Most objects spin faster than the gravity-defined limit; 98 targets require cohesive strengths exceeding that of weak rubble piles, and 22 are compatible only with compact, high-strength interiors. This is the first systematic survey targeting the rotation of such small NEAs, providing the largest homogeneous sample of fast rotators obtained by a single campaign. Our findings demonstrate that fast rotation is the norm for objects smaller than tens of metres, implying that modest cohesive strength is required to prevent their rotational disruption.

The elemental abundances of neon isotopes provide valuable insights into stellar evolution and nucleosynthesis. In this study, we calculate the abundances of the isotopes 18Ne, 19Ne, 20Ne, 21Ne, and 22Ne across the five principal evolutionary phases of the Sun: hydrogen burning, lively old age, onset of rapid growth and red giant, helium burning and helium exhaustion. The calculations were carried out using the open-source nucnet-tools package, developed by the Webnucleo Group at Clemson University. Initial isotope abundances were adopted from standard proto-solar compositions. Their evolution was computed under static hydrostatic burning conditions, assuming constant temperature and density within each phase. The results show that the stable isotopes 20Ne and 22Ne remain dominant throughout the Sun's lifetime, whereas the short-lived isotopes 18Ne and 19Ne decay rapidly during or shortly after the hydrogen-burning phase. The predictions obtained for the helium burning and exhaustion phases provide quantitative neon-isotope abundances that are not extensively reported in the existing literature. These results offer valuable reference values for future studies of solar and stellar evolution, nucleosynthetic pathways, and isotopic modeling.

We investigate the response of space weather events on Earth's upper atmosphere over the polar regions by studying their effect on the drag of the CHAMP and GRACE satellites. Increasing solar activity that results in heating and the expansion of the upper atmosphere threatens low Earth orbit (LEO) satellites. Auroral events are closely related to the stellar energy deposition of solar EUV radiation and precipitating energetic electrons, which influence photochemical processes such as the production of nitric oxide (NO) in the upper atmosphere. To study the production of NO molecules and their influence on the thermospheric structure and satellite drag, we first model Earth's background thermosphere with the 1D upper atmosphere model Kompot by considering the incident X-ray, EUV, and IR radiation during selected space weather events. To investigate the effect of electron precipitation in the production of NO molecules in the polar thermosphere, we apply a Monte Carlo model accounting for the stochastic nature of collisional scattering of auroral electrons in collisions with the surrounding N$_2$-O$_2$ atmosphere, including the production of suprathermal N atoms. The observed effect of the atmospheric drag on CHAMP and GRACE during the two studied events indicates that a sporadic enhancement of NO molecule production in the polar thermosphere and its IR-cooling capability, which counteracts thermospheric expansion and can lead to an ``overcooling'' with decreased density after the space weather event, can have a protective effect on LEO satellites. Their production efficiency, however, is highly dependent on the energy flux of the precipitating electrons. Our results have direct implications for empirical satellite orbit prediction models, as our simulations highlight the need to consider precipitation-induced NO production to improve the predictive power of these models.

We have followed up on the hot-type classical symbiotic outburst reported in YY Her using the Himalayan Chandra telescope. The outburst coincides with the secondary minima of the system. Approximately 12 similar brightening events have been reported between 1890 and 2020, with only the 1993 outburst being studied spectroscopically. In our study, we monitored the system from 2021 to 2023, covering $\sim$1.5 orbital cycles, providing an opportunity to understand the spectral evolution of the outburst over a complete orbital period of YY Her. We found that the temperature and luminosity estimations based on emission line fluxes exhibit orbital phase dependence. The values estimated at phase 0.5, corresponding to the secondary minimum, were the most reliable. The temperature of the hot component is $\approx 1.41 \times 10^5 \, \text{K}$, and the luminosity is $\approx 1020 \, L_{\odot}$ during the outburst, reduced to $\approx 1.3 \times 10^5 \, \text{K}$ and $\approx 830 \, L_{\odot}$ after one orbital cycle at phase 0.5. Temperature estimations during the outbursts suggest that YY Her exhibits both hot-type (2021) and cool-type (1993) behavior, similar to another symbiotic star, AG Dra. Using variations of the Ca II absorption lines, we confirmed the contribution of the ellipsoidal effect in secondary minima in the YY Her light curve.

At low coronal temperatures around or below 1 MK distinct areas in the surroundings of active regions (AR) show emission at a level significantly below the emission coming from the quiet Sun (QS). These areas are referred to as dark halos, dark canopies, or dark moats. To better understand the nature of dark halos we study the connection between the photospheric magnetic field and coronal emission at different temperatures. Combining Solar Orbiter data from the high-resolution Polarimetric and Helioseismic Imager (SO/PHI) and Extreme Ultraviolet Imager (EUI) instruments allows us to identify these areas that are dark in the extreme ultraviolet (EUV) in the immediate vicinity of an AR. We probe both the photospheric magnetic field as well as the coronal intensities as a function of distance to the AR NOAA 12893. The dark halo has an unsigned magnetic flux density similar to the QS, but shows a strong radial dependence with distance from the AR centre. It drops by 38 % from 6.1 G at the inner boundary to 3.8 G at the outer, shifting from above to below QS levels. Coronal emission $\leq$1 MK is $\sim$40 % below QS and shows no dependence on distance to the AR centre. In contrast, at $\geq$1.6 MK, emission exceeds QS levels, but declines outward toward QS values. A few hot loops extend from the AR periphery across the halo, while at lower temperatures no such loops appear and short loops dominate the corona. The reduced unsigned magnetic flux density in the outermost parts of the dark halo, below QS level, suggests that reduced coronal heating due to weak underlying magnetic flux heating could be partially responsible for the reduced emission around 1 MK. Closer to the AR, other mechanisms might lead to reduced heating. The different loops structures detected for hotter and cooler coronal temperatures likely play a crucial role in understanding coronal dark halos.

Using repeat imaging of a galaxy cluster taken over a seventeen-year baseline, we examine the impact that degraded Charge Transfer Efficiency (CTE) has on photometric measurements of extended sources using the ACS/WFC on HST. We examine how measured brightnesses depend on time since ACS installation, source location on the WFC detectors, source brightness, and local background level in individual exposures. We find that global brightness measurements using large apertures are generally reliable within $\sim$0.05 magnitudes across the WFC detectors if exposure backgrounds are above $20e^-/{pixel}$ and sources are brighter than $\sim300e^-$ in a single exposure. However, brightness measurements on smaller scales can suffer deficiencies in excess of 0.1 mags (sometimes, significantly more) in recent data unless sources are very close to the CCD serial registers ($\lesssim 512$ pixels), or brighter than $\sim3000\,e^-$ in a single exposure. We also show how degraded CTE can result in artificial asymmetries in galaxy light distributions, which are largely mitigated if backgrounds are $>20e^-/{pixel}$ and targets are not far ($>1536$ pixels) from the serial registers. As expected, brightness measurements in later epoch data are best when using CTE-corrected images (FLC/DRC), but our results imply that the pixel-based CTE correction algorithm employed by the ACS reduction pipeline does not necessarily place charge back into its original location within extended sources. Based on this study, users are advised to keep backgrounds above the already recommended $30e^-/\mathrm{pixel}$, ensure targets will have at least $\sim 300e^-$ in a single exposure, and place targets close to the serial registers if analysis of their spatially resolved properties is needed.

The boundaries of dark matter and gas in clusters are delineated by the splashback radius and the accretion shock, respectively. Theoretically, both of these boundaries are expected to coincide at the outskirts of halos. However, hydrodynamic cosmological simulations have highlighted significant displacement between them. In this study, we utilise the IllustrisTNG simulation suite to investigate the statistical relationship between the splashback and shock surfaces in a sample of 812 cluster-mass halos. We compute the full angular distribution of both boundaries and examine their relationship, also considering how different moments of this distribution correlate with halo properties. We employ a dispersion-based measure for the splashback boundary and the maximum entropy distance for the shock location. Despite examining various boundary definitions, we consistently observe an offset between the splashback and shock boundaries, with $R_{\rm sh}/R_{\rm sp} \sim 1.3-2$, depending on specific methodological choices. This offset predominantly occurs along void directions. We analyse the redshift evolution of these boundaries for a subset of halos and find that splashback and shock boundaries are not necessarily distinct at earlier times. During mergers, gas dissipates energy and resists contraction via pressure, unlike collisionless dark matter, leading to the observed boundary offset. We also find that the feature in pressure profiles arising from the outer accretion shock is sensitive to the exact method of stacking, which has important implications for observations.

Solar energetic particles (SEPs) are produced in two fundamental ways: at magnetic reconnection sites in solar jets and at collisionless shock waves driven by fast coronal mass ejections (CMEs). "Impulsive" SEP events, on open field lines from jets, have signature abundance enhancements of 3He and of increasingly heavy elements, and their outward streaming electrons drive type-III radio bursts. Similar acceleration for particles trapped on closed loops energizes solar flares. In contrast, fast, wide, CME-driven shocks accelerate seed ions from the ambient corona that grow resonant Alfven waves as they stream outward. These waves can scatter and trap lower-rigidity ions near the shock, limiting outflow, and flattening low-energy spectra upstream at the "streaming limit." Downstream, a spatially-uniform "reservoir" of SEPs is shed by the expanding shock between it and the Sun. These trapped invariant SEP spectra decrease in intensity adiabatically as the volume of the reservoir expands. Ions from the reservoir can seed further acceleration in multi-shock event and energetic proton precipitation can prolong solar gamma-ray emission. Shocks are often additionally seeded by residual impulsive ions which dominate the SEP heavy-ion abundances with their signature enhancements. As samples of the corona, SEP abundances also probe differences with the photosphere that depend upon the first ionization potential (FIP) of the elements.

It has recently been discovered that the time phases (time-varying delays) of the irregular periodicities observed in the Sunspot series, consisting of the approximate 27-day latitude-averaged plasma rotation and approximate 11-yr magnetic pole reversal, share a common pattern or cycle of irregularity. Because the method of measurement contains an unknown constant delay of no practical significance, the average of the measured cycle of delay is subtracted out, and the resultant cycle proceeds from zero delay to maximum delay, back to zero delay, then to maximum advance, and finally back to zero advance/delay. The period of this cycle is estimated to be 164 years, which does not allow for even one full repetition of the pattern in the 205-yr data record. So, it may or may not be a segment of an ongoing (pseudo) periodicity. Nevertheless, the fact that this cycle of irregularity in period is shared by two long-recognized genuine irregular periodicities tied to specific activity from two distinct solar phenomena suggests there is a common phenomenon affecting these periodicities. The physics-based link between the rotation speed of the interior of the Sun and the average over latitude of the rotation period of the surface plasma, approximately 27 days on average, suggests the possibility that there may be a single primary phenomenon that is responsible for irregularity of both rotation of the interior of the Sun and the period, approximately 11 yrs on average, of pole-reversal. This is consistent with recently reported (November 2025) results on hypothesized planetary forcing of the Solar dynamo.

Tidal disruption events which repeat on timescales of months-to-years represent an unambiguous signature of a partial disruption, with the surviving stellar remnant returning to pericentre to be repeatedly stripped by tidal forces. These systems therefore offer the best laboratories to study the differences between partial and full disruptions. One noteworthy observational difference between the two systems is that all known X-ray bright repeating TDEs show rapid transitions between thermal, non-thermal and completely dim states on timescales much shorter than full (non-repeating) TDEs. We argue this can be simply understood as being due to the reduction in fuel supply available to the disk, and that these systems provide evidence that all tidal disruption events undergo a state transition at Eddington ratios $f_{\rm edd} = L_{\rm bol}/L_{\rm edd} \sim 0.01$, similar to X-ray binaries. {As part of this calculation we derive a general expression for the time taken for a TDE disk to fall to a given Eddington fraction, which will be of use to both full and partial TDE science.} Perhaps surprisingly, the late-time optical/UV plateau luminosity observed from these systems is largely unaffected by this reduction in fuel supply, provided the outer disk remains in a thermal state for long enough for this emission to be detected. We then show that collisions between the returning stellar remnant and the disk formed from the last passage will produce potentially observable X-ray flares ($L_{\rm flare} \simeq 10^{42}$ erg/s), but that they are likely to be very difficult to detect as they are generally short-lived ($t_{\rm flare} \simeq 0.1-1$ hr).

Tearing instability, also known as plasmoid instability, is an effective mechanism to speed up magnetic reconnection process, working in a wide range of magnetized plasma systems with different spatial scales, ionization degrees, and collisionality. However, due to observational limitations, observations of {plasma sheet} tearing and the resulting plasmoids are rather scarce. This scarcity significantly hinders our understanding of the role of plasmoids in the reconnection process from an observational perspective. Using high-spatiotemporal multiwavelength observations from the Solar Dynamics Observatory, we traced the entire evolution of a coronal {plasma sheet}. Its formation was driven by the emergence of photospheric magnetic flux, followed by tearing, and eventual decay. The evolution of the {plasma sheet} exhibited two distinct stages. Initially, it rose rapidly, lengthened, and underwent tearing at a low frequency. Subsequently, its ascent slowed, it began to shorten, and the tearing occurred more frequently. Detailed analysis of the reconnecting {plasma sheet} focuses on heating, plasmoid dynamics (formation and ejection), and the resulting reconnection rate change. Two key heating processes are identified: {plasma sheet} tearing and coalescence involving plasmoids and magnetic cusps. More importantly, combining observations with analytical studies suggests that plasmoids are key carriers of magnetic flux fast transferring in the observed torn {plasma sheet}, and their formation and ejection significantly enhance the reconnection rate and facilitate the onset of fast reconnection.

Zeeman observations provide the only direct probe of line-of-sight (LOS) magnetic fields in the interstellar medium. To evaluate their accuracy and limitations, we generate synthetic HI Zeeman spectra from magnetohydrodynamic simulations and idealized cloud models, and analyze the resulting Stokes I and V profiles using two complementary methods. Approach I uses the classical relation between Stokes V and dI/d{\nu} to estimate LOS-averaged magnetic fields, achieving an upper-limit relative error of 16% (half-width of 68.27% confidence interval) for a representative noise level of 0.014 K. Approach II applies Gaussian decomposition to Stokes I and V to estimate component-level magnetic fields, yielding a 13% relative error quantifying the same confidence range, reflecting the intrinsic uncertainty of such Zeeman estimates. Both approaches recover the original fields under uniform-field conditions and remain robust in turbulent environments. Approach I provides a simple and reliable LOS-averaged field estimate, while Approach II, although more complex, offers statistical insight into magnetic field variations along the LOS. We further show that joint fitting of Stokes I and V generally outperforms sequential fitting, particularly in the presence of attenuation. Increasing noise eight-fold produces a more modest rise in uncertainty, doubling to a 26% relative error, while substantial optical depth introduces only a minor additional contribution to the overall uncertainty. Applying these methods to FAST observations of the L1544 star-forming region, we confirm the previously reported LOS magnetic field strength, demonstrating the validity of Zeeman analysis in this benchmark core.

We review our semi-analytic model of stellar black hole (BH) mass growth by gas accretion in gas-rich stellar clusters during their birthstage within the first $\sim 10\,{\rm Myr}$ after the first stellar formation event. Such proto-stellar clusters are massive and compact, with typical masses $\sim 10^6\,{\rm M}_\odot$ and sizes $\sim 1\,{\rm pc}$, suggested by recent James Webb Space Telescope (JWST) observations. We find that the BH masses are shifted by the end of gas depletion to values within and above the BH mass gap, well within the range of components of the recent gravitational-wave (GW) signal GW231123, and up to masses $\sim 10^3\,{\rm M}_\odot$.

The Forbush decrease (FD) of cosmic rays is an important probe of the interplanetary environment disturbed by solar activities. In this work, we study the properties of 8 FDs electrons (including positrons) between 2 GeV and 20 GeV from January, 2016 to March, 2024, with the Dark Matter Particle Explorer. The maximum decrease amplitudes of these events are about 30% - 15%, and the amplitudes reduce with energy. The recovery time of these events shows diverse behaviors of their energy-dependence. Some of them show strong energy-dependence, while some have a nearly constant recovery time. It has been shown that such diverse behaviors could be related with the geometry of the disturbed regions of the interplanetary space by coronal mass ejections (CME), represented by the combined effect of the CME velocity, angular spread, and ejection direction.

We are discussing the capabilities of MASTER network of robotic telescopes to study variable red giant stars. We report the discovery and present archival light curves for three optical transients, connected with long period variable red giants, based on MASTER Robotic Net all-sky survey observations and archival MASTER light curves at Lomonosov database storage. We demonstrate transient detection frames and light curves for three Mira-type variables: MASTER OT $J083717.54-573411.1$, MASTER OT $J190436.33+192828.7$ and MASTER OT $J07010810-6818548$. Spectroscopically validated carbon star MASTER OT $J07010810-6818548$ shows mean magnitude growth but constant pulsation period which may indicate the presence of a long secondary period. As the result of MASTER wide-field images analysis at Lomonosov supercomputer data storage we present the 8-years historical light curves for all three objects with period calculations for all three Mira-type variables. We estimated fundamental astrophysical parameters for all three Mira-type variables: luminosities, radii, and zero age main sequence (ZAMS) masses.

We present a cosmological analysis combining full-shape (FS) clustering measurements from the Dark Energy Spectroscopic Instrument (DESI) DR1 with baryon acoustic oscillation (BAO) measurements from DESI DR2. To achieve a robust combination that accounts for the correlation between the two data releases, we employ the ShapeFit compression method and estimate the joint covariance using EZmocks. This compressed approach inherently mitigates the prior volume effects that have previously dominated Bayesian constraints from DESI data with minimal external priors. Consequently, we obtain--for the first time within a Bayesian framework--reliable DESI-only constraints on extensions to $\Lambda$CDM using only a Big Bang Nucleosynthesis prior on the baryon density and a wide prior on the spectral index. In flat $\Lambda$CDM, we find $\Omega_m = 0.3035 \pm 0.0085$, $h = 0.6876 \pm 0.0059$, and $\sigma_8 = 0.822 \pm 0.034$. For the $w_0 w_a$CDM dynamical dark energy model, we measure $w_0 = 0.49 \pm 0.25$ and $w_a = -1.52 \pm 0.77$, improving constraints by $\sim 30\%$ relative to the analogous DR1 measurement and reducing the discrepancy with $\Lambda$CDM to $1.4\sigma$ when compared to BAO only analyses. We also report competitive limits on the sum of neutrino masses and spatial curvature. This work demonstrates that the ShapeFit compression provides a prior-robust and computationally efficient pathway to constrain beyond-$\Lambda$CDM physics with large-scale structure.

Burst oscillations during thermonuclear X-ray bursts are powered by thermonuclear energy on the neutron star (NS) surface and typically occur close to the spin frequency of the NS. We performed a comprehensive timing analysis of all thermonuclear bursts from the newly discovered millisecond X-ray pulsar SRGA J144459.2$-$604207, observed with NICER, XMM-Newton, and NuSTAR during the 2024 outburst. A total of 39 bursts were detected, allowing for a detailed search for burst oscillations, which had not been previously observed from this source. We report the discovery of burst oscillations at 447.7$-$448.0 Hz from SRGA J144459.2$-$604207 using XMM-Newton and NuSTAR data, consistent with the spin frequency of the NS. The strongest burst oscillation in the XMM-Newton data occurred with a single-trial significance of $5.1\sigma$ and maximum $Z^2$ power of $\sim31$. In the NuSTAR data, the strongest oscillation signal has a significance of $5.2\sigma$ and maximum $Z^2$ power of $\sim32$. The folded pulse profile corresponding to the strongest signal in the 0.5-10 keV band of the XMM-Newton data shows a sinusoidal shape with a fractional rms amplitude of $\sim8.5\%$, while the measurements of the NuSTAR data (3-40 keV range) yield $\sim21\%$. These results represent the first detection of burst oscillations in SRGA J144459.2$-$604207. Additionally, we report the detection of 447.6 Hz oscillations occurring just before a burst onset observed with XMM-Newton. This marks only the second instance in which burst oscillations have been observed before the burst onset.

On Mars, a relatively pure water ice layer lies beneath several centimeters of dry soil at mid-latitudes. Its widespread presence poleward of 60° latitude was detected by remote neutron spectroscopy and confirmed by the Phoenix lander at 68°N. Recent observations of exposed ice indicate that the near-surface ice layer extends to 35° latitude and exhibits pronounced spatial structure. However, previous models did not capture the observed spatial structure of the midlatitude ice layer. Here, on the basis of improved calculations using the Mars Planetary Climate Model, we show that mid-latitude buried ice could be the remnant of a ice layer deposited on the surface when the obliquity was higher than today. Assuming that the ice subsequently sublimated and became buried beneath a sublimation lag, we estimate that surface ice emplaced 630 kyr (4.18 Myr) ago at 35° obliquity (40°), at latitudes of 40-55°N, would today reside at depths of 25-150 (41-255) cm, depending on the regolith and ice properties. For ice emplaced 630 kyr ago, the modeled burial depths align with observations and capture the observed longitudinal depth variations, in contrast to ice emplaced 4.18 Myr ago. We therefore infer that the mid-latitude subsurface ice is younger than 4 Myr.

This work examines the cosmological implications of two functional forms of $f(R,T) = R + \alpha T^n$ gravity: for two different value of $n $ where $n=1$ and $n\neq 1$, and $\alpha$ and $n$ are free parameters. The modified Friedmann equations are derived, and the cosmic evolution of the Hubble parameter $H(z)$ is determined. Cosmological parameters are estimated through $\chi^2$ minimization and MCMC analysis using the emcee algorithm, with model parameters constrained by various observational datasets. Both models reproduce late-time acceleration and remain observationally indistinguishable from $\Lambda$CDM, while allowing small deviations parameterized by $\alpha$ and $n$. The cosmological behavior of the deceleration parameter $q(z)$, the jerk, the snap parameter $s(z)$, and the effective equation of state parameter $\omega$ is also analyzed. The results indicate that the Universe transitioned from deceleration to late-time accelerated expansion, consistent with the $\Lambda$CDM model. Analysis of the energy conditions reveals that the NEC and DEC are satisfied, while the SEC is violated, which explains the transition from a decelerating matter-dominated epoch to the present accelerated phase. These findings indicate that the proposed $f(R,T)$ models are compatible with current observational data and provide a viable alternative to $\Lambda$CDM in describing cosmic acceleration.

Automated characterization of galactic substructure is an essential step in understanding the transformative physical processes driving galaxy evolution. In this study, we investigate the application of deep learning (DL) frameworks to characterize different galactic substructures hosted within parametric light-profile subtracted ``residual'' images of a large sample galaxies from the CANDELS survey. We develop a supervised Convolutional Neural Network (CNN) and unsupervised Convolutional Variational Autoencoder (CvAE) and train it on the single-Sérsic profile fitting based residual images of $10,046$ bright and massive galaxies ($H<24.5\,{\rm mag}$ and $M_{\rm stellar} \geq 10^{9.5}\,M_{\odot}$) spanning $1<z<3$, in conjunction with their visual-based classification labels indicating the nature of residual substructures hosted within them. Using our unique data preprocessing approach, we prepare our residual images such that the inputs to our DL networks comprise only ``galaxy of interest'', and augment them such that our sample span uniformly across different residual characteristics. We assess the latent space of the CNN and CvAE using Principle Component Analysis (PCA) along with independently quantified metrics of residual strength (significant pixel flux $SPF$, Bumpiness, and Residual Flux Fraction). We also employ an unsupervised Gaussian Mixture Modeling (GMM) based clustering scheme with Support Vector Classification (SVC) to identify groupings in PCA space that correspond to similar residual substructure. We find that our supervised CNN latent features in PCA space correlate with the $SPF$ values and distinguish between qualitatively strong and weak residual substructures. While our unsupervised CvAE latent space also correlates with visual and quantitative residual characteristics, but lacks clear discriminatory power when characterizing different residual substructures.

Modern astronomical surveys are producing progressively larger and more complex datasets, making traditional supervised approaches that rely on extensive labelled catalogues increasingly difficult. Consequently, pre-training using self-supervised learning (SSL), which offers a scalable route by extracting structure directly from unlabelled images, is becoming attractive for many downstream applications. In this work we consider the use of coupled self-supervised representation learning approaches for radio galaxy morphology pre-training. In order to account for the more nuanced variations in radio galaxy morphology than are typically included in the augmented views of view-based SSL algorithms, we use a pre-trained Variational Autoencoder (VAE) to generate views for training a larger view-based self-supervised model. To do this, a $\beta$-VAE was trained on the Radio Galaxy Zoo (RGZ) dataset, where moderate regularization ($\beta = 2.3$) was found to provide a good balance between reconstruction quality and disentanglement of generative factors such as source multiplicity and lobe asymmetry. An analysis of the $\beta$-VAE reveals that Fanaroff-Riley class identity manifests as a continuous transition across the latent space, rather than being associated to a single discrete dimension. $\beta$-VAE reconstructions were then incorporated as generative augmentations within a view-based SSL pipeline. Our experiments show that combining these generative views with standard image augmentations improves downstream classification performance, and we present ablation studies clarifying the relative contribution of each augmentation type. These results indicate that generative and contrastive approaches are complementary, and point toward disentanglement-aware self-supervised learning as a promising direction for future radio astronomy surveys.

The ngVLA has adopted a standalone system with a $\sim$~1.5 m diameter antenna and associated ambient temperature receiver and electronics at each antenna as its baseline design for water vapor radiometer (WVR). The WVR is intended to decrease calibration overheads for tropospheric phase correction for the high frequency bands (5 \& 6; 30.5-116 GHz) to levels lower than feasible with fast switched reference gain calibration. In this memo, we present an alternative concept utilizing the standard Band-4 science receivers and the main antennas. This would require the Band-4 feed to be placed between those of Bands 5 \& 6 resulting in a beam offset of 0.8 deg, similar to the beam size of the standalone 1.5 m antenna. The main benefits of the Band-4 WVR concept are: (1) sensing water vapor in a region more representative of the Band 5 \& 6 science beams in the lower layers of the troposphere, in the near-field of the 18 m antenna, than the 1.5 m antenna (2) significantly relaxing the stability requirements due to the lower receiver temperature of the standard Band-4 cryogenic receiver and (3) eliminating the procurement, deployment, and maintenance of a large number of standalone units. Given these advantages, it is recommended that the details and the pros and cons of Band-4 WVR implementation compared a standalone system be carefully considered.

Brown dwarfs and exoplanets are thought to host complex atmospheric phenomena such as clouds, storms, and chemical heterogeneity, akin to weather patterns on Earth. These features can produce pronounced spectral variability. Time-variability monitoring provides a unique window into surface inhomogeneities that cannot be directly resolved with foreseeable imaging technology. Current time-series analysis techniques have provided qualitative constraints on variability mechanisms but lack the ability to quantitatively estimate the extent of variation on atmospheric properties. We present Tempawral, the first data-driven time-resolved atmospheric retrieval framework that quantitatively retrieving variability in atmospheric parameters via an eigen-spectra inversion technique, leveraging the full spectra dataset. We validate this method on simulated time-series spectra of a variable brown dwarf, demonstrating that it successfully recovers key variability drivers, including inhomogeneous cloud coverage, evolving chemical abundances, and changes in temperature structure. We further showcase the utility of Tempawral by applying it to JWST/NIRISS-SOSS time-series observations of a highly variable T2.5 brown dwarf. The observed variability is best explained by a $\sim$300 K temperature perturbation near the 1-bar, accompanied by variations in the abundances of $\rm H_{2}O$, $\rm CO$, $\rm FeH$, as well as changes in the thickness of the iron cloud deck. This work provides a generalized framework for time-resolved atmospheric retrievals in the JWST era, enabling comprehensive interpretations of dynamic atmospheric processes in substellar objects.

We report on a spectroscopic study of the bright, nearby type-IIn supernova SN 2017hcc, and SN 2023usc using data obtained from the Himalayan Chandra Telescope (HCT). SN 2017hcc is well-studied event, and our sampling covers 7 epochs, starting from +14\,d post explosion, and continuing into the nebular stage, at +411\,d. The type-IIn event SN 2023usc was sampled over 5 epochs from +12\,d to +155\,d post explosion. The nearly featureless (except H$\alpha$) late time (+62\,d onward) spectra of SN 2023usc, suggests a novel explosion route for this type-IIn event. Assuming a CSM model created by multi-epoch ejection of material from the pre-explosion progenitor, we present here a comparative study of both events with several other type-IIn / interacting supernovae in progenitors with persistent signatures of a CSM. We find that true narrow lines ($v \ll 1000$\,km\,s$^{-1}$) emerge in the early ($\sim$ +10\,d) spectra only in few events (SN 2017hcc, SN 2023usc and SN 2010jl) initially classified as type-IIn in our sample -- in most cases the line velocity hovers at $\sim 1000$\,km\,s$^{-1}$ even in the very early epochs. CSM line velocity being indicative of its extent and opacity, this suggests that progenitors with a highly extended CSMs, which are also optically transparent in their outer regions may be relatively rare.

We propose that the rapid cooling of the neutron star following its X-ray superburst in MAXI J1752$-$457 over a period of 4 days, observed by two Japanese satellites, MAXI and NinjaSat, is due to enhanced neutrino emission from cycles of electron capture and $\beta^{-}$ decay involving odd-$A$ nuclei (or Urca pairs) in the ocean. Hence, this work provides the first indication of the possible existence of such a ``nuclear Urca process". The observation of MAXI J1752$-$457 implies a hot ignition layer with a maximum temperature of $\sim4~{\rm GK}$, located near the Urca shell in the ocean, such that the nuclear Urca process becomes dominant for up to $\sim2$ days after the superburst. This behavior is distinct from that of normal Type-I X-ray bursts, which are triggered by hydrogen or helium burning in much shallower layers than those of superbursts. Our findings enable probing of superburst ashes through Urca pairs via long-term monitoring of crust cooling on day-long timescales.

We analyze the redshift distribution of high-$z$ galaxies and active galactic nuclei identified in early JWST data, and investigate the presence of periodic structure in the variable $x=\ln(1+z)$. A baseline-corrected unbinned frequency analysis reveals a statistically significant peak corresponding to a spacing $\Delta x \simeq 0.34$, suggesting an approximately log-periodic pattern in the redshift distribution. A periodic structure in $x$ implies a preferred scaling ratio in $(1+z)$, which may be interpreted as a realization of discrete scale invariance. We discuss the possibility that such behavior arises from cyclic condensate dynamics in Bose--Einstein condensate (BEC) cosmology. In the Fukuyama--Morikawa--Tatekawa framework, repeated collapse and re-formation episodes of a self-interacting condensate occur over characteristic timescales of order several $10^8$ years. When mapped into redshift space, this temporal periodicity naturally translates into an approximately constant spacing in $\ln(1+z)$. While the observed frequency is not interpreted as a sharp theoretical prediction, its magnitude is quantitatively consistent with the intrinsic cycle timescale of QCD-axion motivated condensate dynamics. The present analysis therefore provides observational support for cyclic BEC cosmology as a viable dynamical origin of log-periodic structure in the high-redshift universe.

We present a comprehensive study of host galaxies of radio sources within the 1.35$R_{200}$ of the Coma cluster by combining deep 144MHz observations from the LOFAR Two-Metre Sky Survey (LoTSS-DR2) with optical spectroscopy and photometry from DESI and SDSS. We identify 79 spectroscopically confirmed cluster members with reliable radio emission and classify them into compact, extended, and tailed subsamples according to their radio morphologies. By combining their radio and optical properties, we find compact radio sources are predominantly associated with massive, quiescent galaxies driven by AGN activity, while tailed sources are largely hosted by star-forming galaxies, tracing ongoing ram pressure stripping (RPS). Using phase-space analysis and a projected infall time proxy ($d_R$), we find that extended sources are preferentially located in the cluster outskirts ($d_R > 1$), while tailed sources are concentrated in the intermediate infall region ($0.4 < d_R < 1.0$), highlighting the influence of the dense intracluster medium.

The generation of a baryon-antibaryon asymmetry in the Universe via gravitational baryogenesis is investigated for two f(R) modified theories of gravity, the widely used Starobinsky $f(R)=R+R^{2}/M^{2}$ model, and the recently proposed power-law model $f(R)=c_{1}R^{2+k/4}+c_{2}R+c_{3}$ of Odintsov and Oikonomou (2025) that is constructed from the slow-roll inflation parameters, and fits the new high-multipole CMB observations reported by the Planck and ACT collaborations for $k \sim -0.03$. The present investigation is undertaken in the Einstein frame. The motion of the scalaron is studied for the slow roll inflationary era using analytic approximate solutions obtained from its potential and, from these solutions, analytic expressions for the Ricci scalar, its time derivative, the Hubble parameter and the scale factor of the Universe. These expressions for the Starobinsky model were obtained first by Motohashi and Nishizawa (2012) but the expressions for the power-law model of Odintsov and Oikonomou and its generalization are new. The calculated values of the baryon asymmetry factor $\eta $ for the Starobinsky model vary from $(1.05 - 1.46) \times 10^{-11}$, and for the power-law model, from $(1.06 - 1.53) \times 10^{-11}$. The power-law values depend upon the unknown fitting parameters $c_{1}$ and $c_{2}$, and a future fit of the Odintsov and Oikonomou model to data could yield enhanced values. The values for $\eta$ depend upon a mass parameter $M_{\ast}$ which is expected to be of the order of $M_{Pl}=2.435 \times 10^{18}$ GeV. The values of $\eta$ for both models have been calculated for $M_{\ast}=M_{Pl}$ and are quite close to the observed value $\eta = 8.65 \times 10^{-11}$. Since $\eta \propto (M/M_{\ast})^{2}$, reducing $M_{\ast}$ slightly from $M_{Pl}$ to $0.4 M_{Pl}$ would bring the calculated values into agreement with the observed value.

In this paper, we investigate a quintessence field with an exponential potential motivated by the suggestion of time-varying dark energy from the DESI galaxy survey. Assuming a kination epoch in the early Universe, we analytically derive constraints on initial conditions that are consistent with Big Bang Nucleosynthesis and the current dark energy density. Compared to the severe 120-digit fine-tuning required for dark energy to be a cosmological constant, our result suggests that the degree of fine-tuning is naturally relaxed by dozens of orders of magnitude. Furthermore, we discuss the method for testing this model through future observations of the gravitational wave background.

Previous studies based on the latest realisation of the International Celestial Reference Frame (ICRF3) have suggested a correlation between astrometric properties (such as the radio-optical offset) and redshift for active galactic nuclei (AGNs). We extend these investigations by using a large, all-sky sample of approximately 22,000 compact radio sources from the Radio Fundamental catalogue (RFC) to examine this relationship in a systematic and statistically robust manner. We compiled redshifts for about 10,000 RFC sources over the range 0 < z < 5 by combining data from the Dark Energy Spectroscopic Instrument Data Release 1 and the Sloan Digital Sky Survey Data Release 17/19 with additional datasets from the NASA/IPAC Extragalactic Database. Cross-matching with Gaia Data Release 3 yielded a clean sample of 4,068 RFC objects with reliable spectroscopic redshifts and classifications (galaxies and quasi-stellar objects; QSOs). We analysed the redshift dependence of their radio astrometric properties from very long baseline interferometry (VLBI) and their optical astrometric properties from Gaia. We find that the VLBI astrometric properties show no significant dependence on redshift within the achieved level of precision. In contrast, several optical astrometric quantities exhibit clear redshift-dependent behaviour. The median absolute radio-optical offsets decrease markedly over 0 < z < 0.5, where galaxies dominate the sample, decline more gradually over 0.5 < z < 1.3, and exhibit a mild increase at z > 1.3, where QSOs dominate. Similar behaviour is observed for several Gaia astrometric quantities, including astrometric uncertainties, proper motions, and G magnitudes. These behaviours can be largely explained by the dependence of Gaia astrometric performance on G magnitude and by the evolution of the G magnitude with redshift.

RCW 86 is a supernova remnant whose origin has recently been linked to an off-center explosion within a cavity created by its progenitor star. In the southwestern region, the forward shock is thought to have reached the cavity wall, encountering diverse environmental conditions. We report on the spatially resolved X-ray polarimetric observation of RCW 86 with the Imaging X-ray Polarimetry Explorer (IXPE). In the 2--4.5 keV energy band we find no significant detection of polarization. Employing a dedicated background subtraction procedure and Bayesian spectropolarimetric fitting, we derive 99% upper limits on the polarization degree of the synchrotron component: 15% in higher-statistics regions and 30%--40% in lower-statistics regions. These upper limits on the polarization degree in several regions exclude the possibility of a strongly coherent magnetic field down to the subparsec scale, and that of a moderately coherent one on the scale of the synchrotron features as resolved by IXPE. The results indicate that the shocks in the southwestern rim of RCW 86 propagate more slowly than the unshocked ejecta at their locations, yet exceed the measured proper motion speeds. This behavior is consistent with reflected shocks occurring in tenuous regions of the shocked ejecta, distinct from regions that are radio-bright.

Observations of X-ray filaments attached to a couple of powerful pulsars suggest escape of TeV electrons and/or positrons (e$^{\pm}$) from pulsar bow shocks into surrounding large scale magnetic fields. These filaments are usually asymmetric with very weak emission from the other side of the main filaments, and no significant spectral variation has been detected across these filaments, implying inefficient energy loss of emitting particles. We develop a Monte Carlo code to simulate particle transport in a large scale magnetic field and apply the model to PSR B2224+4415 (Guitar). It is shown that, with an injection power of a few tens of percent of the pulsar spin down luminosity, TeV e$^{\pm}$ can explain the observed filament properties with a scattering mean free path along the magnetic field comparable to the length of the observed filament. The model predicts a dim diffuse symmetric X-ray background aligned with the filament on a larger scale, whose flux is proportional to the X-ray emitting e$^{\pm}$ energy loss time for a stable e$^{\pm}$ injection power comparable to the luminosity of this diffuse background. Observations with a large field of view and good sensitivity should be able to detect such a component.

Blazars exhibit variable emission on diverse timescales. Some light curves show signs of quasiperiodic oscillations (QPOs), which may encode clues regarding the physical processes behind the emission or point to supermassive binary black holes. We analyzed five blazars with previously reported high significance year-long QPOs, applying the Lomb-Scargle periodogram and Weighted Wavelet Z-transform methods to Fermi-LAT data up to early 2025. We furthermore examined an additional source (PKS 0139-09), where nascent QPO may be present. As the light curves showed longer term trends, we detrended the data using an STL decomposition, which often revealed a large seasonal component. We find that detrending generally leads to an increase in the strength of the QPO signal. However, except for PG 1553+113, where a clearly persistent QPO signal is present, we detect transience on a timescale of $\lesssim$4000 days. We then forecast the light curves over the following four years, using a traditional statistical method as well as a Transformer-based deep learning model. Applied to a test set, the latter showed significant success in predicting behavior that seems unexpected from simple inspection of the past data. Analyzing the extended time series suggests a markedly weaker QPO signals over the coming years in cases where the transient behavior appears near the end of the observational data. In contrast, in the nascent candidate QPO source (PKS 0139-09) the signal is expected to strengthen significantly. These predictions, which may reflect the physical origin of the QPOs, can be tested against future data.

Bulges and disks are major structural components that define galaxy morphology. The mass ratios of bulges and disks increase statistically with the galaxy mass, with the high-mass end occupied by elliptical galaxies. Although previous theoretical studies have succeeded in reproducing this morphological sequence, it is not yet fully understood why and how this morphological sequence emerged. Galaxy mergers accompanying dark matter halo mergers have been proposed as the major route for bulge formation. On the other hand, it is observationally known that the mass fraction of galaxies (stars plus cold gas) in dark matter halos attains the peak value at $M_{\rm halo} \sim 10^{12} {\rm M}_\odot$ throughout the cosmic time. Using a simple galaxy evolution model including mergers, we show that this feature is the fundamental cause of the morphological sequence. Halos hosting massive galaxies, which stay more massive than this peak mass for long periods during their growth, merge mostly with satellite halos having larger galaxy mass fractions than themselves. Such mergers increase the bulge mass fraction efficiently. In contrast, host halos of low-mass galaxies evolve under unfavorable condition for bulge growth because they stay below the peak mass and merge with satellite halos with smaller galaxy mass fractions. Previous studies suggest that the peak in galaxy mass fraction is created by feedback processes from active galactic nuclei (AGN) and young massive stars including supernovae (SN), which are considered to suppress star formation in high-mass and low-mass galaxies, respectively. This study thus points to a close relationship between the galaxy morphology and feedback processes which have hitherto been considered unrelated and suggests the importance of further investigation into their causal relationship.

For observations of supernovae, kilonovae, tidal disruption events, and other non-repeatable observations, it is important the science data is taken successfully within a specific time window. Part of obtaining that data is often centering objects in the aperture to a higher accuracy than is available from Hubble Space Telescope's (HST's) blind pointing. On HST STIS the sequence of exposures responsible for this centering is the target acquisition or STIS ACQ sequence, and it is most often placed only at the beginning of a visit. Unfortunately, STIS ACQ sequences will fail if the observatory experiences issues locating guide stars in time for the start of the required exposures. If the guide stars are located at a later point in the visit, the remaining science exposures can be taken but the pointing might not be as accurate as is required. This work discusses both the frequency of this issue and the feasibility of placing redundant or "safety" STIS ACQ sequences in a multi-orbit visit to regain the desired pointing accuracy in an affected visit. To do so we select a subset of all 113 STIS ACQ sequences from September 2018 to September 2023 which have experienced this issue. We find that this problem occurs in ~5% of the total STIS ACQ sequences taken during that time period, with a recent increase in the rate to ~9% from March to September 2023. Since the observatory goes through periods of better or worse pointing performance, this recent increased failure rate is not guaranteed to continue. For those failed visits which span multiple orbits, ~39% never obtain a lock on the guide stars and thus take no data. Of the multi-orbit visits that do recover the guide stars, the majority (~78%) do so before the beginning of science exposures in the second orbit. We also provide advice for users on how to make a risk assessment based on the analysis presented here.

Parameter degeneracy in blazar spectral energy distributions (SEDs) is known but rarely quantified. This paper introduces a Fisher Information approach to determine theoretical limits to information extraction in the context of one-zone models. By evaluating the total Fisher Information by varying $\delta$, $B$, $p$, $\gamma_{\rm min}$ and $\gamma_{\rm max}$, we find that EC models encode Fisher information $\gtrsim10^4$ times less than that in SSC models, establishing differences in limits of physical information extraction even in the case of perfect sampling. Moreover, the Fisher information in both SSC and EC models exhibit strong fluctuations across the parameter space, but since the magnitudes are orders of magnitude lower in EC, limits of parameter inference are expected to be worse in FSRQ SEDs than BL Lacs. We also find that the Doppler factor $\delta$ carries at least $10^{2-3}$ more Fisher information than that for $p$ and $B$ in both EC and SSC, making $\delta$ the most constrained SED parameter. Applying our Fisher Information motivated framework to real flaring SEDs of Flat Spectrum Radio Quasars (FSRQs) CTA 102 and 3C 279, we show that mild variations in $\delta$ and $p$ can appreciably produce the flaring SEDs starting from the quiescent model, while two other flares in 3C 279 simple geometric and spectral considerations cannot reproduce the flares, reducing the efficacy of one-zone models. We propose that time-resolved SED models are indispensable to constraining physical parameters in EC-dominated blazars.

Quantifying the maximum energy of fast radio bursts (FRBs) can provide stringent constraints on their emission mechanisms and progenitor models. However, the most energetic bursts are rare, requiring a large sample of FRBs to detect them. In this work, we use the largest available such sample, 2,998 one-off FRBs from the Second CHIME/FRB Catalog, to obtain a lower limit on the maximum energy ($E^{\mathrm{max}}_{\mathrm{iso}}$) of FRBs, assuming isotropic energy distribution from FRB sources. In the absence of known redshifts ($z$) for most sources, we present a framework that uses the dispersion measures (DMs) and fluences of these FRBs, together with the probability distribution of $z$ given DM, to derive the lower limit on $E^{\mathrm{max}}_{\mathrm{iso}}$. We generate simulated FRB samples assuming different parameter values for a log-normal $\mathrm{DM}_{\mathrm{host}}$ distribution and a Schechter function form of the FRB energy function to estimate how many outliers -- FRBs with large DM contributions from the host galaxy or intervening galaxy halos -- could artificially inflate this limit. After accounting for outliers, the lower limit on $E^{\mathrm{max}}_{\mathrm{iso}}$ from Catalog 2 FRBs ranges between $1.2\times10^{41}$ and $1.9\times10^{42}$ erg, with best estimate $1.2\times10^{42}$ erg. This limit is consistent with those derived from much smaller FRB samples. Moreover, inferred energies of hundreds of FRBs appear collectively limited around $\sim10^{42}$ erg, suggesting a physical limit on the energy reservoir of FRB sources. The corresponding isotropic-equivalent FRB source energy is consistent with the total energy available in a magnetar's external dipole magnetic field, supporting magnetars as FRB progenitors.

We present a multitask surrogate for neutron-star equations of state (EoSs) that delivers \emph{distribution-free}, certified uncertainty via split conformal prediction (CP) and its Mondrian variant. The surrogate ingests a six-parameter piecewise-polytropic representation $(\log_{10}p_1,\Gamma_1,\Gamma_2,\Gamma_3,\rho_1,\rho_2)$ -- with fixed transition densities $\rho_1$ and $\rho_2$ -- and jointly performs (i) validity classification under physical/observational constraints and (ii) regression of $M_{\max}$, $R(M_{\max})$, $R_{1.4}$, and $\Lambda_{1.4}$. Trained on a balanced set of $40{,}000$ EoSs, the model attains near-perfect discrimination (AUC $\approx 0.997$) and sub-percent relative errors for masses and radii, with few-percent error for tidal deformability. Across $\alpha\in[0.05,0.25]$, empirical coverages closely track $1-\alpha$ for both Standard and Mondrian CP; in conservative regimes, Mondrian yields narrower average physical widths at comparable coverage. To our knowledge, this is the first application of class-conditioned (Mondrian) conformal calibration to neutron-star EoS surrogates, enabling efficient, reproducible, and uncertainty-aware inference; the framework is readily extensible to functional targets (e.g., full $R(M)$ curves).

Fast sampling photometry is essential for characterising fireballs and their fragmentation episodes which link to the meteoroid internal structure. Accurate measurements remain challenging due to the large required dynamic range of up to 10 stellar magnitudes driving up operational complexity and cost. We developed an all-sky camera system operating at up to 500 frames per second featuring a novel Detection Localised Auto-brightness Control. Custom software manages high data throughput via transient detection and region-of-interest saving with real-time photometry. Two field deployments validate photometric accuracy against conventional 30 frames per second cameras and demonstrate the successful capture of a bright magnitude -15 fireball with minimal saturation. The system achieves an effective dynamic range between apparent magnitudes -3 and -17 capturing minimally saturated light curves for most fireballs. A successful semi-empirical fragmentation analysis verifies its ability to provide data for detailed physical modelling. The primary application for this validated system will be as a core component of the Global Fireball Observatory's next-generation instrumentation. The intention is to deploy it in a hybrid observatory, operating alongside a dedicated high-resolution astrometric camera. This configuration will allow the network to simultaneously capture precise trajectory data for orbit and fall-line calculations and acquire complete, unsaturated high dynamic range light curves at high temporal resolution for detailed physical analysis, combining the strengths of both systems.

The China Space Station Telescope (CSST) will conduct wide-field multiband photometric imaging and slitless spectroscopic surveys, advancing cosmology and galaxy evolution studies. Achieving CSST's cosmological goals requires precise redshifts ($\sigma_{\rm NMAD}\lesssim 0.002-0.005$) from low-resolution ($R\sim200$) and potentially blended slitless spectra. We present BayeSED3, extended for Bayesian full-spectrum analysis, including nebular emission modeling (via \textsc{Cloudy}) and a Bayesian treatment of the model scaling factor, improving reliability over optimization methods for low SNR spectra. Validated on realistic mock data generated with the CESS emulator (median SNR=1.65, including instrumental and self-blending effects), our method achieves excellent redshift precision with three-band (GU+GV+GI) spectroscopy: $\sigma_{\rm NMAD}=0.0008$ ($\sim$80% success) for star-forming and $\sigma_{\rm NMAD}=0.0015$ ($\sim$50% success) for quiescent galaxies. Stellar mass ($\sigma_{\rm NMAD}\approx0.015$ dex for SF, $\approx0.016$ dex for quiescent) and SFR ($\sigma_{\rm NMAD}\approx0.05$ dex for SF, especially at SNR>1) are reliably recovered. Self-blending increases scatter by $\gtrsim30%$, but combining spectroscopy with CSST's seven-band photometry significantly improves accuracy, especially for quiescent galaxies and data-limited cases. Single-band spectroscopy plus photometry yields reasonable redshifts: GU+photometry is limited, GI+photometry gives >60% (SF) and >40% (quiescent) success at $\sigma_{\rm NMAD}\lesssim0.002$, GV+photometry gives >35% (SF) and $\sim$40% (quiescent) at similar precision. The Bayesian framework offers a powerful method for accurate galaxy characterization, enhancing CSST's scientific outcomes despite the challenges of slitless spectroscopy.

We investigate whether the long photometric cycles observed in double-periodic variables (DPVs) can arise from nodal precession of a tilted accretion disk driven by the tidal torque of the companion. Within a simple analytical framework, we derive testable relations linking the long-to-orbital period ratio to the binary mass ratio, the normalized disk size, and the disk tilt angle $\beta$, which itself can be inferred from the long-cycle amplitude, orbital inclination $i$, and disk luminosity fraction. The model naturally reproduces the two observed long-cycle light-curve morphologies -- sinusoidal and double-hump -- distinguished by the geometric criterion $i+\beta \le 90^\circ$ versus $i+\beta>90^\circ$. Applying these relations to a sample of DPVs, we find that the inferred disk sizes are physically reasonable and consistent with independent light-curve modeling for a non-negligible subset of systems. Our results show that tidal nodal precession represents a viable and potentially important contributor to the long-period variability of DPVs and provide a quantitative framework for future observational and theoretical studies.

The Emirates Mission to the Asteroid Belt (EMA) will study the ultra-red asteroid (269)~Justitia. In this work, we present the first detailed investigation of Justitia's surface dynamics using the newly developed 3-D polyhedral shape model with 574 vertices and 1,144 faces. We analyze its geopotential, surface acceleration, escape speeds, slopes, and equilibrium points, and we also search for planar symmetric periodic orbits using an equivalent ellipsoidal approximation. Our results indicate that the lowest geopotential values occur at the poles, which also correspond to regions of maximum surface acceleration. The global slope distribution suggests preferred zones of material accumulation or migration, offering clues to Justitia's long-term morphological evolution. Most slopes remain below $40^\circ$, implying that loose particles may settle stably across large portions of the surface. We identify five equilibrium points consistent with Justitia's estimated density and slow rotational period. Two external points (E$_2$ and E$_4$) exhibit linear stability, and all equilibrium locations lie relatively far from the surface due to the body's slow spin. Additionally, we discover 28 new families of planar symmetric periodic orbits, then classify their topologies and determine their linear stability, providing a dynamical framework relevant to future spacecraft operations near Justitia. Finally, thermal modeling reveals how thermal inertia and heliocentric distance shape Justitia's temperature distribution. The south pole receives more insolation than the north pole, reaching minimum temperatures of about 102~K and 87~K, respectively. These combined dynamical and thermal results offer valuable insights for the EMA mission and for understanding slowly rotating small bodies.

The hydrostatic mass bias ($b_{\mathrm{HSE}}$) is a leading systematic uncertainty in cluster cosmology and a principal source of degeneracy with $\sigma_8$ and $\Omega_m$. We investigate the capability of Stage-IV CMB and optical surveys to calibrate $b_{\mathrm{HSE}}$ using tomographic cross-correlations between the thermal Sunyaev--Zel'dovich (tSZ) effect, galaxy clustering, and weak lensing. We perform a Fisher forecast incorporating realistic survey noise, foreground modeling for clustered CIB and radio sources, and full marginalization over cosmological and astrophysical nuisance parameters, including per-bin galaxy bias perturbations, photometric redshift shifts, intrinsic alignments, and baryonic feedback modeled with HMCode2020. With optimized tomographic binning (10 lens and 5 source bins for LSST; 6 lens and 5 source bins for CSST), we forecast marginalized constraints of $0.98\%$ for SO+LSST, $1.60\%$ for CMB-S4+LSST, and $2.40\%$ for CMB-S4+CSST. Tomography improves $b_{\mathrm{HSE}}$ precision by factors of approximately three relative to non-tomographic analyses, reflecting the role of redshift information in breaking the $b_{\mathrm{HSE}}$--$\sigma_8$ degeneracy. Optical-only probes provide no direct constraint on $b_{\mathrm{HSE}}$, whereas inclusion of tSZ-containing spectra enables percent-level calibration under realistic systematic assumptions. The results demonstrate that multi-probe tomographic analyses with Stage-IV surveys can achieve robust control of hydrostatic mass bias, strengthening cluster-based constraints on structure growth.

Active Galactic Nuclei (AGNs) are characterized by strong temporal flux density variability across the electromagnetic spectrum, offering insights into the complex physical processes governing accretion and plasma outflows. To systematically investigate AGNs flux density variability in radio bands, a long-term program was initiated in late 2016: quasi-Simultaneous Multiwavelength Monitoring of gamma-ray-loud AGNs with the Nanshan 26-m radio telescope (SMMAN). This work presents the first data release of the SMMAN program, spanning over eight years from 2016 to 2024 with observations at 4.8 and 23.6 GHz. The SMMAN sample includes 131 northern ($\delta >\sim0^{\circ}$) sources selected from the Fermi Large Area Telescope third source catalog. The characteristics of variability, spectral index, luminosity, and $\gamma$-ray loudness factor are examined for different AGN classes within the sample. Target sources exhibit stronger variability at 23.6 GHz compared to 4.8 GHz, with BL Lac objects being more variable than flat-spectrum radio quasars (FSRQs). BL Lacs generally have flatter radio spectra, while FSRQs, blazar candidates of uncertain type (BCUs), and radio galaxies (RDGs) span a wider range from flat to steep. FSRQs are more radio-luminous than BL Lacs and other classes, with BCUs intermediate and RDGs generally fainter. FSRQs and BL Lacs have higher $\gamma$-ray loudness factors than RDGs, while BCUs have intermediate values. The SMMAN dataset, incorporated with other historical and ongoing monitoring programs, will provide a unique opportunity to investigate the evolution of spectral energy distributions, search for quasi-periodic oscillations, and analyze supermassive black hole binary systems.

Timing observation of pulsars orbiting around a supermassive black hole (SMBH) can measure the spacetime around the SMBH to a high precision and thus be a novel probe of the gravity theory. Future high-frequency surveys of the Galactic Centre (GC) region to be performed by the next-generation radio telescopes, such as the SKA, may discover pulsars that orbit around Sagittarius A* (Sgr A*), the SMBH dwelling in our GC. In this paper, we present a realistic pulsar-SMBH timing model based on the post-Newtonian equations of motion of the pulsar. Considering the expected timing precision in the future, we take into account several next-to-leading order light propagation time delays in the timing model. For the first time, we include the effects of proper motion of Sgr A*, which were expected to break the spin measurement degeneracy. We forecast the measurement precision of various parameters of Sgr A*, and discuss the data analysis procedure in the presence of red noise, which can be strong if the pulsar is a normal pulsar. The realistic timing model constructed in this study will serve as a useful tool in future searching and timing of pulsar-SMBH systems in the GC.

The Galactic Centre black hole provides a naive environment for understanding unknown matter distribution and new gravitational physics. For this stellar orbits in the nuclear star cluster are reliable probes. We investigate different dark matter mass profiles through pericentre shift of stellar orbits near the black hole. We also study capability of existing and upcoming astrometric facilities to detect dark matter induced precession and to distinguish between several dark matter profiles. Parameters of different dark matter density profiles are estimated by using the most recent upper bound on dark mass near the black hole. These profiles are then used for calculating the gravitational potential and hence the relativistic pericentre shift of both low and high eccentricity orbits of 13 S-stars. We use the recently measured deviation parameter $f_{sp}$ for investigating competition between dark matter and gravitational physics within S2's orbit. The astrometric shift of the pericentres has been calculated and compared with existing and upcoming astrometric capabilities of large and extremely large telescopes. The orbit of S2 is found to be insensitive to dark matter induced precession. Low eccentricity and wider orbits are prominent probes for measuring dark matter induced precession which is accessible to present and upcoming astrometric facilities such as Keck, GRAVITY and TMT. The existing and upcoming facilities can distinguish between different dark matter profiles for some stars and hence they posses the capability to distinguish between possible formation histories of the central region of our Galaxy.

Understanding the properties of strongly interacting matter at extreme densities is a central problem in fundamental physics, but neutron star mergers provide a natural laboratory for probing this regime. However, the complexity of the merger process complicates the interpretation of the associated gravitational-wave and electromagnetic signals. This picture becomes even more complex in the potential scenario in which dark matter accumulates around and in neutron stars, altering their structure and the associated observables. In this work, we study synthetic gravitational-wave observations of binary neutron star mergers with next-generation detectors, investigating their potential to extract both nuclear physics and dark-matter constraints. We also examine how the potential presence of fermionic, non-interacting dark matter inside neutron stars affects the inference of nuclear empirical parameters. We find that combining observations can tighten constraints on nuclear empirical parameters. However, the inferred values remain sensitive to systematic modeling biases and intrinsic degeneracies among the parameters. Conversely, our analysis reveals that even in the presence of dark matter, it will be unlikely to find decisive evidence for dark matter when analyzing gravitational-wave signals. Consequently, systematic biases in nuclear empirical parameter inference potentially resulting from the presence of dark matter are expected to be negligible even for observations with next-generation gravitational-wave detectors.

Dwarf galaxies are important laboratories for studying cosmic magnetism because they can maintain strong magnetic fields via the action of turbulent dynamo despite their low mass and weak gravitational potential. The Magnetic-field Evolution in Dwarf galaxies from Ultra-deep SKA Analysis (MEDUSA) survey is the first SKA-pathfinder programme designed to obtain deep continuum, polarisation, and HI data for dwarf galaxies, enabling a comprehensive study of their radio spectra, magnetic fields, and gas kinematics across a representative population. By analysing the radio continuum spectra and polarisation of the dwarf-dwarf galaxy merger NGC 1487 from the MEDUSA sample, we aim to determine its magnetic field strength and to characterise the large-scale and turbulent components of its magnetic field. We utilise highly sensitive multi-band radio continuum data from MeerKAT L-band (1.28 GHz) and Australia Telescope Compact Array (ATCA) L/S (2.1 GHz), C (5.5 GHz), and X-bands (9 GHz). We analysed the magnetic field configuration using polarisation and rotation measure (RM) synthesis. The integrated spectral energy distribution has a non-thermal spectral index of $\alpha_{\rm nth} = -0.678\pm0.085$. Synchrotron and inverse Compton losses cause a spectral break at $\nu_{\rm b} = 6.2\pm1.3$ GHz. In star-forming regions, the magnetic field exhibits strong small-scale fluctuations in RM, suggesting the action of a small-scale dynamo. Conversely, the field becomes more ordered, aligning with the tidal arms toward the galaxy's outskirts, showing a large-scale magnetic field over $\approx3$ kpc. Observations of the dwarf-dwarf merger NGC 1487 show that even low-mass galaxy mergers, likely the building blocks of larger galaxies in the early Universe, can rapidly amplify and produce coherent large-scale magnetic field structures, highlighting their key contribution in the early magnetisation of galaxies.

Recent studies show that the Milky Way stellar halo is composed of populations of different origins, shaped by multiple accretion events. To better understand the formation of the Milky Way and other spiral galaxies, we characterize the chemical and kinematic properties of halo substructures using GALAH DR4 and Gaia data. We apply wavelet transforms in the space of sqrt(J_r) and azimuthal action (L_z) to identify kinematic overdensities. Stars in the detected structures are analyzed in elemental abundance space to determine their origin. We further assess contamination using the unsupervised machine-learning algorithm t-distributed stochastic neighbor embedding (t-SNE), performing chemical tagging with 15 elemental abundances. We recover five structures: the Galactic disk, the Splash, Gaia-Sausage-Enceladus (GSE), Thamnos1, and Thamnos2. GSE shows two peaks; one at sqrt(J_r) ~ 25 kpc km s^-1 is due to disk contamination, while the other above sqrt(J_r) ~ 40 kpc km s^-1 represents the cleanest GSE population. Thamnos exhibits three peaks linked to Thamnos1 and Thamnos2. Thamnos2 shows higher [alpha/Fe], iron-peak elements are enhanced in the Splash, and halo groups retain a stronger r-process signature. The multiply peaked structures suggest that the splashed disk extends beyond prograde orbits. The distinct chemo-dynamical properties of the halo groups support their extragalactic origin.

The main belt, the region between the orbits of Mars and Jupiter, is home to more than 1 million asteroids. These asteroids form orbital groups, (i.e., asteroid families formed by collisions) and also spectral groups (taxonomies) with different chemical compositions, in particular carbonaceous (C-types) and silicate (S-types). In this paper, we extend the existing main-belt collisional model by finding the appropriate strength-versus-size dependence (also known as the scaling law) for these two groups. We used color indices and geometric albedos of 56 and 72 spectroscopically confirmed C- and S-types (control samples), along with statistical methods on 1 065 034 asteroids, to assign C-, S-, or other types. This allowed us to construct observed size-frequency distributions (SFDs) for several subpopulations constrained by either semimajor axis (inner, middle, outer) or taxonomy (C, S, other). Then we used a Monte Carlo collisional model to compute the long-term collisional evolution (4.5 billion years) and derive synthetic SFDs. Our best-fit scaling laws indicate that S-types must be weaker below approximately 0.2 km than C-types to explain the deficiency of asteroids in the inner part of the main belt near (and below) the observational limit. This may correspond to differences in chemical composition or material porosity. Future research will focus on the scaling laws of asteroids with rare or "extreme" taxonomies (e.g., V, M).

Contemporary exoplanet and brown dwarf atmospheric research relies heavily on retrieval frameworks to recover thermal and chemical properties and perform model comparison in an observational data-driven approach. However, the computational effort required for retrieval modelling has rapidly increased, driven by JWST data that covers large spectral intervals at moderate spectral resolutions, and ground-based, high-resolution spectroscopy. To help tackle the computational burden faced by contemporary retrieval requirements, I present a new sub-stellar atmosphere retrieval modelling framework, Exo Skryer, that utilises the JAX library for Python to enable scalable, computationally efficient forward modelling as well as posterior sampling. I present example retrievals for pre- and current JWST era observations for both transmission and emission spectra, finding consistent results to previous retrieval modelling efforts, apart from a WASP-107b test case. In addition, I present a new method to directly retrieve the real and imaginary optical constants (n, k) of suspected aerosol infrared absorption features. Due to its computational expediency, Exo Skryer will be highly suited for future demanding retrieval efforts that incorporate more spatial dimensionality, complex forward models and high-dimensional parameter sets.

Loop-aligned hydrodynamic modelings help better understand the thermodynamic evolution of flaring plasma confined in solar flare loops. Conventional loop modelings typically assume a uniform loop cross section. With a variation of the cross section taken into account, in this work we carry out both analytical and numerical modelings of the radiative cooling in a solar flare loop. It is found that a cross-sectional expansion with height can efficiently suppress the draining of loop material from the corona while not significantly affecting the decrease of loop temperature. Reflected to the loop energetics, the coronal part of the loop cools more dominantly by radiation, and more importantly, the loop radiative outputs are shifted toward lower temperatures. These findings pose important physical implications for extreme-ultraviolet (EUV) late-phase emissions discovered in some solar flares. The late-phase loops in these flares are believed to bear a more notable cross-sectional expansion owing to their longer lengths. Compared with the main-phase loops, the late-phase loops would emit more effectively at middle temperatures, which could, to a certain degree, mitigate the severe heating requirement for the production of a prominent warm coronal late-phase peak. In addition, the cross-sectional expansion also affects the shape of the emission lights curves, causing a sharper decay after the emission peak. Such an emission pattern has been validated with the observations of an EUV late-phase flare, and could serve as a potential diagnostic tool to judge the degree of loop cross-sectional expansion in an extended flare dataset.

We perform a comprehensive analysis of the chemical and dynamical properties of quasar-damped Lyman-{\alpha} (DLA) galaxies and compare these to the GEAR chemodynamical simulations. Specifically, we aim to constrain the behavior of {\alpha}-element enhancements with metallicity, the dependence of [{\alpha}/Fe] on the specific star formation rate (sSFR), and the absorption-line velocity widths ({\Delta}v90) vs. stellar mass, {\Delta}v90 vs. metallicity, and mass-metallicity relations. For the comparison, we select five galaxies simulated with the chemodynamical Tree-SPH code GEAR with stellar masses in the range of log(M*/Msol) between 6.1 and 10.8, and at six different redshifts between 0.33 and 4.12. We find that the abundance ratios [{\alpha}/Fe] and [M/H] observed in the interstellar medium (ISM) of DLA galaxies overlap with the abundance trends in gas of the simulated galaxies. Our findings corroborate a picture in which DLAs with {\Delta}v90 below and above 100 km/s trace galaxies with masses in the ranges of log M* 6 - 8 and 8 - 11 solar masses, respectively. We suggest that observations should be used with caution when constraining the theoretical [{\alpha}/Fe] vs. sSFR relations because of systematics (if abundances are obtained from emission lines) or differences in the gas properties as probed by a DLA and its counterpart. So far, only the observations in absorption of inner gas of the LMC and SMC are in agreement with the simulated data. We confirm that DLAs detected at large impact parameters most likely probe the gas of satellite or other halo galaxies which are adjacent to the central galaxy. We further find that the velocity widths vs. stellar masses and mass-metallicity relations agree well with observations, while GEAR should be calibrated more carefully to reproduce the {\Delta}v90 vs. metallicity relation.

Triple-double radio galaxies (TDRGs) are amongst the rarest subpopulations of radio galaxies (RGs). They are characterised by three pairs of radio lobes, where each pair of lobes represents an episode of nuclear activity. Such a feature makes them key objects that can be used to constrain the duty cycle of RGs. In this paper, we report the discovery of J022248\m060934, a new TDRG, hosted by a galaxy at a spectroscopic redshift of $z \approx$ 0.94. We have used the MIGHTEE-DR1 data set and MIGHTEE sub-band images as our main data. In total intensity, J022248\m060934 has a bright core and triple-double, edge-brightened-like peaks of radio emission. The polarimetry of the source reveals an inhomogeneous density of the hosting environment which is consistent with the more pronounced bending in its eastern lobes. The spectral index and curvature maps suggest an inverted core and an ultra-steepening of the spectrum towards the outer lobes which reinforce a recurrent nuclear activity. We perform individual spectral age fitting of the components of the source using the JP model and we found a lower limit total age of $\sim$16 Myr. We also derive a short inactive period between the active phases and a rapid duty cycle of 90 per cent for the first cycle of activity. Our spectral ageing analysis suggests that the triple-double structure in TDRGs is not the product of long quiescent periods, as deduced by previous works based on kinematic ages.

Open clusters (OCs) usually evolve gradually as the number of their members changes, which can be manifested in their morphological characteristics. We aim to investigate the morphological stability of 1,490 OCs and further explore the potential change of morphological stability of the OCs at different spatial positions, using the OC catalog from the literature. We define for the first time a new morphological stability parameter Ncore/Nouter, a ratio of member numbers between cluster core and outer areas within tidal radii, which has a significant positive correlation against N, with a slope of 1.140\pm0.039, significantly steeper than the 0.720\pm0.026 measured for Score/Souter. This demonstrates that the stellar density in the core is a more sensitive tracer for morphological stability than geometry. Spatially, the radial sample OCs have larger slopes of Ncore/Nouter and Score/Souter against N, with 1.083\pm0.116 and 0.733\pm0.080, respectively, whereas those in the tangential direction 1.013\pm0.110 and 0.529\pm0.075, respectively, which means that the impact on sample OCs from tidal forces directed toward the Galactic center is possibly stronger than that from the shear force caused by the differential rotation of the Galactic disk. Besides, the sample OCs younger than 30 Myr display a shallow slope of 0.751\pm0.166, with those older than 800 Myr (1.442\pm0.128), reflecting that young OCs likely endure both internal disruptions, such as early dynamical heating weakening core binding and more severe external disturbances, compared to older OCs.

This paper presents new data and analyses of the AT Pyx system, a disk-hosting young star located in a cometary globule in the Gum Nebula. This radiation-driven structure is an unusual environment for observations of planet formation and differs greatly from the low-mass star-forming regions disks are most commonly observed in. Aided by a collection of visual and spectroscopic data available for this system, this paper aims to infer the possibility of embedded planets existing within the disk and how the system's environment may affect its disk morphology. Using data from the VLT's instruments XSHOOTER, ESPRESSO and - most prominently - SPHERE along with data from ALMA, we make a variety of measurements (geometric, photometric and otherwise) to characterise the observed disk features and attributes such as spiral arms and eccentricity. Mapping of the velocity components is also undertaken using the ALMA gas line data to characterise disk orientation and determine the likelihood that the system is experiencing a late-stage infall event. The disk is found to be eccentric when deprojected. Under the assumption that the formation of a planet is wholly responsible for the primary and secondary spiral arms, we find the mass of such a planet can range between 0.004 and 3 Jupiter masses. Measurement of the velocities associated with nearby globule cloud material returns reasonable velocities for a late-stage infall event. We estimate far-ultraviolet field strength at AT Pyx's location to be low in comparison to other surveyed disks. We also find that AT Pyx is possibly a binary system. AT Pyx is the first disk within a cometary globule to be spatially resolved, and is now the first such disk to be investigated to this extent. The work of this paper could potentially be a first step into the further study of disks in the moderate-FUV environment of the Gum Nebula and its globules.

The use of silicon photomultipliers (SiPMs) alongside conventional photomultiplier tubes (PMTs) is a remarkable technological development in modern ground-based very high energy gamma-ray astronomy. SiPMs exhibit comparable or even higher photon detection efficiencies (PDEs) than PMTs. The sensitivity of a PMT matches well the spectral shape of Cherenkov radiation from extended air showers. In contrast to a PMT, the sensitivity of a SiPM is shifted toward longer wavelengths, where the intensity of light of night sky (LoNS), considered as unwanted noise, increases significantly. It is obvious that a SiPM with a higher PDE will indeed measure more Cherenkov light than a PMT, but it will also detect significantly higher LoNS noise; the question is which factor will predominate in the signal-to-noise-ratio (SNR). To compare the performance of a PMT with that of a SiPM, we built SiPM-based modules and installed these and operated in parallel in the imaging camera of the 17 m diameter MAGIC telescope. Our long-term studies show that SiPM, despite their higher PDE, can deliver only a comparable to PMT performance. As already the name SiPM suggests, we use these semiconductor sensors analogously to classical PMTs: We amplify their small signals, digitize, and calibrate the converted amplitudes. Although SiPM is essentially a digital sensor, its common-anode design does not allow one to directly profit from it. Numerous arrays of single-photon avalanche diodes (SPADs) are being developed in various laboratories worldwide. Unlike SiPM, SPAD arrays digitize the incident photons from the outset and count their number. We will dwell on the potential further developments of SPADs.

Primordial black holes formed in the early universe are compelling candidates for dark matter. We investigate their production in Rastall gravity, a modification of general relativity that introduces a non-minimal coupling between matter and geometry through the non-conservation of the energy-momentum tensor. Analyz- ing cosmological perturbations during radiation domination, we demonstrate that the Rastall parameter fundamentally alters the collapse dynamics, modifying the growth of density fluctuations, the critical threshold for black hole formation, and the fluctuation amplitude at horizon crossing. Current cosmological constraints from Big Bang Nucleosynthesis, the cosmic microwave background, and large-scale structure restrict the Rastall parameter to small values, yet within this allowed range PBH production can be altered by orders of magnitude compared to general relativity. Our results establish primordial black holes as novel probes of modified gravity.

We investigate how the magnetic loop length of solar flares relates to the timing between their thermal and nonthermal emission signatures. Our study analyzes a sample of 96 C-, M-, and X-class flares observed between 2013 and 2015 with soft X-rays, hard X-rays, and extreme UV. For each event, we determine the time delay {\Delta}t between the hard X-ray and soft X-ray peak, and estimate the flare loop length L from UV footpoints assuming a semicircular geometry. In every case, longer flare loops are consistently associated with larger timing delays. Across the full sample, we find a strong correlation, R = 0.88 between L and {\Delta}t. We also quantify how closely each flare follows the Neupert effect using a coefficient RN, defined as the Pearson correlation between the time derivative of the soft X-ray flux and the hard X-ray light curve. Applying correlation thresholds of RN > 0.5 and RN > 0.8 yields subsets of 87 and 46 events, respectively. In both cases, the linear relationship between loop length and peak delay remains clearly expressed. For the RN > 0.5 subset, the correlation is R = 0.87, while the more selective subset with RN > 0.8 displays an even stronger correlation of R = 0.91. These results show that the overall trend persists across increasingly stringent correlation thresholds. The results provide direct observational confirmation that magnetic loop geometry plays a key role in governing the temporal evolution of energy transport in solar flares.

We provide context for Apophis' 2029 Earth passage by analyzing its possible source populations, in particular, the Flora family, which has a similar composition, corresponding to LL chondrite meteorites. Out of ${\sim}3380$ NEOs larger or equal than Apophis (${\ge}420\,{\rm m}$), $610\pm 140$ are LL-like NEOs from Flora. Their mean encounter probability is $p = 86\times 10^{-18}\,{\rm km}^{-2}\,{\rm y}^{-1}$, corresponding to once per 13000 y frequency of encounters closer than 38000 km. However, this does not apply to Apophis alone, for which the specific encounter probability is higher, $p' = 1603\times 10^{-18}\,{\rm km}^{-2}\,{\rm y}^{-1}$, but the frequency is lower, only once per 430000 y, when we consider it as a single object. Our simulation of the Flora family over $\sim$1 billion years indicates that Apophis-like bodies have orbits that are particularly persistent in near-Earth space. The temporal distribution of encounter probabilities exhibits peaks (up to ${>}10^4$ in the same units) and the specific value for Apophis is not unusual (occurring ${\sim}70\%$ of time). In other words, there is always at least one Apophis-like body among NEOs. We find that such persistence also creates favorable opportunities for temporary capture as Earth coorbitals. Apophis-like bodies are ultimately removed from the inner solar system by approaching the Sun or by impact into one of the terrestrial planets, where the relative split between these outcomes is $(45\pm 2)\,\%$ and $(50\pm 2)\,\%$. While our current knowledge of Apophis' orbit guarantees no threat from Apophis in the next few centuries, we cannot predict any specific outcome for Apophis in the coming thousands or millions of years. Evaluating this statistically over the long term, we find that objects in Apophis-like orbits have a $(19\pm 2)\,\%$ chance of Earth impact over their lifetime of ${\sim}30\,{\rm My}$.

The Milky Way (MW) stellar halo, containing debris from past accretion events, serves as a fossil record of hierarchical mass assembly. Due to their distinct properties, RR Lyrae stars (RRLs) serve as excellent tracers for identifying and characterising the halo's substructures. We analysed a sample of 4933 RRLs, for which we calculated the integrals of motion and orbital parameters. We applied the domain-informed novelty detection CLustering in Multiphase Boundaries (CLiMB) framework to identify RRL membership in the MW substructures. We analysed the metallicity distributions of RRLs in major accreted system remnants as a snapshot of their chemical evolutionary status during early epochs. We calculated the weighted mean metallicity ([Fe/H]) and the corresponding standard deviation for Gaia Sausage/Enceladus ([Fe/H] = $-1.57 \pm 0.25$ dex), Sequoia ([Fe/H] =$ -1.64\pm0.26$ dex), and the Helmi streams ([Fe/H] = $-1.66\pm0.19$ dex). The metallicity distribution of RRLs in Thamnos was found to be bimodal, with the metal-poor peak likely representing the genuine accreted Thamnos population ([Fe/H] = $-1.94\pm0.20$ dex), in agreement with recent works based on spectroscopic abundances. Our analysis shows that the substructures ED-1 and L-RL3 are highly contaminated by thick disc stars. However, the metal-poor tails in their metallicity distributions may be signatures of remnants from small accreted systems. We also identify over-densities of RRLs in correspondence with the recently reported substructures Shiva and Shakti, which we suggest are of in-situ origin. Finally, we applied the RRL-based mass-metallicity relation of galaxies to test the nature of the identified dynamical substructures.

The radio detection of extensive air showers provides a powerful method for studying the origin of high-energy cosmic rays. The Low-Frequency Array (LOFAR) offers unprecedentedly detailed measurements of the radio emission footprint. However, fully exploiting this information requires advanced reconstruction techniques. In this paper, we introduce a novel framework for air shower reconstruction based on Bayesian inference and Information Field Theory (IFT). Our method is built on a fully differentiable forward model of the radio signal, which incorporates a physical emission parameterization and a precise wavefront model. Additionally, we augment this physical model with Gaussian processes to account for systematic uncertainties in both the signal fluence and arrival timing. By leveraging gradient information, our approach enables efficient (three orders of magnitude acceleration w.r.t.\ the legacy method) and robust inference of the underlying physical shower parameters, such as primary energy and the depth of shower maximum, $X_\text{max}$. This work provides not only point estimates but also a rigorous quantification of uncertainties. We achieve a resolution in $X_\text{max}$ of $25\,\mathrm{g/cm^2}$ and a radiation energy resolution of $12\%$ on simulations for LOFAR.

NGC 4051 is a nearby (16.7 Mpc), Narrow Line Seyfert 1 galaxy (NLS1), which has a low black hole mass of $10^6$ M$_{\odot}$. It is also known for its rapid X-ray variability, on timescales of kilo-seconds and has a complex, multi component wind in both the soft X-ray and Fe K bands. Here we present the first high resolution XRISM Resolve spectrum of NGC 4051, which was captured in a historically bright state for a 150 ks exposure. XRISM resolves two blue-shifted Fe K shell absorption troughs in the mean spectrum, which can be ascribed to H-like iron and arises from two outflow components with outflow velocities of 0.025c and 0.04c. A time dependent spectral analysis shows that the iron K absorption is variable on timescales of less than a day, increasing in velocity over the duration of the observation. The velocity changes may be explained either by the passage of two separate transiting absorbers, of different velocities, or by a single accelerating outflow of approximately constant column density. In the latter case, the wind acceleration is likely to be too large to be caused by radiation pressure and instead magnetic driving is favored to accelerate the wind up to 0.04c. The outflow can originate from an accretion disk wind, whose kinetic power is sub-Eddington in contrast to recent examples of winds from powerful, luminous quasars observed by XRISM.

There is still much debate around the inner workings of the GRB prompt emission phase with many questions still left unanswered. Polarization signatures offer a promising new avenue to discriminate between the various GRB prompt emission models. The aim of this study is to estimate energy and time resolved polarization signatures resulting from Inverse Compton (IC) scattering for two specific GRB prompt emission models, namely the backscatter-dominated cork model by \citet{Vyas2} and a Compton drag model by \citet{Lazzati2}. In order to achieve this we apply an IC polarization Monte Carlo (MC) algorithm to those two GRB models in order to estimate the expected polarization signatures. For the backscatter-dominated cork model we find polarization signatures below $\sim 10$~\%, likely below the detection limits of current or near-future X-ray and $\gamma$-ray polarimeters. Our results for the Compton drag model indicates polarization results consistent with that found by \citet{Lazzati3}. Furthermore, we find some energy and time dependence of the estimated polarization, particularly for the polarization angle.

We present the evolution of optical spectra and lightcurves of the fast nova V1674 Herculis during 150 days past its eruption. Using the post-eruption AAVSO light curve, we have calculated the orbital period of V1674 Her to be 0.153 days. There is no unambiguous white dwarf spin period in our data. The optical spectra show that the ionisation increases with time. A morpho-kinematic analysis of the H$\alpha$ line profile indicates a bipolar morphology with polar blobs and an equatorial ring. Lyman beta fluorescence is found to be the dominant mechanism for the excitation of neutral oxygen. On day 19.87, [Ne III] & [Ne V] lines are present, indicating the presence of the ONe white dwarf. On day 147.66, the nebular lines are still present, implying that the nova had not gone into quiescence yet; this spectrum is accretion-dominated.

We present deep, high-resolution ($\sim$100 mas) Karl G. Jansky Very Large Array (VLA) Ka-band (9.1 mm) observations of the disk around MWC 480, and infer dust properties through a combined analysis with archival Atacama Large Millimeter/submillimeter Array (ALMA) data at 0.87, 1.17, 1.33, and 3.0 mm. The prominent dust ring at 95 au (B95) is detected at 9.1 mm for the first time, while the faint outer ring at 160 au is not revealed. Through non-parametric visibility modeling, we identified two new annular features: a plateau within 20-50 au across all wavelengths, and a shoulder exterior to the B95 ring at 0.87, 1.17 and 1.33 mm, consistent with signatures of planet-disk interaction. We find that the width of the B95 ring remains constant across wavelengths, suggesting that fragmentation dominates over radial diffusion or that unresolved substructure is present within the ring. Resolved spectral modeling yields two families of dust solutions that reproduce the observations equally well: compact grains or highly porous (90\%) grains, with carbonaceous components dominated by refractory organics or amorphous carbon, respectively. The inferred maximum grain sizes peak at the locations of the two rings and reach centimeter within the B95 ring. The total dust masses are $860^{+95}_{-78}\rm~M_\oplus$/$1500^{+440}_{-330}\rm~M_\oplus$ (large/small-grain solution in inner disk) and $230^{+14}_{-13}\rm~M_\oplus$ for the two dust mixtures. The B95 ring alone contains $100^{+5}_{-5}\rm~M_\oplus$ and $43^{+2}_{-2}\rm~M_\oplus$, respectively, sufficient to assemble the cores of giant planets. Finally, we highlight the power of broadband, multi-wavelength observations in placing better constraints on dust composition and porosity in protoplanetary disks.

How the solar atmosphere is heated from a temperature of about $5,000-6,000$\,K in the lower atmosphere to about $1-2$\,MK in the corona has challenged the astrophysical community for about 80 years. The same puzzle exists for the stellar coronae heating as well. In this study, we present a series of findings on solar spicules and their subsequent impact on the corona within a coronal hole environment, characterized by locally open magnetic field lines, combining insights from MHD simulations with observations. We find that the convective and turbulent motions around the solar surface cause plenty of shocks and small-scale magnetic reconnection in the lower atmosphere. The combined effects of shock compression and reconnection outflows then drive the formation of groups of spicules with a quasi-period of about $300$\,s and width of $\sim 200-500$\,km. The spicule upflows provide an averaged mass flux above $10^{-9}$\,kg\,m$^{-2}$\,s$^{-1}$ in the lower corona to sustain the solar wind in coronal holes, and they continuously trigger further new local slow-mode waves and shocks. These waves supply an energy flux of $10-100$\,W\,m$^{-2}$ in the lower corona, and they are dissipated by heat conduction and compression heating to sustain the corona temperature of about $1$\,MK. The results also indicate that the upward propagating disturbances (PDs) observed in extreme ultraviolet (EUV) passbands are caused by both spicule upflows and slow-mode waves and shocks. Our findings help to understand the long standing problem of coronal heating and the origin of solar winds in coronal hole regions.

This study presents a Normal Behavior Model (NBM) developed to forecast monitoring time-series data from the ASTRI-Horn Cherenkov telescope under normal operating conditions. The analysis focused on 15 physical variables acquired by the Telescope Control Unit between September 2022 and July 2024, representing sensor measurements from the Azimuth and Elevation motors. After data cleaning, resampling, feature selection, and correlation analysis, the dataset was segmented into fixed-length intervals, in which the first I samples represented the input sequence provided to the model, while the forecast length, T, indicated the number of future time steps to be predicted. A sliding-window technique was then applied to increase the number of intervals. A Multi-Layer Perceptron (MLP) was trained to perform multivariate forecasting across all features simultaneously. Model performance was evaluated using the Mean Squared Error (MSE) and the Normalized Median Absolute Deviation (NMAD), and it was also benchmarked against a Long Short-Term Memory (LSTM) network. The MLP model demonstrated consistent results across different features and I-T configurations, and matched the performance of the LSTM while converging faster. It achieved an MSE of 0.019+/-0.003 and an NMAD of 0.032+/-0.009 on the test set under its best configuration (4 hidden layers, 720 units per layer, and I-T lengths of 300 samples each, corresponding to 5 hours at 1-minute resolution). Extending the forecast horizon up to 6.5 hours-the maximum allowed by this configuration-did not degrade performance, confirming the model's effectiveness in providing reliable hour-scale predictions. The proposed NBM provides a powerful tool for enabling early anomaly detection in online ASTRI-Horn monitoring time series, offering a basis for the future development of a prognostics and health management system that supports predictive maintenance.

We present the low-resolution optical transmission spectrum of the inflated hot Saturn HAT-P-44b. The planet is a close sibling in radius (1.24 $\mathrm{R_{Jup}}$), temperature (1100 K), and mass (0.35 $\mathrm{M_{Jup}}$) to the exceedingly well-characterized WASP-39b. Using the ACAM instrument on the William Herschel Telescope (WHT), we obtain a transmission spectrum with sub-scale height precision of 246 ppm, with a wavelength range of 495 -- 874 nm and a 20 nm resolution, despite a relatively faint host star ($V\mathrm{_{mag} = 13.2}$). We detect absorption due to sodium with 3.9$\sigma$ confidence. Atmospheric retrieval of the transmission spectrum also reveals evidence for \ch{H2O} absorption and Rayleigh scattering from \ch{H2} gas consistent with a cool 800 K atmosphere and a super-solar metallicity of 7$\substack{+16 \\ -5}$$\times$solar. Comparison of retrieval models disfavour the inclusion of a super-Rayleigh scattering slope or high-altitude clouds (at $<1$ mbar) while being agnostic towards the presence of mid-altitude clouds. Our transmission spectrum of HAT-P-44b shows strong similarity to that of its sibling WASP-39b.} This is the tenth planet in the LRG-BEASTS (Low-Resolution Ground-Based Exoplanet Atmosphere Survey using Transmission Spectroscopy) survey.

Bars are among the most prominent galactic structures, yet their formation mechanisms remain incompletely understood. They can form either internally, via dynamical instabilities, or externally, triggered by interactions with other galaxies. The impact of mergers on bar formation and survival, however, has not been thoroughly investigated. To explore the influence of mergers on bars, we construct a suite of \textit{N}-body merger pairs where a Gaia-Sausage-Enceladus-like radially biased satellite disk galaxy merges with a central disk galaxy during its bar formation. With the central galaxy fixed, the satellite varies in merger parameters: the mass ratio $m/M$ relative to the central galaxy, the impact parameter $b$, and the orbital inclination angle $\theta_i$ relative to the central disk. We find that the bar survival probability decreases with increasing $m/M$. Mergers with $m/M\lesssim1/10$ generally preserve the forming bar, whereas those with ${m/M}\geq1/2$ tend to destroy it, producing more early-type-like remnants. For intermediate mass ratios ($1/5 \leq m/M \leq 1/3$), several models yield ``weakening bars'', in which the bar survives the merger but gradually decays during subsequent secular evolution, possibly due to interactions between nested double bars formed from merger debris. In contrast to $m/M$, $b$ and $\theta_i$ have only secondary and stochastic effects on bar survival. The different influences of these three merger parameters can be naturally explained by the tidal force exerted by the satellite on the forming bar, which tends to weaken the bar when the satellite crosses it nearly perpendicular to its major axis.

We analyse CMB constraints on stable networks of cosmic strings and domain walls using for the first time full Planck 2018 data together with BICEP/Keck 2018 $B$-mode measurements. The defect-induced anisotropies are computed using the Unconnected Segment Model for Nambu-Goto and Abelian-Higgs strings, as well as for stable domain walls, and included in a full Markov Chain Monte Carlo analysis jointly varying all $\Lambda$CDM parameters, the tensor-to-scalar ratio, and the string/domain wall tension. No statistically significant evidence for defects is found, although we observe a mild preference for non-zero cosmic string tension. Our results improve previous constraints on the defect power spectrum by up to a factor of two. In the particular case of strings, the improvement is driven by the $B$-mode data, and is especially pronounced for Abelian-Higgs strings. We also present forecasts for upcoming Simons Observatory data and find that, with the baseline noise configuration, the constraints on the string tension could improve by about a factor of three. Finally, we assess the impact of Nambu-Goto string loops on CMB anisotropies in light of both current and future observations.

The fundamental parameters of a low-mass star can potentially be determined from its photometry and astrometry. This is complicated by the fact that 10-20 percent of low-mass stars are predicted to be equal-mass binaries. These unresolved systems appear more luminous compared to single stars with the same fundamental parameters. We present a method to differentiate binary stars from single-star main sequence K and M dwarfs using their \textit{Gaia} DR3 XP spectra. We assemble a training set of stars which have pristine astrometry and photometry, are located within 100pc of the Sun, and exclude stars with \textit{Gaia} DR3 flags suggesting they may be unequal mass systems, thereby leaving stars that are predominantly either single- or equal-mass binaries. We then iteratively train Random Forest Regression (RFR) models to predict absolute magnitude and color given the RP spectral coefficients of a star. After each model, we remove the stars that have absolute magnitudes significantly brighter than their predicted values. This method converges on a model trained only on single stars. We then use this model to identify the ``overluminous'' K and M stars in \textit{Gaia} DR3 within 100 parsecs, with some quality cuts. We find that $\sim13\%$ of the sample is significantly overluminous and assume these to be unresolved binaries. We aggregate several multiplicity surveys across different projected separations and incorporate our overluminous binaries to create a general \textit{Catalog of Systems} within 100 pc. We use this \textit{Catalog} to provide lower limits on the multiplicity fraction for stars between $0.1$ and $0.7~M_{\odot}$.

HelioSpectrotron~5000 (HS5000) is an interactive, multi-resolution solar spectral atlas designed to facilitate direct comparison between high-resolution reference spectra and observations obtained with a wide range of ground-based instruments. Based on the Hamburg FTS atlas, the HS5000 provides both absolute intensity and continuum-normalized spectra at arbitrary spectral ranges and resolutions, as well as curated line identifications and optional telluric contamination. This framework enables rapid wavelength calibration, line identification, and context image generation.

We argue that all dynamical astrophysical black hole merger channels are expected to result in a common eccentricity distribution at gravitational wave (GW) frequencies relevant for LIGO/Virgo/KAGRA (LVK) in the high eccentricity limit. This follows from the large separation of scales between the GW regime required for creating eccentric mergers in LVK, and the underlying astrophysical formation environment. Our analytical solution shows exceptional agreement with numerical studies. This finding has important implications for both theoretical studies and ongoing searches for eccentric GW sources.

High-energy multimessenger background analyses over the past decade have provided evidence for a population of hidden neutrino sources that are opaque to GeV-TeV gamma rays, a picture bolstered by recent observations of the nearby active galaxy NGC 1068. The coronal regions in the hearts of active galactic nuclei (AGNs) have been proposed as the most promising sites for such hidden nonthermal particle production, and NGC 1068 is expected to be the most neutrino-active galaxy for IceCube. We demonstrate that the latest all-sky neutrino spectrum, exhibiting a spectral bend around 3-30 TeV, is consistent with predictions of the magnetically powered corona scenario, and the models for the all-sky neutrino flux can simultaneously explain the multimessenger data from NGC 1068 within observational and modeling uncertainties. We further show, in a largely model-independent way, that the contribution from NGC 1068-like sources does not overshoot the observed medium-energy neutrino flux. Finally, we highlight the key role of the Eddington ratio, which can drive substantial variations in the predicted neutrino fluxes of nearby AGNs, and we encourage systematic multimessenger searches for the neutrino-brightest AGNs.

Among potentially habitable worlds, rocky planets orbiting M dwarfs offer the most favorable prospects for atmospheric characterization, yet their climates may differ substantially from those of Earth analogs. In the tidally locked limit, the nightside's tendency to radiatively cool and potentially trap volatiles as permanent ice introduces a strong dependence of habitability on the planet's surface and atmospheric boundary conditions. We perform a suite of synchronously rotating experiments spanning a wide range of topographic and orographic realizations with different mean elevations and landmass distributions. Across a grid of $p_{\mathrm{N2}} = 0.5$-$8~\mathrm{bar}$ and $F_{\star} = 1200$-$1700~\mathrm{W\,m^{-2}}$, we find that surface relief breaks the flow symmetry, replacing the circumpolar vortices with mechanically forced stationary waves. Steep orography produces standing Rossby gyres that strengthen the cross-terminator jet and align vertical uplift with the day--night boundary. These new circulation regimes enhance moisture transport, increasing the infrared optical depth and promoting additional nightside cloud formation, which produces a stronger cloud-greenhouse feedback and lower the critical fluxes required for global planetary deglaciation. Broad, elevated plateaus drive a similarly fragmented but slightly weaker circulation, yielding less effective moisture transport. These results show that the relief and spatial distribution of landmasses, parameters unconstrained for most exoplanets, can exert strong controls on the climatic bifurcations of tidally locked M-dwarf exoplanets.

We report an updated constraint on the Solar System quadrupole parameter $Q_2$, which encodes the external field effect predicted by modified gravity versions of the MOND paradigm. Using the dataset employed to compute the DE440 planetary ephemerides, and estimating it simultaneously with other parameters included in the planetary ephemerides, we find $Q_2 = (1.6 \pm 1.8) \times 10^{-27}\,\mathrm{s}^{-2}$ (1-$\sigma$), representing an improvement of 40% over previous estimates. We also show explicitly that the contribution to the MOND prediction of $Q_2$ from the Solar System's largest planet, Jupiter, is at the 0.05% level, validating the approximation of retaining only the Sun in theoretical calculations. With this new constraint on $Q_2$, we update previously acknowledged tensions with external galaxy rotation curves, now leading to discrepancies at the $3$-$15\sigma$ level depending on the detailed mass modeling or the subset of galaxies considered. Within the Milky Way itself, the $Q_2$ constraint imposes an upper bound of only 2% (at 95% confidence) on the MOND boost to the galactic radial acceleration (i.e., the ratio of the observed over baryonic Newtonian acceleration) at the position of the Sun, in strong tension with current observational limits. The updated $Q_2$ posterior finally confirms that Solar System measurements provide stronger constraints than current wide-binary data on classical modified gravity versions of MOND.

A new contextualised reading of Fritz Zwicky's 1933 article ''The redshift of extragalactic nebulae'' about the virial analysis of the velocity dispersion of galaxies in the Coma cluster leads to a reconsideration of the traditional discourse on the introduction of dark matter. We argue that this component of matter was not only already on the stage of the scientific debates of the time, but also, in a more concealed form, played a central role in Zwicky's epistemic context. We thus reject the narration that dark matter is the result of a ``na{ï}ve'' astrophysical observation and emphasise the cosmological motivations that prompted Zwicky to presciently search for it. Moreover, with regard to its abundance, we argue that the discrepancy between the observed amount of luminous matter in the Coma Cluster and Zwicky's higher mass estimate derived from virial analysis was not, in fact, astonishing. What Zwicky described as a surprising excess of dark matter was of precisely the order of magnitude he had set out to identify. Consequently, we challenge the widespread view that dark matter was merely an ad hoc hypothesis introduced to rescue Newtonian theory. Instead, we suggest it may represent one of the earliest cosmological indications supporting a new emerging theory of gravitation: General Relativity. This reinterpretation contributes to ongoing debates in the philosophy of science concerning the epistemic status of ad hoc hypotheses.

Long-lived particles (LLPs) provide an interesting window into physics beyond the Standard Model, offering characteristic signatures at colliders and in cosmology. In this work, we investigate LLPs decays into dark matter. If the lifetime of LLPs are longer than $10^4$s, the decay products can disrupt the synthesis of light nuclei in the early universe and alter Big Bang Nucleosynthesis (BBN) predictions. If the LLP is much heavier than the dark matter particle, the decay contributes to the number of effective neutrino species, $N_{eff}$. We describe these decays via dimension-six effective operators and outline the parameter space in which such decays obey cosmological bounds stemming from BBN, structure formation, Cosmic Microwave Background, and Baryon Acoustic Oscillation data.

We use the Velocity-dependent One Scale Model for topological defect evolution to explore and classify the possible scaling solutions for string networks with time-varying tension, in cosmological and non-cosmological settings and under two different phenomenological assumptions for the behavior of these variations, which rely on different stretching and damping contributions to the string dynamics. We discuss how these assumptions impact the standard scaling solutions, as well as the evolution of the string network density. In addition to simple power-law cosmological epochs solutions, we also discuss the behavior of the network during the radiation-to-matter and matter-to-acceleration transitions. Overall, our results show that for the same amount of tension variation, a change in the stretching length scale tends to have a more significant impact on the network than a change in the damping length.

We consider glueball dark matter (DM) in a Yang-Mills dark sector confined at $\Lambda_D$ scale and coupled to the Standard Model through electrically and dark-color charged vector-like fermion portals, with the mass scale $m_\psi$. In a simple case with two lightest mass-degenerate vector-like fermions with opposite electric charges the effective amplitudes with one $C$-odd glueball (oddball) and odd number of photons vanish, rendering the lightest $C$-odd spin-1 state with mass $m_\chi$ a viable DM candidate provided that $m_\psi\gtrsim 5.5 \Lambda_D$. We develop a controlled effective field theory framework with non-perturbative information supported by QCD phenomenology leading to a quantitative prediction for coherent elastic glueball scattering off nuclei. We find a steep scaling of the spin-independent cross section $\sigma_{\rm SI}\propto \Lambda_D^{2.15} m_\psi^{-8}$. This implies that the sensitivity of the current and next-generation xenon experiments in the range of $\sigma_{\rm SI} \sim 10^{-46} - 10^{-48}$ cm$^2$ corresponds to $m_\psi \simeq 3-30$ GeV, respectively, for $\Lambda_D\simeq 0.55-5.5$ GeV. We provide a minimal UV completion of the portal sector compatible with collider phenomenology. Our results pave a quantitative foundation for testing glueball DM in direct-detection experiments.

In this work, we perform a relativistic test of the spacetime geometry of Sagittarius A* (Sgr A*) using the orbit of the S2 star. We consider a broad class of compact object models, including Schwarzschild, Reissner-Nordström, Bardeen, Hayward, and Simpson-Visser black holes, as well as the Janis-Newman-Winicour naked singularity spacetime. For each geometry, we integrate the timelike geodesic equations and consistently project the resulting trajectories onto astrometric and spectroscopic observables, incorporating Rømer time delay and relativistic redshift effects. The theoretical predictions are tested with current Very Large Telescope (VLT) observations of the S2 star, while simultaneously imposing constraints from the Event Horizon Telescope shadow size. We find that several spacetimes that are degenerate at the level of shadow imaging, most notably Schwarzschild, Reissner-Nordström, and Bardeen regular black hole geometries, remain statistically indistinguishable when tested against present S2 data. We further carry out a statistical model comparison based on the Akaike and Bayesian information criteria (AIC and BIC) to evaluate the relative performance of the alternative spacetime models. Our analysis also constrains the generalized charge like parameter $q/M$ in non-Schwarzschild spacetimes based on current S2 star observations, and identifies specific black hole and horizonless geometries that can be further tested with forthcoming high precision astrometric observations from the VLT and Keck telescopes.

Observations of supermassive black holes by the Event Horizon Telescope reveal significant inhomogeneities, most likely related to density and magnetic field perturbations. To model these features, we conduct high-resolution 2D general-relativistic magnetohydrodynamics (GRMHD) simulations of a Fishbone-Moncrief torus around a Kerr black hole using the Black Hole Accretion Code $\texttt{BHAC}$. We compare unperturbed accretion with a case featuring plasma density bubbles with pressure balanced magnetic islands of different amplitudes. Power spectrum analysis of accretion time series, performed via the Blackman-Tukey method, shows that the perturbed case exhibits (1) steeper spectral indices compared to the unperturbed case, deviating from the characteristic $1/\omega$ noise spectrum, and (2) increased correlation times, providing evidence for absorption of macro-structures at the event horizon. Spatial auto-correlation analysis of near-horizon turbulence confirms larger energy-containing coherent structures in the perturbed case altering the accretion rate. These results provide new insights for interpreting observations of supermassive black hole environments, where near-horizon turbulence may play a key role in the accretion process.

We investigate the impact of nuclear pasta on crustal structure and torsional oscillations using a Bayesian ensemble of unified neutron-star equations of state based on relativistic mean-field models constrained by nuclear experiments, empirical saturation properties, chiral effective field theory, and multimessenger observations. For each posterior sample, we compute the pasta sequence within a compressible liquid-drop model and quantify the onset density, thickness, and mass fraction of the pasta layers. We show that the appearance and extent of nuclear pasta are primarily controlled by the symmetry-energy slope parameter $L$. While spherical and rod-like pasta configurations are present for all equations of state, only a small fraction of the posterior supports slab, tube, or bubble geometries. The transition from spherical nuclei to rods is tightly constrained to occur at a density of $\rho_{\rm sr} = 0.0588^{+0.0045}_{-0.0065}\,\mathrm{fm^{-3}}$. We further predict that the nuclear pasta layer occupies a relative radial thickness of $\Delta R_{\rm pasta}/\Delta R_{\rm c} = 0.140^{+0.025}_{-0.036}$ and contributes a relative mass fraction of $\Delta M_{\rm pasta}/\Delta M_{\rm c} = 0.475^{+0.071}_{-0.113}$. Using the resulting crust models, we present the first quasi-periodic oscillations (QPOs) analysis based on a Bayesian posterior ensemble of neutron-star equations of state and systematically assess their compatibility with observed low-frequency quasi-periodic oscillations. We find that the predicted QPO frequencies are strongly correlated with the curvature of the symmetry energy evaluated at sub-saturation density, $K_{\rm sym}(\rho_0/2)$, and that uncertainties in the equation of state translate into a range of angular indices $\ell$ consistent with the observed frequencies.

A fundamental challenge in low-frequency gravitational-wave detectors is the limited signal-to-noise ratio (SNR) of calibration lines, particularly in torsion-bar systems where the response is governed by rotational dynamics. In this work, we resolve this issue by optimizing the geometrical configuration of a torque-coupled gravity field calibrator (GCal), achieving an improvement in calibration-line SNR by more than an order of magnitude compared to conventional layouts. For the Cryogenic sub-Hz cROss torsion-bar detector with quantum NOn-demolition Speed-meter (CHRONOS), the calibration signal appears as a monochromatic line within the $0.1$--$10~\mathrm{Hz}$ band. At $1~\mathrm{Hz}$, the strain-equivalent calibration amplitude reaches $|h_{\rm GCal}| = 1.18 \times 10^{-14}$, corresponding to an SNR density of $|h_{\rm GCal}|/S_h = 4.25 \times 10^{3}$. This demonstrates for the first time that a high-SNR calibration line can be directly injected into the sub-Hz band of a torsion-bar detector. A first-order perturbative error propagation analysis yields a total fractional systematic uncertainty of $\delta h_{\rm GCal}/h_{\rm GCal} = 0.24\%$, dominated by geometric alignment uncertainties, while contributions from mass uncertainties and the gravitational constant remain subdominant. The corresponding absolute systematic uncertainty is $\delta h_{\rm GCal} \sim 10^{-17}$ at $1~\mathrm{Hz}$. These results establish torque-coupled gravitational calibration as a practical solution to the longstanding low-SNR problem in sub-Hz torsion-bar detectors and provide a robust pathway toward precision absolute calibration in the low-frequency regime.

Extreme mass-ratio inspirals (EMRIs) are among the most powerful probes of strong-field gravity and of the environments surrounding supermassive compact objects. Motivated by the expected presence of dark matter near galactic centers, we investigate the emergence and gravitational-wave imprints of chaotic dynamics in EMRIs evolving in non-vacuum spacetimes. Within a unified dynamical framework, we analyze test-particle motion in a broad class of dark-matter-embedded geometries, including singular black holes, regular black holes, naked singularities, and Einstein-cluster configurations. We show that environmental perturbations generically break integrability in the strong-field regime, giving rise to chaotic motion whose onset, duration, and termination depend sensitively on horizon structure, core regularization, and matter distribution. Using the numerical Kludge approach, we demonstrate that chaotic trajectories produce systematic qualitative modifications of the emitted gravitational radiation, such as irregular amplitude modulation and loss of phase coherence, in contrast to the smooth, quasi-periodic waveforms generated by regular motion. Our results establish the robustness of chaos in environmentally perturbed EMRIs and provide a clear conceptual link between nonlinear orbital dynamics, spacetime structure, and observable gravitational-wave signatures.

We have developed a compressed neighbor list for short-range particle-particle interaction based on a space- filling curve (SFC) memory layout and particle clusters. The neighbor list can be constructed efficiently on GPUs, supporting NVIDIA and AMD hardware, and has a memory footprint of only 4 bytes per particle to store approximately 200 neighbors. Compared to the highly-optimized domain-specific neighbor list implementation of GROMACS, a molecular dynamics code, it has a comparable cluster overhead and delivers similar performance in a neighborhood pass. Thanks to the SFC-based data layout and the support for varying interaction radii per particle, our neighbor list performs well for systems with high density contrasts, such as those encountered in many astrophysical and cosmological applications. Due to the close relation between SFCs and octrees, our neighbor list seamlessly integrates with octree-based domain decomposition and multipole-based methods for long-range gravitational or electrostatic interactions. To demonstrate the coupling between long- and short-range forces, we simulate an Evrard collapse, a standard test case for the coupling between hydrodynamical and gravitational forces, on up to 1024 GPUs, and compare our results to the analytical solution.

Leptophilic dark matter (DM) is a well-motivated thermal WIMP framework that can evade stringent nuclear-recoil searches while remaining testable via DM-induced heating of neutron stars (NS). In this work, we study leptophilic scalar DM in a $\mathrm{U(1)}_{L_\mu-L_\tau}$ gauge extension of the Standard Model, which provides a common leptophilic portal for all scenarios considered. To reproduce the observed relic abundance while suppressing direct-detection signals, we investigate three benchmark realizations: (i) a secluded DM scenario in which the relic density is set by annihilation into $\mathrm{U(1)}_{L_\mu-L_\tau}$ gauge bosons, and two pseudo-Nambu-Goldstone boson (pNGB) DM models based on (ii) an SO(4) symmetry and (iii) an SO(3) symmetry. In the SO(4) pNGB model, the DM mass arises at tree level from a soft breaking term, while the elastic scattering amplitude is suppressed by a symmetry-protected cancellation. In the SO(3) pNGB model, the DM mass is generated radiatively at one loop via the $\mathrm{U(1)}_{L_\mu-L_\tau}$ gauge interaction, and we show that this gauging preserves the same cancellation mechanism, maintaining compatibility with direct-detection null results. We perform a systematic parameter scan imposing relic density, direct- and indirect-detection, and neutrino trident constraints, and identify viable sub-TeV to TeV DM candidates. Assuming maximal capture in NSs, we find that the remaining parameter space can be tested by near-infrared observations of NSs, providing sensitivity complementary to terrestrial searches in regions that are currently weakly constrained.

Mapping gravitational-wave background (GWB) anisotropy with pulsar timing arrays (PTAs) is affected by harmonic-space mode suppression and mode coupling arising from an array's nonuniform sky response. Spherical harmonic expansions must be truncated at finite multipole l_max^rec, often set to l_max^N_pair$\equiv {\rm int}\left[\sqrt{\text{N_pair}}-1\right]$, where N_pair is the number of distinct pulsar pairs in an array. This choice is motivated by the counting argument that cross-correlations provide at most N_pair independent constraints. We obtain the multipole l_max^res corresponding to the maximum informative angular scale of a PTA. It is defined such that expansions to l_max^res (approximately) span the space of "observable skies" encoded in the N_pair eigenmaps of the Fisher information matrix, and therefore depends on the array configuration. We explicitly show that GWB power contained in multipoles l$\gtrsim$l_max^res do not significantly affect analyses that use expansions out to l_max^res, because the PTA response acts as a low-pass filter. In contrast, truncating at l_max^rec< l_max^res leads to leakage of small-scale angular power from l_max^rec<l$\leq$l_max^res. Even choosing l_max^rec=l_max^res, the standard frequentist estimator of the angular power spectrum C_l remains biased by the modes unobservable by the array. Although we can (partially) debias the standard estimator -- improving its agreement with an injected spectrum -- this reduction in bias comes at the expense of an increase in variance, particularly for poorly constrained modes with l$\gg$l_eff. We therefore recommend: (i) using l_max^res for PTA analyses involving spherical harmonic expansions, and (ii) using the debiased standard estimator for C_l recovery, but only out to multipoles l<l_eff ($\ll$l_max^res) corresponding to sufficiently constrained modes.

The aim of this work is to develop high precision Time-of-Flight (TOF) devices based on high refractive index solid Cherenkov radiators read out by silicon photomultipliers (SiPMs). Cherenkov light is prompt and therefore ideal for reaching the intrinsic timing limits of TOF systems. By utilizing a thin, high-refractive-index radiator a nearly instantaneous signal is generated by particles exceeding the Cherenkov threshold. In order to achieve the ultimate time resolution, we carried out a rigorous optimization of the radiator material and geometry, alongside the efficiency of the optical coupling to the SiPM sensors. The key factors limiting the time resolution were characterized by comprehensive Monte Carlo simulations, subsequently validated against experimental beam test data. We assembled small-scale prototypes instrumented with various Hamamatsu SiPM arrays sensors with pitches ranging from 1.3 to 3 mm coupled with various window materials, such as fused silica and MgF2, featuring various thickness values. The prototypes were successfully tested in beam test campaigns at the CERN-PS T10 beam line. The data were collected with a complete chain of front-end and readout electronics based on either the Petiroc 2A or the Radioroc 2 interfaced to a picoTDC to measure charges and times. By comparing the time measurements with two SiPM arrays we were able to measure a time resolution better than 33.2 ps at the full system level with a charged particle detection efficiency of 100%. Our results demonstrate the expected performance benchmarks for the charged particle detection efficiency and time resolution and highlight the potential of the developed Cherenkov-based TOF detectors for next-generation particle identification systems.