Papers for Wednesday, Jan 21 2026

We investigate the orbital dynamics of four-planet systems consisting of Earth-mass planets on initially-circular, coplanar orbits around a star of one solar mass. In our simulations, the innermost planet's semimajor axis is set at 1 AU, with subsequent semimajor axes spaced equally in terms of planets' mutual Hill radii. Several sets of initial planetary longitudes are investigated, with integrations continuing for up to $10^{10}$ orbits of the innermost planet, stopping if a pair of planets pass within 0.01 AU of each other or if a planet is ejected from the system. We find that the simulated lifetimes of four-planet systems follow the general trend of increasing exponentially with planetary spacing, as seen by previous studies of closely-spaced planets. Four-planet system lifetimes are intermediate between those of three- and five-planet systems and more similar to the latter. Moreover, as with five-planet systems, but in marked contrast to the three-planet case, no initial semimajor axes spacings are found to yield systems that survive several orders of magnitude longer than other similar spacings. First- and second-order mean-motion resonances (MMRs) between planets correlate with reductions in system lifetimes. Additionally, we find that third-order MMRs between planets on neighboring orbits also have a substantial, though smaller, destabilizing effect on systems very near resonance that otherwise would be very long-lived. Local extrema of system lifetimes as a function of planetary spacing occur at slightly smaller initial orbital separation for systems with planets initially at conjunction relative to those in which the planets begin on widely-separated longitudes. This shift is produced by the asymmetric mutual planetary perturbations as the planets separate in longitude from the initial aligned configuration that cause orbits to spread out in semimajor axis.

Direct observations of exoplanets probe the demographics and atmospheric composition of young self-luminous companions, yielding insight into their formation and early evolution history. In the near future, Gaia will reveal hundreds of nearby young exoplanets amenable to direct follow-up observations. Long-baseline interferometry with current and future facilities is best capable of exploiting this unique synergy which is poised to deliver a statistical sample of benchmark planets with precise dynamical masses and in-depth atmospheric characterization. This will enable tackling the longstanding question of how giant planets form from multiple angles simultaneously, shining light on the complex physical processes underlying planet formation.

We present a new free-form lens model for the $z=0.348$ galaxy cluster AS1063, based on previously spectroscopically confirmed lensed galaxies and new images from the GLIMPSE program. We use the ultra-deep data to identify new counterimages for previously confirmed (spectroscopically) lensed systems. We use the full set of spectroscopically confirmed systems to derive a new lens model, which is later used to confirm many of the previous lensed system candidates and discover new lensed system candidates in the images. We compute the geometric redshifts, time delays, and magnification for all counterimages (confirmed and not confirmed). Among the new systems we find a peculiar multiply lensed galaxy with a strong emission line at $\approx 4\, \mu$m that likely corresponds to H$-\beta$ and/or OIII at $z\approx 7.5$. This galaxy could be a little-red-dot or an extreme emission line galaxy. We also identify a yellow supergiant lensed star candidate at $z\approx 3.1$. This star shows some similarities with previous Kaiju stars and we nickname it "Hedorah", in honor of the famous yellow-eyed Kaiju. Previous lensed stars at $z>0.1$ are either blue supergiants or red supergiants, making Hedorah the first yellow supergiant discovered beyond $z=0.1$ and confirming that, despite their rarity, they can also be found at these redshifts. Since many Cepheid stars are yellow supergiants, we consider the possibility that Hedorah could also be the first Cepheid discovered at cosmological distances, but we conclude that Hedorah is more likely a hypergiant yellow star approaching the end of its life. Alternatively, Hedorah could be a small group of stars, although this is less likely based on Hedorah's peculiar colors and additionally may require the more exotic fuzzy dark matter to help explain the lack of counterimage.

Photometric light curves suffer from fundamental degeneracies that limit surface information recovery. We demonstrate that astrometry enables access to complementary information through photocentre variations induced by rotating surface features. The forthcoming commissioning of microarcsecond-precision astrometric missions presents an opportunity to improve stellar surface mapping. This paper extends a previous theoretical framework for stellar surface mapping, along three primary directions: (1) we derive analytical selection rules showing that astrometry is sensitive to spherical harmonic modes not detectable via photometry, particularly odd-$\ell$ modes that encode north-south asymmetries; (2) we quantify the information content of combined photometric and astrometric observations, showing that the rank of observable modes grows faster for combined observations than for either technique alone, though the fraction of recoverable modes still decreases asymptotically with increasing spatial resolution; and (3) we reframe astrometric jitter-traditionally treated as noise in exoplanet studies-as a signal encoding stellar surface structure. Given the limited proposed target lists of high-precision astrometric missions, this capability is particularly valuable: understanding host star surfaces is crucial for both removing stellar signals from exoplanet detections and characterising star-planet interactions. We show that while Sun-like stars require sub-microarcsecond precision, evolved stars with angular diameter and larger spots present immediate opportunities with current technology, such as the Gaia mission.

The outer regions of planetary systems host dusty debris discs analogous to the Kuiper belt (exoKuiper belts), which provide crucial constraints on planet formation and evolution processes. ALMA dust observations have revealed a great diversity, and that some belts contain CO gas, whose origin and implications are uncertain. Most of this progress, however, has been limited by low-resolution observations. We conducted the first ALMA large programme dedicated to debris discs: the ALMA survey to Resolve exoKuiper belt Substructures (ARKS). We selected the 24 most promising belts to constrain their detailed radial and vertical structure, and to characterise the gas content. We constrained the radial and vertical distribution of dust, as well as the presence of asymmetries. For a subset of six belts with CO gas, we constrained the gas distribution and kinematics. To interpret these observations, we used a wide range of dynamical models. The first ARKS results are presented as a series of ten papers. We discovered that up to 33% of our sample exhibits multiple dusty rings. For highly inclined belts, we found that non-Gaussian vertical distributions are common and are indicative of multiple dynamical populations. We also found that 10 of the 24 belts present asymmetries. We find that the CO gas is radially broader than the dust, but this could be an effect of optical depth. At least one system shows non-Keplerian kinematics due to strong pressure gradients, which may have triggered a vortex that trapped dust in an arc. Finally, we find evidence that the micron-sized grains may be affected by gas drag in gas rich systems. ARKS has revealed a great diversity of structures in exoKuiper belts that may arise when they are formed in protoplanetary discs or subsequently via interactions with planets and/or gas. We encourage the community to explore the reduced data and data products.

Debris discs were long considered to be largely gas-free environments governed by collisional fragmentation, gravitational stirring, and radiative forces. Recent CO detections show that gas is present, but its abundance and origin remain uncertain. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) revealed a narrow gas and dust ring in the disc HD 121617 with an asymmetric arc 40% brighter than the rest of the ring. We aim to constrain the total gas mass in HD 121617 assuming the dust arc is produced by hydrodynamical gas-dust interactions. We used the Dusty FARGO-ADSG code, modelling dust as Lagrangian particles, including radiation pressure and dust feedback, and varying the total gas mass. Simulations were compared to observations using radiative transfer. An unstable gas ring creates a size-dependent radial and azimuthal dust trap whose efficiency depends on gas mass. Two models, with 50 and 5 Earth masses of gas, reproduce both the ALMA band 7 arc and the outward offset of the VLT/SPHERE scattered-light ring via gas drag and radiation pressure. We infer a conservative gas-mass range of 2.5 to 250 Earth masses. If the ALMA asymmetry is caused by gas drag, the required gas mass compared with the observed CO implies substantial H2, consistent with primordial gas. HD 121617 would then be a hybrid disc between protoplanetary and debris stages. Since a planet could also create an arc, future observations are needed to distinguish these scenarios.

ExoKuiper belts around young A-type stars often host CO gas, whose origin is still unclear. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) includes 6 of these gas-bearing belts, to characterise their dust and gas distributions and investigate the gas origin. As part of ARKS, we observed the gas-rich system HD121617 and discovered an arc of enhanced dust density. In this paper, we analyse in detail the dust and gas distributions and the gas kinematics of this system. We extracted radial and azimuthal profiles of the dust (in the millimetre and near-infrared) and gas emission ($^{12}$CO and $^{13}$CO) from reconstructed images. To constrain the morphology of the arc, we fitted an asymmetric model to the dust emission. To characterise the gas kinematics, we fitted a Keplerian model to the velocity map and extracted the azimuthal velocity profile by deprojecting the data. We find that the dust arc is narrow (1-5 au wide at a radius of 75 au), azimuthally extended, and asymmetric; the emission is more azimuthally compact in the direction of the system's rotation, and represents 13% of the total dust mass (0.2$M_\oplus$). The arc is much less pronounced or absent for small grains and gas. Finally, we find strong non-Keplerian azimuthal velocities at the inner and outer wings of the ring, as was expected due to strong pressure gradients. The dust arc resembles the asymmetries found in protoplanetary discs, often interpreted as the result of dust trapping in vortices. If the gas disc mass is high enough ($\gtrsim20M_\oplus$, requiring a primordial gas origin), both the radial confinement of the ring and the azimuthal arc may result from dust grains responding to gas drag. Alternatively, it could result from planet-disc interactions via mean motion resonances. Further studies should test these hypotheses and may provide a dynamical gas mass estimate in this CO-rich exoKuiper belt.

CO gas is detected in a significant number of debris discs, but its origin and evolution remains unclear. Key constraints are its mass and spectro-spatial distribution, which are coupled through optical depth and have only been analysed at low to moderate resolution so far. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is the first ALMA large program to target debris discs at high spectro-spatial resolution. We used $^{12}$CO and $^{13}$CO J=3-2 line data of 18 ARKS debris belts, 5 of which were already known to host gas, to analyse the spectro-spatial distribution of CO, constrain the gas masses, and to search for gas in the remaining systems. We developed a line-imaging pipeline and produced line cubes for each disc, with a spatial resolution down to $\sim$70 mas and spectral resolution of 26 m s$^{-1}$. Using spectro-spatial shifting and stacking, we produced high signal-to-noise maps, and radial and spectral profiles that reveal the distribution and kinematics of gas in 5 gas-bearing discs. For these discs, we constrained the inner radius of the $^{12}$CO, and found the radial brightness profile of CO peaked interior to the dust ring, but that CO was more radially extended than the dust. We present the first radially resolved $^{12}$CO/$^{13}$CO isotopologue flux ratios in gas-bearing debris discs, which are constant with radius for the majority of systems, indicating $^{12}$CO and $^{13}$CO are both optically thick or thin throughout the discs. We report CO line fluxes/upper limits for all systems and optical depth dependant masses for the 5 gas-bearing systems. Finally, we analysed the $^{12}$CO J=3-2 line luminosities for the ARKS debris discs and discs from the literature. We confirm that gas is mostly detected in young systems. However, the high scatter seen in young/high fractional luminosity systems indicates no trend within the systems with detected gas.

We present a catalog of 643 diffuse galaxies identified through a citizen science search of the Fornax cluster, of which we estimate 21.8% are nucleated (139/637; 6 inconclusive). This marks the first crowd-sourced effort to construct a cluster-scale census of diffuse galaxies. These objects were visually identified using a combination of the Fornax Deep Survey and Dark Energy Camera Legacy Survey imaging across 26 deg$^2$. Over 1,400 volunteers cataloged the candidates within this sky area at a rate of 1.15 days/deg$^2$. Our catalog is highly complete relative to existing dwarf catalogs of Fornax ($> 80\%$ of objects recovered) down to an effective radius $r_{\mathrm{eff}} = 5^{\prime \prime}$, the minimum size we suggested volunteers classify, and to an effective r-band surface brightness as faint as $\langle \mu_r \rangle \simeq26$ mag arcsec$^{-2}$. We detect 97 candidates that existing automated searches of Fornax did not find and three candidates not found by any prior search, automated or visual. The stellar mass distribution of our sample is consistent with similar dwarf studies of Fornax, with the nucleated fraction peaking at 80% for a host galaxy mass of $\sim$10$^{8.5}M_{\odot}$. The efficiency and completeness of our catalog thus establishes citizen science as a valuable tool for mapping diffuse galaxy populations in future sky surveys, such as the Legacy Survey of Space and Time.

We present results from simultaneous X-ray polarimetric and spectroscopic observations of the bright neutron star low-mass X-ray binary Cyg X-2, performed by the Imaging X-ray Polarimetry Explorer (IXPE) and the Nuclear Spectroscopic Telescope Array (NuSTAR). IXPE detected significant polarization (15 sigma) from the source in the 2-8 keV energy band with an average polarization degree (PD) of 4.5% +/- 0.3% and a polarization angle (PA) of 128 +/- 2 degrees as the source moved along the horizontal branch of its Z-track. The PD increases with energy reaching 9.9% +/- 2.8% in the 7-8 keV band, with no evidence for energy-dependent variation in the PA. The PA is roughly consistent with previous measurements obtained during the normal and flaring branches and also with the known radio jet axis. From spectropolarimetric analysis, the main contribution to the polarized radiation is due to Comptonized photons, but the polarization is higher than predicted in typical spreading layer geometries. The observed high polarization may be due to a combination of a highly polarized reflected component and a moderately polarized spreading layer on the neutron star surface or produced by electron scattering in an equatorial wind.

Dusty discs detected around main-sequence stars are thought to be signs of planetesimal belts in which the dust distribution is shaped by collisional and dynamical processes, including interactions with gas if present. The debris disc around the young A-type star HD 131835 is composed of two dust rings at ~65 au and ~100 au, a third unconstrained innermost component, and a gaseous component centred at ~65 au. New ALMA observations show that the inner of the two dust rings is brighter than the outer one, in contrast with previous observations in scattered light. We explore two scenarios that could explain these observations: the two dust rings might represent distinct planetesimal belts with different collisional properties, or only the inner ring might contain planetesimals while the outer ring consists entirely of dust that has migrated outwards due to gas drag. To explore the first scenario, we employed a state-of-the-art collisional evolution code. To test the second scenario, we used a simple dynamical model of dust grain evolution in an optically thin gaseous disc. Collisional models of two planetesimal belts cannot fully reproduce the observations by only varying their dynamical excitation, and matching the data through a different material strength requires an extreme difference in dust composition. The gas-driven scenario can reproduce the location of the outer ring and the brightness ratio of the two rings from scattered light observations, but the resulting outer ring is too faint overall in both scattered light and sub-millimetre emission. The dust rings in HD 131835 could be produced from two planetesimal belts, although how these belts would attain the required extremely different properties needs to be explained. The dust-gas interaction is a plausible alternative explanation and deserves further study using a more comprehensive model.

Astrometric jitter noise arises when starspots on a rotating stellar surface move in and out of view, shifting the photocenter. This noise may limit our ability to detect and weigh small, sub-Neptune-sized planets around active stars. By deriving a linear forward model for the astrometric jitter signal of a rotating star in a spherical-harmonic coordinate system, we show that jitter noise can be used to reconstruct surface-brightness maps and, in principle, disentangle jitter from stellar reflex motion due to an orbiting planet. Furthermore, we show that astrometry and photometry probe complementary surface information: photometry measures even-degree spherical harmonic surfaces that are symmetric about the equator, while astrometry measures odd-degree modes. Their joint use, therefore, breaks degeneracies in surface mapping. Our model further quantifies the variation in the astrometric signal with inclination angle, which is foundational for studies of worst-case configurations of astrometric star-spot noise. For example, we show that pole-on stellar inclinations lead to poorly constrained inversions, as any stellar surface produces a purely circular astrometric jitter signal. We characterize the degeneracy in jointly identifying the stellar surface and stellar inclination, and develop a surface estimation approach. Using this approach, we present example simulations and reconstructions that demonstrate the use of astrometry data alongside light-curve data to improve stellar surface mapping and localize spot positions in latitude and longitude. With forthcoming high-precision Gaia astrometry, astrometric surface mapping provides a promising new approach to probe stellar activity.

We present Green's function solutions of the Newtonian time-dependent thin disk equations in the presence of outflows, showing that simple and exact analytical expressions exist in various natural limits of the problem. These Green's functions are mathematically very similar to the classical Lynden-Bell & Pringle solutions in the absence of outflows, but differ strongly in their precise physical details and observational implications. Solutions are presented for phenomenological radius-dependent outflows which both do and do not torque the local accretion flow, and for outflows which are launched proportional to the local accretion rate. Generically, outflows lead to a more rapid decay of the bolometric luminosity of the disk, flatten the radial dependence of the disk temperature, and suppress variability in the accretion rate at small radii and low frequencies (on long timescales). Observational implications of these four results are discussed in detail.

The chemistry of shielded molecular gas is primarily driven by energetic, charged particles dubbed cosmic rays (CRs), in particular those with energies under 1 GeV. CRs ionize molecular hydrogen and helium, the latter of which contributes greatly to the destruction of molecules. CR ionization initiates a wide range of gas-phase chemistry, including pathways important for the so-called "carbon cycle", C$^+$/C/CO. Therefore, the CR ionization rate, $\zeta$, is fundamental in theoretical and observational astrochemistry. Although observational methods show a wide range of ionization rates -- varying with the environment, especially decreasing into dense clouds -- astrochemical models often assume a constant rate. To address this limitation, we employ a post-processed gas-phase chemical model of a simulated dense molecular cloud that incorporates CR energy losses within the cloud. This approach allows us to investigate changes in abundance profiles of important chemical tracers and gas temperature. Furthermore, we analyze analytical calibrators for estimating $\zeta$ in dense molecular gas that are robust when tested against a full chemical network. Additionally, we provide improved estimations of the electron fraction in dense gas for better consistency with observational data and theoretical calibrations for UV-shielded regions.

Asymmetries in debris discs provide unique clues to understand the evolution and architecture of planetary systems.** The aim of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is to expand our understanding of radial and vertical dust structures, as well as gas distributions and kinematics, in debris discs.** Here, in ARKS~VI, we present a systematic analysis of the asymmetries and stellocentric offsets present in the ALMA continuum data for the ARKS survey. Our aims are to identify asymmetries in debris disc dust distributions, quantify debris disc asymmetry properties, and discuss the potential origins of debris disc asymmetries.** We utilised empirical methods to identify emission asymmetries** and the presence of offset emission between disc centres and the locations of the host stars, via an analysis of their calibration procedures and disc properties. We associated observational asymmetry types** and plausible physical classes** associated with each source. We show that there are ten systems, almost half of the ARKS sample, that host either a continuum emission asymmetry or offset emission. Three systems host offsets (HD15115, HD32297, and HD109573 (HR4796)), four host azimuthal asymmetries (HD9672 (49Ceti), HD92945, HD107146, and HD121617), two host an asymmetry in their major axis (HD10647 (q$^1$ Eri), and HD39060 ($\beta$ Pic)), and one hosts an asymmetry in their minor axis (HD61005). We attribute the offset asymmetries to non-zero eccentricities, and three of the azimuthal asymmetries to arcs. The presence of an asymmetry or offset in the ARKS sample appears to be correlated with the fractional luminosity of cold dust.** Conclusions: This study demonstrates that debris disc asymmetries in the ARKS sample are common, and plausibly so in the wider population of debris discs at (sub)-millimetre wavelengths.** ** = ABRIDGED FOR ARXIV: FULL ABSTRACT IN PAPER

Long-period comets, which are often considered to be representative of material in the protoplanetary disk that formed the Solar System, are ideal to investigate the question of chemical inheritance in astronomy. Determining the chemistry of comets, both individually and as a population, has become of great importance in comparative studies against sources representative of evolutionary precursors to planetary systems. Contemporaneous observations of long-period comet C/2022 E3 (ZTF) were obtained with the JWST and the Atacama Large Millimeter/submillimeter Array (ALMA) in early 2023 March. This work focuses on \ce{CH3OH} measurements from both ALMA and JWST as well as \ce{H2O} measurements from JWST. Radiative transfer modeling of \ce{CH3OH} and \ce{H2O} was performed to investigate spatial variations in rotational temperature, column density, and production rates, as well as a comparison of derived values between the two telescopes. Most of the spatial distributions of the modeled values are centrally peaked, and the modeled values from JWST are all within the error bars of the average values from ALMA. C/2022 E3 (ZTF) also displays an enhancement in modeled rotational temperature in the anti-Sunward direction that is shown to be statistically significant. Based on non-LTE radiative transfer modeling, the declining \ce{H2O} rotational temperatures as a function of nucleocentric distance observed by JWST can be explained primarily as a result of rotational line cooling. The values derived in this work are in general agreement with single-dish millimeter-wave observations.

EP J171159.4-333253 is a new neutron-star low-mass X-ray binary discovered in outburst by the Einstein Probe (EP) on 2025 June 23, exhibiting clocked type-I X-ray bursts, eclipses and dips. In this paper, we report on the results of the X-ray spectral and timing analyses for EP J171159.4-333253 using data collected by EP and NuSTAR during the first 21 days of the outburst. The X-ray burst recurrence time can be characterized over a subset of nine bursts spanning 1.6 days around the NuSTAR observation, and the result is $t_{\rm rec}=8196 \pm 177\,$s with indications of a possible decreasing trend. From the X-ray eclipse events, the binary orbital period and the eclipse duration are estimated to be $P_{\rm orb}=6.48301 \pm 0.00003\,$hr and $D_{\star,X} = 1245.5^{+6.9}_{-6.5}\,$s, respectively. These enable an estimate of the mass and radius of the companion star and the binary inclination, which are $M_2\approx0.6-0.8\,M_\odot$, $R_2\approx0.7-0.8\,R_\odot$ and $i\approx73-75^\circ$, respectively. We also report on joint ULTRACAM and EP observations on 2025 July 21--22, detecting the source optical counterpart and covering an eclipse in both X-ray and optical bands. The optical eclipse is wavelength-dependent and broader than in X-rays, indicating that part of the optical emission arises from an extended region in the accretion flow. Despite a moderate variation in the source flux, the properties of the persistent X-ray emission are typical of a hard spectral state. We further evaluated the ratio of the accretion energy to the thermonuclear energy to be 120--130, implying helium bursts with the accreted hydrogen being depleted in-between bursts.

Recent research raises the possibility that 3 and 7 million years ago, the Sun encountered massive clouds that shrank the heliosphere--the solar cocoon protecting our solar system--exposing Earth to its interstellar environment, in agreement with geological evidence from 60Fe and 244Pu isotopes. Here we show that during such encounters Earth was exposed to increased radiation in the form of high-energy particles. During periods of Earth's immersion in the heliosphere, it received particle radiation that we name Heliospheric Energetic Particles (HEPs). The intensity of < 10 MeV protons was at least an order of magnitude more intense than today's most extreme solar energetic particle (SEP) events. SEPs today last minutes to hours, but HEP exposure then lasted for extensive periods of several months, making it a prolonged external driver. During Earth's excursion outside the heliosphere, it was exposed to a galactic cosmic ray radiation with the intensity of < 1 GeV protons at least an order of magnitude more intense than today. Therefore, the space surrounding Earth was permeated by a variable high-energy radiation. We discuss the implications for Earth's climate and biodiversity.

The Habitable Worlds Observatory is the first telescope ever designed to search for life and will be a powerhouse of discovery across topics in astrophysics. The observatory was the top recommendation of the Astro2020 Decadal Survey for large missions and a new HWO Technology Maturation Project Office was formed in August 2024 to mature the architecture, science and technology. In this paper we review the overall approach taken to mature the mission concept. We show progress on architecture development, integrated modeling, science cases, and technology roadmaps consistent with pre-formulation studies. We discuss plans for instrument studies and international engagement and science engagement including a Community Science and Instrument Team. Finally, we describe the plan forward to the Mission Concept Review.

CO gas has been detected in $\sim$20 debris discs. We present ALMA observations of the CO-rich HD 121617 debris disc from the ARKS survey. Using high-resolution Band 7 observations of $^{12}CO \ J=3-2$, we analyse local CO line profiles to investigate optical depth, CO mass, and temperature. Spectra are aligned and stacked in concentric annuli to produce local line profiles. The resulting profiles are Gaussian-shaped and broadened by Keplerian shear. The line profiles are modelled using both a simplified toy model and a RADMC-3D model including projection effects and Keplerian shear. Fitting the RADMC-3D model to the $^{13}$CO data, we find that an optically thick model with a temperature of 38 K and a CO mass of $2 \times 10^{-3}$ M$_{\oplus}$ reproduces the observations. The model reproduces the enhanced emission at orbital azimuths of $\sim \pm45^{\circ}$ and $\pm135^{\circ}$, forming an X-shaped structure in the velocity-integrated intensity map, as well as the broader $^{12}$CO linewidth relative to $^{13}$CO. Scaling the model by the ISM abundance ratio ($\sim$77) also reproduces the $^{12}$CO data, though high optical depths and model assumptions limit mass constraints. We find that azimuthally averaged local line profiles appear Gaussian regardless of optical depth, cautioning against their use for distinguishing optically thin and thick emission. We constrain the mean molecular weight to $12.6_{-1.1}^{+1.3}$, dependent on model assumptions. Our $^{13}$CO results suggest that C$^{18}$O may also be optically thick in CO-rich debris discs and that the mean molecular weight is significantly higher than if H$_2$ were the dominant gas species, suggesting a non-primordial composition.

Observations of molecular structures on scales of $\sim 0.1-50$ pc show that the specific angular momentum ($j$) scales with radius ($R$) as $j\sim R^{3/2}$. We study the effects of turbulence, gravity, and the magnetic field in shaping this scaling, by measuring clump size and specific angular momentum in three SPH simulations of the formation of giant molecular clouds, progressively adding these three ingredients. In each simulation, we define ``full'' and ``reduced'' clump samples, the latter restricted to aspect ratios $A<3$. We find that, in the non-magnetic runs, elongated clumps deviate the most from the \jR\ relation, which is best reproduced by the reduced sample in the gravity+turbulence run. In the purely hydrodynamic case, no dense elongated structures form, suggesting that turbulence alone is insufficient to generate dense filaments, although clumps have $j$ magnitudes consistent with observations. In the gravity+turbulence+magnetic field run, most of the clumps are filamentary, yet the full sample appears to follow the observed \jR\ relation. This result, rather than being a real trend, could be the combination of the increase in $j$ by the filamentary geometry, and its reduction by turbulence inhibition by the magnetic field. Finally, we measure the gravitational, magnetic, pressure-gradient, and hydrodynamic torques (which involve turbulent viscosity) in our clump samples. We find that, in magnitude, the hydrodynamic torques tend to be larger than the rest. This result is consistent with our previous work, where we proposed that gravity drives cloud formation and contraction, while turbulence redistributes angular momentum through fluid-parcel exchanges.

Gravitational-wave (GW) observatories have used template-based search to detect hundreds of compact binary coalescences (CBCs). However, template-based search cannot detect astrophysical sources that lack accurate waveform models, including core-collapse supernovae, neutron star glitches, and cosmic strings. Here, we present a novel approach for template-free search using coincident anomaly detection (CoAD). CoAD requires neither labeled training examples nor background-only training sets, instead exploiting the coincidence of events across spatially separated detectors as the training loss itself: two neural networks independently analyze data from each detector and are trained to maximize coincident predictions. Additionally, we show that integrated gradient analysis can localize GW signals from the neural-network weights, providing a path toward data-driven template construction of unmodeled sources, and further improving precision by frequency matching. Using the CodaBench dataset of real LIGO backgrounds with injected simulated CBCs and sine-Gaussian low-frequency bursts, CoAD achieves recall up to 0.91 and 0.85 respectively at a false-alarm rate of one event per year, and achieves recall above 0.5 at signal-to-noise ratios below 10. The fully-unsupervised nature of CoAD makes it especially well-suited for next-generation detectors with greater sensitivity and associated increases in GW event rates.

As the only known unambiguous star in a Maunder minimum-like chromospheric activity state, the properties of HD 166620 can provide valuable insight into the behaviour of the Sun during the historic extended low-states of its activity cycle. The coronal X-ray activity of HD 166620 has so far only been probed with a ROSAT/HRI observation in 1996, near the chromospheric activity maximum before the star entered its grand minimum around 2004. We conducted a deep {\it XMM-Newton} observation of HD 166620 during its chromospheric Ca II H&K activity grand minimum to achieve a better understanding of its magnetic activity. We detected HD 166620 with an X-ray luminosity of ${{\rm log}\,L_{\rm X}\,\rm{(erg\,s^{-1})}=26.56^{+0.10}_{-0.12}}$, corresponding to ${{\rm log}\,(L_{\rm X}/L_{\rm bol}) = -6.58^{+0.10}_{-0.12}}$ and an X-ray surface flux of log Fx (erg/cm^2/s) = 3.97+0.10-0.12. With respect to the earlier ROSAT observation, the X-ray brightness of HD 166620 has decreased by a factor of 2.5 during its Maunder minimum-like state. To place its X-ray properties into context, we constructed an X-ray sample of late-type stars within 10 pc of the Sun. The activity of HD 166620 is below the levels of all other K dwarfs in the 10 pc sample. The corona of HD 166620 during its grand minimum emits at the level of the solar background corona, which implies that it has no large active magnetic structures. Along with long-term Ca II H&K monitoring of HD 166620, this result provides evidence that the solar activity during the Maunder minimum was not reduced significantly below the levels seen during its present-day cycle minima. The similar X-ray surface flux of HD 166620 and the modern quiet Sun, and also their Rossby number near the critical value of spin-down models, suggest a connection between the regime of weakened magnetic braking and the occurrence of Maunder minimum states.

Escape trajectories from the Earth-Moon system play an important role in interplanetary transfer. This paper focuses on the escape trajectories from a 167 km circular Earth orbit in the Earth-Moon planar circular three-body problem and the Sun-Earth/Moon planar bicircular four-body problem and is denoted to providing a comprehensive analysis on these escape trajectories. To achieve these purposes, the global maps of escape trajectories are constructed, and escape trajectories with one lunar gravity assist are pre-filtered. Then, an effective method to identification escape families is proposed based on dynamical analysis and data mining techniques. Once the escape families are identified, the corresponding characteristics are analyzed to provide insights into the construction of escape trajectories. Based on these escape families, the effects of the solar gravity perturbation on the number of escape trajectories, the emergence and disappearance of escape families, variation in generalized energy, and transfer characteristics are further summarized, providing insights into the model selection in the escape trajectory construction. This paper establishes an analysis methodology of escape trajectories from a perspective of escape families, deepening the understanding of escape dynamics.

Fast radio bursts (FRBs) are mysterious radio transients with uncertain origins and environments. Recent studies suggest that some active FRBs may originate from compact objects in binary systems. In this work, we develop a unified theoretical framework to model the multi-wavelength afterglows of FRBs resided in binary systems and apply it to two representative repeaters, FRB 20200120E and FRB 20201124A. By solving the dynamics and radiation processes of FRB ejecta interacting with the surrounding medium, we compute afterglow light curves in the radio, optical, and X-ray bands. Our results show that radio afterglows offer the best prospects for detection, with their brightness highly sensitive to ejecta kinetic energy and ambient density. Future high-sensitivity radio telescopes, such as the Square Kilometre Array (SKA), could detect these signals. Optical afterglows, though short-lived and challenging to observe, may be significantly enhanced in dense environments, potentially making them detectable with facilities like the Large Synoptic Survey Telescope (LSST). In contrast, X-ray afterglows are predicted to be too faint for detection with current instruments. Our study highlights the potential of multi-wavelength afterglows as probes of FRB progenitors and their surrounding environments, offering crucial insights into the nature of these mysterious transients.

Fast Radio Bursts (FRBs) are millisecond-duration radio transients of mysterious origin, with growing evidence linking at least some of them to magnetars. While FRBs are primarily observed in the radio band, their potential multi-wavelength afterglows remain largely unexplored. We investigate the possible afterglow of FRB 20171020A, a rare nearby and bright FRB localized in a galaxy at only 37 Mpc. Assuming that this source produces a future bright burst, we model the expected afterglow emission in the radio, optical, and X-ray bands under both uniform and wind-like ambient media, within the framework of the magnetar model. Our results show that the optical afterglow is the most promising for detection, but it fades rapidly and requires follow-up within a few hundred seconds post-burst. The radio afterglow may be detectable under favorable conditions in a dense stellar wind, whereas the X-ray counterpart is too faint for current telescopes. These findings suggest that rapid optical follow-up offers the best opportunity to detect the afterglow of the next bright burst from FRB 20171020A, providing unique insights into the progenitor and its environment. To assess observational feasibility, we estimate the event rate of nearby FRBs with sufficient energy to power detectable afterglows, finding a rate of $\sim$0.3 per year for CHIME surveys. Although this rate is low and the optical detection timescale is short, coordinated fast-response strategies using global telescope networks could significantly improve the chance of success. As more nearby FRBs are discovered, multi-wavelength observations will be essential in unveiling the physical nature of these enigmatic events.

The neutron-star X-ray source 1E 161348-5055, associated with the supernova remnant RCW 103, exhibits clear intensity variations with a period of 6.67 hr. To clarify the nature of this object and its long periodicity, detailed timing studies were applied to its archival X-ray data, taken with ASCA (in 1993), XMM-Newton (in 2001, 2005, and 2016), and NuSTAR (2016 and 2017). It was assumed that the 6.67 hr period arises due to the beat between the rotation and free precession periods of the star that is slightly aspherical. By removing timing perturbations to be caused by this long periodicity, the six data sets consistently yielded evidence for pulsations at periods of P~1.01 s, to be interpreted as the objects' spin period, although the optimum energy range differed among the data sets. The measured six periods accurately line up on a linear spin-down trend of dP/dt = 1.097x 10^{-12} s/s. The object is implied to have a characteristic age of 14.7 kyr, a spin-down luminosity of 4.2x10^{34} erg/s, which is insufficient to power the X-ray luminosity, a dipole magnetic field of ~4.6x10^{13} G, and a toroidal field of ~7 x10^{15} G. Its similarity and dissimilarity to magnetars are discussed. An emission geometry, which crudely explain these results, is presented.

A small body orbiting around an accreting massive object and periodically crossing its accretion disk is a common configuration in astrophysics. In this work, we investigate the secular evolution of extreme mass-ratio inspirals (EMRIs), in which a stellar-mass object (SMO), e.g., a star or a stellar-mass black hole (sBH), collides with the accretion disk of a central supermassive black hole (SMBH), within a fully relativistic framework. We find (1) the disk always tends to align the SMO no matter what the initial orbital inclination $\iota$ relative to the disk is, (2) the final orbital eccentricity of the SMO captured by the disk is always low though the orbital eccentricity may temporarily grow when the orbital inclination $\iota$ is large and the SMO is an sBH, and (3) via collisions with the accretion disk only, only a small fraction of sBHs that are initially close to the SMBH and close to the disk can be captured by the disk within typical disk lifetime of active galactic nuclei. Two-body scatterings between SMOs in the nuclear stellar cluster play an essential role in randomly kicking sBHs towards the disk and boosting the capture rate.

At the base of the Sun's convective zone, a narrow shear layer called the tachocline separates strong latitudinal differential rotation above from nearly rigid rotation in the radiative zone below. The observed thinness of the tachocline is a long-standing dynamical puzzle because the tachocline should have spread significantly due to inward-burrowing meridional circulation, also called "radiative spreading." We recently presented the first pair of global simulations to reveal a statistically stationary tachocline confined against radiative spreading by the Maxwell stresses from the nonaxisymmetric modes of a dynamo, which penetrated into and below the tachocline through a novel magnetic skin effect. In the work presented here, we systematically examine how this "dynamo confinement scenario" works against radiative spreading in a suite of simulations as the governing parameters trend in the direction of the true solar regime. We find that as the stable stratification of the radiative zone is made progressively stronger, the dynamo cycles get longer, the magnetic field consequently penetrates deeper due to the skin effect, and the tachocline becomes more confined. Furthermore, these results have interesting consequences for solar spin-down. In all of our radiatively spreading simulations, the tachocline region spins down due to the burrowing circulation. Below the tachocline, the Maxwell stresses transmit this spin-down further to rigidify the deeper radiative zone. We thus speculate that, in addition to confining the tachocline, the dynamo may provide a pathway to communicate spin-down from the near-surface layers to the deep interior.

A transient (LPA J0108+13) with repeated bursts was detected on the Large Phased Array (LPA) radio telescope at a central frequency of 110.4 MHz in the direction of the radio galaxy 3C 33. The flux density of bursts ranges from tens to hundreds of Jy, and the duration of the bursts is \approx 1^m - 4^m. In daily observations, the total duration of which at the location of the transient exceeds 200 hours in the observation interval 2013-2025, 6 bursts were found. The nature of the source could not be determined. We believe that a new type of transients has been discovered.

Herbig Ae/Be stars (HAeBes) are critical tracers of intermediate- and high-mass star formation, yet their census remains incomplete compared to low-mass young stellar objects like T-Tauri stars. To expand the known population, we systematically searched for HAeBes in LAMOST DR7 low-resolution spectra. Following Sun et al., we applied Uniform Manifold Approximation and Projection (UMAP) for dimensionality reduction and Support Vector Machine (SVM) classification, identifying $\sim$240,000 spectra with potential H$\alpha$ emission. After removing contaminants (non-stellar objects, extragalactic sources, CVs, and Algol systems) and restricting to B/A-type stars, we obtained 1,835 candidates through 2MASS/WISE visual inspection. Spectral energy distribution analysis confirmed 143 sources with infrared excess ($J$-band or longer wavelengths), including 92 known HAeBes. From the remaining 51 candidates, we classified 26 with strong infrared excess as new HAeBes. Color-index analysis of confirmed HAeBes and classical Ae/Be stars (CAeBes) revealed that the $(K-W1)_0$ vs. $(W2-W3)_0$ diagram effectively separates these populations: CAeBes predominantly occupy $(K-W1)_0 \leq 0.5$ and $(W2-W3)_0 \leq 1.1$, while other regions trace transition disks ($(K-W1)_0 < 0.5$ and $(W2-W3)_0 > 1.1$), globally depleted disks ($(K-W1)_0 > 0.5$ and $(W2-W3)_0 < 1.1$), and Class I/Flat/II HAeBes ($(K-W1)_0 > 0.5$ and $(W2-W3)_0 > 1.1$). More importantly, the HAeBes exhibit a clear evolutionary gradient on this diagram, with those in the Class III, Class II, Flat-SED, and Class I evolutionary stages being effectively distinguished by concentric ellipses that are roughly centered at (0,0) with semi-major axes of $a$=1.5, $a$=3.0, and $a$=4.0, and a semi-major to semi-minor axis ratio of 1.6:1.

This study examines the dynamics of the third body in an elliptic restricted three-body problem (ERTBP) framework, taking into account perturbations from radiation pressure, oblateness, and elongation of the primary bodies, as well as disk-like structures. The objectives are to determine the positions and stability of the equilibrium points, asses how these points shift under the influence of perturbations, and evaluate the dependence of their stability on the orbital eccentricity and perturbation parameters. The ERTBP model is modified to include a radiating, oblate primary body and an elongated secondary body modeled as a finite straight segment, alongside perturbations from a surrounding disk. The system's equations of motion are numerically solved using parameters from perturbed and classical cases. Equilibrium positions are computed over a range of eccentricities and perturbation values, and stability is analyzed using linearized equations and eigenvalue methods. In all cases, we have found three collinear ($L_1$, $L_2$, $L_3$) and two non-collinear ($L_4$, $L_5$) equilibrium points solutions. The inclusion of radiations, oblateness, elongation using a finite straight segment, and disk perturbation systematically displaces each equilibrium point from its classical location, with the magnitude and direction of the displacement varying with the perturbation parameter. Stability analysis confirms that the collinear points remain linearly stable under all tested conditions. Meanwhile, non-collinear points are stable under a specific condition. We investigate the stability boundary of these points as a function of orbital eccentricity and we found there is a critical range of eccentricity values within which stability is preserved.

In the core accretion scenario, forming planets start to acquire gaseous envelopes while accreting solids. Conventional one-dimensional models assume envelopes to be static and isolated. However, recent three-dimensional simulations demonstrate dynamic gas exchange from the envelope to the surrounding disk. This process is controlled by the balance between heating, through the accretion of solids, and cooling, which is regulated by poorly-known opacities. In this work, we systemically investigate a wide range of cooling and heating rates, using three-dimensional hydrodynamical simulations. We identify three distinct cooling regimes. Fast-cooling envelopes ($\beta \lesssim 1$, with $\beta$ the cooling time in units of orbital time) are nearly isothermal and have inner radiative layers that are shielded from recycling flows. In contrast, slow cooling envelopes ($\beta\gtrsim10^3$) become fully convective. In the intermediate regime ($1\lesssim\beta\lesssim300$), envelopes are characterized by a three-layer structure, comprising an inner convective, a middle radiative, and an outer recycling layer. The development of this radiative layer traps small dust and vapour released from sublimated species. In contrast, fully convective envelopes efficiently exchange material from inner to outer envelope. Such fully convective envelopes are likely to emerge in the inner parts of protoplanetary disks ($\lesssim$ 1 au) where cooling times are long, implying that inner-disk super-Earths may see their growth stalled and be volatile depleted.

Supernovae are observed to occur approximately 1-2 times per century in a galaxy like the Milky Way. Based on historical records, however, the last core-collapse galactic supernova observed by humans occurred almost 1,000 years ago. Luckily, we are well positioned to catch the next one with the advent of new neutrino detectors and astronomical observatories. Neutrino observatories can provide unprecedented triggers for a galactic supernova event as they are likely to see a supernova neutrino signal anywhere from minutes to days before the shock breakout causes the supernova to brighten in optical wavelengths. Given its large etendue, the Vera C. Rubin Observatory is ideally positioned to rapidly localize the optical counterpart based on the neutrino trigger. In this paper we simulate events to study the efficiency with which supernovae are optimally localized by the Vera C. Rubin Observatory. We find that the observatory is ideal for initial localization of nearly all observable supernova triggers and has a 57-97% chance of catching any supernova based on theoretical stellar mass density predictions and observations. We provide an analysis of optimal filter selection and exposure times and discuss observational caveats.

Coherent oscillations at kilohertz frequencies have recently been detected in a small number of gamma-ray bursts (GRBs), suggesting quasi-periodic dynamics in their central engines. A prominent example is GRB~230307A, which exhibited a brief, highly coherent, energy-dependent periodic signal interpreted as the possible spin signature of a nascent millisecond magnetar formed after a compact binary merger. Motivated by these developments, we conducted a comprehensive search for similar signals, accounting for both temporal and spectral dependencies, in 532 short GRBs with time-tagged event data recorded by the Burst and Transient Source Experiment (BATSE) onboard the \textit{Compton Gamma-Ray Observatory}. Within this sample, we identify a single statistically significant case: GRB~960616 (BATSE trigger~5502), in which the $\sim$30~ms main emission episode is coherently modulated at 1100~Hz, with the strongest modulation above 320~keV and a fractional amplitude of $\sim$47\%. Assuming the presence of a coherent periodic modulation, we use data-driven Monte Carlo simulations to place an upper limit of $\sim$8\% on the fraction of the total radiated energy that can be modulated by the QPO. This event, exhibiting a periodicity at $\sim$0.91~ms, further supports the possibility that millisecond periodicities can arise during GRBs in merger-driven scenarios.

Intermediate polars are a class of cataclysmic variables in which a white dwarf accretes material from a companion star. The intermediate polar nature confirmation usually derives from the detection of two periods in both $X$-ray and optical photometry. In this respect, the high energy signal is often characterized by modulations on the white dwarf spin and the orbital period. However, noting that the periodograms may be characterized by strong features also at the synodic period and/or other sidebands, the timing analysis of the $X$-ray signal may offer the unique possibility to firmly discover an intermediate polar candidate. Here, we concentrate on a sample of five cataclysmic variable {binary} candidates: {i.e. SAXJ1748.2-2808, 1RXS J211336.1+542226, CXOGC J174622.7-285218, CXOGC J174517.4-290650, and V381 Vel, listed in the IPHome catalogue. Our main aim is to confirm if they belong to the intermediate polar class or not. The results of our analysis show that we can safely assess the intermediate polar nature of all the considered sources, apart for the case of V381 Vel which instead behaves like a cataclysmic variable of the polar subclass. Moreover, the source SAXJ1748.2-2808, previously classified as a HMXB, appears to be, most likely, an intermediate polar variable.

Pulsar glitches are generally viewed as stochastic events driven by sudden angular momentum transfer from the neutron star's superfluid interior to its crust. Except two peculiar pulsars with quasi-periodic glitches, this stochastic view has prevailed. Here, by clustering temporally proximate small glitches of the Crab pulsar, we uncover clear evidence of an underlying quasi-periodic modulation, challenging the paradigm of purely random behavior. Furthermore, our correlation analyses reveal a strong positive relationship between glitch cluster size and waiting time since the preceding clusters. These findings demonstrate the effect of angular momentum reservoir operating over long-term scales and enable the predictions of next glitching window. Remarkably, two minor glitches detected in July and August 2025, which align with our initial prediction made in June, should be confirmed as the onset of this predicted activity. Inspired by the initial success, we forecast the occurrence of a major glitch from now until August 2026, with possible glitch size up to a relative change in rotational frequency of $697.2 \times 10^{-9}$. Physically, the observed long-term quasi-periodicity and cluster size-waiting time correlations imply that each glitch event releases only a fraction of the stored superfluid angular momentum. This partial-release mechanism provides a unified framework for both stochastic and quasi-periodic glitch behaviors across different pulsars, underscoring the universality of the superfluid angular momentum reservoir effect. As the most intensively monitored object, the Crab pulsar serves as a natural laboratory for studying angular momentum inside neutron stars.

In this Letter, we report the detection of soft X-ray time lags-i.e. variability in the softer photons lagging behind that in the harder photons-in seven XMM-Newton observations of the tidal disruption event (TDE) candidate AT2021ehb. We find correlated variability between the soft (0.3-0.7 keV) and hard (0.9-10 keV) bands on about 10^4 s time-scales, and measure a soft lag of about 500 s. This behaviour is broadly consistent with the disk-corona reverberation scenario established in active galactic nuclei (AGNs). Together with the previously reported strong hard X-ray emission and broad Fe K line, our results suggest the presence of a compact corona and prominent relativistic disk reflection in AT2021ehb. The unusually high blackbody temperature (peaking at about 200 eV) is difficult to reconcile with thermal emission from a standard accretion disk around a about 10^7 Msun black hole, and may instead be analogous to the soft excess commonly observed in AGNs, whose physical origin remains debated. Finally, the measured lags offer a possible explanation for the rapid X-ray flux decline that occurred only three days after the peak, pointing to a scenario in which the corona cools following a sudden loss of the magnetic support required to sustain it.

Millisecond pulsars (MSPs) have been proposed as evolutionary products of low-mass X-ray binaries (LMXBs) through a stage in which they are spider pulsars (i.e., redbacks and black widows). However, recent work has found that the systemic kicks of observed MSPs are significantly lower than the kicks of LMXBs and spiders, which appears to be in tension with this evolutionary model. We argue that this tension can be relieved, at least to some degree, by considering the fact that the observed MSPs are located at relatively short distances, whereas spider pulsars are located at greater distances and LMXBs are situated even further away. We model the distance-dependent kinematic bias for dynamically old objects, which favors observing objects that have received low kicks at short distances and correct the observed systemic kicks for this bias. We find that this kinematic bias can be big enough to close the gap between the MSP and LMXB kicks, although the spider pulsars appear to come from a slightly different systemic kick distribution, but this difference is not necessarily physical. All corrected systemic kick distributions are consistent with predictions from binary population synthesis for progenitor systems with a post-supernova orbital period of $P_{\text{orb}}\leq10\,$d and a companion mass of $M_{c}\leq1\,M_{\odot}$, where the natal kicks are calibrated to the velocities of young isolated pulsars. We conclude that the difference in observed systemic kicks is not necessarily in tension with a common origin for MSPs, spider pulsars, and LMXBs.

We present a comprehensive analysis of the dark matter (DM) content and its structural dependence in star-forming disk-like galaxies at intermediate redshifts ($0.1 \lesssim z \lesssim 0.85$), utilizing spatially resolved kinematic data from the MAGPI survey. We report the following: (1) Low stellar mass galaxies ($M_{\rm star} < 10^{9.5}\, M_\odot$) are strongly DM dominated across all radii, with average $\langle f_{_{\rm DM}} \rangle \sim 0.85$, while high-mass ($M_{\rm star} > 10^{10.5}\, M_\odot$) systems exhibit relatively low DM fractions in their inner regions ($\langle f_{_{\rm DM}} \rangle \sim 0.47$) which is equivalent to local massive disk galaxies (e.g., Milky Way and Andromeda). This suggests a mass-dependent structural dichotomy, most-likely governed by a combination of internal galactic processes and environmental influences. (2) A tight inverse correlation between $f_{_{\rm DM}}$ and baryon mass surface density ($\Sigma_{\rm bar}$), with intrinsic scatter of $\sim 0.11$ dex. This is consistent with an inside-out baryon assembly scenario and suggests that the fundamental structural correlations of galaxies were already established by $z\sim 0.85$. (3) No significant evolution in $f_{_{\rm DM}}$ with redshift across the MAGPI window, and when combined with higher-redshift ($0.6 \leq z \leq 1.5$) data from Sharma et al. 2025, we quantitatively show that the reported decline in $f_{_{\rm DM}}(z)$ is most-likely due to observational biases against low-mass systems at $z > 1$. These results offer empirical evidence for a scenario in which disk-like galaxies evolve through a co-regulated build-up of baryonic and DM components, preserving internal structural regularities (such as the total mass distribution and rotation-curve shape) throughout cosmic time.

We investigate the circumstellar environment of the O-type bloated star candidate IRAS 19520+2759 (I19520) using high-resolution observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Expanded Very Large Array (EVLA). Radio continuum emission traced by the EVLA (C, K, and Q bands) exhibits a spectral index of 0.5, consistent with a thermal jet. ALMA 1.3 mm continuum map reveals a compact source coincident with the optical counterpart of I19520, likely tracing the dense core hosting the central massive young stellar object. A prominent molecular outflow in the east-west direction, along with a possible secondary outflow oriented northeast-southwest, is identified in the $^{13}\mathrm{CO}$ emission. A hot molecular core and a Keplerian disk are detected in several $\mathrm{SO}_2$ transitions. Assuming an edge-on disk geometry, the dynamical mass of the central object is estimated to be in the range of $10$-$15~M_\odot$.

The radius valley,a deficit of planets near 2 $\mathrm{R_{\oplus}}$, was observed among exoplanets of radius $\lesssim$ 5 $\mathrm{R_{\oplus}}$ with periods $<$ 100 days by NASA's $Kepler$ mission. It separates super-Earths (rocky, $\lesssim 1.9$ $\mathrm{R_{\oplus}}$) from sub-Neptunes (volatile-rich, $\gtrsim 2$ $\mathrm{R_{\oplus}}$) and may arise from formation conditions or atmospheric loss. Disentangling these mechanisms has led to numerous studies of population-level trends, although the resulting interpretations remain sensitive to sample selection and the robustness of host-star parameters. We re-examine its existence, depth, and dependence on period, flux, stellar mass, and age. Using SWEET-Cat and MAISTEP tool, we derived stellar parameters for 1,221 main-sequence stars (1,405 planets), with effective temperatures 4400--7500 K and radii 0.62--2.75 $\mathrm{R_{\odot}}$, achieving 2\% precision in radius and mass. Planetary radii were recomputed from radius ratios, yielding 5\% median uncertainty. The valley is partially filled near 2 $\mathrm{R_{\oplus}}$ and depends on period, flux, and stellar mass, with slopes $-0.12^{+0.02}_{-0.01}$, $0.10^{+0.02}_{-0.03}$, and $0.19^{+0.09}_{-0.07}$. Sub-Neptunes show a stronger stellar mass-dependent trend than super-Earths ($0.17^{+0.04}_{-0.04}$ vs $0.11^{+0.05}_{-0.05}$). With stellar age, the super-Earth/sub-Neptune ratio rises from $0.51^{+0.11}_{-0.08}$ ($<3$ Gyr) to $0.64^{+0.11}_{-0.11}$ ($\gtrsim3$ Gyr), and the valley becomes shallower and shifts to larger radii. A 4D fit shows consistent slopes with 2D analyses and a weaker age trend ($0.07^{+0.03}_{-0.04}$). These results suggest prolonged atmospheric loss, which is consistent with a core-powered mass loss scenario and emphasize the need for improved determinations, a goal expected to be achieved by future missions like PLATO.

Recent work has demonstrated that widely used strong-line oxygen abundance indicators, such as O3N2, $\rm R23$, and $\widehat{\rm R}$, suffer from large uncertainties when applied to high-redshift galaxies. We show that this loss of precision primarily arises because, at fixed \Oabund, galaxies span a wide dynamic range in ionization parameter and nitrogen enrichment. Here we develop updated indicators that explicitly incorporate both effects via the proxies O32 and N2O2. We define ${\rm R}_{\rm u}\equiv \rm R23+\alpha_1 O32+\alpha_2 N2O2$, $\widehat{\rm R}_{\rm u}\equiv \rm \widehat{R}+\beta_1 O32+\beta_2 N2O2$, and ${\rm O}_{\rm u}\equiv \rm O3N2+\gamma_1 O32+\gamma_2 N2O2$, and calibrate \Oabund~as low-order polynomials in each composite indicator. Applied to a JWST sample with $T_{\rm e}$-method abundances, the updated indicators substantially tighten the correlations with \Oabund, boosting adjusted coefficients of determination from $\mathbb{R}^2\lesssim 0$ (classical indicators) to $\mathbb{R}^2\gtrsim 0.5$ for the full sample and to $\sim 0.7$ at $z>2$. The residuals reveal a redshift evolution in the mapping between \Oabund, strong lines, ionization, and nitrogen enrichment, with a pivotal turning point near the cosmic noon ($z\sim 2$). Our calibrations provide a practical, physically grounded path to precise metallicity measurements in the JWST era and a firmer basis for quantifying early chemical enrichment and feedback.

We study high-energy (HE) neutrino production from interactions between supernova (SN) ejecta and the surrounding circumstellar material (CSM), focusing on regular Type~II and Type~IIn SNe. Using observationally inferred CSM density distributions, we calculate the resulting neutrino fluxes and examine their dependence on key parameters, including the CSM density normalization $D_*$, outer radius $R_{\rm csm}$, proton acceleration efficiency $\epsilon_p$, and magnetic energy fraction $\epsilon_B$. Detection prospects are assessed with a binned likelihood analysis for IceCube, indicating that nearby SNe with moderately dense, confined CSM can produce detectable signals, with a typical detection horizon of $\sim 0.1$ - 1 Mpc. For a Galactic SN at $\sim 10$ kpc, high-statistics neutrino data with detailed temporal and spectral information can constrain $D_*$, $R_{\rm csm}$, and $\epsilon_p$ to within a factor of $\sim 10$ or to a precision of $\sim 20\%$, depending on the assumed values of $D_*$ and $R_{\rm csm}$. These neutrino signals thus provide a complementary probe of the CSM profile and shock acceleration, alongside traditional electromagnetic observations.

We present a method to test the impact of ad-hoc modifications of some of the generic parameters of hadronic interactions -- cross section, elasticity, and multiplicity -- on any observable quantity using full 3-dimensional simulations of extensive air showers induced by ultra-high-energy cosmic rays. Our approach not only extends the existing 1-dimensional tools to three dimensions, but also introduces more flexible features to better respond to the needs of both theory and experiment. We first thoroughly validate the \conexD framework for the simulation of both longitudinal and lateral features of air showers, in particular for a non-standard configuration of the framework in which different energy thresholds for modifications are applied. Moreover, we show that the implementations of the ad-hoc modifications in this configuration are consistent with the previous one-dimensional simulations. Lastly, we discuss the importance of studying the interaction modifications in three dimensions and the effects of parallel modifications of multiple parameters.

The Flat Spectrum Radio Quasar OP 313 is a high-redshift (z = 0.997) blazar that entered an intense gamma-ray active phase from November 2023 to March 2024, as observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. We present a multiwavelength analysis covering 15 years of data, from August 2008 to March 2024, to contextualize this period of extreme gamma-ray activity within the long-term emission of the source. We analyzed a long-term, comprehensive, multiwavelength dataset from different facilities and projects from radio to gamma-rays. We identified the 7 most intense gamma-ray flaring periods and performed a kinematic analysis of Very Long Baseline Array (VLBA) data to determine whether new jet components emerged before or during these flares. For 2 of these flaring periods, we performed the modeling of the spectral energy distribution (SED). The VLBA-BU-BLAZAR and MOJAVE datasets reveal a new jet component appearing in both visibility datasets prior to the onset of one of the strongest gamma-ray flares. By comparing the timing of the VLBA-BU-BLAZAR knots ejection with the gamma-ray flaring periods, we constrained the setup of the SED modeling. We also found that the first gamma-ray flaring period is less Compton-dominated than the others. Our results suggest that the recent activity of OP 313 is triggered by new jet components emerging from the core and interacting with a standing shock. The {\gamma}-ray emission likely arises from dusty torus photons upscattered via Inverse Compton (IC) by relativistic jet electrons. The SED modeling indicates that this component is less dominant during the first {\gamma}-ray flaring period than the later ones.

Intervals of enhanced turbulent fluctuations are typically less frequent within the magnetic cloud region of an interplanetary coronal mass ejection (ICME). We investigate two such intervals inside an ICME observed by the \textit{Wind} spacecraft on 8--9 June 2000 and characterize their associated wave populations. We focus on spectral analysis and plasma instability analysis, using ion-scale normalized magnetic helicity and polarization properties with respect to the background magnetic field $B_0$. In the first interval, the ion-scale normalized magnetic helicity shows a left-handed circularly polarized signature. In the second interval, the left-handed signature persists and an additional high-frequency right-handed population appears. The propagation is approximately parallel to $B_0$. The left-handed fluctuations are compatible with Alfvén ion-cyclotron (AIC) waves, while the right-handed fluctuations are consistent with fast magnetosonic/whistler (FM/W) waves. The ICME plasma accesses resonance conditions that support multiple ion-scale wave modes. Evolving anisotropies in the plasma and the approach to marginal stability allow the coexistence of AIC-like and fast-magnetosonic/whistler-like fluctuations, with enhanced electron heating favoring the growth of the FM/W contribution and strengthening the density--magnetic-field magnitude correlation.

The tidal disruption of planets by their host stars represents a growing area of interest in transient astronomy, offering insights into the final stages of planetary system evolution. We model the hydrodynamic evolution and predict the multi-wavelength observational signatures of planetary TDEs around a solar-mass host, focusing on Jupiter-like and Neptune-like progenitors and examining how different eccentricities of the planet's pre-disruption orbit shape the morphology and emission of the tidal this http URL perform 2D hydrodynamic simulations using the FARGO3D code to model the formation and viscous evolution of the resulting debris disk. We employ a viscous alpha-disk prescription and include a time-dependent energy equation to compute the disk's effective temperature and subsequently derive the bolometric and multi-band photometric light this http URL simulations show that planetary TDEs produce a diverse range of luminous transients. A Jupiter-like planet disrupted from a circular orbit at the Roche limit generates a transient peaking at $L_{bol} \approx 10^{38}$ erg s$^{-1}$ after a 12-day rise. In contrast, the same planet on an eccentric orbit (e=0.5) produces a transient of comparable peak luminosity but on a much shorter timescale, peaking in only 1 day and followed by a highly volatile light curve. We find that the effect of eccentricity is not universal, as it accelerates the event for Jupiter but delays it for Neptune. A robust "bluer-when-brighter" colour evolution is a common feature as the disk cools over its multi-year lifetime. The strong dependence of light curve morphology on the initial orbit and progenitor mass makes these events powerful diagnostics. This framework is crucial for identifying planetary TDEs in time-domain surveys.

Reactions between atomic nuclei are measured in great detail in terrestrial laboratory experiments; transferring and extrapolating this knowledge to how the same reactions act within cosmic environments presents major challenges. Cross-disciplinary efforts are needed in view of the many nuclear reactions that govern the chemical evolution of the universe, and occur in a broad range of stellar plasma conditions that require astrophysical exploration. Since the early identification of 'processes' of nucleosynthesis, new insights have been obtained on the complexity of nuclear reaction mechanisms. We use 12C induced capture and fusion processes to illustrate the challenge of low-energy measurements and of using theoretical methods to extrapolate measurements towards energy regimes within cosmic sources. Particle beam experiments at accelerator facilities above and deep underground simulate stellar reactions, new experimental facilities and methods complement these, and this is further complemented by improved theoretical tools to calculate the quantum effects of nuclear reactions at the various cosmic conditions. Astronomical signatures of cosmic nuclear reactions are deduced from light curves characterizing cosmic explosions through gamma-ray lines and presolar grains to the detection of rare neutrino particles from our Sun to distant cosmic events. High resolution spectroscopy of stars has been expanded to objects measured in the X-ray and the gamma energy range of the electromagnetic spectrum. Astro-seismology and isotopic analysis of meteoritic inclusions provide new tools. Chemical-evolution models describe the complex dynamics during the evolution of galaxies. This article summarizes the experimental and theoretical work, and the broad range of observational tools that test the experimental data and the theoretical interpretation of nuclear processes in the cosmos.

Dark galaxies - systems rich in neutral hydrogen (HI) gas but with no stars - are a common prediction of numerous theoretical models and cosmological simulations. However, the unequivocal identification of such sources in current HI surveys has proven challenging. In this work, we present interferometric follow-up observations with the VLA of a former dark galaxy candidate J0139+4328, originally detected with the single-dish FAST telescope. The improved spatial resolution of the VLA data allow us to identify a faint optical counterpart and characterize the galaxy. Located at a distance of about 31 Mpc, J0139+4328 has a stellar mass of 3 x 10^6 M_Sun and a relatively high gas richness of M_HI/M_star = 18. Despite its high ratio, the galaxy is consistent, within the scatter, with the stellar-to-HI mass relation of HI-selected samples in the literature and with the baryonic Tully-Fisher relation (BTFR), although its kinematic measurement is subject to large uncertainties. This case highlights the potential of modern high-sensitivity HI surveys for detecting low surface brightness, gas-rich galaxies, but underscores the need for careful interpretation of low-resolution HI data, with potentially large centroid errors, and for sufficiently deep optical imaging to ensure robust identification.

Despite the importance of magnetic fields in massive stars, their origin is widely debated and still not well understood. With the mounting evidence for the importance of studying magnetic fields in interacting massive binary and multiple systems, it appears necessary to investigate the presence of magnetic fields in semi-detached systems with ongoing mass transfer, and in contact systems where mass is actively being exchanged. We present an analysis of 53 high-resolution HARPSpol spectropolarimetric observations of a sample of 14 massive binary and multiple systems using the least-squares deconvolution technique. The majority of the studied systems are classified as semi-detached or contact binaries. Definite detections of the presence of a magnetic field are achieved in all studied systems apart from the rather faint system SV Cen, for which only a marginal detection was obtained. The fact that the presence of magnetic fields is detected in all but one of the studied systems strongly suggests that interaction between the system components plays a definite role in the generation of magnetic fields in massive stars. The measured mean longitudinal magnetic field strength for all targets is of the order of a few hundred Gauss to a few kiloGauss. The strongest longitudinal magnetic fields of 4 to 5kG are discovered in the massive O-type triple system MY Ser in both components of the contact binary. kiloGauss-order magnetic fields are also detected in two other systems, V1294 Sco and V606 Cen. It is possible that there is an implication of some system characteristics, such as multiplicity, the mass ratio between the components, and a large fillout factor, on the measured magnetic field strength. Our results for the magnetic field measurements in interacting binaries present the first assessment of the occurrence rate of magnetic fields in a representative sample of such systems.

We present a new first-principles analytic approach to interpreting eclipses and phase curves of rocky planets. Observations with JWST have reported nondetections of atmospheres around the majority of hot rocky planets orbiting M dwarfs. However, most of these "bare rock" inferences are based on models that are ill-suited to many currently observable planets, as they were developed for use on cooler, slower-rotating bodies. In particular, these models rely on the weak temperature gradient assumption, in which rotation is neglected and temperature gradients can be simply related to wind speeds. We find that this assumption may not be valid for over 40% of terrestrials observed with JWST, including TRAPPIST-1b, GJ 367b, and TOI-2445b. Our simple new four-box model does not rely on this assumption, and instead allows the heat transport efficiency to be specified or follow scalings derived herein. This method is fast, interpretable, physically motivated, reproduces previous general circulation model results, and can be used as a starting point for more detailed modeling. We observe that the longitudinal temperature structure of tidally locked terrestrials depends strongly on the atmospheric circulation. Considering the applicable range of atmospheric dynamical regimes, we find that a given planet's nightside temperature can plausibly vary by 100s of Kelvin (from detectable to undetectable). Furthermore, a planet's dayside energy balance can display complex behavior, with degeneracies between surface pressure and dayside temperature. Illustrating an application to observations, we find that assumptions about atmospheric dynamics and longitudinal temperature structure can bias atmospheric constraints at the order-of-magnitude level.

The James Webb Space Telescope (JWST) offers exceptional spectral resolution and wavelength coverage, which are essential for studying the coldest brown dwarfs, particularly Y dwarfs. These objects are at the cold end of the sub-stellar sequence and exhibit atmospheric phenomena such as cloud formation, chemical disequilibrium, and radiative-convective coupling. We examine a curated sample of 22 late-T to Y dwarfs through Bayesian atmospheric retrieval (nested sampling) and supervised machine learning (random forests). Bayesian model comparison indicates that cloud-free models are generally favored for the hottest objects in the sample (T6-T8). Conversely, later-type dwarfs exhibit varying preferences, with both gray-cloud and cloud-free models providing comparable fits. The atmospheric parameters retrieved are consistent across the applied methodologies. Evidence of vertical mixing and disequilibrium chemistry is found in several objects; notably, the Y1 dwarf WISEPAJ1541-22 favors a gray cloud model and shows elevated abundances of both CO and CO2 compared to equilibrium chemistry calculations. As anticipated, the abundances of H2O, CH4, and NH3 increase with decreasing effective temperature over the T-Y sequence.

Knowing the masses and sizes of protoplanetary disks is of fundamental importance for the contemporary theories of planet formation. However, their measurements are associated with large uncertainties. In this proof of concept study, we focus on the very early stages of disk evolution, concurrent with the formation of the protostellar seed, because it is then that the initial conditions for subsequent planet formation are likely established. Using three-dimensional hydrodynamic simulations of a protoplanetary disk followed by radiation transfer postprocessing, we constructed synthetic disk images at millimeter wavelengths. We then calculated the synthetic disk radii and masses using an algorithm that is often applied to observations of protoplanetary disks with ALMA, and compared the resulting values with the actual disk mass and size derived directly from hydrodynamic modeling. We paid specific attention to the effects of dust growth on the discrepancy between synthetic and intrinsic disk masses and radii. We find that the dust mass is likely underestimated in Band 6 by factors of 1.4-4.2 when Ossenkopf & Henning opacities and typical dust temperatures are used, but the discrepancy reduces in Band~3, where the dust mass can be even overestimated. Dust growth affects both disk mass and size estimates via the dust-size-dependent opacity, and extremely low values of dust temperature (~ several Kelvin) are required to recover the intrinsic dust mass when dust has grown to mm-sized grains and its opacity has increased. Dust mass estimates are weakly sensitive to the distance to the source, while disk radii may be seriously affected. We conclude that the accuracy of measuring the dust mass and disk radius during the formation of a protoplanetary disk also depends on the progress in dust growth. (Abridged)

Debris discs are analogues to our own Kuiper belt around main-sequence stars and are therefore referred to as exoKuiper belts. They have been resolved at high angular resolution at wavelengths spanning the optical to the submillimetre-millimetre regime. Short wavelengths probe the light scattered by such discs, which is dominated by micron-sized dust particles, while millimetre wavelengths probe the thermal emission of millimetre-sized particles. Determining differences in the dust distribution between millimetre- and micron-sized dust is fundamental to revealing the dynamical processes affecting the dust in debris discs. We aim to compare the scattered light from the discs of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) with the thermal emission probed by ALMA. We focus on the radial distribution of the dust. We used high-contrast scattered light observations obtained with VLT/SPHERE, GPI, and the HST to uniformly study the dust distribution in those systems and compare it to the dust distribution extracted from the ALMA observations carried out in the course of the ARKS project. We also set constraints on the presence of planets by using these high-contrast images combined with exoplanet evolutionary models. 15 of the 24 discs comprising the ARKS sample are detected in scattered light, with TYC9340-437-1 being imaged for the first time at near-infrared wavelengths. For 6 of those 15 discs, the dust surface density seen in scattered light peaks farther out compared to that observed with ALMA. These 6 discs except one are known to also host cold CO gas. Conversely, the systems without significant offsets are not known to host gas, except one. This observational study suggests that the presence of gas in debris discs may affect the small and large grains differently, pushing the small dust to greater distances where the gas is less abundant.

We present the first systematic spectroscopic observations of extragalactic transients from the Dark Energy Spectroscopic Instrument (DESI), as part of the DESI Transients Survey program. With 5,000 fibers and an ${\sim} 8$ deg$^2$ field of view, we exploit DESI as a machine for the discovery and classification of transients. We present transient classifications from archival DESI data in Data Releases 1 and 2, relying on a combination of a secondary target program and serendipitous observations. We also present observations from the first 6 months of the DESI spare fiber program dedicated to transients. The program is run in coordination with a dedicated DECam time-domain survey, serving as a pathfinder for what we will be able to achieve in conjunction with the Rubin Observatory Legacy Survey of Space and Time (LSST). We classify over 250 transients, of which the majority were previously unclassified. The sample comprises thermonuclear and core-collapse supernovae and tidal disruption events (TDEs), including a TDE observed before its discovery in imaging. We demonstrate DESI's ability to classify a population of faint transients down to $r\sim 22.5$ mag during main survey operations, with negligible impacts on DESI's main observations. can you make this plain text for arxiv abstarct

Binaries contain rich physical information, and the study of binaries has always been a hot topic in stellar physics research. The stars LAMOST J1132 and LAMOST J0529 have not yet been recorded in the SIMBAD astronomical database. We have investigated their physical properties via methods such as spectral analysis, photometric analysis, and light curve analysis. Based on comprehensive analysis, we conclude that they are two newly discovered binary systems, each consisting of a hot white dwarf and a flaring companion star. Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectra indicate that both stars contain hot white dwarfs. The spectral fitting yields $T_{eff}$=53728$\pm$2467\,K, log$g$=7.98$\pm$0.08 for LAMOST J1132, and $T_{eff}$=47381$\pm$494\,K, log$g$=7.84$\pm$0.05 for LAMOST J0529. The weak neutral metal lines in the LAMOST spectra and the discrepancy between the Global Astrometric Interferometer for Astrophysics (GAIA) and LAMOST spectra both indicate that these two sources are likely binary systems. The relatively high flux values for both sources in the near-infrared and mid-infrared bands support our preliminary judgment. The color index in the near-infrared bands suggests that the companion star is K or M type for LAMOST J1132 and M type for LAMOST J0529. Light curve data from the Zwicky Transient Facility (ZTF) indicate that the companion stars of both sources are stars exhibiting flare activity. The eclipse probability is very low, indicating that these two sources are non-eclipsing binary systems. The physics of binaries is fascinating, and future data from LAMOST Medium Resolution Spectra are expected to enable the detection of magnetic fields in these two hot white dwarfs.

PyIRD is a Python-based pipeline for reducing spectroscopic data obtained with IRD (InfraRed Doppler; Kotani et al. (2018)) and REACH (Rigorous Exoplanetary Atmosphere Characterization with High dispersion coronagraphy; Kotani et al. (2020)) on the Subaru Telescope. It is designed to process raw images into one-dimensional spectra in a semi-automatic manner. Unlike traditional methods, it does not rely on IRAF (Tody, 1986; 1993), a software used for astronomical data reduction. This approach simplifies the workflow while maintaining efficiency and accuracy. Additionally, the pipeline includes an updated method for removing readout noise patterns from raw images, enabling efficient extraction of spectra even for faint targets such as brown dwarfs. The code is open source and available at this https URL .

Statistics on tidal disruption events (TDEs) may be contaminated by repeating TDEs (rTDEs), which have been extensively discovered recently. However, the origin of rTDEs remains unclear. In addition, no statistical research on rTDEs with time intervals $>5$ years has been made yet. In this work, we searched for rTDEs with time intervals of 5--19 years using CRTS data in a sample of 16 ZTF BTS TDEs at $z<0.05$. We found 2 rTDE candidates, AT 2019azh and AT 2024pvu, with time intervals of 13.2 and 17.1 years, respectively. The peak luminosities of CRTS flares are close to those of ZTF flares. For the CRTS flare of AT 2024pvu, using GALEX UV observations near the peak, we measured a blackbody temperature of $\sim19500$ K, consistent with TDEs and higher than SNe. Moreover, we estimated the expected number of SNe in the sample to be $\lesssim0.08$, and hence the probability that both CRTS flares are SNe is only 0.3\%. Therefore, the possibility that both CRTS flares are SNe can be ruled out, and it is likely that both are TDEs. Using the two rTDEs, we inferred that the TDE rate is 2--3 orders of magnitude higher than the average over 5--19 years prior to TDE detection. Considering another two rTDEs with intervals of $\sim$2 years in the sample and possible rTDEs missed by CRTS, rTDEs with intervals of $<20$ years may account for 25\%--60\% of the TDE sample. We prefer to explain rTDEs as repeating partial TDEs, but the possibility of independent TDEs cannot be ruled out and requires future observational tests.

Solar white-light flares (WLFs) are solar flares exhibiting enhanced emission in the optical continuum. They are critical for understanding energy release and transport mechanisms in solar flares and for conducting comparative studies with stellar WLFs. However, the scarcity of accurately and reliably measured optical continuum light curves for solar WLFs significantly hampers related studies. Based on the optimized solar WLF identification method, we construct a dataset of optical continuum light curves for 70 solar WLFs using 6173 Å continuum intensity images from the Solar Dynamics Observatory. Moreover, for each solar WLF event, we also provide the location of the white-light emission enhancement signals and key parameters including bolometric energies and durations derived from both the traditional fixed-temperature blackbody model and the refined variable-temperature blackbody model. This dataset will serve as a valuable resource for future statistical investigations of solar WLFs and for comparative studies between solar and stellar flares.

We report the discovery of multiple compact molecular features exhibiting extremely broad velocity widths toward the W44 molecular cloud. ALMA CO $J$=3--2 data reveal eight ``Petit--Bullets'' surrounding the previously known ``Bullet.'' Each Petit--Bullet shows a distinct V-shaped structure in position--velocity space, reminiscent of the Y-shaped morphology of the Bullet, suggesting a common origin. These features are interpreted as the result of high-velocity plunges of compact gravitational objects into dense molecular gas. The spatial and kinematic properties of the Petit--Bullets suggest that the plunging material was not a single object but rather a small cluster of compact bodies. A virial mass of $1.0\!\times\! 10^{5}\, M_\odot$ inferred from their velocity dispersion is comparable to that of typical globular clusters. Momentum analysis further implies that the main Bullet likely formed by an isolated black hole. These findings provide new evidence for dynamical interactions between halo clusters and disk molecular gas.

The systematic analysis of non-adiabatic effect on convective mode has been conducted using wave energy relation. In the adiabatic analysis, the "propagation diagram" for convective mode is proposed as a useful tool to see its behavior. In the non-adiabatic analysis, it is found that for strongly non-adiabatic case, a monotonically growing convective mode becomes oscillatory. In this phase, the radial displacement and the distribution of wave energy show only one bump, in which the distribution of entropy energy eS almost overlaps with the distribution of gravity energy eg. Entropy energy eS seems to act as potential energy of oscillatory convection. In addition to this, this change occurs not gradually, but abruptly with change of non-adiabatic indicator.

This work develops low-energy spacecraft (SC) trajectories using Venus gravity assists to study asteroids during heliocentric transfer segments between planetary encounters. The study focuses on potentially hazardous asteroids (PHAs) as primary exploration targets. This paper proposes a method for calculating SC trajectories that enable asteroid flybys after a Venus gravity assist. The method involves formulating and solving an optimization problem to design trajectories incorporating flybys of selected asteroids and Venus. Trajectories are calculated using two-body dynamics by solving the Lambert problem. A preliminary search for candidate asteroids uses an algorithm to narrow the search space of the optimization problem. This algorithm uses the V-infinity globe technique to connect planetary gravity assists with resonant orbits. The resonant orbit in this case serves as an initial approximation for the SC's trajectory between two successive planetary flybys. Four flight schemes were analyzed, including multiple flybys of Venus and asteroids, with the possibility of an SC returning to Earth. The proposed solutions reduce flight time between asteroid approaches, increase gravity assist frequency, and enhance mission design flexibility. The use of Venus gravity assists and resonant orbits ensures a close encounter with at least one asteroid during the SC's trajectory between two consecutive flybys of Venus, and demonstrates the feasibility of periodic Venus gravity assists and encounters with PHAs. The developed method was applied to construct trajectories that allow an SC to approach both Venus-resonant asteroids and PHAs via multiple Venus gravity assists. An additional study was carried out to identify asteroids accessible during the Earth-Venus segment in launch windows between 2029 and 2050.

Binaries in which both stars are pulsating are rare but extremely valuable. We present the first study of an asteroseismic binary system consisting of a core helium-burning red clump (RC) star and a red giant branch (RGB) star. The Kepler target KIC 10841730 is a wide binary (period $2917 \pm 8$ d) that provides ideal conditions to test the accuracy of RC models. While prior studies of RC stars have revealed discrepancies in modelling the period spacings of mixed modes, other model parameters remain largely untested. We perform a detailed modelling analysis using individual mode frequencies and cover a large parameter space in mass, metallicity, He-abundance, mixing length, overshooting, and mass-loss, and we also explore different methods to correct for surface effects. We find two possible results for the red clump models. One solution requires introducing an unexpected offset of the phase shift in the red clump model, yielding an age consistent with the companion star and current masses of $1.01 \pm 0.06$ and $1.08 \pm 0.06$ M$_\odot$ for the RC and RGB star, respectively. Alternatively, we find that excluding the identification of two questionable radial modes resolves the phase-shift offset issue but results in a higher mass and thus a much younger age for the red clump star, contradicting the age obtained from its companion. We conclude that uncertainties in red clump models affect not only the g-mode period spacings but also the properties of the p modes. We show the power of asteroseismic binaries in validating and constraining stellar models and highlight the need for refining red-clump models.

Classical Be (CBe) stars are rapidly rotating B-type stars with Balmer emission lines that originated from the decretion disks surrounding them in their spectra. Accounting for $\sim$20% of all B-type stars, most CBe stars are thought to form through mass and angular momentum transfer from their companions. It follows that in most close CBe star binaries, the companions are expected to be post-main-sequence stars rather than main-sequence (MS) stars. Hitherto, $\sim$100 CBe star binaries have been identified, the majority of which are Be/X-ray binaries. As expected, none of the others have indeed been confirmed as CBe+MS binary stars. To further study and verify the origin of CBe stars, identifying additional CBe star binaries is indispensable. In this study, we report 504 CBe stars identified using data from Data Release 12 of the Large sky Area Multi-Object fiber Spectroscopic Telescope. Among these, 141 are newly identified and 14 exhibiting radial velocity variations are identified as CBe star binaries. Besides, 60 CBe stars with high normalized unit weight error (RUWE) but not confirmed by dynamics are proposed as potential CBe star binaries. We also find that 34 CBe stars are potential cluster members. By calculating peculiar velocities, 37 runaway stars are identified with peculiar velocities ranging from $\sim$40 km s$^{-1}$ to $\sim$101 km s$^{-1}$.

We present a method for identifying false positive double-lined spectroscopic binary (SB2) candidates by analysing medium-resolution survey spectra from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR10. Specifically, we focus on spectra contaminated by moonlight, which exhibit near-zero radial velocity (RV) and solar-like spectral lines from the secondary component. By applying strict constraints on the contamination parameters and fitting the contaminated spectra, we ultimately confirmed that 126 false positive binaries are single stars contaminated by moonlight. Additionally, we identify several key factors contributing to moonlight contamination: the lunar phase during observation, the G-band magnitude of the star, and the angular distance between the star and the moon. Notably, artificial satellites in low-orbital can also introduce contamination from solar-like spectral components, but they typically display significantly higher transverse velocities. In a follow-up study, we will expand our analysis to identify additional false positive SB2 systems and systematically classify them according to their contamination sources.

Context. Understanding the survival and orbital evolution of brown dwarf (BD) companions during the post-main-sequence (MS) evolution of their host stars is increasingly important, especially with recent discoveries of many substellar companions around white dwarfs (WDs). Aims. We investigate the long-term evolution and final outcomes of BDs orbiting low-mass MS stars as these evolve into WDs. By comparing forward-modeling populations with observed WD-BD binaries, we test evolutionary models and predict the existence of yet-undetected systems. Methods. We employ the COMPAS binary population synthesis code to evolve observed MS-BD systems through the post-MS phases of their host stars into the WD stage, tracking orbital changes driven by mass loss, tides, and common-envelope (CE) evolution. Results. Our simulations reproduce a period gap in the distribution of detached WD-BD binaries, consistent with observations. We also identify a boundary separating detached and semi-detached systems on the period-mass diagram, located at orbital periods of $\sim$1-2 hours depending on the BD mass. Conclusions. We predict that a subset of currently known MS-BD binaries will survive post-MS evolution and emerge as detached WD-BD systems, while others will undergo CE evolution and potentially form cataclysmic variables with BD donors. Our results reproduce the observed period gap in WD-BD binaries and provide quantitative predictions for the role of CE efficiency in shaping their distribution. This work predicts that many WD-BD systems remain undetected, motivating targeted searches with microlensing and high-contrast imaging techniques using next-generation large telescopes.

The long-term dynamical evolution of asteroid families is governed by the interplay between orbital and rotational evolution driven by thermal forces and collision. We aim to observationally trace the rotational evolution of main-belt asteroid families over Gyr timescales. We analyzed rotational properties of 8739 asteroids with spin period measurements and 3794 asteroids with obliquity determinations across 28 asteroid families spanning ages from 14~Myrs to 3~Gyrs. We introduced a dimensionless timescale that normalizes each asteroid's family age by its classical YORP timescale, enabling direct comparison of rotational states across different evolutionary stages. We examined two key observables: the fraction of slow rotators (periods greater than or equal to 30 hours) and the polarization fraction (the degree to which asteroid spin poles align correctly with their position in the family's V-shape distribution according to the Yarkovsky theory). Evolution of both quantities were fitted to identify characteristic transition timescales. We discovered that the slow-rotator fraction increases steeply with $t$ and saturates at $f_{\rm slow} \simeq 0.25$ around a breakpoint $t_{\rm bp} \simeq 20$. This implies a stochastic YORP timescale $\tau_{\rm YORP,stoc} \simeq 10\,\tau_{\rm YORP}$ by comparison with rotational evolution models that include tumbling and weakened YORP torques. The polarization fraction reaches a maximum of $\simeq 0.8$ at $t \simeq 16$ and then decays toward the random limit $f_{\rm pol} \rightarrow 0.5$ for $t \gtrsim 20$, indicating an increasing dominance of collisional spin reorientation over time. The rotation properties within different asteroid families offer crucial clues to rotation evolution and can serve as a new dimension for age estimation of asteroid families with more data in the LSST era.

We present a joint analysis of galaxy clustering and galaxy--galaxy lensing measurements from BOSS galaxies using a simulation-based emulation method combined with a halo occupation distribution model. Our emulators are constructed with the Aemulus $\nu$ simulations, a suite of $w\nu$CDM $N$-body simulations with massive neutrinos as independent particle species. We combine small-scale analysis of clustering from $0.1h^{-1}$Mpc to $60.2~h^{-1}$Mpc and lensing from $1.7h^{-1}$Mpc to $60.2~h^{-1}$Mpc to perform cosmological constraints. We split the BOSS galaxies into three redshift bins to measure their clustering and employ galaxies from Dark Energy Camera Legacy Survey and Hyper Suprime-Cam as source galaxies to measure lensing separately. We find that the addition of lensing significantly improves the constraining power on $S_{8}=\sigma_8(\Omega_m/0.3)^{0.5}$, with a weak improvement for $f\sigma_{8}$. Our results of $f\sigma_{8}$ indicate tensions of around $1\sim4\sigma$ below the results of CMB observations of Planck. For $S_{8}$, our results are also lower than Planck, and the tension can be mitigated when considering possible systematics in lensing measurement. As a byproduct, our analysis prefers a non-zero neutrino mass but without strong significance, with the constraining power dominated by the clustering. Given the accuracy and precision of our model and the observational data, it is anticipated that larger and higher-quality spectroscopic datasets will improve the constraints on this fundamental property in the near future.

Pulsar timing arrays (PTAs) can detect the low-frequency stochastic gravitational-wave background (GWB) generated by an ensemble of supermassive black hole binaries (BHBs). Accurate determination of BHB merger timescales is essential for interpreting GWBs and constraining key astrophysical quantities such as black hole (BH) occupation fractions and galaxy coalescence rates. High-accuracy $N$-body codes such as \texttt{Griffin} can resolve sub-pc BHB dynamics but are too costly to explore a wide range of initial conditions, motivating the need for surrogate models that emulate their long-term evolution at much lower computational cost. We investigate neural ordinary differential equations (NODEs) as surrogates for the secular orbital evolution of BHBs. Our primary contribution is a parameterised NODE (PNODE) trained on an ensemble of $N$-body simulations of galaxy mergers spanning a two-dimensional parameter space defined by the initial orbital eccentricity and particle resolution $(e_i, N)$, with the learned vector field explicitly conditioned on these parameters. A single PNODE thereby learns a simulation-parameter-conditioned dynamical model for the coupled evolution of the BH pair's orbital state across the ensemble, yielding smooth trajectories from which stable hardening and eccentricity growth rates can be extracted. The PNODE accurately reproduces the secular evolution of the specific orbital energy and angular momentum, and the corresponding Keplerian orbital elements, for held-out trajectories, with modest generalisation to a partially unseen high-resolution case. Combining PNODE predictions with semi-analytical prescriptions for stellar hardening and gravitational-wave emission yields BHB merger timescales consistent with those obtained from direct $N$-body inputs within current theoretical uncertainties.

We present a study of the transition of Neutral Hydrogen (HI) gas from the end of the Epoch of Reionization (EoR) to late-time large-scale structure. We examine the signature of the transition as traced through the redshifted 21-cm line with SKA-Low at $3 < z < 7$. To do so, we use the semi-numerical simulation \textsc{21cmFAST} to model the HI during the EoR and add a HI-halo based post-processing model of the late-time HI. This approach gives a robust estimate of the amplitude of the HI temperature field and predicts the observable power spectrum during the transition period. We find that our simulation pipeline reproduces the expected power spectrum trends from existing observations and theory, in addition to replicating current observational constraints on $\Omega_{\text{HI}}$. Our simulations predict a drop in power of four orders of magnitude between $4 < z < 7$. Assuming an inhomogeneous recombination model, we find a flattening of the power due to lingering neutral islands masking the late-time HI signal for $5 < z< 6.5$. Using SKA-Low deep survey parameters, we find HI power spectrum detectability at scales $k \leq 1$ $h$ Mpc$^{-1}$ for redshifts $3< z < 7$, even when using the horizon limit to mitigate foregrounds. Our results suggest a sufficient SNR of the HI power spectrum tracing the underlying halos $z < 5$, which can be used for late-time cosmology. Our results suggest that the resulting $\Omega_{\rm HI}$ constraints can trace different reionization scenarios such as a decreased escape fraction. This study implies that deep SKA-Low observations for $3< z< 7$ will be an important probe to constrain reionization parameters as well as cosmological models.

Context. In the Local Group, dwarf spheroidal galaxies (dSphs) and ultra-faint dwarf galaxies (UFDs) exhibit large velocity dispersions. These values are generally attributed to the presence of substantial amounts of dark matter (DM), in line with the predictions of the standard model of galaxy formation. However, alternative, more conservative explanations exist, such as non-virialized dynamical states induced by tidal interactions, the presence of stellar streams, and artificial inflation of the velocity dispersion caused by binary-star orbital motion. Aims. We study the dynamical evolution of UFDs using purely stellar ("dry") dynamics, without invoking DM. We dynamically evolve our systems up to a Hubble time and compare our results with observational studies and previous theoretical work. Methods. We employ direct high precision NBODY simulations performed with the NBODY6++GPU code. We explore the role of binaries in inflating the velocity dispersion of low-mass host galaxies. We also present both the stellar and dynamical evolution of the stellar population, which is necessary to properly interpret our results. Results. We find that, in all our models, the UFD remains globally quasi-stationary for approximately 3000 Myr. Subsequently, the system undergoes mass segregation and experiences a phase resembling core collapse. Red giants and white dwarfs (WD) are found to play significant, but distinct, roles. Red giants provide the dominant contribution to the luminosity, whereas WDs constitute the largest fraction of the non-luminous component, accounting for approximately 13% of the total stellar population. Finally, if not taken into account properly, velocity dispersion measurements can be strongly biased by the presence of a significant binary population, which can lead to substantial overestimates of velocity dispersion in UFDs

Inflationary cosmology successfully accounts for the observed properties of primordial fluctuations using quantum field theory in an expanding background. However, the quantum nature of these fluctuations has not been experimentally established, since classical stochastic models could reproduce the observed two-point statistics by construction. Existing approaches to testing primordial quantumness focus primarily on Bell inequalities, which provide a sharp conceptual criterion but are difficult to implement with cosmological observables. In this work we adopt a falsification-based approach. We define a precise classical hypothesis for the origin of primordial perturbations (local stochastic fields admitting a positive probability distribution) and identify inequality constraints that must be satisfied within this class. We show how violations of these classicality inequalities can be probed using realistic cosmological observables, without invoking Bell tests or non-commuting measurement settings. We further identify symmetry-protected spectator sectors in which quantum coherence is parametrically preserved during inflation, allowing violations of observable magnitude to survive decoherence. Our results show that large-scale structure and future 21 cm surveys provide a viable and quantitative route to falsifying classical stochastic descriptions of primordial fluctuations.

This study extends the 3D magnetohydrodynamic (MHD) simulation of a jet emerging from a binary neutron star (BNS) merger presented in Pavan et al. (2023), in which an incipient jet was manually injected into the realistic environment imported from a previous general-relativistic MHD simulation of a merging BNS system. The jet evolution is followed up to almost 10 seconds without loss of resolution. Our results reveal that the jet faces challenges in penetrating the dense surroundings, leading to a barely successful outflow that exhibits structural asymmetries and low Lorentz factors. By the end of the extended simulation, 98% of the jet energy is converted to kinetic form and its angular structure is stabilized. The physical quantities inferred thus provide reliable inputs for afterglow emission calculations. This work demonstrates a method for simulating jets in 3D up to nearly ballistic regimes that is general and ready to be applied to any jet in a BNS merger context.

We present the JWST NIRSpec/G395H transmission spectrum of the misaligned hot Jupiter HAT-P-30b from 2.8--5.2 $\mu$m as part of the BOWIE-ALIGN survey, a comparative survey designed to probe the link between planet formation and atmospheric composition in samples of misaligned and aligned hot Jupiters orbiting F-type stars. Through independent data reductions and retrieval analyses, we find evidence for absorption features of H$_2$O and CO$_2$ in the atmosphere of HAT-P-30b. Our retrieved abundances are consistent with equilibrium chemistry, from which we infer a sub-solar C/O ratio (0.16--0.45), and sub-solar and sub-stellar metallicity (0.2--0.8$\times$solar, compared to a stellar metallicity of 1.1--1.6$\times$solar), with muted spectral features. This composition challenges formation models of continuous migration and accretion within a steady disc of stellar metallicity, and could be the result of low C/O ratio gas accretion within the water ice line, low metallicity accretion due to the trapping of volatiles further out in the disc, or the combined accretion of low metallicity gas and carbon-poor solids.

When an object is scattered near a supermassive black hole (SMBH), tidal oscillations excited within it reduce its orbital energy, leading to capture by the SMBH. This process, called tidal capture, can also occur when the object approaches even closer to the SMBH, resulting in a partial tidal disruption event (pTDE). Previous studies on pTDEs of main-sequence stars have shown that as the disruption intensifies, dynamical effects dominate over tidal oscillations, causing the remnant material to acquire a kick velocity instead of being captured by the SMBH. In this work, we performed hydrodynamic numerical simulations of pTDEs involving giant stars. We found that for weaker disruptions, the dynamical behavior of the remnant material resembles that of main-sequence stars. However, as the disruptions deepen, the remnant material transitions from gaining energy to losing energy, leading to capture by the SMBH. This behavior markedly differs from that of main-sequence stars, demonstrating that the presence of a compact core significantly influences the dynamical processes in pTDEs. Our simulations reveal that the energy change of the remnant material strongly correlates with asymmetric mass -- lossspecifically, the difference in mass outflow between the Lagrange points L1 and L2. This suggests that the energy change stems from asymmetric mass loss, consistent with conclusions from previous studies on main-sequence stars. However, quantitative analysis contradicts earlier models, indicating that the dynamical model of pTDEs requires further refinement. Finally, we discuss the characteristics of repeating pTDEs produced by this process and their potential observability, as well as the implications for the long-term orbital evolution of high eccentricity extreme mass ratio inspiral systems.

We investigate the impact of flickering variability in jet power on the luminosity and morphology of radio galaxies. We use a Lagrangian particle method together with relativistic hydrodynamics simulations using the PLUTO code to track the evolution of electron spectra through particle acceleration at shocks and cooling processes. We introduce an adapted version of this method which improves tracking of adiabatic cooling in regimes where low density jet material mixes with high density from the ambient medium in the lobes. We find that rapid increases in jet power can lead to large increases in hotspot luminosity due to the interaction of a travelling shock structure with the pre-existing shock structure at the jet head. We show that in some cases it may be possible to identify a bright region of emission corresponding to a shock travelling along the jet axis. We find that the time-averaged radiative efficiency of variable jets is similar to their steady counterparts, but find significant departures from this on an instantaneous basis. We suggest that, together with environmental effects and differences in the average powers of jets, variable jet powers may have a significant impact on how we understand the diversity of radio jets seen in observations and have significant implications for interpretations of jet powers, energy budgets and luminosity-linear size diagrams.

Carbon monoxide is thought to be the second most abundant molecule in the Universe. This makes observation of both its parent isotopologue ($^{12}$C$^{16}$O) and its stable isotopologues, $^{13}$C$^{16}$O, $^{12}$C$^{18}$O, $^{12}$C$^{17}$O, $^{13}$C$^{18}$O and $^{13}$C$^{17}$O, important in variety of objects. Here the MARVEL (Measured Active Rotational-Vibrational Energy Levels) algorithm is used to determine precise rotational vibrational energy levels for the five minor isotopologues of carbon monoxide in their electronic ground state. A review of 27 literature sources yields 3716, 1454, 89, 728 and 57 validated transitions for $^{13}$C$^{16}$O, $^{12}$C$^{18}$O, $^{12}$C$^{17}$O, $^{13}$C$^{18}$O and $^{13}$C$^{17}$O, respectively, giving 863, 499, 33, 345 and 45 empirically determined, rotation vibration energy levels, respectively.

Using a recently suggested magneto-morphological classification (MMC, Abramenko, 2021, MNRAS Vol 507) of solar active regions (ARs), we explored 3048 ARs, observed from12 May 1996 to 27 December 2021. Magnetograms were acquired with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) and with the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). ARs were separated between three classes: class A - regular ARs (bipoles which follow the empirical rules compatible with the mean field dynamo theory); class B - irregular ARs (``wrong'' bipoles and multipolars); class U - unipolar sunspots. An aim of the present study is to explore time variations of a typical unsigned magnetic flux of ARs of different classes. The typical flux was acquired as the mean flux over all ARs of a given class observed during one solar rotation. The time profiles of the mean fluxes for different classes were compared. We found that, except for periods of deep solar minima, the mean flux of B-class ARs always dominate that of A-class ARs, and, what is the most important, the time profile of B-class ARs is highly intermittent versus the rather smooth and quazi-constant A-class profile. Intermittency implies a direct involvement of turbulence. We conclude that, through the entire active phase, the Sun is capable of producing regular moderate ARs at a quazi-constant rate along with the production of large and complex irregular ARs in the very intermittent manner. The result is the first observational evidence for the long-standing speculative assumption on the involvement of the convection zone turbulence into the regular global dynamo-process on a stage of the active regions formation.

The bulk compositions of small planets ($R_p< 2 \mathrm{R}_\oplus$) are directly linked to their formation histories, making reliable compositional constraints imperative for testing models of planet formation and evolution. Because exoplanet interiors cannot be directly observed, their make-up must be inferred from mass-radius-composition models that link assumed stellar abundances to the direct observables: planetary mass and radius. There are a variety of such models in the literature, each adopting different equations of state (EOS) to describe the materials' properties at depth and varying assumptions about the minerals present within the planets. These EOS+mineral suites provide the foundations for compositional inferences, but they have not yet been systematically compared. In this work, we review several suites, with a detailed description of the basic structure, mineral physics, and materials within standard small planet models. We show that EOS+mineral suites predict planet densities whose differences are comparable to current observational uncertainties, which present a challenge for robustly interpreting and classifying small planets. We apply a powerful small-planet characterization framework, which illustrates that variations among EOS+mineral suites lead to inconsistent conclusions for both individual planets and sample-level demographics. Our results demonstrate the need for more careful considerations of the materials and EOS used in mass-radius-composition models, especially given the current focus on finding and characterizing potentially habitable rocky planets. We conclude with recommendations for best practices so that future interpretations of small planets and their formation are accurate and consistent.

We present a comprehensive temporal and spectral study of the flat-spectrum radio quasar PKS~0402$-$362 using \textit{Fermi}-LAT/Swift-XRT/UVOT observations spanning from MJD 54686-60321. The $\gamma$-ray light curve exhibits multiple phases of enhanced activity, with the fractional variability parameter ($F_{\mathrm{var}}$) showing larger amplitudes at longer timescales, consistent with variability trends observed in other FSRQs. Statistical analysis of the flux and spectral index distributions using the Anderson--Darling test and histogram fitting reveals that both distributions deviate from a single log-normal form and are better represented by a double log-normal profile, indicating two distinct flux states. A search for quasi-periodic oscillations in the $\gamma$-ray emission using the Lomb--Scargle periodogram identified a significant periodic signal at $\sim$413~days with a confidence level exceeding $3\sigma$. However the proximity of the timescale to one year and limited number of observed cycles prevents a definitive interpretation. Broadband spectral energy distributions for six flux states were modeled using a one-zone leptonic framework incorporating synchrotron, synchrotron self-Compton (SSC), and external Compton (EC) components. The SEDs are well reproduced with physically reasonable parameters: high-flux states exhibit harder electron spectra and lower magnetic field strengths ($B \sim 0.2--0.6\,\mathrm{G}$), while low-flux states show softer spectra and stronger magnetic fields ($B \sim 1.3\,\mathrm{G}$). The fitted break energy decreases during high-flux states, suggesting enhanced radiative cooling and a transition toward a particle- or kinetic-energy-dominated jet. These trends are consistent with the ``harder-when-brighter'' behavior commonly observed in blazars.

The evolution of protoplanetary discs is a function of their internal processes and of their environment. It is unclear if angular momentum is mainly removed viscously or by magnetic winds, or by a combination of the two. While external photoevaporation is expected to influence disc evolution and dispersal, there are observational limitations towards highly irradiated discs. The interplay between these ingredients and their effect on the gas and dust distributions are poorly understood. We investigate the evolution of both the gaseous and solid components of viscous, MHD-wind or hybrid discs, in combination with external FUV-driven mass loss. We test which combinations of parameters protect discs from external irradiation, allowing the solid component to live long enough to allow planet formation to succeed. We run a suite of 1D simulations of smooth discs with varying initial sizes, levels of viscous and MHD-wind stresses modeled via an $\alpha$ parametrisation, and strengths of the external FUV environment. We track disc radii, various lifetime diagnostics, and the amount of dust removed by the photoevaporative wind, as a function of the underlying parameters. The biggest role in determining the fate of discs is played by a combination of its ability to spread radially outwards and the strength of FUV-driven erosion. While MHD wind-driven discs experience less FUV erosion due to the lack of spread, they do not live for longer compared to viscously evolving discs, especially at low-to-moderate FUV fluxes, while higher fluxes yield disc lifetimes that are insensitive to the disc's angular momentum transport mechanism. For the solid component, the biggest role is played by a combination of inward drift and removal by FUV winds. This points to the importance of other physical ingredients, such as disc substructures, even in highly-irradiated disc regions, in order to retain solids.

The DragonCam Microscopic Camera is an instrument being developed for NASA's Dragonfly mission [1] to Saturn's moon Titan. The Microscopic Camera will be body-fixed to the Dragonfly vehicle and will image the surface at a distance of about one meter (98.6 cm nominal) with a pixel scale of better than 60 microns/pixel and a nominal 52 degree angle to the Titan surface. With the 4.8 um pixel pitch of the sensor, this is a focal length of about 77.5 mm. To accommodate range variations due to vehicle pose and surface topography, the Microscopic Camera has a focus mechanism to give it a depth of field (DOF) of about 130mm. Since the Microscopic Camera's boresight is tilted by 52° off the vertical, the optical configuration has a compensating tilted focal plane, taking advantage of the Scheimpflug imaging principle. The optics are all-refractive with nine elements, a six-element stationary group and a three-element moving group. A plano-plano window seals the optics from the environment and also serves as the substrate for a bandpass filter. The optomechanical system is derived from the Mars Hand Lens Imager [11]; the moving group is mounted to a linear slide which is translated via a cam follower by the rotation of a cam driven by a stepper motor. The Microscopic Camera is designed to survive at temperatures as low as -130C without power. The camera is enclosed in a cavity in the foam insulation covering the spacecraft and looking through a single-pane window. Prior to imaging, the camera will be heated to operating temperature (nominal -30C) for proper actuation of the mechanism. STOP analysis has been performed to demonstrate that optical performance is maintained after heating. Software focus merging will be performed in the onboard camera control electronics to minimize image data downlink requirements.

Covariance matrices are essential cosmological probes of fundamental physics, providing information on numerous fundamental physical parameters and varying with any change in the underlying cosmology. However, this cosmology dependence, while providing excellent information, also makes them computationally intensive to compute, as a new covariance matrix must explicitly be calculated for every variation in cosmology before comparisons to observational data can be made. In this paper, we develop an efficient model for estimating the parameter dependence of the covariance matrix of the matter power spectrum by Taylor expanding around a known value of the parameter space. This method allows us to use a relatively small number of input cosmologies, specifically one fiducial cosmology and two further cosmologies for each parameter. We explicitly calculate the covariance matrices for these cosmologies and then develop a new model that allows us to interpolate from these the form of the covariance matrix with a cosmology that is located elsewhere in that given parameter space without explicit perturbation theory calculations. This method speeds up covariance matrix calculations in new cosmologies by orders of magnitude compared to explicit perturbation theory calculations at each point in a given parameter space. Using different approximations, we develop three versions of our interpolated covariance matrix and validate the model by recreating all of our input cosmologies using all three forms, both with and without super-sample covariance corrections in each case, and show that the models provide robust recreations of the original results, with the different approximations being valid in certain regimes.

The microlensing event OGLE-2011-BLG-0950 exhibits the well-known ``Planet/Binary'' degeneracy, in which distinct lens configurations produce similar light curves but imply substantially different mass ratios between the lens components. A previous study suggested that high-resolution imaging could break this degeneracy through differences in the lens-source relative proper motion. In this work, we identify a new planetary model for this event that arises from a newly identified degeneracy, simultaneously reproducing the observed light curve and remaining consistent with the relative proper motion measured from high-resolution imaging. By combining constraints from the light-curve modeling and high-resolution observations, we infer a lens system consisting of a $\sim 1~M_{\odot}$ host star orbited by a $\sim 1.5~M_{\rm Jup}$ planet, with a projected separation of about 2 or 8 au, subject to the ``Close/Wide'' degeneracy. Our reanalysis of the color-magnitude diagram further indicates that the source star has unresolved companions that contribute non-negligible blended light, highlighting the importance of carefully accounting for source and lens companions in future Roman microlensing analyses. Finally, we show that adopting a single mass--luminosity relation significantly underestimates the uncertainties in the inferred lens properties for host masses $\gtrsim 1~M_{\odot}$.

We quantify the utility of HCN and HNC to characterize gas conditions in the nearby starburst galaxy NGC 253. We use measurements from the Atacama Large Millimeter/Submillimeter Array (ALMA) Large Program ALCHEMI: the ALMA Comprehensive High-resolution Molecular Inventory. Using different subsets of the eight total HCN and HNC transitions measured by ALCHEMI, we test the number and combinations of transitions necessary for constraining the temperature, H$_2$ volume and column densities, cosmic-ray ionization rate, and beam-filling factor in three representative regions within NGC 253. We use these combinations of HCN and HNC transitions to constrain chemical and radiative transfer models and infer the gas conditions using a Bayesian nested sampling algorithm combined with neural network models for increased efficiency. By comparing the shapes of the resulting posterior distributions, as well as the medians and uncertainties for each gas parameter, from each test case to what we obtain with the full set of eight transitions (the control), we quantify how well each test reproduces the control. We find that multiple transitions each of both molecules are required to obtain a median parameter value within a factor of 2 of the control with an uncertainty less than 2-3 times that of the control. We also find that transition combinations that feature a range of upper-state energies are most effective. We show that single transitions, such as HCN J = 1-0 or 3-2, are among the worst-performing combinations and result in parameter values up to an order of magnitude different than the control.

The matter distribution in the Universe exhibits a rich variety of structures forming the cosmic web. These structures arise from the anisotropic gravitational collapse of primordial density fluctuations and define the pathways along which galaxies flow from voids to high-density clusters. Local density variations within these structures play a fundamental role in driving the environmental evolution of galaxies. To characterise filament boundaries, we analysed galaxy overdensity profiles around filaments in two redshift ranges: $0.05 < z < 0.1$ and $0.1 < z < 0.3$. Perpendicular and parallel profiles were derived by averaging galaxy overdensity as a function of distance. Characteristic scales and central overdensities were then analysed by fitting analytical models, specifically exponential and power-law families. We also introduced normalised density profiles to account for survey incompleteness. The perpendicular overdensity profiles show a nearly constant value in the central regions $D_{fila} < 1$ Mpc, decreasing at distances up to $\approx 10$ Mpc. The mean physical widths (scale radii) at $0.05 < z < 0.1$ and $0.1 < z < 0.3$ are $2.39 \pm 0.69$ and $5.56 \pm 2.29$ Mpc, respectively. This scale difference between redshift ranges is also evident in the normalised profiles. Conversely, profiles along filaments remain constant at distances larger than $\approx 20$ Mpc from the nearest intersection. Our results show that the influence zone of cosmic filaments extends up to $\sim 10$ Mpc from their spines. Furthermore, a mild evolution in structural parameters is observed over the past $\sim 4$ Gyr, suggesting that filaments undergo measurable changes even at relatively low redshifts.

Between the groundbreaking detections of stellar-mass black holes by LIGO/Virgo/KAGRA and JWST's revelation of a surprisingly abundant population of supermassive black holes, one crucial missing link remains: the elusive intermediate-mass black holes (IMBHs). IMBHs represent a key phase in the hierarchical growth of black holes, yet they have persistently evaded detection. Traditional methods, effective for both actively accreting and quiescent black holes, have largely failed to uncover this hidden population. Here, we argue that novel observational strategies--particularly time-domain variability studies of active galactic nuclei (AGN) and tidal disruption events--provide a promising path forward. Finding IMBHs will resolve critical gaps in our understanding of black hole formation and the various mechanisms driving their subsequent growth. The upcoming Vera C. Rubin Observatory, with its unprecedented capacity to monitor the dynamic sky, stands to revolutionize our ability to detect these long-sought IMBHs, shedding new light on the assembly history of black holes across cosmic time.

Our understanding of the $\gamma$-ray sky has greatly advanced, yet studying the unresolved $\gamma$-ray background (UGRB) can unveil the nature of the faintest $\gamma$-ray source populations in the Universe. Statistical cross-correlations between the UGRB and tracers of large-scale cosmic structure allow us to infer which sources contribute the most to this emission. In this work, we examine the angular correlation between the UGRB and the matter distribution traced by galaxies, using twelve years of Fermi Large Area Telescope (LAT) observations along with three years of Dark Energy Survey (DES) data. We detect a correlation with a signal-to-noise ratio of 7.96, primarily driven by large angular scales. We then perform a multi-tracer analysis that combines this measurement with the cross-correlation between $\gamma$ rays and DES weak lensing. The two single-tracer results are mutually consistent, and their combination yields a total significance of 8.6, firmly establishing the extragalactic origin of the UGRB. Intriguingly, the properties inferred for the sources contributing to the UGRB show departures from those of the resolved {\gamma}-ray population, suggesting that the faint end of the $\gamma$-ray sky is not a simple extrapolation of currently resolved sources.

Radio halos are megaparsec-scale diffuse, non-thermal radio sources located at the centers of galaxy clusters, tracing relativistic particles and magnetic fields in the intra-cluster medium. Their origin is generally attributed to cluster mergers that generate turbulence and re-accelerate aged electrons. We study the diffuse radio emission, spectral properties, and the connection between thermal and non-thermal components in the massive galaxy clusters Abell 773 and Abell 1351 ($M_{500} \sim 7 \times 10^{14}\,M_{\odot}$), both of which are dynamically disturbed. We combine LOFAR LoTSS-DR2 observations at 144 MHz with uGMRT observations at 650 MHz, supplemented by archival XMM-Newton X-ray imaging. We confirm that both clusters host radio halos extending up to a largest linear size of $\sim 2$ Mpc. We measure an integrated spectral index $\alpha_{144}^{650} \sim -1.0$ for both clusters. The radio halo in Abell 773 resembles a classical halo and follows a sublinear radio--X-ray surface brightness relation. In contrast, Abell 1351 shows a more complex and asymmetric morphology, influenced by embedded radio sources including the brightest cluster galaxy, a tailed radio galaxy, and a ridge-like feature. These contaminating sources lead to deviations from the sublinear trend in the point-to-point radio--X-ray analysis of Abell 1351.

We present the Meerkat Asteroid Guard, an imminent impactor warning service developed and maintained by the European Space Agency's Near-Earth Object Coordination Centre (NEOCC). The software uses the method of systematic ranging to perform orbit determination on tracklets in the Near-Earth Object Confirmation Page (NEOCP), which typically have short observational arcs. Fitted orbits are propagated to determine the likelihood of an impact with Earth. In addition, magnitude fitting and Monte Carlo sampling are performed to estimate the object's size, possible impact locations and times, and suggest a best telescope pointing for object follow-up. A set of object scores are produced from computed posterior probabilities across the grid, giving a statistical description of the object's orbital and physical characteristics. The scores are packaged with several informative plots in an email alert, which is sent to Meerkat subscribers in the event of a significant impact probability, close approach, or other scientifically interesting event. The highlights of the five years of Meerkat's operational service are presented, including the successful warnings for all of the past six imminent impactors discovered before impact and several interesting close approaches.

We analyse a sample of 225 star-forming regions from the DESIRED-E project, each with simultaneous determinations of the electron temperature from ionized nitrogen and oxygen, $T_{\rm e}$([NII]) and $T_{\rm e}$([OIII]), respectively. We derive new empirical relations connecting the gas-phase metallicity to the global electron temperature, $T_{\rm e}$(H$^+$), as determined via radio observations. We establish two calibrations: one assuming a homogeneous temperature distribution ($t^2 = 0$, the ``direct method''), and another accounting for internal temperature fluctuations ($t^2 > 0$). Applying these calibrations to 460 radio observations of Galactic HII~regions spanning Galactocentric distances from $\sim0.1$ to 16 kpc, we determine the radial O/H gradient in the Milky Way under both assumptions. We further compare these nebular gradients to independent metallicity estimates from young O- and B-type stars and Cepheid variables. We find that the $t^2 > 0$ calibration yields a gradient in excellent agreement with stellar-based determinations, whereas the $t^2 = 0$ method underestimates metallicities by up to $\sim$0.3 dex. This discrepancy cannot be reconciled by invoking oxygen depletion onto dust grains or nucleosynthetic processing via the CNO cycle in massive stars. We also find that one widely used relation in the literature, assuming $t^2 = 0$, produces an excessively steep gradient -- likely due to the use of outdated atomic data and pre-CCD observations. Finally, we explore potential azimuthal variations in the Galactic metallicity distribution driven by the presence of the spiral arms, finding no evidence for variations larger than $\sim$0.1 dex with respect to the general radial gradient.

We present the longest and the densest quasi-simultaneous radio, X-ray and optical campaign of the black hole low mass X-ray binary GX339-4, covering five years of weekly GX339-4 monitoring with MeerKAT, Swift-XRT and MeerLICHT, respectively. Complementary high frequency radio data with the Australia Telescope Compact Array are presented to track in more detail the evolution of GX339-4 and its transient ejecta. During the five years, GX339-4 has been through two "hard-only" outbursts and two "full" outbursts, allowing us to densely sample the rise, quenching and re-activation of the compact jets. Strong radio flares were also observed close to the transition between the hard and the soft states. Following the radio flare, a transient optically thin ejection was spatially resolved during the 2020 outburst, and was observed for a month. We also discuss the radio/X-ray correlation of GX339-4 during this five year period, which covers several states in detail from the rising phase to the quiescent state. This campaign allowed us to follow ejection events and provide information on the jet proper motion and its intrinsic velocity. With this work we publicly release the weekly MeerKAT L-band radio maps from data taken between September 2018 and October 2023.

We investigate the environmental parameters and gas fraction (f$_{gas}$) properties of type~2 AGN and non-AGN galaxies, utilizing a large sample of galaxies from SDSS DR7 with z $\le$ 0.3. We find that the environment affects type~2 AGN and non-AGN galaxies in similar ways and does not impact the strength of AGN-driven outflows. The f$_{gas}$ of type~2 AGN and non-AGN host galaxies show no variation between group and isolated environments, suggesting that host galaxy gas content is largely independent of large-scale environment. We find that type~2 AGN host galaxies possess systematically lower f$_{gas}$ than their non-AGN counterparts when matched in stellar mass and star formation rate (SFR). This suggests that AGN activity plays a significant role in regulating the molecular gas reservoir and, consequently, the star formation processes within galaxies. We find that Type~2 AGNs exhibiting strong outflows are associated with higher gas fractions, higher star-formation rates, and younger stellar populations than those with weak or no outflows. This may indicate either concurrent star formation in gas-rich systems hosting powerful outflows, or a time delay between AGN activity and its effect on star formation consistent with a delayed AGN feedback scenario.

Ultraluminous X-ray sources (ULXs) are point-like sources that exhibit apparent X-ray luminosities exceeding the Eddington limit for stellar-mass compact objects. A widely accepted interpretation is that these systems are X-ray binaries accreting matter possibly at super-Eddington rates. In this regime, photon trapping inflates the accretion disk, making it geometrically and optically thick. Radiation-driven winds launched from the supercritical disk form funnel-shaped walls along the symmetry axis. While the apparent X-ray luminosity can exceed the Eddington limit due to geometrical beaming within this funnel, a misalignment with the observer's line of sight strongly suppresses the X-ray emission, rendering the ULX electromagnetically obscured. This work explores the potential for high-energy neutrino production in black hole-hosting ULXs. We model proton acceleration via magnetic reconnection in the region above the super-accreting black hole. Although electromagnetic emission is efficiently absorbed by the dense wind and radiation fields, neutrinos generated from photomeson interactions can escape. Our model self-consistently accounts for energy losses of pions and muons in this environment. The results indicate that misaligned, electromagnetically obscured Galactic ULXs could produce a neutrino flux detectable by instruments like KM3NeT and IceCube within several years of observation.

We present a study of the gas kinematics within the Hestia project, a state-of-the-art set of simulations of the Local Group, with a particular focus on the velocity patterns of different ions and the large-scale motion of gas and galaxies towards the Local Group barycentre. Using two high-resolution Hestia runs, we examine the distribution and velocities of H I, C IV, Si III, O VI, O VII, and O VIII and their imprints on sightlines observed from the Sun's location in different reference frames. To mimic observational strategies, we assess the contribution of rotating disc gas, assuming simple kinematic and geometrical considerations. Our results indicate that local absorption features in observed sightlines most likely trace material in the circumgalactic medium of the Milky Way. Some sightlines, however, show that intragroup material could be more easily observed towards the barycentre, which defines a preferred direction in the sky. In particular, H I, Si III, and C IV roughly trace cold gas inside the Milky Way and Andromeda haloes, as most of their mass flux occurs within the virial region of each galaxy, while oxygen high ions mostly trace hot halo and intragroup gas, with comparable mass fluxes in the Local Group outskirts and the circumgalactic medium of the two main galaxies. Additionally, we find that pressures traced by different ionic species outside the Milky Way halo show systematically higher values towards the barycentre direction in contrast to its antipode in the sky. Kinematic imprints of the global motion towards the barycentre can be seen at larger distances for all ionic species as the Milky Way rams into material in the direction of Andromeda, with gas towards the anti-barycentre lagging behind.

Identifying the mechanisms responsible for angular momentum transport in protoplanetary disks, and the extent to which those mechanisms produce turbulence, is a crucial problem in understanding planet formation. The bulk of the gas in protoplanetary disks is weakly ionized, which leads to the emergence of three non-ideal effects, Ohmic diffusion, ambipolar diffusion, and the Hall effect. These low-ionization processes can in some cases suppress turbulence driven by the magnetorotational instability (MRI). However, it has recently been shown that these non-ideal terms can also affect the dynamics of the gas in fundamentally different ways than simple diffusion. In order to further study the role of low-ionization on disk gas dynamics, we carry out a 3D local shearing box simulation with both Ohmic diffusion and ambipolar diffusion and an additional simulation with the Hall effect included. The strength of each non-ideal term, when present, is representative of gas at a radius of 5 AU in a realistic protoplanetary disk. We find the Hall effect increases the saturation strength of the magnetic field, but does not necessarily drive turbulence, consistent with previous work. However, interactions between ambipolar diffusion and the Keplerian shear lead to the ambipolar diffusion shear instability (ADSI), which can drive the initial growth, not damping, of magnetic perturbations. To our knowledge, this is the first work that explicitly demonstrates the viability of the ADSI in the non-linear regime within protoplanetary disks. At later times in the disk, the MRI (reduced in strength by ambipolar-diffusion), may also be present in regions of weak magnetic field between strong concentrations of vertical magnetic flux and sustain turbulence locally in protoplanetary disks.

We report the detection of the polycyclic aromatic hydrocarbon (PAH) 1H-cyclopent[cd]indene (c-C11H8) in TMC-1 with the QUI- JOTE line survey. We detected 22 independent lines corresponding to 88 rotational transitions with quantum numbers ranging from J=19 up to J=24 and Ka <= 5 in the Q-band range. The identification of this new PAH was based on the agreement between the rotational parameters derived from the analysis of the lines and those obtained by quantum chemical calculations. The column density derived for 1H-cyclopent[cd]indene is (6.0 +- 0.5) x 10^12 cm-2, with a rotational temperature of 9 K. Its abundance is high, as is that of the rest of the PAHs, but it is the lowest of all those detected to date in TMC-1, being 2.66 times less abundant than indene and 4.66 times less than phenalene. This result will help us to better understand the growth of five- and six-membered rings in dark clouds. Chemical models explaining their formation through the bottom-up model are still very incomplete and require further experimental and theoretical effort. Even so, the most likely formation reactions would occur between the smallest rings with small hydrocarbons; the most probable reaction for the formation of cyclopentindene is that between indene and C2H, C2H3, and/or their cation.

We study the local evolution of the Universe around galaxy formation sites in the EAGLE50 large-volume reference simulation. Using the reduced inertia tensor (r-TOI), we followed the anisotropic evolution of initially spherical Lagrangian volumes (LVs) centred at galaxy formation sites, both in dark matter (DM) and in cold baryons (CB), from very high redshift $z=15$ onward. We describe LV deformation in terms of the r-TOI eigen-directions, principal axes, their derived shape parameters, and the timescales for the freezing-out of these principal directions and axes. Of particular interest are the age of the Universe, $t_{\rm U}$, when the local Cosmic Web (CW) spine emerges, and that when anisotropic DM mass arrangements (i.e., migrant mass flows) cease. We find that the shapes LVs acquire along their evolution affect the halo and stellar mass of their central galaxy: prolate-shaped LVs show a tendency to host low-mass galaxies at $z=0$, while massive galaxies tend to form within triaxial or oblate LVs. Also, the local CW spine tends to set in earlier on in LVs that are to host massive galaxies than in those harbouring less massive galaxies. In addition, anisotropic DM-mass rearrangements stop late on average, at $t_{\rm U}\sim 10.5\,$Gyr, and even slightly later for CB. Interestingly, $z=0$ LVs with either flattened configurations in CB or those that are highly prolate in DM, are more likely to host rotation-dominated galaxies. This effect increases from $z=1$ to $z=0$. Finally, the CB spine of LVs that are more likely to host rotation-dominated galaxies emerges at later times.

We present a novel tunnel architecture for the Einstein Telescope that departs from the traditional large-cavern approach and reduces the excavated volume by an order of magnitude. In the proposed design, all seismic isolation systems are housed in raise-bore wells drilled upward from the main tunnel toward an upper service tunnel. The pre-isolators for the most sensitive optics are located in the service tunnel, seating directly on strong and compact rock, while the other filters are distributed along the wells within compact, side-access vacuum chambers. Shorter, separate wells accommodate the seismic isolation systems for less demanding optics. This configuration provides substantial advantages: easier lock acquisition and improved robustness of the interferometers, lower-frequency pendulum stages, reduced congestion around the test masses, simplified installation and maintenance, improved vacuum partitioning, strong physical decoupling between the high- and low-frequency interferometers, and enhanced compatibility with future advances of Newtonian-noise cancellation. A novel technique for real-time, precision monitoring of rock motion and tilt provides a new signal for Newtonian noise cancellation and enables correction of seismic disturbances even during earthquakes, offering unique geophysical measurement capabilities.

We discuss the most recent observations made with the Southern Connecticut Stellar Interferometer (SCSI), which is a three-station stellar intensity interferometer located on the campus of Southern Connecticut State University, in New Haven, Connecticut. Two different kinds of observations are presented. We first analyze observations of Vega taken in a three-telescope mode. (Previously, the instrument had only two operational stations.) We show that, while the efficiency remains nearly identical to that reported in our last paper, the addition of the third station allows more photon data to be recorded simultaneously, and therefore we can build up the photon-bunching peak in the data stream in fewer hours on sky for an unresolved source. In the second part of the paper, we report our observations to date of the nearby red giant star, Arcturus, most of which occurred in the first half of 2025. These show that, as a partially resolved source at the baselines we used, we detect fewercorrelations in the photon-bunching peak than for an unresolved source of comparable brightness. Combining the data with speckle imaging observations taken at Apache Point Observatory, we derive a new measurement of Arcturus' diameter that extends the time baseline of interferometric observations of the star and is consistent with previous analyses made by other investigators.

We present the analysis of four microlensing events, KMT-2020-BLG-0202, KMT-2022-BLG-1551, KMT-2023-BLG-0466, and KMT-2025-BLG-0121, which exhibit extended and complex anomalies in their light curves. These events were identified through a systematic reanalysis of KMTNet data aimed at detecting planetary signals that deviate from the typical short-term anomaly morphology. Detailed modeling indicates that all four anomalies were produced by planetary companions to low-mass stellar hosts. The events have mass ratios of $q \sim (5$--$14)\times10^{-3}$ and Einstein timescales of $t_{\rm E} \sim 20$--$43$ days. Bayesian analyses based on Galactic models show that the companions are super-Jupiters with masses of a few to approximately 10 $M_{\rm J}$, orbiting sub-solar-mass hosts located at distances of $D_{\rm L} \sim 4$--$7$~kpc. All planets lie well beyond the snow line of their hosts, placing them in the regime of cold giant planets. These detections demonstrate that extended and complex microlensing anomalies, which are often challenging to recognize as planetary in origin, can nonetheless contain planetary signals. This work underscores the unique sensitivity of microlensing to cold, massive planets beyond the snow line and highlights the importance of systematic reanalyses of survey data for achieving a more complete and unbiased census of exoplanets in the Galaxy.

We present results from a coordinated, multi-epoch near-infrared and centimeter radio survey of young stellar objects (YSOs) in the Coronet, aimed at probing the connection between mass accretion and ionised mass loss. Using VLT-KMOS, we detect Br$\gamma$ emission in 5 of the 26 targets, which also exhibit 3.3-cm continuum emission in VLA images, consistent with partially ionised jets. For seven additional sources, stringent flux upper limits were obtained. The derived accretion and ionised mass-loss rates for class I and class II YSOs follow a sublinear correlation $\dot{M}_{\mathrm{ion}} \propto \dot{M}_{\mathrm{acc}}^{0.3}$, consistent with previous results for class II YSOs but extended here to earlier stages. Multi-epoch observations reveal modest variability in both tracers but no clear temporal correlation between accretion and ejection within timescales of a few months. The ratio $\dot{M}_{\mathrm{ion}}/\dot{M}_{\mathrm{acc}}$ shows an anti-correlation with $\dot{M}_{\mathrm{acc}}$, increasing with time from class I YSOs to class II YSOs, suggesting an increase in jet-launching efficiency or ionisation fraction with evolution. These findings support a direct connection between accretion and outflow across the $\sim$ Myr timescale of YSO evolution, while highlighting the complexity of their short-term interplay.

Magnetic fields and turbulence may play a key role in the evolution of protoclusters, influencing the formation of dense cores and stars. Here, we examine the morphology of the magnetic fields in the G327.29 protocluster using both the velocity gradient technique (VGT) extracted from molecular line emissions and linear polarization in the dust continuum emission. The VGT analysis is performed using four molecular tracers: DCN (3-2), C18O (2-1), HN13C (3-2), and H13CO+ (3-2) - which probe gas across different density regimes, observed with the ALMA 12 m array. Owing to its sensitivity to gas dynamics, a comparison between VGT and dust polarization provides a powerful probe of the evolutionary processes in massive star-forming regions. From our analysis we reveal a complex magnetic-field structure, shaped by the combined influence of turbulence and gravity. In addition, it also appears that there is a large-scale (beyond the core scale) gravitational infall from the surrounding medium on to the filament and the central densest region. Furthermore, we observe that cores are dominated by a mix of turbulence and gravity. Overall, this work presents, likely for the first time, the application of VGT to a massive protocluster, G327.29, using high-resolution ALMA observations.

The Kepler-observed distribution of planet sizes have revealed two distinct patterns: (1) a radius valley separating super-Earths and sub-Neptunes and (2) a preference for intra-system size similarity. We present a new model for the exoplanet population observed by Kepler, which is a "hybrid" of a clustered multi-planet model in which the orbital architectures are set by the angular momentum deficit (AMD) stability (He et al. 2020; arXiv:2007.14473) and a joint mass-radius-period model involving envelope mass-loss driven by photoevaporation (Neil & Rogers 2020; arXiv:1911.03582). We find that the models that produce the deepest radius valleys have a primordial population of planets with initial radii peaking at $\sim 2.1 R_\oplus$, which is subsequently sculpted by photoevaporation into a bimodal distribution of final planet radii. The hybrid model requires strongly clustered initial planet masses in order to match the distributions of the size similarity metrics. Thus, the preference for intra-system radius similarity is well explained by a clustering in the primordial mass distribution. The hybrid model also naturally reproduces the observed radius cliff (steep drop-off beyond $\sim 2.5 R_\oplus$). Our hybrid model is the latest installment of the SysSim forward models, and is the first multi-planet model capable of simultaneously reproducing the observed radius valley and the intra-system size similarity patterns. We compute occurrence rates and fractions of stars with planets for a variety of planet types, and find that the occurrence of Venus and Earth-like planets drops by a factor of $\sim 2$-4 for the hybrid models compared to previous clustered models in which there is no envelope mass-loss.

We investigate the redshift 1.98 galaxy cluster XLSSC 122 using the Hubble Space Telescope (HST) from the core of the cluster out to 3 Mpc, a scale equivalent to 10 times the R500 = 295 kpc radius. We present an expanded photometric and spectroscopic catalogue of the cluster, bringing the total number of spectroscopically classified member galaxies to 74, with 35 new member galaxies added in the outer regions of the cluster. We compute the radial galaxy number density profile in the cluster, and observe no clear evidence of infalling groups or cosmic filaments. We observe a clear bimodal colour relation in member galaxies, with red fraction increasing towards the cluster centre. This rapid increase of red fraction upon infall is indicative of a fast quenching mechanism, such as ram pressure stripping, as galaxies enter the cluster centre. We fit a luminosity function to the cluster members, finding a similar low mass slope but fainter scale magnitude than z = 1 clusters of similar temperature, implying a similar galaxy evolution rate to clusters at lower redshift.

Deuterium fractionation becomes highly efficient in cold, dense cores where CO is frozen out. Cyclopropenylidene (c-C3H2), an early-formed carbon ring, and its deuterated isotopologues trace gas-phase deuteration in these environments. We present a statistical study of c-C3H2 deuteration in starless and prestellar cores of the Perseus Molecular Cloud using observations of c-C3H2, c-C3HD and c-C3D2 obtained with the Yebes 40 m, ARO 12 m and IRAM 30 m telescopes towards 16 cores. Gaussian fits and RADEX modeling yield column densities for the detected species. c-C3H2 is detected in 14/15 covered cores, c-C3HD in 15/16, and c-C3D2 in 9/16. Derived column densities range from 0.5-8.1 x 10^{13} cm^{-2} for c-C3H2, 0.2-2.1 x 10^{12} cm^{-2} for c-C3HD, and 0.6-1.6 x 10^{11} cm^{-2} for c-C3D2. The ortho-to-para ratio of c-C3H2 is obtained for all but one core, with a median value of 3.5\pm0.4. Statistically corrected D/H ratios span 0.5-9.2% (median 1.5\pm0.2%), and D2/D ratios 9-55% (median 25.9\pm4.3%). No trend is found between the c-C3H2 ortho-to-para ratio and core evolutionary stage traced by n(H2). The median D/H ratio in Perseus appears lower than values reported for Taurus and Chamaeleon, while the D2/D ratio agrees with Taurus within uncertainties. A positive correlation between D/H and n(H2) supports the use of D/H as an evolutionary tracer. D2/D does not correlate with n(H2), but shows a positive correlation with T_{kin}, suggesting that its formation is influenced by a mildly endothermic pathway.

One of the primary goals of next-generation gravitational lensing surveys is to measure the large-scale distribution of dark matter, which requires accurate mass inversion to convert weak-lensing shear maps into convergence (kappa) fields. This work develops a mass inversion method tailored for upcoming space missions such as CSST and Euclid, aiming to recover both the mass distribution and the convergence power spectrum with high fidelity. We introduce MIU2Net, a versatile deep-learning framework for kappa-map reconstruction based on the U2-Net architecture. A new loss function is constructed to jointly estimate the convergence field and its frequency-domain energy distribution, effectively balancing optimal mean squared error and optimal power-spectrum recovery. The method incorporates realistic observational effects into shear fields, including shape noise, reduced shear, and complex masks. Under noise levels anticipated for future space-based lensing surveys, MIU2Net recovers the convergence power spectrum with 4% uncertainties up to l approximately 500, significantly outperforming Wiener filtering and MCALens. Beyond two-point statistics, the method accurately reconstructs the convergence distribution, peak centroid, and peak amplitude. Compared to other learning-based approaches such as DeepMass, MIU2Net reduces the root-mean-square error by 5% without smoothing and by 38% with a 1-arcmin smoothing scale. MIU2Net represents a substantial advancement in mass inversion methodology, offering improved accuracy in both RMSE and power-spectrum reconstruction. It provides a promising tool for mapping dark matter environments and large-scale structures in the era of next-generation space lensing surveys.

Correlations between gamma-ray and X-ray spectral properties of pulsars are investigated in order to provide observational hints on physics involved in pulsars' high-energy emissions. Using a sample of 43 pulsars detected in both X-ray and gamma-ray bands, we find that pulsars' gamma-ray luminosity, $L_\gamma$, clearly correlates with the luminosity of non-thermal X-ray emission, $L_{\rm p}$, and anti-correlates with non-thermal X-ray photon index. Other gamma-ray spectral parameters show weaker or negligible correlations. The found relation that $L_\gamma \propto L_{\rm p}^{0.49\pm 0.05}$ implies a certain connection between radiation mechanisms and energy distributions of radiating particles for these high-energy emissions. Pulsars with and without detected thermal emissions seem to show different dependencies in those correlations, suggesting the possible existence of two different kinds of pulsars. The ones without detected thermal emissions may represent a population of pulsars with low surface temperature. The origin and energetics of high-energy emitting electron-positron pairs for this group of pulsars probably do not depend on their surface thermal emissions, while that of the other group do. The low surface temperatures might be evidence for the working of some exotic processes of neutron-star cooling. Similar to $L_{\rm p}$, some tempting relationships are found among each gamma-ray spectral parameter, surface temperature and thermally radiating area radius. It again strengthens the connection between gamma-ray and X-ray emissions from pulsars.

China's Tianwen-2 mission, launched on 29 May 2025, targets the near-Earth object (469219) Kamo`oalewa, an Earth quasi-satellite trapped in a 1:1 mean-motion resonance with our planet. Determining the origin of Kamo`oalewa is central to understanding the formation pathways and dynamical evolution of Earth's quasi-satellite population. Here we show a strong possibility of main-belt origin for Kamo`oalewa using long-term dynamical simulations. We examine three candidate source regions: the $\nu_6$ secular resonance ($\nu_6$), the 3:1 mean-motion resonance with Jupiter (3:1J MMR), and the Flora family. A total of 42,825 test particles were integrated over 100 Myr. We find that asteroids from all three regions can be transported onto Kamo`oalewa-like orbits, albeit with markedly different efficiencies. Particles originating near the $\nu_6$ show the highest transfer probability (3.31%), followed by the Flora family (2.54%) and the 3:1J MMR (0.39%). We further identify representative dynamical pathways linking these source regions to Earth quasi-satellite orbits. The Tianwen-2 spacecraft is expected to rendezvous with Kamo`oalewa in 2026, performing close-proximity operations and returning samples. The mission will provide decisive observational constraints on the asteroid's composition and physical properties, offering a critical test of its proposed origin.

Using new Institut de Radioastronomie Millimétrique (IRAM) 30m telescope $\rm N_2H^+$, $\rm C^{18}O$ $J$=1-0 and Atacama Pathfinder Experiment (APEX) telescope $\rm ^{13}CO$ and $\rm C^{18}O$ $J$=2-1 maps together with archival far-infrared continuum data, and $\rm ^{12}CO$, and $\rm ^{13}CO$ $J$=1-0 data, we present a comprehensive analysis of the massive filament CFG024.00$+$0.48 (G24) across clump-to-cloud scales. Our results show that G24 is an $\sim$80 pc giant filament with a total mass of $\sim$$10^5$ M$_{\odot}$. In the different tracers the filament width is measured to be about $\sim$2 times the beam size of the observations, as expected for power-law density distributions, giving beam-deconvolved widths in the range from 0.8 to 2.8 pc. We determine a line-of-sight thickness of $\sim$2.2 pc demonstrating that G24 is not an edge-on, flatten structure. The virial parameter obtained from line mass ($\alpha_{\rm line,vir}=M_{\rm line,vir}/M_{\rm line}$) from the $\rm C^{18}O$ (1-0) data is 0.85, and that obtained from $Herschel$-based H$_2$ column density is 0.52, suggesting G24 is globally close to virial equilibrium. The distribution of the 40 dust clumps appears to have a ''two-tier'' fragmentation pattern. For the clump groups, the separation, with a mean/median of 3.68/3.46 pc, is very close to expected length associated with the maximum fragmentation growth rate of $\lambda_{\rm max}=3.55 \pm0.32$ pc estimated for the dust. However, the longitudinal centroid velocity profiles of $\rm C^{18}O$ and $\rm N_2H^+$ show oscillation patterns with wavelengths of 9.8$\pm$0.1 pc and 9.9$\pm$0.1 pc, respectively. This is $\sim$2 times larger than the corresponding values of $\lambda_{\rm max}$ of 4.96$\pm$0.63 pc and 4.65$\pm$1.34 pc, respectively. This suggests that the velocity structure is not dominated by flows directly associated with the fragmentation seen in the dust emission.

With the increasing number of high-precision astronomical observations, physical quantities that were previously inaccessible to accurate calculations, such as cosmic birefringence, have once again become a focal point of interest. Such phenomena induce a nonvanishing cross-correlation between the $E$- and $B$-mode polarizations of the cosmic microwave background (CMB), thereby providing a direct observational signature of parity violation. The Chern-Simons coupling between the scalar field in early dark energy (EDE) models and CMB photons is regarded as a plausible mechanism for generating cosmic birefringence. Recent data from the Atacama Cosmology Telescope (ACT) deliver $EB$ measurements at higher multipole moments than those previously achieved by {Planck}, while DESI and PantheonPlus datasets provide new and stringent constraints on the late-time expansion history. Using a joint analysis of {Planck}, DESI DR1, Pantheon+, and ACT data, we perform a full-parameter constraint on the cosmic birefringence effects induced by the EDE-CMB photon coupling. Our results favor a higher Hubble constant, $H_0 = 76.9^{+2.9}_{-2.5}\,\rm km\,s^{-1}\,Mpc^{-1}$, and a relatively large EDE fraction, $f_{\mathrm{EDE}} = 0.232^{+0.074}_{-0.047}$. By comparing the cosmological evolution of this model across different data combinations, we find that the ACT-$EB$ data combined with {Planck} + DESI + PantheonPlus provide good constraints to both early- and late-Universe observations.

We present the first comprehensive photometric and spectroscopic analysis of the RS CVn system HD 26172, robustly determining the previously debated evolutionary state of its primary star. Since this system is a single-lined spectroscopic binary with spot-induced light curve modulations, we derived its physical parameters by combining the TESS light curves, the radial velocity curve from our observations, and the primary-star mass estimates based on three complementary this http URL results reveal that HD 26172 is a detached binary system composed of a $1.25 \pm 0.32 M_{\odot}$ subgiant and a $0.63 \pm 0.11 M_{\odot}$ main-sequence star. The conclusion of subgiant primary is also supported by the absence of lithium absorption and no observed infrared excess. Using long-term photometry from the KWS survey, we detected a tentative stellar activity cycle of 5635 days with an amplitude of 0.04 mag in HD 26172. Additionally, we identified ten optical flare events exhibiting temporally clustered outburst behavior. The presence of a long-term activity cycle, pronounced starspot activity, and frequent optical flares makes HD 26172 a valuable laboratory for studying magnetic activity in subgiants within close binary systems.

We present imaging and spatially resolved spectral analyses of eight Chandra data taken for the Circinus galaxy in $\approx$ 22 years to reveal neutral iron-K$\alpha$ emission on a circumnuclear scale ($\sim$ 10--100 pc) and search for time variability in the emission. By simulating and taking account of point-source emission from the active galactic nucleus (AGN), we detect iron-line emission $\sim$ 20--60 pc away from the nucleus, particularly in the eastern and western regions. In the two regions, possible time variability in the line flux was also detected. Our spectral analysis then finds that the observed equivalent widths can reach $\sim$ 2 keV and the slopes of underlying continua are rather inverted with $\Gamma < 0$. These are consistent with a scenario in which the iron emission originates from clouds illuminated by AGN X-rays; our result could provide the first extragalactic example of AGN X-ray echoes. In this scenario, we estimated the physical sizes of the illuminated clouds based on the timescale of variability to be less than 6 pc. Furthermore, we compared the iron emission distribution with the cold molecular distribution inferred by Atacama Large Millimeter/submillimeter Array (ALMA) observation of CO($J$=3--2), revealing that in the region of bright iron-line emission, the molecular emission seems to be weak. This might suggest that the AGN X-ray emission affects the chemical composition in the form of AGN feedback.

This study aims to investigate the relationship between the spatial-temporal evolution of shock properties and the longitudinal dependence of SEP intensities and spectra. The shock parameters, including the normal speed, oblique angles, compression ratio, and Alfven Mach number, were derived by combining a steady-state solar-wind simulation with the three-dimensional (3D) reconstruction of the shock surface based on multi-view observations. We compared the local shock parameters at the magnetic connecting points with in situ proton intensities and peak spectra to establish the link between shock evolution and SEP characteristics. The shock nose consistently exhibited higher particle-acceleration efficiency with the largest normal speed, compression ratio, and supercritical Alfven Mach number, while the flanks showed delayed transition to supercritical Alfven Mach number with weaker efficiency. The earliest and most rapid proton enhancement of STEREO-B correlated with efficient shock acceleration and prompt magnetic connectivity to the shock. Spectral analysis revealed that proton energy spectra were consistent with the relativistic diffusive shock acceleration (DSA) estimations. The initial shock acceleration began at about 1.4-5 Rsun and caused the widespread longitudinal SEP distribution. The longitudinal dependence of SEP intensity and spectral variations arise from the combined influence of 3D shock properties, magnetic connectivity, and particle transport processes. The agreement between in situ proton indices and relativistic DSA estimations supports DSA in this SEP event and provides insights into the early-stage acceleration at the source region.

We report the third observation of the Seyfert-1.9 active galactic nucleus (AGN) MCG-05-23-16 with the Imaging X-ray Polarimetry Explorer (\textit{IXPE}), together with optical spectro-polarimetry obtained at the Very Large Telescope (VLT), and combined with archival near-ultraviolet, optical and near-infrared polarimetric data. No X-ray polarization was detected in the 2-8 keV band, with a 99\% confidence upper limit of $\leq$2.9\%, further reduced to $\leq$2.5\% when combined with the two past IXPE observations of the same target. Monte Carlo simulations suggest that equatorial coronal models are disfavored if the AGN is indeed a type-1.9/2 AGN, but coronae coplanar to the accretion disk remain consistent if the source is less inclined than previously assumed. \textit{VLT}/FORS2 data reveal a typical type-2 spectrum in total flux, a broad H$\alpha$ line in polarized flux, and strongly wavelength dependent polarization degree and angle, rotating by nearly 70$^\circ$ across the optical band. Comparison with historical measurements confirms long-term stability of the polarization spectrum and a $\sim$90$^\circ$ rotation in the near-ultraviolet. Interpreting the multi-wavelength polarization relative to the AGN ionization axis indicates that the main obscurer is not a compact circumnuclear torus, but a distant kpc-scale dust lane crossing the galaxy. This result implies that MCG-05-23-16 is in fact a type-1 AGN seen through foreground dust. The low X-ray column density becomes consistent with the absence of polarization, provided that the nuclear inclination is low.

We present upgraded Giant Metrewave Radio Telescope (uGMRT) observations of 0731+438, an \mbox{FR II} radio galaxy at a redshift of 2.429 with two lobes separated by 82 kpc. A blueshifted, faint and broad \mbox{H{\sc i}} 21 cm absorption line with velocity full width at half maximum (FWHM) $\sim 600\,\rm km\,s^{-1}$ is detected against the southern radio lobe that is 47 kpc from radio core, indicating a neutral hydrogen outflow associated with jet-gas interaction beyond the host galaxy. The outflow has a mass outflow rate of $\sim\,0.4T_{\rm s}\Omega\rm\, M_\odot\,{\rm yr}^{-1}$, which could increase to $\sim\,4.0T_{\rm s}\Omega\rm\,M_\odot\,{\rm yr}^{-1}$, corresponding to an energy outflow rate of $2.4T_{\rm s}\Omega\times10^{40}$ -- $1.5T_{\rm s}\Omega\times10^{41}\,\rm erg\,s^{-1}$, where $T_{\rm s}$ is the spin temperature and $\Omega$ is the solid angle of the outflow. Previous optical observations identified an extended emission line region aligned with the radio axis, ionized by the central Active Galactic Nucleus (AGN). Within this region, a warm and ionized outflow with a mass outflow rate of $\sim\,50\rm\, M_\odot\,{\rm yr}^{-1}$ and an energy outflow rate of $\sim1.7\times10^{43}\,\rm erg\,s^{-1}$ was detected. We propose that both the extended emission line region and the optical outflow are results of synergistic effect between jet and AGN radiation. The AGN likely exerts negative feedback on the host galaxy, as evidenced by the gas expulsion by the jet and the high velocity dispersion of ionized gas observed optically. So far, detections of jet-driven neutral hydrogen outflows remain rare. The high redshift, large outflow radii, substantial mass outflow rate and energy outflow rate of the neutral hydrogen outflow in 0731+438 expand the known parameter space of such outflows.

Dwarf galaxies have a wide variety of structures, such as dark matter (DM) distribution, stellar-to-halo mass ratio, and stellar density. Recent high-resolution simulations have shown a variety of stellar-to-halo mass ratios for dwarf galaxies with a DM halo mass of $\sim 10^9 M_{\odot}$ at $z=0$. In this study, we performed cosmological $N$-body/smoothed-particle hydrodynamic zoom-in simulations of dwarf galaxies with the highest gas and DM particle mass resolutions of 2.37 $M_{\odot}$ and 12.8 $M_{\odot}$, respectively. The stellar-to-DM halo mass ratio of one of our simulated dwarf galaxies was $\sim 10^{-4}$, typical for satellites of the Milky Way. The stellar mass ($10^5 M_{\odot}$) and half-mass radius (68 pc) were also similar to those of the satellites of the Milky Way. The power-law slope of the DM halo was $\alpha = -1.1$. On the other hand, the other simulated galaxy exhibited a stellar-to-halo mass ratio of $\sim 10^{-3}$ and a steeper power-law slope ($\alpha=-1.9$) than the other; the presence of baryonic matter deepened the cusp. The mass of $>10^6 M_{\odot}$ and a half-mass radius of $\sim 36$ pc of this galaxy were similar to those of ultra-compact dwarf galaxies rather than the satellites of the Milky Way. This DM halo grew in mass earlier than the former one, and the central DM density was higher than that of the other even in the DM-only simulations.

Type Ia supernovae are a cornerstone of modern cosmology, providing first evidence for cosmic acceleration and new tests of dark energy. Son et al. 2025 (S25) claim a strong redshift evolution in standardized supernova luminosities driven by supernova progenitor age, with dramatic cosmological implications: rapidly evolving dark energy, decelerating expansion, and a $9\sigma$ tension with $\Lambda$CDM. We show that the underpinning evidence required for this conclusion -- the supernova progenitor-age dependence, the redshift-dependent age difference, and their combined impact -- is either negligible or relies on effects already corrected for in modern supernova analyses. First, the S25 analysis omits the standard host-galaxy stellar mass correction that captures known environmental dependencies that also correlate with stellar age. Applying this correction to the S25 sample, we find no dependence of standardized supernova brightness on host age. Independent data also show no significant difference at low-redshift in standardized brightness between star-forming galaxies and several Gyr older quiescent galaxies of the same stellar mass. Second, the S25 scenario predicts strong redshift evolution of the host-mass effect. Data from the Dark Energy Survey supernova survey measure evolution of $-0.028 \pm 0.034~\mathrm{mag}\,z^{-1}$, consistent with zero and altering the dark-energy equation-of-state measurement ($w$) by $<$0.01 if included. Third, we demonstrate that the claimed $\sim5$~Gyr progenitor age difference between nearby and distant supernovae is overstated by factors of three to five largely due to a conflation of host galaxy age with supernova progenitor age. We conclude that type~Ia supernova cosmology remains robust for current measurements of dark energy.

Doppler lensing, a relativistic effect resulting from the peculiar velocities of galaxies along the line of sight, provides insight into the large-scale structure of the Universe. Relativistic simulations are essential for modeling Doppler lensing because they incorporate gravity and motion in spacetime. We compare two relativistic $N$-body simulation frameworks, $\texttt{GEVOLUTION}$ and $\texttt{SCREENING}$, to calculate Doppler lensing convergence in cosmic voids of different sizes and halos of different masses. Our analysis reveals scale-dependent performance: $\texttt{SCREENING}$ shows larger differences in small voids (radius range: 15--25 Mpc/h) with a mean absolute relative difference of 38.5\%, due to linearized dynamics failing in nonlinear regimes. Medium voids (25--35 Mpc/h) show better agreement (9.5\% mean difference). For large voids (35--45 Mpc/h), $\texttt{SCREENING}$ exhibits intermediate differences (16.9\% mean difference) with central instabilities. Moreover, our Doppler convergence analysis with massive halos ($10^{11.5}$--$10^{14} {~h^{-1}\mathrm{M}_\odot}$) demonstrates excellent consistency (1.6--3.6\% mean difference). These findings provide clear guidance for simulation choice: $\texttt{GEVOLUTION}$ is recommended for precision studies critical to $\Lambda$CDM or modified gravity tests, while $\texttt{SCREENING}$ offers a computationally efficient alternative for relativistic treatments with large catalogs of voids and halos, assisting future astrophysical surveys.

In this study, we revisit the well-known cosmic microwave background (CMB) anomaly referred to as Hemispherical Power Asymmetry (HPA), using CMB temperature maps from the Planck mission public release 4 (PR4) and the WMAP nine-year data release. Employing the Local Variance Estimator (LVE) method, we systematically reexamine the properties of HPA to investigate possible frequency dependence as well as scale dependence in its amplitude and direction. We model the HPA as a scale-dependent dipole modulation following a power-law form, rather than assuming a scale-invariant case. Our analysis incorporates seven cleaned frequency-specific CMB temperature maps from both the Planck and WMAP missions to test the robustness of the observed asymmetry across instruments and frequency channels. We find that the dipolar modulation characteristic of HPA is present in all cases examined, with consistent estimates of the preferred direction and scale-dependent variation in dipole amplitudes. These results support the conclusion that the observed asymmetry is unlikely to arise from instrumental artifacts or data-processing effects, and instead points toward a persistent large-scale feature of the CMB sky with a possible cosmological origin.

Exocomets are small bodies that evaporate when they approach their host star. They are detected through variability of non-photospheric features with spectroscopy and/or asymmetric transits in time-series photometry. In the past four decades ~30 systems have shown such variations, and were therefore classified as exocomet host stars. However, some publications have pointed out mechanisms that might mimic exocometary features, and therefore, careful monitoring is needed to confirm the origin of the observed variability. With this paper we aim to investigate the exocomet nature of the non-photospheric variable features observed in the exocomet candidate stars HD 36546, HD 42111 and HD 85905. All of them have shown some degree of variability, particularly in their Ca II K line. We analysed the non-photospheric Ca II K line features from high-resolution spectra obtained using new NOT/FIES and Mercator/HERMES, and some additional archival spectra of the target stars. The variability was quantified through the changes in the equivalent widths of those features, which are assumed to be of circumstellar origin. Column densities were also estimated for each variable feature. Strong variability was found for HD 85905, consistent with a potential link to exocometary activity. However, the binarity of the system, which we confirmed through interferometric VLTI/PIONIER observations, complicates the interpretation of these signatures and prevents us from drawing definitive conclusions. The remaining two sources do not show any significant variability, but due to the sporadic nature of the exocometary events, we cannot discard the exocomet hypothesis. Further monitoring of the stars will be necessary to carry out a robust determination of the variability patterns and timescales that would completely rule out other scenarios.

Galaxies in the nearby Universe, particularly dwarf systems, exhibit inner mass profiles of dark matter haloes that systematically depart from canonical cold dark matter expectations, signalling an interplay between baryonic feedback and the collisionless halo. We update an analytical cusp-core transition model by incorporating the effect of supernova-driven mass loss. Adapting this model to SPARC galaxies, we measure the energy conversion efficiency epsilon, defined as the fraction of supernova feedback energy that is used to change the central dark-matter potential. We find epsilon ~ 0.01 for nearby SPARC galaxies. Building on these measurements, we compare the dynamical energy required for a cusp-core transformation with the feedback energy available over burst cycles and identify a cusp-core transition forbidden region on the halo-stellar mass plane where transformation cannot occur. Galaxies with halo masses from 10^8 to 10^11 M_sun lie outside the forbidden region, whereas ultra-faint dwarf galaxies < 10^8 M_sun, galaxy groups and clusters > 10^11 M_sun fall within it, consistent with their high central densities and the inefficiency of core formation at very low and very high masses. This approach also explains the observed diversity of inner density profiles in low-mass systems, showing that both the star formation rate and the energy conversion efficiency govern them, with the latter emerging as a key parameter setting the strength of the cusp-core transition. Beyond the cusp-core problem, our observationally inferred energy conversion efficiency provides a model independent benchmark that strongly constrains galaxy formation models.

We investigated the properties of boxy/peanut-shaped (B/PS) bulges in a sample of 71 galaxies from the Edge-on Galaxies in the Pan-STARRS Survey (EGIPS) and 20 simulated galaxies from Illustris TNG50 using multicomponent photometric decomposition. For each real and simulated galaxy, we obtained a suitable photometric model in which the B/PS bulge was represented by a dedicated 2D photometric function. For real galaxies, we found that more flattened X-structures are generally residing in larger B/PS bulges. When tested against the galaxy masses, we verified that both larger bulges and more flattened X-structures are typically found in more massive galaxies. Since large bars are also known to reside in more massive galaxies, we conclude that the flatness of X-structures in larger B/PS bulges has a physical origin, rather than being solely a result of projection effects due to differences in observed bar viewing angles. When comparing the properties of B/PS bulges between EGIPS galaxies and TNG50 galaxies, with bars rotated for different viewing angles, we found that B/PS bulges in TNG50 are considerably smaller and less luminous in terms of total intensity. This is consistent with previous studies of bar properties in TNG50, indicating the B/PS bulges in TNG50 differ from those in real galaxies, as do their parent bars.

We present a novel application of Vector quantised variational autoencoders (VQ-VAEs) to deep 1.28 GHz radio continuum images taken from the MeerKAT Galaxy Cluster Legacy Survey (MGCLS).VQ-VAEs are deep learning models widely used in modern computer vision applications and pipelines. Designed for image generation, VQ-VAEs are trained to reconstruct the input dataset via a low-dimensional discrete embedding. VQ-VAEs effectively learn the distribution of training data, where samples that do not fit the distribution well yield the highest reconstruction errors. This property makes VQ-VAEs a good candidate for the task of anomaly detection. In this work, we examine the effectiveness of VQ-VAEs in identifying radio continuum sources with anomalous structures in the image-plane domain. We find VQ-VAEs to be useful as part of a solution for searching such large datasets. We observe that they are able to remove a majority of the typical sources in such data, even when trained in an unsupervised manner on unlabelled data. We also provide our testing set of a large sample of manually labelled radio sources, in particular radio galaxies, taken from the MGCLS at 1.28 GHz. Automated approaches to searching through high volumes of data are key in extracting the full scientific potential of the Square Kilometre Array and its pathfinders.

In a gravitational lensing system, the relative transverse velocities of the lens, source, and observer induce a frequency shift in the observed radiation. While this shift is typically negligible in most astrophysical contexts, strategies for its detection have been proposed for both electromagnetic and gravitational waves. This paper provides a rigorous theoretical treatment of the effect, deriving general expressions for the frequency shift within a lensing system embedded in a cosmological spacetime. Our formulation remains valid for arbitrary distances and velocities - including highly relativistic regimes - under any Friedmann-Lemaître-Robertson-Walker metric. Expanding upon previous papers on moving lenses, we provide a detailed analysis of frequency effects induced by lenses moving at relativistic speeds. Furthermore, we extend standard lensing theory by deriving an exact formula for the variation in the source's viewing direction. This result is of interest for strongly anisotropic emitters, such as compact binary systems emitting gravitational waves. Finally, we quantify the apparent misalignment between the lens and the source's two images produced by time-delay effects in lens systems moving with high velocity.

Primitive asteroids include C-, P-, and D-classes, known to be dark and having spectra mostly featureless. They differ in the spectral slope, which ranges from moderate values for C-types, and progressively increases in P- and D-types, the latter being the reddest. While C- and P-types are commonly observed in the asteroid main belt, D-types are commonly found further from the Sun, in the Cybele, Hilda, and Jupiter Trojans regions, and very few are reported in the main belt. This study aims at characterizing the abundance of primordial and red asteroids, belonging to the P-, D-, and Z-classes in the Mahlke et al. (2022) taxonomy, in the 2-5.2 AU region using the third data release by the Gaia mission spectral catalog, which includes more than 60000 spectrophotometric data of asteroids. We have applied the following methodology to identify primordial asteroids in the catalog: 1) selection of objects with signal to noise ratio greater than 20; 2) albedo value less than 12%; 3) chi-squared fit to automatically identify potential D-, Z-, and P-types using Bus-DeMeo and Mahlke taxonomy; 4) visual inspection of every spectrum to confirm the taxonomic classification. Referring to Mahlke taxonomy, we have found 318 new D-types across the main belt, as well as 124 Z-types, and is in agreement with theoretical estimations. We computed the spectral slope in the visible range (0.55 - 0.81 \mu m). We also have identified 265 P-types in the main belt. For the Cybele and Hilda asteroids, we characterize the taxonomic class of all the bodies with SNR higher than 20 in the Gaia catalog, for a total sample of 193 and 180 asteroids, respectively.

PSR J2021+4026 is a remarkable $\gamma$-ray pulsar exhibiting repeated transitions between high $\gamma$-ray flux (HGF) and low $\gamma$-ray flux (LGF) states. With 17-yr Fermi-LAT monitoring, we reveal persistent secular evolution and enhanced spin-down rate variability within individual emission states -- beneath the quasi-periodic state transitions. After removing discrete jumps, the jump-corrected flux $\delta F_\gamma$ shows a three-phase evolution: rise ($+2.02^{+0.17}_{-0.15}\%~\mathrm{yr}^{-1}$), decline ($-3.72^{+0.34}_{-0.47}\%~\mathrm{yr}^{-1}$), and rapid rise ($+14.9^{+6.4}_{-4.4}\%~\mathrm{yr}^{-1}$), with all rates quoted relative to the long-term mean flux $\langle F_\gamma \rangle=7.8\times 10^{-10}\,\mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$. Moreover, the flux of the LGF state is gradually approaching the stable HGF level at a rate of $+0.72 \pm 0.11\%~\mathrm{yr}^{-1}$. These results demonstrate that secular flux evolution in PSR J2021+4026 operates largely independently of discrete state transitions, yet jointly with them drives the system toward a stable high-flux equilibrium.

We present the first systematic study of rigidly rotating protoquark stars based on isentropic equations of state (EOS) within the density-dependent quark mass (DDQM) framework. Using a quasi-static equilibrium approach, we follow the Kelvin--Helmholtz evolution from hot, lepton-rich matter to a cold, catalyzed quark star. Rotation substantially enhances the maximum stable mass (by up to $\sim 40\%$), equatorial radius, and key rotational observables, with the ratio of rotational kinetic to gravitational potential energy, $T_{\rm kin}/|W|$, reaching $0.18$--$0.19$ near the Keplerian limit, indicating a heightened susceptibility to gravitational-wave--emitting instabilities. Thermal evolution introduces a clear ordering: all stellar properties peak during the lepton-rich stages and decrease monotonically as the star cools. Compared to hadronic stars, rotating protoquark stars exhibit larger radii, higher moments of inertia, and stronger quadrupolar deformation, producing a distinct signature in the mass--radius--spin plane that can accommodate objects such as HESS~J1731--347 and PSR~J0740+6620. These results demonstrate that future multimessenger observations must account for both thermal history and rotation to robustly identify quark matter in compact stars.

Pulsar timing arrays (PTAs) have recently reported compelling evidence for the presence of a gravitational-wave background signal. Mapping the gravitational-wave background is key to understanding how it is formed, since anisotropy is a tracer for, for example, a supermassive black hole binary origin. In this work we refine the frequentist regularised gravitational-wave mapping analysis developed in our previous work (as part of the MeerKAT PTA 4.5-year data release). We derive a point-spread function describing the angular resolution of a PTA. We investigate how the point spread function changes for different PTA constellations and determine the best possible angular resolution achievable within our framework. Using simulated data, we demonstrate that previous methods do not capture the actual resolution - especially in regions of the sky with a high density of pulsars. We propose an improved scheme that accounts for a variable local resolution and test it using realistic simulations of the latest MeerKAT dataset. We demonstrate that we are able to identify a continuous gravitational wave signal in a region with good pulsar sky coverage with approximately a factor of two increase in significance compared to our previous method.

The luminous accreting pulsar SMC X-1 is an appropriate target to explore the accretion dynamics. SMC X-1 shows unique quasi-periodic flux variability of 40-65$\,$days known as superorbital modulation. To constrain the accretion structure of SMC X-1 based on timing and spectral study, we have analysed X-ray data of SMC X-1 observed by Suzaku and NuSTAR at various epochs between 2011 and 2022. The spectral analysis shows that the hydrogen column density ($N_\mathrm{H}$) increases from $1.1 \times 10^{22}\,\mathrm{cm^{-2}}$ to $1.24 \times 10^{23}\,\mathrm{cm^{-2}}$ as the flux decreases with the superorbital modulation. The neutral iron K$\alpha$ line at 6.4$\,$keV has a broad width of 0.3$\,$keV, and its equivalent width increases as toward superorbital low states. The line broadening is consistent with Keplerian motion at the inner disc rather than the stellar wind velocity of the donor star. These findings support that the superorbital modulation is a consequence of X-ray attenuation by the warped accretion disc. To test this interpretation, we have conducted photon transport simulations of a system consisting of a neutron star, a warped disc, and optically-thin disc atmosphere. Occultation of the central source by the disc successfully reproduces the observed variations in the equivalent width of neutral iron K$\alpha$ line, pulse profiles, and flux in hard X-rays. Notably, a disc precession angle of approximately $30^\circ$ can account for the observational features. For the radiation pattern of the photon source, the preferred beam width corresponds to a standard deviation of $30^\circ$.

In this work, we explore the possibility of using quantum computers provided for usage in cloud by big companies (such as IBM, IonQ, IQM Quantum Computers, etc.) to run our quantum neural network (QNN) developed for data analysis in the context of LISA Space Mission, developed with the Qiskit library in Python. Our previous work demonstrated that our QNN learns patterns in gravitational wave (GW) data much faster than a classical neural network, making it suitable for fast GW signal detection in future LISA data streams. Analyzing the fees from hardware providers like IBM Quantum, Amazon Braket and Microsoft Azure, we found that the fees for running the first segment of our QNN sum up to \$2000, \$60000, and \$1000000 respectively. Using free plans, we succeed to run the 3-qubit feature map of the QNN for one random data sample on {\fontfamily{qcr} \selectfont ibm\_kyoto} and {\fontfamily{qcr}\selectfont IQM Quantum Computers\_Garnet} quantum computers, obtaining a fidelity of 99\%; we could also run the first prediction segment of our QNN on {\fontfamily{qcr} \selectfont ibm\_kyoto}, implemented for 4 qubits, and obtained a prediction accuracy of 20\%. We queried providers such as IBM Quantum, Amazon Braket, Pasqal, and Munich Quantum Valley to obtain access to their plans, but, with the exception of Amazon Braket, our applications remain unanswered to this day. Other major setbacks in using the quantum computers we had access to included Qiskit library version issues (as in the cases of IBM Quantum and IQM Quantum Computers) and the frequent unavailability of the devices, as was the case with the Microsoft Azure provider. All the results presented in this paper were accumulated in 2024.

We present the discovery and characterisation of two transiting exoplanets orbiting the metal-poor, solar-type star TOI-5788. From our analysis of six \textit{TESS} sectors and a dedicated \textit{CHEOPS} programme, we identify an inner planet (TOI-5788~b; $P = 6.340758\pm0.000030\,\si{\day}$) with radius $1.528\pm0.075\,\mathrm{R_\oplus}$ and an outer planet (TOI-5788~c; $P = 16.213362\pm0.000026\,\si{\day}$) with radius $2.272\pm0.039\,\mathrm{R_\oplus}$. We obtained 125 radial-velocity spectra from HARPS-N and constrain the masses of TOI-5788~b and~c as $3.72\pm0.94\,\mathrm{M_\oplus}$ and $6.4\pm1.2\,\mathrm{M_\oplus}$, respectively. Although dynamical analyses indicate that a third planet could exist in a stable orbit between 8 and 14 days, we find no evidence of additional planets. Since the TOI-5788 system is one of the few systems with planets straddling the radius gap, and noting that there are even fewer such systems around metal poor stars, it is a promising system to constrain planet formation theories. We therefore model the interior structures of both planets. We find that TOI-5788~b is consistent with being a rocky planet with almost no envelope, or having an atmosphere of a high mean molecular weight. We find that TOI-5788~c is consistent with both gas-dwarf and water-world hypotheses of mini-Neptune formation. We model the atmospheric evolution history of both planets. Whilst both scenarios are consistent with the atmospheric evolution of TOI-5788~c, the gas-dwarf model is marginally preferred. The results of the atmospheric evolution analysis are not strongly dependent on stellar evolution. This makes the system a promising target to test internal structure and atmospheric evolution models.

Fast Radio Bursts (FRBs) are millisecond-duration radio transients of extragalactic origin. Classifying repeating FRBs is essential for understanding their emission mechanisms, but remains challenging due to their short durations, high variability, and increasing data volume. Traditional methods often rely on subjective criteria and struggle with high-dimensional data. In this study, we apply an unsupervised machine learning framework that combines Uniform Manifold Approximation and Projection (UMAP) and Hierarchical Density-Based Spatial Clustering of Applications with Noise (HDBSCAN) to eight observed parameters from FRB 20220912A. Our analysis reveals three distinct clusters of bursts with varying spectral and fluence properties. Comparisons with clustering studies on other repeaters show that some of our clusters share similar features with sources such as FRB 20201124A and FRB 121102, suggesting possible common emission mechanisms. We also provide qualitative interpretations for each cluster, highlighting the spectral diversity within a single source. Notably, one cluster shows broadband emission and high fluence, which are typically seen in non-repeating FRBs. This raises the possibility that some non-repeaters may be misclassified repeaters due to observational limitations. Our results demonstrate the utility of machine learning in uncovering intrinsic diversity in FRB emission and provide a foundation for future classification studies.

Detecting small planets via the radial velocity method remains challenged by signals induced by stellar variability, versus the effects of the planet(s). Here, we explore using Gaussian Process (GP) regression with Transiting Exoplanet Survey Satellite (TESS) photometry in modeling radial velocities (RVs) to help to mitigate stellar jitter from oscillations and granulation for exoplanet detection. We applied GP regression to simultaneous TESS photometric and RV data of HD 5562, a G-type subgiant ($M_\star=1.09M_{\odot}$, $R_\star=1.88R_{\odot}$) with a V magnitude of 7.17, using photometry to inform the priors for RV fitting. The RV data is obtained by the Magellan Planet Finder Spectrograph (PFS). The photometry-informed GP regression reduced the RV scatter of HD~5562 from 2.03 to 0.51 m/s. We performed injection and recovery tests to evaluate the potential of GPs for discovering small exoplanets around evolved stars, which demonstrate that the GP provides comparable noise reduction to the binning method. We also found that the necessity of photometric data depends on the quality of the RV dataset. For long baseline and high-cadence RV observations, GP regression can effectively mitigate stellar jitter without photometric data. However, for intermittent RV observations, incorporating photometric data improves GP fitting and enhances detection capabilities.

This is a collection of papers presented by the JEM-EUSO Collaboration at the 39th International Cosmic Ray Conference (ICRC 2025) (Geneva, Switzerland, July 14--24, 2025).

We investigate accretion disk dynamics and the formation of quasi-periodic oscillations (QPOs) in the infrared limit around Kerr-like black holes in asymptotically safe gravity. Relativistic hydrodynamic solutions of Bondi-Hoyle-Lyttleton (BHL) accretion reveal that quantum corrections significantly modify the structure of the shock cone formed around the black hole. The black hole spin controls the asymmetric of the shock cone through frame-dragging effects, whereas the quantum correction parameter softens the effective gravitational potential, resulting in a wider shock opening angle, weaker post-shock compression, and reduced density concentration within the cone. Time-dependent mass accretion rates reveal oscillation modes trapped within the shock cone. The power spectral density (PSD) investigations suggest that these modes naturally generate low-frequency QPOs, whose amplitudes, coherence, and harmonic structure depend on both the spin and the quantum correction parameter. The PSD analyses performed at different radial locations reveal that identical QPO frequencies are obtained in all cases. The numerically detected frequencies result from the excitation of global oscillation modes trapped within the post-shock region. The resulting global modes are found to consist of fundamental frequencies, their associated harmonic overtones, and near-commensurate frequency ratios such as 2:1 and 3:2. Coherent oscillations are enhanced and near-commensurate frequency ratios are produced when moderate rotation and moderate quantum corrections are coupled. Large quantum correction parameters, on the other hand, wash out unique spectral peaks and suppress oscillation amplitudes.

We present multi-wavelength observations of the first recorded low state of the intermediate polar BG CMi. Optical monitoring of the source by members of the American Association of Variable Star Observers reveals a decrease of ~0.5 mag that lasted ~50 d in early 2025. During the low state the optical timing properties imply that BG CMi underwent a change in the accretion mode, as power at the spin frequency $\omega$ dramatically dropped. An XMM-Newton observation revealed a substantial decrease in intrinsic absorption and a slight increase in intrinsic X-ray luminosity, compared to archival Suzaku data. Timing analysis of the X-ray light curves shows that power shifted from the orbital frequency $\Omega$ (prominent in Suzaku data) to $2\Omega$ in the low state XMM-Newton data, along with the strengthening of certain orbital sidebands. We suggest that BG CMi transitioned to disk-overflow accretion, where the white dwarf accreted matter via both a disk and a stream, the latter becoming more dominant during the low state due to a decrease in the mass and size of the disk.

The unprecedented wavelength coverage and sensitivity of the James Webb Space Telescope (JWST) permits to measure the absorption features of a wide range of condensate species from Silicates to Titan tholins. Atmospheric retrievals are uniquely suited to analyse these datasets and characterize the aerosols present in exoplanet atmospheres. However, including the optical properties of condensed particles within retrieval frameworks remains computationally expensive, limiting our ability to fully exploit JWST observations. In this work, we improve the computational efficiency and scaling behavior of aerosol models in atmospheric retrievals, enabling in-depth studies including multiple condensate species within practical time scales. Rather than computing the aerosol Mie coefficients for each sampled model, we pre-compute extinction efficiency (Qext), scattering efficiency (Qscat) and asymmetry parameter (g) grids for seven condensate species relevant in exoplanet atmospheres (Mg2SiO4 amorph sol - gel, MgSiO3 amorph glass, MgSiO3 amorph sol - gel, SiO2 alpha, SiO2 amorph, SiO and Titan tholins). The pre-computed Qext grids significantly reduce computation time between 1.4 and 17 times with negligible differences on the retrieved parameters. They also scale effortlessly with the number of aerosol species while maintaining the accuracy of cloud models. Thereby enabling more complex retrievals as well as broader population studies without increasing the overall error budget. The Qext, Qscat and g grids are freely available on Zenodo as well as a public TauREx plugin -TauREx-PCQ- that utilize them.

Asteroseismology is widely used for precise determining of masses of solar-like oscillating stars by performing individual-frequency modeling or applying homological scaling relations. However, these methods lack dynamical validation on the main sequence due to the absence of eclipsing double-lined binary system (SB2) as benchmark objects. By providing the orbital inclination, astrometric binary systems from ESA Gaia DR3 offer an abundant alternative for eclipsing systems. We present KIC693187 as the first SB2, hosting a solar-like oscillating post-main-sequence star with dynamical masses. By combining Gaia astrometry with spectroscopic obtained with the Las Cumbres Observatory network (LCO), we find $M_1^\mathrm{dyn}$=0.99$\pm$0.05$M_\odot$ and $M_2^\mathrm{dyn}$=0.89$\pm$0.04$M_\odot$ for the primary and secondary, respectively. Asteroseismic parameters were extracted from photometry of the NASA \Kepler satellite. The mass from individual frequency modeling is $M_1^\mathrm{IF}$=0.92$\pm$0.01$M_\odot$. Taking into account the systematic uncertainty of 0.04$M_\odot$ for best fit models from individual frequency fitting, we find an agreement within 1.2$\sigma$. From scaling relations we obtain a mass range of 0.93 to 0.98$M_\odot$ by using the observed large frequency separations (\dnu) in the scaling relations for the primary. By using standard corrections for departures from the asymptotic regime of \dnu, we obtained a mass range of 0.83 to 1.03$M_\odot$. The upper ends of both ranges agree well with the dynamical mass of the primary. This approach provides the first empirical validation for main-sequence solar-like oscillators and opens a new window for validating asteroseismology. Through a dedicatded program targeting astrometric SB2 binary systems, ESA's PLATO space mission will provide will enlarge the benchmark sample substantially.

Binary stars are as common as single stars. The number of detected planets orbiting binaries is rapidly increasing thanks to the synergy between transit surveys, Gaia and high-resolution direct imaging campaigns. However, global planet formation models around binary stars are still underdeveloped, which limits the theoretical understanding of planets orbiting binary star systems. Hereby we introduce the PAIRS project, which aims at building a global planet formation model for planets in binaries, and to produce planet populations synthesis to statistically compare theory and observations. In this first paper, we present the adaptation of the circumstellar disc to simulate the formation of S-type planets. The presence of a secondary star tidally truncates and heats the outer part of the circumprimary disc (and vice-versa for the circumsecondary disc), limiting the material to form planets. We implement and quantify this effect for a range of binary parameters by adapting the Bern Model of planet formation in its pebble-based form and for in-situ planet growth. We find that the disc truncation has a strong impact on reducing the pebble supply for core growth, steadily suppressing planet formation for binary separations below 160 au, when considering all the formed planets more massive than Mars. We find as well that S-type planets tend to form close to the central star with respect to the binary separation and disc truncation radius. Our newly developed model will be the basis of future S-type planet population synthesis studies.

Roughly half of Sun-like stars have at least one stellar companion, whereas it is widely assumed that most known exoplanets orbit single stars, largely due to observational biases. However, astrometric surveys, direct imaging, and speckle interferometry are steadily increasing the number of confirmed exoplanets in binaries. A stellar companion introduces additional effects, such as circumstellar disk truncation and gravitational perturbations, which can strongly impact planet formation. While global planet formation models, for example the Bern model, have been broadly applied to single stars, modeling S-type binaries requires key modifications to capture these effects. This study extends the Bern model by incorporating the gravitational influence of a stellar companion into its N-body integrator, allowing us to quantify how this perturbation affects planetary formation and final system architecture across a range of binary configurations. By comparing binary and single-star systems under identical initial conditions, we can assess the specific impact of binary-induced dynamics. We ran three sets of simulations: (i) a grid of in situ single-embryo cases to quantify gravitational effects; (ii) formation simulations with and without migration to compare outcomes with single-star analogs; and (iii) multi-embryo runs to evaluate impacts on multi-planetary systems. Planets forming beyond half the host star's Hill radius are much more likely to become unbound especially in systems with high binary eccentricity. Even within stable zones, growth is suppressed by both reduced material availability and increased eccentricity from stellar perturbations. Both disk truncation and stellar perturbations must be included to model planet formation in S-type binaries accurately. Neglecting either one will end up misrepresenting planetary growth and survival.

The forbidden lines of [Fe XIII] at 10,747 Angstrom and 10,798 Angsrtom are among the most prominent lines in the near-infrared spectrum of the solar corona. They have been used routinely, both outside and during eclipses, as sensitive probes of the electron density and polarization in the solar corona. Many novae pass through a coronal phase, wherein the highly ionized nova ejecta have physical conditions that are remarkably similar to those of the solar corona. Many of the coronal emission lines that are seen are common to the spectra of both the Sun and novae. Yet, it appears that no robust detection of the [Fe XIII] lines has been made in a nova. Here we report the detection of these two infrared [Fe XIII]lines in the spectrum of the recurrent nova V3890 Sgr, taken 23.43 and 31.35 days after its August 2019 outburst. From their line strengths, we derive values of 10^10 per cubic cm and 10^[8.5-9] per cubic cm for the electron density on the two. The decrease in density between epochs can be explained if the density decreased with a power law n ~ r**alpha with a alpha inferred to be -3. The average temperature of the coronal gas is estimated to be T = (2.51\pm0.06) x 10^6~K. We find that recurrent novae with giant secondaries, including T CrB whose eruption is imminent, are the most suitable sources for further detections of the [Fe XIII] lines. epochs.

The evolution of radio jet morphology and its energetics is significantly influenced by the environment in which the host galaxy resides. As giant radio galaxies (GRGs) often extend to the scale of entire galaxy clusters ($\sim$Mpc) and beyond, they are a suitable class of objects for studying jet--intracluster medium interactions. This paper presents a multiwavelength study of a GRG, J1007+3540, using the LOFAR Two-metre Sky Survey second data release (LoTSS DR2) at 144 MHz and the upgraded Giant Metrewave Radio Telescope (uGMRT) at 400 MHz. The source has a projected linear extension of 1.45 Mpc and is hosted by MaxBCG J151.77665+35.67813, within the WHL 100706.4+354041 cluster. At both frequencies, the source exhibits clear signatures of recurrent jet activity, a one-sided, extended, tail-like diffuse structure with a morphological break in the tail. The estimated radiative ages of the inner lobes and outer north lobe are $\sim$140 Myr and $\sim$240 Myr, respectively. In addition to the radio analysis, we performed optical--to--infrared spectral energy distribution modelling. The host galaxy is an evolved elliptical system with a stellar mass of $\log_{10}(M_\star/M_\odot) = 11.0$ and an old stellar population age of $\sim$12 Gyr. The high infrared-derived star formation rate ($\sim106~M_\odot$~yr$^{-1}$) of the source implies significant dust-obscured star formation, potentially linked to merger-driven gas inflows. J1007+3540 presents a rare combination of a restarted jet, a detached tail-like structure, and unusual spectral flattening beyond the tail break, which is very rare to report together in a GRG. This rare and remarkable system offers a unique laboratory for probing the interplay between active galactic nucleus activity, star formation, and environmental effects in cluster-surrounded GRGs.

The persistent Hubble tension and the lack of a fundamental explanation for dark energy motivate the exploration of alternative mechanisms capable of reproducing late-time cosmic acceleration. In this work, we revisit gravitationally induced particle creation as a phenomenological non-equilibrium process that can effectively mimic a dynamical dark-energy component. Within the thermodynamic framework of open systems, we model the production of an unspecified particle species with constant intrinsic equation-of-state parameter and consider four phenomenological parametrisations of the particle-creation rate. The modified continuity and Friedmann equations lead to an effective negative pressure and a redshift-dependent effective equation of state, which we constrain using Cosmic Chronometers, Pantheon+ supernovae, DESI DR2 BAO, a compressed CMB likelihood, and SH0ES data. Using the full dataset combination, we find that particle-creation models provide fits comparable to $\Lambda$CDM, yielding $H_0 \simeq 69.3\,\mathrm{km\,s^{-1}\,Mpc^{-1}}$ and present-day effective dark-energy equation-of-state values close to $w^{\rm eff}_{\rm DE}(0)\simeq -1$, with all models predicting an accelerating Universe ($q_0\simeq -0.55$). When the Hubble tension is assessed using early- and late-time dataset splits, particle-creation scenarios reduce its statistical significance to the $\simeq 2.4\sigma$--$3\sigma$ level, compared to the $4.3\sigma$ discrepancy obtained in $\Lambda$CDM. Although deviations from $\Lambda$CDM remain mild and Bayesian model comparison indicates no statistical preference between models, gravitationally induced particle creation emerges as a viable late-time extension of the standard cosmological model and provides a consistent phenomenological framework for exploring departures from $\Lambda$CDM.

We reconsider back-reaction from large amplitude, short-scale perturbations onto a long wavelength adiabatic mode. In a loop expansion of the long-mode power spectrum, this back-reaction appears first at 1-loop. Due to the separation between the long and short scales, the separate universe method provides a simple and efficient framework for this computation. In this paper, building on our earlier work, we employ a $\delta N$ formula for the long mode, which captures the effect of short scales. We show that back-reaction at 1-loop is due to either (i) non-linearity of the $\delta N$ formula, or (ii) 1-loop corrections to the initial conditions. We argue that contributions of type (ii) cannot themselves be described within the separate universe framework, but their properties can be constrained using soft theorems and a ''multi-point propagator'' expansion. When applied to a band of enhanced short-scale perturbations that crossed the horizon during inflation, our result shows that the loop correction decouples from their detailed properties. Furthermore, the back-reaction we obtain is scale-invariant. Its magnitude is model-dependent, but is degenerate with effects from modes that were still sub-horizon at the end of inflation. In this scenario (but not necessarily in all scenarios), we conclude that the effect is not observable.

An increasing number of discoveries of isolated and quenched dwarf galaxies are challenging the idea that the present-day local environment of low-mass systems is the main determinant of their quenching. We present new Hubble Space Telescope (HST) data of one such system, the dwarf galaxy Canes Venatici C (CVn C). CVn C is a low-mass (3.4(+4.2-2.6)*10^6 M_sun) galaxy with a Tip of the Red Giant Branch distance of 8.43(+0.47-0.32) Mpc determined from the resolved stars in the HST imaging, which we also use to derive CVn C's structural parameters. CVn C's distance places CVn C in the Local Volume and in an isolated environment with the most tidally influential L* galaxy > 5Rvir away. Additional constraints from the HST color-magnitude diagram, archival Far-Ultraviolet (FUV), and neutral hydrogen (HI) data show that CVn C is quenched, with no evidence of star formation in the last 100 Myr and no detectable gas (MHI < 1.5*10^6 M_sun). Circumstantial evidence suggests that CVn C may have quenched via past interactions with the L* galaxy NGC 4631 (L_K = 10^10.4 L_sun), and was possibly sent on an extreme backsplash orbit by the tidal dissolution of a subhalo group. However, other quenching mechanisms-such as stripping via the cosmic web-cannot be ruled out. CVn C adds to the growing number of quenched dwarf galaxies in under-dense environments, a population that will be critical to defining the mass and environment regimes in which different quenching mechanisms operate.

The internal composition of sub-Neptunes remains a prominent unresolved question in exoplanetary science. We present a technique to place constraints on envelope mean molecular weight that utilises observations of escaping hydrogen or helium exospheres. This method is based on a simple timescale argument, which states that sub-Neptunes require a sufficiently large hydrogen or helium reservoir to explain on-going escape at their observed rates. This then naturally leads to an upper limit on atmospheric mean molecular weight. We apply this technique to archetypal sub-Neptunes, namely GJ-436 b, TOI-776 b and TOI-776 c, which have all been observed to be losing significant hydrogen content as well as relatively featureless transit spectra when observed with JWST. Combining constraints from atmospheric escape and transit spectroscopy in the case of TOI-776 c allows us to tentatively rule out the high mean molecular weight scenario, pointing towards a low mean molecular weight atmosphere with high-altitude aerosols muting spectral features in the infra-red. Finally, we reframe our analysis to the hycean candidate K2-18 b, which has also been shown to host a tentative escaping hydrogen exosphere. If such a detection is robust, we infer a hydrogen-rich envelope mass fraction of $\log f_\text{env} = -1.67\pm0.78$, which is inconsistent with the hycean scenario at the $\sim 4\sigma$ level. This latter result requires further observational follow-up to confirm.

GRB~221009A has been the brightest gamma-ray burst (GRB) observed to date, and its afterglow has been characterised with unprecedented detail at TeV energies by LHAASO. Quite puzzlingly, it is also the most energetic GRB known. Among the riddles posed by this mysterious source, however, the sheer energetics are hardly the most intriguing: an unprecedented, narrow, luminous emission line at around 10 MeV has been uncovered by a detailed spectral analysis of \textit{Fermi}/GBM data immediately following the brightest peak in the GRB prompt emission and the peak of the TeV afterglow. As noted in the discovery article, the temporal evolution of the line properties can be explained as being due to high-latitude emission from a geometrically thin, relativistically expanding shell where annihilation of a large number of electron-positron pairs took place. We show that this interpretation yields stringent constraints on the properties of such shell, that point to a process that happens at radii typical of external shocks. We then demonstrate that the shell could have been the blastwave associated with the GRB precursor, with the line arising after pair loading of such blastwave as it was illuminated by the bright and hard radiation of the GRB main event. The scenario, which also explains the abrupt initial rise of the LHAASO afterglow, requires the progenitor of the GRB to have been surrounded by a circum-stellar medium (CSM) extending out to a few $10^{15}\,\mathrm{cm}$, with a density $n_\mathrm{ext}\sim 10^{8}-10^{9}\,\mathrm{cm^{-3}}$ reminiscent of those found from studies of Type IIn supernovae. This provides a precious clue to the nature of the progenitor of this peculiar GRB, which could also be present in other bursts that feature a long quiescence followed by a bright emission episode with a hard spectrum.

We refine the constraints on the lifetime of decaying super-heavy dark matter particles (SHDM), with masses ranging from $10^7$ to $10^{15}$ GeV, by analyzing ultra-high-energy (UHE) gamma-ray data. Our approach involves an accurate comparison of the primary gamma-ray emissions resulting from prompt SHDM decays in the galactic halo with the most recent upper limits on isotropic UHE gamma-ray fluxes provided by various extensive air shower experiments. We demonstrate that a precise consideration of the field of view and the geometric acceptance of different UHE gamma-ray observatories has significant implications for the inferred limits of dark matter lifetime. In addition, we examine the influence of uncertainties linked to the current models of the galactic dark matter distribution, employing diverse halo density profiles while varying both their radial extent and the local dark matter density. Our findings indicate that the newly established UHE gamma-ray constraints are marginally less stringent than earlier evaluations, thereby revisiting the SHDM parameter space and allowing for observable neutrino fluxes.

Scalar-tensor theories predict fundamental scalar fields of considerable interest in astrophysics and cosmology. We investigate the superradiant instability of scalar clouds around Kerr black holes, showing that stimulated decay generates detectable electromagnetic signals. The growth of the superradiant scalar cloud differs from that of other bosonic fields and depends sensitively on the matter distribution surrounding the black hole, which originates from the scalar-matter coupling realized by the chameleon mechanism in modified gravity theories. In non-uniform matter distributions, stimulated emission from scalar clouds offers an observational signature that distinguishes fundamental scalars from other light bosonic fields.

We derive the affine tensor associated with the energy and momentum densities of both gravitational and matter fields, the complex pseudo-tensor, for $f(Q)$ non-metric gravity, the straightforward extension of Symmetric Teleparallel Equivalent of General Relativity (STEGR), characterized by a flat, torsion-free, non-metric connection. The local conservation of energy-momentum complex on-shell is satisfied through a continuity equation. An important analogy is pointed out between gravitational pseudo-tensor of teleparallel $f(T)$ gravity, in the Weitzenböck gauge, and the same object of symmetric teleparallel $f(Q)$ gravity, in the coincident gauge. Furthermore, we perturb the gravitational pseudo-tensor $\tau^{\alpha}_{\phantom{\alpha}\lambda}$ in the coincident gauge up to the second order in the metric perturbation, obtaining a useful expression for the power carried by the related gravitational waves. We also present an application of the gravitational pseudotensor, determining the gravitational energy density of a Schwarzschild spacetime in STEGR gravity, adopting the concident gauge. Finally, analyzing the conserved quantities on manifolds, the Stokes theorem can be formulated for generic affine connections

We propose a simple yet testable framework for light fermion dark matter (DM) with mass in the MeV--GeV range, charged under a dark $U(1)_D$ gauge symmetry. The $U(1)_D$ is spontaneously broken by a scalar field $\Phi$, giving mass to the dark gauge boson $X_D$. The dominant DM annihilation proceeds via a forbidden channel, where the DM pair annihilates into slightly heavier dark gauge bosons and scalars after the dark-sector phase transition. Once the dark-sector phase transition occurs, the induced mass gap activates the forbidden annihilation channel, which in turn determines the DM relic abundance and naturally suppresses late-time annihilation. As a result, the scenario avoids stringent cosmic microwave background and indirect detection constraints that typically exclude thermal light DM. Moreover, the same symmetry-breaking phase transition is strongly first-order, producing a stochastic gravitational wave background that could be probed by upcoming space-based interferometers and pulsar timing arrays. We demonstrate that achieving the observed DM abundance tightly correlates the DM mass with the nucleation temperature of the phase transition. Thus, this setup links the DM relic abundance, dark-sector dynamics, and gravitational wave signals, offering complementary paths for discovery in both terrestrial and cosmological observations.

The primordial 7Li abundance predicted by standard Big Bang Nucleosynthesis (BBN) exceeds that inferred from old, metal-poor stars by a factor of about 3-4. In standard BBN, most primordial 7Li is produced as 7Be in the early Universe and later converted by electron capture. Additional production or destruction channels of 7Be, such as proton capture or antineutrino capture during BBN, may therefore affect the final lithium yield. We quantify the depletion of 7Be due to in-situ electron capture, including the associated antineutrino channel, positron decay from nuclear excited states, and proton capture through the radiative 7Be(p,gamma)8B reaction. We also investigate stimulated emission induced by the dense photon background during the nuclear statistical equilibrium epoch, as well as a three-body Auger-like variant transferring the capture energy to a continuum electron. Decay rates are computed using first-order perturbation theory, modelling weak interactions with a Fermi contact term and factorising hadronic and leptonic currents. Thermally averaged rates are obtained by folding cross-sections with Maxwell-Boltzmann distributions and accounting for particle densities in the temperature range 10-100 keV. We find that the electron-capture rate decreases rapidly with temperature and is significantly enhanced by the inclusion of the antineutrino channel. Stimulated emission and plasma screening increase the radiative proton-capture rate by only 1-3 percent at temperatures around 87 keV. The Auger-like channel contributes at the level of a few thousandths of a percent and becomes negligible at lower temperatures. Overall, our total rate revises previous estimates by nearly an order of magnitude. Electron capture, proton capture, and positron decay provide corrections to the dominant depletion channel 7Be(n,p)7Li.

We investigate the nonlinear formation of plasmoids in 2D low-beta current sheets through the interplay between the Kelvin-Helmholtz instability (KHI) and the lower-hybrid drift instability (LHDI). Using a hybrid kinetic-gyrokinetic model-based Super Simple Vlasov (ssV) code with fully kinetic ions and drift-kinetic electrons, we simulate Harris-type current sheets and velocity shear layers with strong cross-field density gradients. Our central hypothesis is that steep density gradients drive LHDI, which can grow faster than KHI and initiate an inverse cascade from kinetic to fluid scales, potentially suppressing KHI. Our simulations confirm that, in thin current sheets, LHDI develops rapidly at the sheet edges and nonlinearly merges into larger-scale magnetic islands before KHI can evolve. These LHDI-driven structures distort the velocity shear and suppress classical KH vortices. In contrast, for thicker current sheets or weaker density gradients, KHI dominates and produces the expected rolled-up vortices and associated plasmoids. These findings demonstrate that LHDI-induced turbulence can act as both a seed and a regulator of plasmoid-generating instabilities, mediating cross-scale energy transfer. This mechanism is relevant to thin boundary layers in space plasmas, such as the solar wind magnetosphere interface, and suggests that microturbulence can govern large-scale magnetic topology during collisionless reconnection.

Hokkaido University and Tohoku University have been developing and operating a constellation of 50-cm-class microsatellites for Earth observation. DIWATA-1, launched in 2016, was deployed into a circular orbit at an altitude of approximately 400 km from the International Space Station (ISS). For the subsequent satellite developed in 2021, the structural design and vibration test campaign were optimized to meet a strict one-year development schedule. This paper summarizes how the structural design of the previous satellite was reviewed and updated, and how the vibration test was successfully completed in a single trial to minimize schedule and technical risks. These lessons learned provide valuable insights, as there are only a limited number of reported cases of 50-kg-class microsatellites deployed from the ISS.

Collective excitations are fundamental in quantum many-body physics, yet their spectra have traditionally been studied within Newtonian dynamics. In this Letter, we investigate collective excitations in Bose gases under Modified Newtonian Dynamics (MOND). We derive an anisotropic excitation spectrum in the MOND regime. This anisotropy arises directly from the intrinsic nonlinear structure of the MOND Poisson equation, forming a distinctive signature of the modified gravitational response. We then analyze the Jeans instability, obtaining analytic expressions for the direction-dependent critical wavelength and mass. These results advance our understanding of collective behavior in quantum systems under modified dynamics and establish clear theoretical signatures for testing MOND-like effects in quantum simulators.

This paper describes the system design methodology derived from the development and evaluation tests of deployable solar panels to be mounted on a 3U CubeSat. The study mainly includes structural analysis, thermal analysis, and a review of vibration test results. Hokkaido University is developing the 3U CubeSat HOKUSHIN-1 in collaboration with Tohoku University and Muroran Institute of Technology. Deployable solar panels are a key technology for future planned lunar exploration missions, as they enable power-intensive communication and propulsion required for orbit control. The satellite also demonstrates a newly developed compact and efficient propulsion system. The satellite has dimensions of approximately 10x10x34 cm, a mass of 3.99 kg, and will be deployed into a circular orbit at an altitude of about 400 km with an orbital inclination of 51.6 degrees from the International Space Station.

Since the launch of its first satellite in 2009, Tohoku University has continuously developed and operated Earth observation satellites and engineering demonstration satellites in the 50cm-class and CubeSat-class (up to 3U). The 50cm-class satellite launched into operation in 2021 enabled efficient operations through cloud-based management functions for both the satellite and ground stations, including automatic command generation. By 2022, up to eight operational satellites were simultaneously managed on a daily basis using three ground stations (Sendai, Hakodate, and Sweden). This paper presents the operational achievements to date and introduces the system that supports efficient satellite operations

Within the framework of the extra-dimensional Randall-Sundrum set-up, we investigate the freeze-in production of Standard Model (SM) gauge-singlet scalar, fermionic, and massive vector dark matter (DM). Assuming that both the DM and SM fields reside on the IR brane and interact solely through the graviton and radion, we demonstrate that the observed DM relic abundance measured by Planck can be achieved across a wide range of reheating temperatures, all while naturally addressing the hierarchy problem, satisfying constraints from collider, early Universe cosmology including $\Delta{N}_{\rm eff}$. We further show that the same set-up can accommodate TeV-scale leptogenesis capable of generating the observed baryon asymmetry of the Universe. Remarkably, we find that current graviton searches at the Large Hadron Collider (LHC) already impose strong constraints on the reheating temperature in this scenario.

The AMACA project (Astronomy education with a Multi-sensory, Accessible, and Circular Approach) develops multi-sensory activities for accessible education and engagement in astronomy. Despite promising innovations, existing resources are often poorly documented, designed for one-time events, expensive, and lack interdisciplinary collaboration, user testing, and broad dissemination. AMACA addresses these challenges by creating multi-sensory activities for education and outreach, with a particular focus on accessibility for people with sensory disabilities. A circular approach informs its educational structure: (1) a PhD course on multi-sensory astronomy outreach develops hands-on activities with the support of astronomers, psychologists, and organizations for the visually impaired and the deaf; (2) PhD candidates teach High School (HS) students how to deliver the activities; (3) HS students lead the activities at the Astronomy Festival "The Universe in All Senses"; (4) HS students train teachers to implement the activities in their classrooms. AMACA also develops tools to guide project development and track participants' learning. Key findings show improved communication and accessibility awareness among PhD candidates, increased emotional engagement with astronomy among HS students, enhanced public engagement with research and accessibility awareness, and high teacher satisfaction with the flipped-roles, hands-on approach. Overall, AMACA enhances accessibility and engagement in astronomy education across audiences.

The trace anomaly of dense matter, $\Delta \equiv 1/3 - P/\varepsilon$, defined in terms of the ratio of pressure $P$ to energy density $\varepsilon$, quantifies deviations from conformal symmetry and plays a central role in both the hydrodynamic response and gravitational equilibrium. While $\Delta(\varepsilon)$ has recently been inferred from neutron star observations, we report the first Bayesian extraction of the trace anomaly from collective flow observables in intermediate-energy heavy-ion collisions. By employing transport-model simulations that explicitly decouple the cold-matter mean-field potential from thermal effects, we directly constrain the cold dense-matter equation of state (EOS). Remarkably, the trace anomaly inferred from laboratory flow data agrees quantitatively, within $68\%$ credible intervals, with independent astrophysical posterior bands. This nontrivial agreement demonstrates that heavy-ion collisions and neutron star observations probe the same universal macroscopic properties of dense matter, establishing the trace anomaly as a composition-insensitive descriptor of dense matter across widely different physical environments.

Within the standard framework in which electroweak sphaleron processes relate lepton and baryon number, we derive an upper limit on the amplitude of a chiral gravitational wave background produced prior to the electroweak epoch. This bound is independent of the production time of chiral GWs for superhorizon modes, while it becomes sensitive to the production time for subhorizon modes. For sufficiently high reheating temperatures, the bound becomes significantly more stringent than the conventional big bang nucleosynthesis constraints at frequencies above the MHz scale, thereby providing a powerful and \emph{model-independent} probe of parity-violating physics in the early Universe.

This paper presents AAFIYA (Antenna Analysis in Frequency-domain for Impedance and Yield Assessment), a modular Python toolkit for automated characterization of radio-frequency antennas using measurement and simulation data. The toolkit provides a unified workflow for processing S-parameters, impedance, realized gain, beam patterns, polarization metrics, and calibration-based yield estimation, with support for standard Touchstone files and beam pattern data. AAFIYA is validated using measurements from an electromagnetic anechoic chamber involving Log Periodic Dipole Array (LPDA) reference antennas and Askaryan Radio Array (ARA) Bottom Vertically Polarized antennas over 100-850 MHz. Extracted metrics, including impedance matching, realized gain patterns, vector effective lengths, and cross-polarization ratio, are compared against full-wave simulations from HFSS and WIPL-D, showing good agreement across frequency and angle. The results demonstrate that AAFIYA enables accurate, reproducible, and publication-ready antenna analysis, and provides a flexible foundation for future extensions, including automated optimization and data-driven antenna design.

Black holes contain more entropy than any other component of the observable universe. Gravitational-wave observations from LIGO and Virgo have shown evidence of a previously unknown black hole mass range, which provides new information to update the entropy budget. Increases in entropy due to binary black hole mergers, as implied in the second law of thermodynamics, should also be added to the budget. In this study, we update the cosmological entropy budget for black holes in the stellar to lite-intermediate-mass range $(5-300~M_\odot)$, originating from either supernovae or binary mergers, by utilizing a suite of population synthesis models and phenomenological fits derived from numerical relativity. We report three new insights: Firstly, the cumulative entropy from merging black holes surpasses the total entropy from cosmic microwave background photons around the onset of the Over-massive Black Hole Galaxy phase at $z\sim 12$, suggesting that mergers played a more significant role in shaping the thermodynamic state of the early universe than relic radiation. Secondly, if primordial black holes constitute a nonzero fraction of dark matter, their early binary mergers establish an ``entropy floor" in the Dark Ages and can dominate the cumulative merger-generated entropy history even for small abundances. Thirdly, by computing the cosmological density parameters, we highlight the thermodynamic asymmetry in black hole mergers, where the production of gravitational-wave energy is inefficient compared to the immense generation of Bekenstein-Hawking entropy.

Simulation-based inference provides a powerful framework for Bayesian inference when the likelihood is analytically intractable or computationally prohibitive. By leveraging machine-learning techniques and neural density estimators, it enables flexible likelihood or posterior modeling directly from simulations. We introduce Residual Neural Likelihood Estimation (RNLE), a modification of Neural Likelihood Estimation (NLE) that learns the likelihood of non-Gaussian noise in gravitational-wave detector data. Exploiting the additive structure of the signal and noise generation processes, RNLE directly models the noise distribution, substantially reducing the number of simulations required for accurate parameter estimation and improving robustness to realistic noise artifacts. The performance of RNLE is demonstrated using a toy model, simulated gravitational-wave signals, and real detector noise from ground based interferometers. Even in the presence of loud non-Gaussian transients, glitches, we show that RNLE can achieve reliable parameter recovery when trained on appropriately constructed datasets. We further assess the stability of the method by quantifying the variability introduced by retraining the conditional density estimator on statistically identical datasets with different optimization seeds, referred to as training noise. This variability can be mitigated through an ensemble approach that combines multiple RNLE models using evidence-based weighting. An implementation of RNLE is publicly available in the sbilby package, enabling its deployment within gravitational-wave astronomy and a broad range of scientific applications requiring flexible, simulation-based likelihood estimation.

In recent years, significant progress has been made in the velocity-moment-based quasi-linear (QL) theory of waves and instabilities in plasmas with nonequilibrium velocity distributions (VDs) of the Kappa (or $\kappa$) type. However, the temporal variation of the parameter $\kappa$, which quantifies the presence of suprathermal particles, is not fully captured by such a QL analysis, and typically $\kappa$ remains constant during plasma dynamics. We propose a new QL modeling that goes beyond the limits of a previous approach, realistically assuming that the quasithermal core cannot evolve independently of energetic suprathermals. The case study is done on the electron-cyclotron (EMEC) instability generated by anisotropic bi-Kappa electrons with $A=T_\perp/T_\parallel > 1$ ($\parallel, \perp$ denoting directions with respect to the background magnetic field). The parameter $\kappa$ self-consistently varies through the QL equation of kurtosis (fourth-order moment) coupled with temporal variations of the temperature components, relaxing the constraint on the independence of the low-energy (core) electrons and suprathermal high-energy tails of VDs. The results refine and extend previous approaches. A clear distinction is made between regimes that lead to a decrease or an increase in the $\kappa$ parameter with saturation of the instability. What predominates is a decrease in $\kappa$, i.e., an excess of suprathermalization, which energizes suprathermal electrons due to self-generated wave fluctuations. Additionally, we found that VDs can evolve toward a quasi-Maxwellian shape (as $\kappa$ increases) primarily in regimes with low beta and initial kappa values greater than five. Instability-driven relaxation only partially resolves temperature anisotropy in bi-Kappa electron VDs, as wave fluctuations generally act to further energize suprathermal electrons.

We construct posterior distributions of the equation of state (EoS) for matter beyond the inner crust of neutron stars by incorporating finite nuclei (FN) constraints within relativistic mean field models. These constraints are implemented in three complementary ways: (i) through theoretical bounds on the EoS, (ii) implicitly via nuclear matter parameters, and (iii) explicitly by enforcing consistency with experimental binding energies and charge radii of selected nuclei. The resulting low-density nucleonic EoSs are subsequently matched to a model-agnostic speed-of-sound parametrization, constrained by astrophysical observations, including NICER mass-radius measurements, tidal deformability limits from GW170817, and lower bounds on the maximum neutron-star mass inferred from radio pulsar observations. We find that the admissible range of the transition density is strongly sensitive to the choice of the low-density EoS. In particular, the inclusion of explicit FN constraints significantly reduces the allowed parameter space of the nucleonic EoS at low densities, narrowing the transition-density range by nearly a factor of two. Consequently, neutron-star properties inferred from EoSs with explicit FN constraints differ substantially, with especially pronounced effects for low-mass neutron stars and their correlations with nuclear matter parameters. A quantitative comparison, using metrics based on Mahalanobis distance, shows consistency of the explicit constraints with PSRs J0740+6620, J0030+0451, and J0437-4715, but suggest a possible tension with PSR J0614-3329. These findings underscore the critical importance of a consistent treatment of finite-nuclei properties for reliably inferring the behavior of high-density matter and the presence of possible phase transitions from astrophysical observations.