Papers for Thursday, Mar 12 2026

Two independent electron temperature diagnostics applied to the quiet solar corona yield systematically different results. Radio brightness temperatures from the Nancay Radioheliograph indicate T_e ~ 0.6 MK, while hydrostatic scale-height modeling of the same plasma requires T_e ~ 1.5 MK (Mercier & Chambe 2015). Both diagnostics probe electrons; they disagree by a factor of R = 2.4 +/- 0.3. This discrepancy persists across an eight-year dataset spanning solar minimum and is consistent with LOFAR observations at lower frequencies (Vocks et al. 2018). We consider the propagation alternative (turbulent scattering of radio emission), which operates in the correct direction to suppress the apparent brightness temperature, but the ratio R is invariant over the solar cycle despite expected variations in turbulence levels. We propose that the residual, cycle-invariant discrepancy reflects non-Maxwellian electron velocity distributions. Radio bremsstrahlung samples the distribution core, while ionization rates and scale heights are dominated by the suprathermal tail. For kappa distributions, the predicted ratio is kappa/(kappa - 3/2); the observed R = 2.4 implies kappa ~ 2-3. This is consistent with spectroscopic measurements in active regions but in tension with perturbative theoretical predictions of kappa ~ 10-25. We make falsifiable predictions: Active Region cores should show a collapsed ratio (R <= 1.5) as collisionality restores thermal equilibrium. Applying fluid transport equations (Spitzer-Harm conductivity) to plasmas with kappa ~ 2-3 is physically invalid, but we do not compute the resulting heat flux, which remains an open problem.

Cryogenic sub-Hz cROss torsion-bar detector with quantum NOn-demolition Speed meter(CHRONOS) is a proposed next-generation ground-based gravitational-wave observatory designed to explore the sub-Hz frequency band with unprecedented sensitivity. Utilizing a cryogenic torsion-bar interferometric configuration with quantum non-demolition speed-meter readout, CHRONOS targets a frequency window that bridges space-based missions and current high-frequency ground-based detectors, opening a new frontier in gravitational-wave astronomy. The observatory will enable long-duration tracking of compact binary inspirals well before merger, significantly improving source localization, parameter estimation, and tests of general relativity. In addition to transient signals, CHRONOS is optimized to probe the stochastic gravitational-wave background (SGWB) in the sub-Hz regime, providing powerful constraints on primordial gravitational waves, inflationary tensor spectra with red or blue tilts, first-order phase transitions, cosmic strings, and other relics of high-energy physics. By connecting gravitational-wave measurements across cosmological frequency scales-from cosmic microwave background polarization to pulsar timing arrays and high-frequency interferometers-CHRONOS will contribute to a coherent reconstruction of the gravitational-wave spectrum over more than twenty orders of magnitude. Crossing critical sensitivity thresholds in the sub-Hz band, CHRONOS will establish a new pillar of gravitational-wave astronomy and cosmology, enabling transformative advances in astrophysics and fundamental physics.

Studying exoplanet flow and variability requires solving atmospheric dynamics equations accurately. Here we use the shallow-water equations to evaluate and employ Dedalus3, a spectral method-based software package for solving differential equations. A well-known jet instability test is used for the evaluation; then, the package is used to investigate the nonlinear evolution of observed, Jupiter's zonal (east--west) jets; finally, the package is used to compare hot-Jupiter flows with different initial conditions. Our results indicate that Dedalus3 can be a useful tool for investigating planetary flow dynamics, but careful testing and execution are necessary for each problem.

We present new observations from JWST NIRCam that reveal a striking kpc-wide cavity in the stellar distribution of the central galaxy in the cluster Abell402. Supporting data from HST allow us to rule out extinction due to dust as an explanation and, instead, suggest that this is a localized depression in the stellar density field corresponding to ~2x10^9 Msun in missing stars within a volume of 0.5kpc^3. On larger scales, both the JWST and HST data show evidence for a 2.2kpc flattened core in the stellar distribution (on which the smaller-scale cavity is superimposed), which implies the presence of a central ultra-massive black hole with M_BH = 6 +/- 4 x10^10 Msun. We report evidence for a mid-IR-bright point source at one edge of the cavity, suggesting that this black hole is actively accreting. MUSE spectroscopy reveal that this source is a LINER AGN and that there is a second candidate AGN on the opposite side of the cavity with a relative velocity of 370km/s -- if real, this implies the presence of a kpc-separation dual AGN with a total binary mass of 6 +/- 2 x10^10 Msun, which would make this the most massive binary black hole system discovered to date. We propose that this unique stellar cavity is the result of a short-lived dynamical interaction between at least one supermassive black hole and the background stellar density field, caused either by three-body scattering during binary hardening or the induction of a dipole instability in the stellar density field.

Dark matter (DM) could be subject to non-gravitational self-interactions which is relevant to resolve potential problems of cold DM on small scales. Their impact on astrophysical objects such as galaxies and galaxy clusters allows for constraining the strength of this scattering and eventually further properties of the cross-section. To model self-interacting dark matter (SIDM), N-body simulations are a crucial tool widely employed by the SIDM community. In this paper, we describe the SIDM implementation in the cosmological hydrodynamical N-body code OpenGadget3 and release it to the public. It is capable of simulating elastic scattering for various differential cross-sections, including strongly anisotropic cross-sections. Beyond single-species models, the code also allows simulating a two-species model with cross-species interactions. In addition to describing the numerical schemes for modelling various flavours of SIDM, we discuss the technical challenges of implementing them. Moreover, we demonstrate through several test problems that OpenGadget3 can accurately simulate DM self-interactions. Furthermore, we assess the performance of the code and provide scaling tests. Lastly, we highlight remaining challenges in the context of SIDM and describe directions for improving the current state of the art.

We present the results of a detailed kinematic analysis of a significant fraction of the known population of Galactic star clusters aimed at constraining the physical mechanisms driving the onset and evolution of cluster rotation. Our study reveals for the very first time the presence of rotation in clusters at any age, with about $25\%-30\%$ of systems in the sample showing significant evidence of rotation. This result increases by a factor of $\sim5$ the number of clusters identified as rotators so far and it finally enables an observational reading of cluster rotation as a function of time. Young ($<500$ Myr) clusters show a larger range of rotation velocities than older systems. In addition, at young ages we observe a significantly larger fraction ($50\%-60\%$) of rotating systems than at older ones ($\sim 15\%$). These purely empirical results are compatible with rotation being imprinted during the very early stages of cluster formation and early evolution and then being progressively erased by the long-term effects of dynamical evolution. For the sub-sample of clusters for which we were able to perform a full 3D analysis, we calculated the angle between the internal rotation axis and that of the cluster orbital motion. Interestingly, while for clusters with an age smaller than their orbital period we observe similar fractions of prograde and retrograde systems, more evolved clusters appear to be preferentially prograde. We argue that such a behavior is in qualitative agreement with the expectations for the evolution of systems in which primordial rotation was imprinted by the parent molecular cloud and/or by the following hierarchical cluster assembly processes, and in which internal cluster dynamics and interactions with the Galactic field have induced a torque-driven alignment between cluster rotation and orbital motion.

We report on a pair of X-ray Multi-Mirror Mission (XMM-Newton) observations of the Supergiant X-ray binary 4U 1909+07, which were performed on 2021 October 3 and 2021 October 8, respectively. We measure the neutron star rotation period in each observation to be $\sim$602.62 s. This continues a long spin-up trend that has persisted since 2001 where the neutron star spin period was found to be $\sim$604.66 s. In our timing analysis, we observe strong variations in the amplitude of the 1--10 keV pulse profile as a function of time, and for the first time we find a low flux interval extending for a single pulse period in which pulsations are no longer detected. We interpret this low flux interval as a pulse dropout similar to those observed in Vela X-1 and GX 301-2, which were each explained by a low-density cavity in the wind driving the propeller effect. In our time-resolved spectral analysis, we observed the spectral continuum, which can be described as an absorbed power law modified by a high-energy cutoff, to significantly soften during the pulse-dropout phase. No evidence of an increasing absorption column density was found. The observed softening in 4U 1909+07 also supports an interpretation that the observed pulse dropout may be driven by the propeller effect, but the quasi-spherical settling accretion regime cannot be ruled out.

We investigate the cosmological constraining power of combined weak galaxy lensing and galaxy clustering probes, i.e. $3\times2$-point analyses, assuming flexible models for redshift uncertainty, and Lagrangian perturbation theory and hybrid effective field theory models for galaxy intrinsic alignments, galaxy bias and baryonic physics. In this context, we provide a detailed accounting of the limiting systematics on $3\times2$-point analyses. Our main finding is that in the presence of current levels of uncertainty on baryonic physics, the information content of weak lensing analyses saturates on quasi-linear scales, allowing the use of source galaxy samples that are significantly less dense, e.g. with number densities of $5\rm \, arcmin^{-2}$, without sacrificing constraining power, provided that redshift distributions can be calibrated at the $\sigma(\langle z\rangle)=0.005$ level. We show that for sufficiently narrow lens and source redshift distributions, intrinsic alignment contributions can be largely self-calibrated, though sufficient flexibility must be given to the redshift and scale dependence of this signal. The near optimality of such relatively sparse source galaxy samples opens the possibility to directly calibrate the redshift distributions and intrinsic alignment contamination of such a sample using a spectroscopic instrument like DESI, thus mitigating the dominant systematics in weak lensing analyses.

The LMC's Globular Clusters (GCs) bring a novel opportunity to understand the LMC's assembly history and dark matter (DM) properties, provided the kinematically outlying GCs can be reliably identified. However, traditional diagnostics like the Energy-Angular Momentum space fail because of large uncertainties on the GC velocities. In this work, we develop a new, robust statistical framework for identifying kinematically outlying LMC GCs, by using their Gaia-DR3 Proper Motions (PMs) combined with previous Line-of-Sight (LoS) velocity measurements. We use the difference between a GC's velocity vector and the average velocity vector of the surrounding red clump stars as a metric for quantifying a GC's kinematic peculiarity. We account for both the velocity measurement uncertainties and the LMC's intrinsic velocity dispersion. We find 5 LMC GCs to be kinematically outlying based on PM differences alone, and additional 6 GCs if LoS velocity information is also used. Majority of the GCs with outlying PMs are clustered at a distance of 3-4 kpc from the LMC center. The inclusion of outlying LMC GCs introduces a bias of upto 30% in the LMC's enclosed mass estimates derived using GCs as dynamical tracers; caution must be exercised in choosing the GC sample for precisely determining the LMC's DM content. We discuss the possibility that the kinematically outlying LMC GCs may have been accreted from external galaxies, and motivate future spectroscopic follow-up of the GC population to better understand the assembly history of massive satellite galaxies of Milky Way like hosts.

We provide a complete study of the factors influencing gravitational-wave signal localization using pulsar timing arrays. We derive analytical expressions for the Cramér-Rao sky localization precision that delineate the impact of the angular proximity, $\xi$, between the pulsar and the gravitational wave source, and the precision, $\sigma_L$, with which pulsar distances are known. Interference between the Earth and pulsar terms creates rapid angular oscillations for sky-coordinate Fisher matrix elements that aids localization, which is complemented by more broadly varying antenna response gradient information. The relative importance of these factors depends on whether pulsar distances are known precisely [i.e., $\sigma_L\leq\lambda_\mathrm{GW}/(1-\cos\xi)$] or imprecisely, respectively. If the former, tightening pulsar distance precisions improves signal localization according to $\Delta\Omega_\mathrm{sky}\propto\sigma_L^2$ until the Earth-pulsar system reaches its diffraction limit. If the latter, localization precision is degraded, but more pulsars in close proximity to the source is the best means of improving. With $\alpha$ indexing pulsars, this scales as $\Delta\Omega_\mathrm{sky}~\propto~(\sum_\alpha \mathrm{SNR}_\alpha^2/\xi_\alpha^2)^{-1}$ in the small-angle limit of the unmarginalized Fisher matrix, and we derive the analytic generalization to any angle between a pulsar and the source. Finally, we study a scenario where pulsar-term phases are treated as nuisance variables that are unconnected to binary or PTA properties. This phase-decoupled scenario, which is how all PTA continuous wave searches are currently conducted, delivers localization performance similar to the antenna-response--driven case, and does not exhibit significant improvement as pulsar distance precisions are tightened.

We study the transport of baryons in magnetically arrested accretion flows and relativistic jets using general relativistic magnetohydrodynamic simulations that incorporate a passive Eulerian tracer. The tracer allows us to reconstruct a proxy for the physical baryon density supplied by the accretion disk while excluding the mass injected numerically to maintain stability in highly magnetized, low-density regions. Applying this method to axisymmetric black hole simulations with varying spin, we show that baryon loading of the jet is intrinsically episodic and regulated by magnetic flux eruption cycles occurring in the inner accretion flow. Each eruption evacuates baryons from the innermost equatorial region, drives reconnection in extended current sheets, and expels moderately magnetized disk material along the funnel wall, establishing a recurrent mass-loading channel. In spinning black holes, shear-driven waves along the jet boundary further enhance baryon entrainment, whereas this mechanism is suppressed in the non-spinning case. For parameters representative of the black hole accretion flow in M87, we map the global structure and time evolution of the Goldreich-Julian screening boundary, defined as the surface separating regions where the plasma density is sufficient to supply the charges required to screen electric fields parallel to the magnetic field from regions that are charge starved. For spinning black holes, we find that the electromagnetic power of the jet is predominantly carried by baryon-poor plasma, with extended time intervals of charge starvation. Our results provide a framework for diagnosing jet composition, charge starvation, and reconnection-driven mass loading in magnetically arrested black hole systems, with direct implications for particle acceleration and non-thermal emission in low-luminosity accretion flows.

The tidal field of a black hole can turn a star into a gas stream whose orbit can precess, especially if the a black hole is rapidly spinning. In this work, we investigate the impact of precession on the light curves of tidal disruption events (TDE). To do so, we perform two-dimensional radiation-hydrodynamic simulations of the interaction of the TDE wind and luminosity with the precessed stream wrapped around the black hole. Our results show that in events with black holes of $\sim10^6~\text{M}_{\odot}$ and no orbit-spin inclination, the line of sight has little effect on the light curves, since the stream covers a small fraction of the solid angle as the precession is confined to the orbital plane. In the case of black holes of $\gtrsim10^7~\text{M}_{\odot}$ and high inclination ($i\sim90^{\circ}$), the light curve peaks can be delayed by $\sim$100 days due to presence of the precessed stream blocking the radiation in the early phase of the event. We also discuss our efforts to model self-consistently the hydrodynamic evolution of a tidal stellar stream on curved spacetimes by the presence of a massive black hole.

Photometric selection of luminous high-redshift ($z\gtrsim4$) quasars is plagued by contamination from numerous low-mass Galactic stars, reddened lower-redshift quasars, as well as compact luminous red galaxies. Confirmation of these rare objects thus requires extensive spectroscopic campaigns on 4 and 8-meter-class telescopes with relatively low success rates. Here we demonstrate the utility of SPHEREx spectrophotometric survey data for quasar confirmation with no ground-based follow-up required, "from the couch," applied to candidates from a purposefully simplistic photometric and astrometric Gaia+WISE selection down to low Galactic latitudes ($|b|\geq8^\circ$). Primarily from the detection of their strong broad H$\alpha$ emission lines, we discover 87 new luminous $4.0 < z < 5.7$ quasars with median $M_\text{1450} = -27.5$, including 19 quasars at $z>5$, and recover 219 previously published quasars at $z>4$. We validate our SPHEREx selection with a 100% confirmation rate in ground-based spectroscopic follow-up of 29 of our new $z>4$ quasars, including 11 unpublished archival spectra. We also discover 203 additional lower-redshift quasars at $0.3 < z < 4$, consisting primarily of relatively rare highly-reddened and strong broad-absorption-line objects that are likely missed by traditional quasar surveys. Finally, we show that the Ly$\alpha$ absorption breaks and H$\alpha$ lines of luminous quasars are already detectable at redshifts $5.7\lesssim z\lesssim6.5$ after the completion of only the first of four all-sky surveys to be performed by SPHEREx during its planned two-year mission.

Accreting white dwarfs (AWDs) are among the best natural laboratories for understanding disk accretion. Their proximity, brightness, and purely classical nature make them ideal systems in which to probe the fundamental physics that governs the transport of angular momentum, the generation of outflows, and the coupling between disks, magnetospheres, and accretors. Yet despite decades of study, many critical questions remain unresolved. In this ``unreview'', we therefore focus not on what is known, but on what is unknown. What drives viscosity and sustains accretion in largely neutral disks? How are powerful winds launched, and how do they feed back on the disk and binary evolution? Why do so many systems show persistent retrograde precession, and what drives bursts in magnetic AWDs? By identifying these open problems -- and suggesting ways to resolve them -- we aim to motivate new observational, numerical, and theoretical efforts that will advance our understanding of accretion physics across all mass scales, from white dwarfs to black holes.

We investigate primordial magnetogenesis by assuming the gauge field is prepared in a thermal state during inflation rather than the standard Bunch-Davies vacuum. The temperature $\mathcal{T}$ introduces a physical scale that breaks conformal invariance at the level of the state while preserving the standard Maxwell action. This modification results in a {\it dissipative boost} that alters the magnetic energy density scaling from $a^{-4}$ to $a^{-3}$, resulting in a present-day magnetic field $B_0$ enhancement that can potentially range from about $10^{8}$ to $10^{16}$ on galactic scales. While this toy model alone does not satisfy observational lower bounds, it demonstrates that thermal initial conditions can significantly mitigate the conformal obstruction. Our results suggest that embedding this mechanism within a fully dynamical warm inflation framework, where dissipation continuously maintains the thermal bath, provides a highly promising path towards successfully realizing a minimal model of inflationary magnetogenesis without the need to invoke non-minimal couplings, anomalous background dynamics or nonlinear extensions of electrodynamics.

We investigate the origin of the observed X-ray weakness in high z LRDs and other JWST-selected broad line AGN by comparing their X-ray and optical properties with those of a diverse sample of low z AGN, including super-Eddington accreting massive black holes (SEAMBHs), NLS1s, and type I AGN from large surveys. We examine the relations between X-ray luminosity, broad H{\alpha} line luminosity, Eddington ratio, bolometric luminosity and X-ray-to-bolometric luminosity correction, and explore whether high z sources may represent analogues of local highly accreting systems. While a few LRDs and JWST-selected AGN are consistent with the SEAMBHs population in the $L_x/L_{H{\alpha}}$ versus $\lambda_{Edd}$ plane, most lie below it, suggesting either more extreme accretion conditions, suppressed coronal emission or heavy obscuration. We identify an anti-correlation between $L_x/L_{H\alpha}$ and $\lambda_{Edd}$ in the low z, high accreting subsample, consistent with theoretical expectations of slim-disc accretion. We further show that, for SEAMBHs, $H\alpha$-based bolometric luminosities underestimate SED-based values even after dust correction. We find that SEAMBHs, LRDs, and JWST-selected AGN occupy a similar high-$\kappa_{bol,x}$ regime, indicating that the relative deficit of X-ray emission compared to the bolometric output could potentially support the view that suppression of the hot corona emission is a common feature of highly accreting systems across cosmic time. Our results are consistent with the idea that the observed X-ray weakness of LRDs and JWST-selected AGN may be linked to the physics of highly accreting SMBHs. Moreover, observational limitations at high z, including instrumental sensitivity and the steep X-ray spectra expected for highly accreting systems, likely further suppress the detected X-ray signal.

In this paper I review recent results on high-energy neutrino astronomy and what they can reveal about some of the most extreme cosmic accelerators. I discuss recent measurements of the diffuse TeV-PeV cosmic neutrino spectrum by the IceCube observatory and the current flux limits in the ultra-high-energy regime, contextualizing the recent detection of an ultra-high-energy neutrino by the KM3NeT observatory. I review the recent emergence of a TeV signal from nearby Seyfert galaxies such as NGC 1068, the potential of $\gamma$-ray blazars as neutrino sources above the PeV regime, and the current status of tidal disruption events and other transient classes as possible neutrino sources. For each of these topics, I discuss ongoing developments in source models and their current limitations. I argue for the indispensable role of next-generation multi-messenger facilities, such as IceCube-Gen2, in solidifying current source associations, probing the ultra-high-energy regime, and resolving vast transient populations that remain unidentified with current statistics.

In globular clusters, hierarchical mergers are among the most promising pathways to forming massive black holes such as GW231123. A key factor determining whether a merger-remnant black hole will be retained in these environments and thus participate in subsequent hierarchical mergers is the recoil kick velocity. Analytic models for the recoil velocity are currently employed in nearly all population-synthesis frameworks. We instead use a state-of-the-art recoil-kick model gwModel_flow_prec developed from a combination of numerical-relativity and black-hole perturbation-theory data, together with data-driven techniques such as normalizing flows and the post-Newtonian structure of the kick. Employing both back-of-the-envelope estimates and detailed N-body as well as semi-analytical cluster simulations, we show that gwModel_flow_prec leads to a noticeable increase in the retention probability of hierarchical-merger remnants compared to the previously used analytic model and changes the mass and spin distribution of the black holes formed through hierarchical mergers. Additionally, we discuss the implications of our results in the context of massive binaries such as GW231123.

We tested the effectiveness of radio selection to discover heavily obscured AGNs, particularly at high-z, and we measured their abundance for the first time from a radio perspective. We consider the radio sources detected in the J1030 field, which is one of the fields with the deepest combination of 1.4 GHz radio and X-ray observations. We defined a radio excess parameter as the ratio between the star formation rate (SFR) that would correspond to the observed radio luminosity and the one directly derived from the spectral energy distribution (SED) fitting, $\rm REX=SFR_{1.4GHz}/SFR^{corr}_{SED}$. We then select as radio excess AGN those sources with $\rm REX>8.5$, corresponding to a $3\sigma$ excess above the median value. In this way, we find 145 radio-excess sources falling into the \textit{Chandra} X-ray image footprint but without X-ray detection. From the deep X-ray upper limits, we estimated a lower limit to the obscuration of each radio-excess AGN, finding on average $\log (N_H/\rm{cm^{-2}})>23.7$. A CTK AGN scenario is also supported by the results of the X-ray stacking analysis performed on sources at $z>1.5$, which revealed X-ray luminosities and hardness ratios compatible with very highly obscured AGN. Finally, we computed the number density of these radio-selected CTK AGN. While at $z\sim 2$ the radio number density agrees well with the CTK AGN predictions of different population synthesis models, at $z\sim3$ the radio selection returns a CTK AGN number density $\sim 2-3$ times larger than what is predicted by the X-ray models and observations. This result supports the effectiveness of radio emission in selecting the most obscured sources, unraveling a population of AGN potentially missed by X-rays surveys at $z>3$, paving the way to a synergistic use of the future radio and X-ray facilities such as the \textit{SKAO} and \textit{NewAthena}.

Accurate modelling of redshift-space distortions (RSD) in galaxy clustering is essential for extracting cosmological information from current and forthcoming large-scale structure surveys. While perturbation theory is reliable on large scales, much of the constraining power lies at intermediate and small separations, where nonlinear dynamics within and between dark matter haloes dominate. We present a halo streaming model for nonlinear galaxy clustering in redshift space that is accurate and physically interpretable. Our framework combines the streaming model for RSD with a halo-model decomposition of the galaxy clustering into central/satellite and one-/two-halo contributions. We build dedicated emulators for the key physical ingredients, trained on a suite of $N$-body simulations: halo mass functions, real-space halo two-point correlation functions, and pairwise velocity moments. By emulating these modular building blocks rather than the final redshift-space observable, this approach preserves physical transparency, enables targeted optimisation for each ingredient, and remains flexible to changes in tracer populations and galaxy-halo connection models. The resulting halo streaming model reproduces the simulated nonlinear anisotropic clustering signal down to highly nonlinear scales, while achieving the computational efficiency required for cosmological parameter inference. This framework is designed to support full-shape RSD analyses for surveys such as DESI and \textit{Euclid}, facilitating precision measurements of structure growth and tests of gravity. All codes and trained emulators are publicly available in the \href{this https URL}{\texttt{freyja}} repository.

We present a synopsis of the project to establish thirty-two new faint ($ 16.5 \leq V \leq 19.8 $) DA white dwarfs as spectrophotometric standards distributed over the whole sky. Our results validate the use of fully radiative pure hydrogen model fluxes for hot DA white dwarfs to predict the observed broadband fluxes from near ultraviolet through the near infrared to accuracies of a few parts per thousand. After fitting the line of sight reddenings simultaneously with the model spectral energy distributions of these stars against spectroscopic and multi-band photometric observations, we have shown that residuals have an rms of typically 0.4 percent. This indicates that the complications from interstellar dust extinction have been adequately mitigated. Our stars supplement the three brighter DA white dwarfs that define the flux scale of CALSPEC. The consequent photometric accuracy, their all sky coverage, and their brightness range that matches the dynamic range of large telescopes, constitutes an unprecedented ensemble of standard stars for both ground as well as space based use. This paper targets readers who may wish to use these as standard stars, and provides for them the essential content to understand their strengths and limitations, without traversing the technical details of analysis that are already captured in a series of papers since 2016. The narrative here describes the motivation, justification, and evolution of the analysis methods; the input data that constrain the modeling; as well as the stability of our results in the face of future improvements in models.

Radio bridges are diffuse synchrotron emission observed between merging galaxy clusters. Recent radio observations have reported both detections and non-detections of radio bridges between clusters. The detections imply the presence of cosmic rays (CRs) and magnetic fields permeating the cosmic web that produce synchrotron emission observable with current facilities, whereas the non-detections suggest that specific physical conditions are required for their formation. We study the CR reacceleration by solenoidal turbulence in the filament connecting two massive clusters at an early stage of the merger. Our aim is to test whether this mechanism can generate diffuse emission in the inter-cluster region. We perform a cosmological magneto-hydrodynamical (MHD) simulation using the Enzo code. We improved a run-time Lagrangian tracer method implemented in Enzo, and follow the trajectories of baryonic matter using $N=\mathcal{O}(10^7)$ tracer particles. In post-processing, we conduct a parallel computation of the Fokker-Planck (FP) equation for all tracers, with cooling and reacceleration efficiencies evaluated from the local quantities recorded along each tracer trajectory. Our simulation generate a Mpc-sized radio bridge in the early stage of the cluster merger. Within a reasonable parameter range, the reacceleration model produces a broad variety of spectra. In our fiducial model, the simulated bridge matches several properties of the one found between Abell 399 and Abell 401, such as its spectral shape, intensity profile, and pixel-by-pixel correlation between radio and X-ray intensities. The inter-cluster region is filled with turbulence induced by infalling mass clumps and subsequently amplified by the approaching motion of the clusters. The CR reacceleration by the turbulence is a viable mechanism to power a Mpc-sized synchrotron emission observed as radio bridges.

Cygnus~X-3 is a Galactic X-ray binary with a 4.8-hr orbital period operating in the ultraluminous regime. Although the system is viewed at relatively low inclination ($i\approx28^\circ$), it exhibits a deep orbital modulation. Recent IXPE observations show strong linear polarization orthogonal to the radio jet, indicating that the X-ray emission is dominated by reflection from the inner walls of a supercritical outflow funnel. We propose a Hybrid Roche-lobe overflow (RLOF) scenario in which a massive Wolf-Rayet donor effectively fills its Roche lobe with a focused wind driving a super-Eddington accretion stream. Using a numerical synthesis of the folded light curve, we show that the modulation is reproduced when the central funnel is periodically occulted by a vertically extended, shock-heated Turbulent Wall formed by stream impact on the outer disk rim. This produces a phase lag ($\Delta\phi\approx0.11$) between X-ray minimum and binary conjunction, with extended attenuation by the WR wind defining a broader Suppression Region. This geometry explains the enhanced iron-line equivalent width during X-ray minimum via a coronagraphic effect. The large radial-velocity amplitude of FeXXVI measured by XRISM ($K_{\rm obs}\approx430$ km s$^{-1}$) and its zero-crossing at $\phi_X=0.0$ arise naturally in the stream-impact region rather than from orbital motion of the compact object. Finally, we show that the observed secular orbital expansion ($\dot P>0$) follows directly from highly non-conservative mass transfer with inner-disk mass loss, indicating that Cygnus~X-3 is a stable, long-lived system in a supercritical accretion regime.

Tidal disruption events (TDEs) occur when stars are destroyed by supermassive black holes and are among the brightest nuclear transients. It has been thought that strong relativistic effects rapidly dissipate orbital energy and produce prompt disk formation when the stellar pericenter is smaller than $\sim 10$ gravitational radii. Using a general relativistic hydrodynamic simulation of a strongly relativistic TDE around a $10^{6}\,M_{\odot}$ black hole, we find instead that the overall evolution is similar to weakly relativistic TDEs: the debris remains highly eccentric, with most of the returned mass residing near the orbital apocenter ($\sim 250\times$ the initial pericenter distance), and shocks, rather than accretion, power the event. The simulation starts from the initial stellar approach and follows the debris evolution up to $35\,\text{days}$ after the peak mass-return time ($\simeq$ $23\,\text{days}$). Although early shocks driven by strong relativistic apsidal precession and pericenter nozzle compression dissipate orbital energy efficiently, they last only about a week ($\sim 0.3$ of the peak mass-return time). Stream self-interactions increase the incoming stream's angular momentum, thereby expanding its pericenter distance, weakening precession and shocks, and reducing dissipation. These results, along with previous work on weakly relativistic TDEs, suggest that circularization may be intrinsically slow regardless of the strength of relativistic effects, and the flow remains highly eccentric and extended during the peak of optical/UV luminosity.

M-dwarf stars are the most favorable hosts for detecting small transiting planets, yet mid-to-late M-dwarfs that acquired sufficient TESS multi-sector coverage only through recent Cycle 6+ observations represent a newly accessible discovery space. This paper presents a systematic transit survey of 121 M3-M6 dwarfs (Teff = 2700-3400 K) selected as "newly enabled" targets -- stars with <=2 archival TESS sectors that only recently crossed the multi-sector detection threshold, covering P = 0.5-100 d. The sample was selected from 498,312 TIC M-dwarfs via a 9-step funnel. The pipeline combines TLS with a signal validation cascade, TRICERATOPS vetting, Gaia DR3 verification, and three empirical signal reliability tests. Pipeline validation achieved 100% recovery (16/16 planets) on 10 known systems with zero false positives. The survey identifies 20 transit-like signals across 16 systems, none with prior TOI designations. The reliability framework classifies 2 as Tier 1 (High Robustness), 7 as Tier 2 (Moderate), and 10 as Tier 3 (Noise-Susceptible); one monotransit is excluded. No candidate SDE significantly exceeds its host star's noise floor. The global false alarm rate is 17.4% (21/121; Wilson 95% CI: [11.6%, 25.1%]).The 2 Tier 1 candidates are priorities for RV confirmation. The 10 Tier 3 candidates require additional TESS sectors to establish signal persistence; 9 systems need high-resolution imaging.

Jets are common eruptive phenomena in the solar atmosphere which may occur repeatedly. Many studies of their fine dynamics have been conducted. However, the fine dynamics of persistent interactions among various features that drive recurrent coronal jets have not been studied in detail. In this paper, we use observations from the Solar Orbiter to report persistent interactions between rising loops/mini-filaments and a fan-spine-like structure, which produced more than 22 ejections. Many loops and mini-filaments under the fan-spine-like structure rose with speeds of $8 - 58\,\mathrm{km~s^{-1}}$ and an average of $27\,\mathrm{km~s^{-1}}$. These rising loops and mini-filaments interacted with the fan-spine-like structure successively, producing ejections with speeds ranging from 25 to 186\,\rm km\,s$^{-1}$ and an average at $80\,\rm km~s^{-1}$. We observed the fine dynamics of the drivers of these recurrent jets in detail, including partial eruption of mini-filaments, formation of a new mini-filament by contraction of remaining threads from the partially-erupted mini-filament, and interaction between rising loops (or mini-filaments) and the fan-spine-like structure. Brightenings appeared near the footpoint of these rising structures, followed by the formation of current sheets. Some arcades at the outflow region contracted with speeds of around $10\,\rm km~s^{-1}$, and the outflow region moved at around $8\,\rm km~s^{-1}$ toward the opposite direction. Bright blobs were observed in the current sheets, and they propagated at speeds averaging at $21\,\rm km\,s^{-1}$ and had an average width of 296\,km. We emphasize the vital roles of persistent rising loops and/or mini-filaments in producing recurrent jets by interacting with the fan-like structure, and show their detailed dynamics with unprecedentedly-high-resolution observations.

We propose a variable ADAF disk model for X-ray binary systems. In this model, the accretion flow consists of an outer thin disk and an inner optically thick advection-dominated accretion flow (ADAF) torus. The size of the turbulent ADAF is variable. A complete cycle of ADAF contraction, transition to a thin disk, and subsequent re-expansion corresponds to the rapid rise, peak, and decay phases observed in the X-ray outbursts of black hole binaries. This cycle also tracks the canonical evolution through the low-hard, high-soft, and back to the low-hard state in the hardness-intensity diagram. The model unifies the presence of near-ISCO Fe emission lines with the truncated disk paradigm, as observed in the black hole system GX 339-4. It explains the 35-day period in the neutron star system Her X-1 more effectively through variable ADAF sizes than through a precessing disk. This variable ADAF framework may be extended to explain similar phenomena in active galactic nuclei.

Galaxy clusters are the largest gravitationally bound systems, and they continue their growth through mergers in a hierarchical {\Lambda}CDM Universe. Therefore, we can describe the merger stage of a cluster as the dynamical state of clusters. Previous studies have investigated this phenomenon, but several limitations remain, including reliance on dichotomous classifications, constraints on the number of indicators used, absence of reliability, and incompatibility of methods between observation and simulation studies. To overcome this, we developed an enhanced and observation-applicable cluster dynamical state classification method using the Bayesian classifier with the class-conditional Gaussian mixture distribution model using the N-cluster Run simulation data. The Bayesian classifier was designed for two merger stages (merger and relaxed) as well as three merger stages (recent merger, ancient merger, and relaxed) to provide a more detailed interpretation of the merger processes. In the results, using a larger number of indicators yields better results, with their order of importance being: magnitude difference, center offset, sparsity, Kuiper V statistic, and mirror asymmetry. Additionally, our analyses show that a projected classifier (built on the 6D space, but evaluated on lower dimensional projections) consistently produces better outcomes than non-projected classifiers (i.e., classifiers built directly on the corresponding low dimensional spaces), which means limited observation data can be used to classify with enhanced performance. Furthermore, the new classification method outperforms our previous research. This new method can suggest a way of overcoming previous limitations and provides new insights by providing the reliability of dynamical state classification results.

Planetary atmospheres cannot remain hydrostatic at all altitudes because they approach finite density at infinite radius, implying infinite mass. Classical treatments address this in two directions: either retain a hydrostatic structure while allowing particles in the high-velocity tail to decouple and escape in a Jeans-type manner, or promote the gas to a continuum outflow to obtain a transonic Parker-type solution. The usual criterion compares the local mean free path to the sonic point radius. If the mean free path is shorter, the atmosphere is hydrostatic with an imposed Jeans escape flux; if it is longer, the gas is hydrodynamic with Jeans escape neglected. Here, we show that hydrogen-rich atmospheres do not separate cleanly into hydrodynamic and Jeans-escape regimes. At any radius, some particles still collide and behave as a fluid, while others have already experienced their last collision and move collisionlessly on ballistic trajectories. The relative importance of these two behaviors changes smoothly with radius rather than switching at a single boundary. The hydrodynamic channel accelerates and passes through a sonic point, whereas the collisionless channel decelerates under gravity and grows with altitude, removing mass and momentum from the collisional flow. As the collisionless component grows, the bulk flow speed reaches a maximum and then decelerates thereafter, producing profiles similar to Parker breeze solutions even though escape is carried by the collisionless channel. This two-channel framework provides a first step toward a self-consistent treatment that unifies hydrodynamics and kinetics in atmospheric loss models.

Ultraluminous X-ray sources represent extreme super-Eddington accretion regimes, and a subset is now known to host highly magnetized neutron stars. However, direct observational probes of their surface magnetic fields remain scarce. In this Letter, we report the detection of a narrow X-ray absorption feature at $3.29\pm0.02$ keV in the XMM$-$Newton/EPIC-pn spectrum of NGC 4656 ULX-1. The source exhibits a hard-ultraluminous state, while our timing analysis reveals a candidate pulsation at $\sim$0.9736 Hz, with a local significance of $5.5\sigma$ and a pulsed fraction of $\sim11\%$. The feature is robust against changes in continuum modeling and data-selection criteria, retaining a statistical significance of $\gtrsim3\sigma$ in Monte Carlo simulations. Interpreting the absorption as a proton cyclotron resonant scattering feature implies a local magnetic field of $B\sim(6-7)\times10^{14}$ G in the line-forming region. This value is consistent with strong magnetic fields anchored near the neutron star surface, even if the large-scale dipole is substantially weaker. Although we discuss electron cyclotron features and atomic transitions as possible alternatives, they appear less consistent with the observed phenomenology.

Extreme solar particle events (ESPEs) are caused by rare, enormously strong solar eruptions and can produce globally detectable spikes in tree-ring radiocarbon 14C, known as Miyake events, which serve as precise chronological tie-points and indicators of extreme solar activity. After production, radiocarbon is subjected to the complex carbon cycle, including large-scale atmospheric transport, which is crucially important for fast and strong Miyake events with highly inhomogeneous 14C production. A new 3D dynamical model, SOCOL:14C-Ex, of the radiocarbon atmospheric production and transport is presented here, which can model fast changes in the 14C atmospheric concentrations with high temporal and spatial resolution. Precise response curves of $\Delta^{14}$C to a reference ESPE (100xGLE#69) were computed for various event dates. They can be directly applied to analyse Miyake events under different conditions. Seven strong events over the past 14 millennia (AD 993, AD 774, 664 BC, 5260 BC, 5411 BC, 7177 BC, and 12351 BC) were analysed by fitting the reference curves to the available annual D14C data, identifying the most probable values and confidence intervals of their parameters -- strength, event's date and background level. By applying corrections for the geomagnetic and atmospheric (CO2) factors, the strengths of the corresponding ESPEs were assessed. The strongest ESPE is confirmed to be that of 12351 BC, while that of AD 774 remains the strongest event during the Holocene. To conclude, a new tool, based on the radiocarbon atmospheric transport model SOCOL:14C-Ex, is presented to analyse fast changes in the $^{14}$C production.

We propose a pulsar candidate cross matching algorithm to sift radio pulsar search candidates from repeated observations of the same sky location such as globular clusters, high energy sources, or supernova remnants. Our method uses both the candidate spin period ($P$) and dispersion measure (DM) value; if two or more candidates from different observations have similar spin periods to within 1\%, and dispersion measure values within 10\%, they are likely to correspond to the same candidate detection. We have demonstrated the effectiveness of our method through the discovery of the pulsar M12B with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). This pulsar has a spin period of 2.76\,ms and a dispersion measure of $42.70 \pm 0.05\,\mathrm{cm}^{-3}~\mathrm{pc}$. This pulsar has a profile with three peaks, being faint, showing scintillation. It is in an approximately 0.53-day orbit. Our discovery indicates that more pulsars might be effectively discovered if the algorithm is applied to the search results from other archival globular cluster observations.

Intra-pixel sensitivity variations (IPSVs) in charge-coupled devices (CCDs) and complementary metal-oxide-semiconductor (CMOS) detectors constitute a significant source of astrometric error for undersampled stellar observations. Since laboratory-based IPSV measurements suffer from limited applicability, we propose a computational method to directly infer IPSV from stellar images and validate it with simulated data. By minimizing the flux residuals between theoretical and observed stellar models through least-squares fitting, we can successfully recover the IPSV, which is treated as nearly identical across pixels. Simulations demonstrate that the reconstructed IPSV achieves high accuracy, and the instrumental point spread function (IPSF) restored using this IPSV improves stellar centroiding by nearly 30$\times$, effectively eliminating periodic pixel-phase errors. The method remains robust under different morphologies of IPSV and varying sampling conditions. Additionally, the framework can be extended to an iterative IPSF-IPSV closed-loop scheme that updates both components simultaneously, providing a practical pathway for continuous detector calibration in future space-based astronomical surveys.

We present a multi-wavelength study of the Flat Spectrum Radio Quasar PKS 1222+216, analyzing its long-term variability of radio data obtained in 2013-2020 from the iMOGABA, MOJAVE, and VLBA-BU-BLAZAR programs, along with $\gamma$-ray data from Fermi-LAT. We found that the radio flux densities at 15, 22, 43, and 86 GHz declined exponentially by 37%-56% over a year following a $\gamma$-ray flare in November 2014. We estimated jet physical parameters through Gaussian model fitting of VLBA 43 GHz data, identifying 10 jet components. The cooling timescales of the jet emission regions, i.e., newly ejected components C9, C10, and C11, range from 43 to 222 days, with the estimated jet viewing angles of approximately 8 degrees and magnetic field strengths of 77-134 mG in the jet emission regions. Additionally, by determining the magnetic field strength at different frequencies, we found that the magnetic field scales as $B\propto r^{-0.3\pm0.04}$, indicating a non-equipartition condition ($k_\text{r}\gtrsim 1$) or a slow decline in magnetic field strength profile ($m<1$). By analyzing component ejection times, we discovered that the $\gamma$-ray flare in 2014 coincided with the interaction between the moving component C9 and the stationary feature A1. We estimated that the $\gamma$-ray emission region is located at $9.2\pm1.0$ pc from the central engine, beyond the BLR and dusty torus, suggesting that the seed photons for inverse Compton scattering originate from the jet itself, external CMB radiation, or a surrounding sheath. Our results favor a scenario where $\gamma$-ray emission originates further downstream from the central engine through interactions between moving and stationary components. Additionally, our study presents an alternative method for estimating magnetic field strengths in AGNs undergoing long-term synchrotron cooling based on the associated timescale.

Infrared (IR) excess observed around white dwarfs (WDs) is typically attributed to companions or debris disks. These systems are interesting because they offer a unique opportunity to study the late stages of stellar evolution and the interactions between WDs and surrounding material. The 11th data release (DR11) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) - one of the largest spectroscopic surveys to date - has recently provided spectra for 3092 WDs, many of which have yet to be systematically investigated for IR excess. In this study, we cross-correlated the LAMOST DR11 WD catalog with optical and IR surveys, including the Sloan Digital Sky Survey (SDSS), Two Micron All-Sky Survey (2MASS), UKIRT Infrared Deep Sky Survey (UKIDSS), and Wide-field Infrared Survey Explorer (WISE). We performed spectral energy distribution fitting using the VOSA tool for 1818 WDs and identified 167 IR excess WD candidates. After excluding 23 sources with potential contamination within 6" and five additional sources identified through WISE ccf flag analysis, we identified 139 objects with candidate IR excess. These include 30 candidate WD + M dwarf binaries (18 new systems), 19 candidate WD + brown dwarf (BD) binaries (eight new systems), 66 candidate WD + dust disks (38 new systems), and 24 candidate either WD + BD or WD + dust disks (19 new systems). Given the limited spatial resolution of WISE, all candidate systems require follow-up IR observations for confirmation, such as high spatial resolution imaging or IR spectroscopy. This will help expand the parameter space of dust disks, allowing us to explore a broader range of possibilities.

Transit timing variations (TTVs) provide a powerful means to detect and characterise additional bodies in known planetary systems, even when they do not transit their host stars. We investigate the dynamical architecture of the HD 332231 system by analysing the TTVs of its close-in gas giant, HD 332231 b. Our goal is to assess whether the observed deviations from a linear ephemeris can be explained by the presence of an additional planetary companion. We refine the transit ephemeris of HD 332231 b using high-precision TESS photometry and complementary ground-based observations. We extract individual transit mid-times, construct an O-C diagram for transit timing data, and model the observed TTV signal through an extensive suite of N-body integrations covering a broad range of possible companion masses and orbital configurations. We detect a coherent TTV pattern with a period of approximately 6.7 years and an amplitude of about 45 minutes. Although numerous orbital configurations reproduce the observed TTVs, the combination of current radial velocity and photometric constraints yields a modest improvement in likelihood for solutions with an external planet on an orbit longer than 60 days, likely near a high-order mean-motion resonance and with moderate to high eccentricity. Our results suggest that HD~332231 b is part of a dynamically interacting multi-planet system. Continued transit monitoring and radial velocity follow-up will be essential to confirm the perturber's nature and refine the system's dynamical architecture.

Matter ejected from the magnetar crust during giant flares (GFs) may undergo $r$-process nucleosynthesis, producing short-lived optical transients termed "novae breves". Although intrinsically much fainter than kilonovae from compact binary mergers, novae breves may occur within or near the Galaxy, making them promising observational targets. We aim to investigate how the neutron star (NS) equation of state (EOS) and the mass of the central magnetar affect the ejecta properties following GFs and the resulting nova brevis emission. We employ a semi-analytical ejecta model combined with nuclear reaction network calculations to compute nucleosynthesis yields and multi-band light curves for different EOSs and magnetar masses, and assess their detectability with current and future facilities. We find that variations in the EOS and magnetar mass modify the ejecta mass and its density and velocity distributions, etc., leading to observable differences in nova brevis light curves. In particular, both the peak luminosity and the characteristic peak timescale are EOS-dependent. Assuming a fixed Galactic magnetar mass of 1.4 solar mass and taking the $u$ band as an example, we find that the minimum apparent AB magnitudes range from 7 mag (H4 EOS) to 8.5 mag (WFF EOS) with peak timescales of 100-1000 s. A more massive magnetar produces fainter emission with a shorter peak timescale. For a magnetar mass of 1.4 solar mass, novae breves associated with known magnetars may reach peak luminosities of 1e37-1e39 erg/s, enabling targeted searches, particularly following high-energy GF alerts. Moreover, a detection horizon of 10 Mpc or beyond is achievable with current and future facilities, allowing searches for novae breves from previously unknown magnetars in the Local Volume. Although challenging, detection of such rapidly evolving transients is feasible.

This study presents an investigation of nearly two dozen candidate Luminous Blue Variables (cLBVs) in the galaxies M31 and M33. Eight stars have been studied in detail, while an additional sixteen objects are briefly mentioned. Multi-epoch spectra of confirmed cLBVs from LAMOST and previous literature show broad hydrogen, He I lines, abundant Fe II and [Fe II] emission lines, and discernible spectral variability, consistent with the characteristics of known LBVs. Low outflow velocities inferred from P Cygni profiles are also incorporated into the classification criteria. Moreover, key stellar properties, including temperature and luminosity, are determined using the Spectral Energy Distribution (SED) fitting and spectral modeling. By comparison with stellar evolutionary tracks on the temperature luminosity diagram, the initial masses are estimated to be in the range of approximately 32 to 60 $M_{\odot}$. Except for J013401 and J013411, other stars locate within the typical LBV region between the S Doradus instability strip and their outburst phase. More importantly, our sample, except for the binary system, are all positioned in the LBVs region rather than that of B[e]SGs in the near-infrared color-color diagram. Based on all available information, one of the eight sources is confirmed as an LBV, four stars are designated as high-probability cLBVs, and the remaining three stars await further photometric observations to secure their classification. Given the current scarcity of known cLBVs, our study has the potential to make a significant increase in the number of LBVs in M31 and M33.

Competing processes shape Europa's surface: geological activity replenishes material through resurfacing, while bombardment by charged particles alters surface chemical composition. Each process leaves distinct spectral signatures. We present a novel data-driven analysis of JWST NIRSpec-IFU observations of Europa's leading hemisphere across three observing geometries, targeting nine spectral bands sensitive to water ice, radiolytic products, and volatiles. Through spectral factorization, we isolate the dominant components of spectral variability and reconstruct their spatial distributions. We find that CO2 enrichment extends beyond Tara Regio, and covers multiple chaos units in a lens-like pattern. These CO2-enriched areas co-occur with anomalous ice-texture signatures. Together, these findings suggest that enrichment in volatiles on Europa may reflect retention-favorable near-surface microphysics as well as emplacement, refining how they are interpreted in the context of surface--interior exchange. This has implications for interpreting the sources and supply rates of extant carbon-bearing species and, ultimately, for assessing Europa's habitability.

We aim to investigate the blowout jet-like prominence eruption, which occurred on October 6$^{th}$, 2023, with the help of imaging and spectroscopic observations. Firstly, the prominence rises slowly with a speed of 33 km/s, followed by a fast rise (i.e., 338 km/s). Later, the northern leg breaks completely, and the eruption forms the blowout jet. The jet consists of different plasma threads, which show a range of upflow (i.e., 125 to 593 km/s) and downflow velocities (i.e., 43 to 158 km/s). The jet plasma column exhibits transverse oscillations, and this motion (untwisting motion) propagate at the speed of 267 km/s, are consistent with being Alfev{é}n waves. The transverse motion has the time period, amplitude, and transverse velocity of 1332 s, 26.19 Mm, and 126.18$\pm$7.27 km/s, respectively, and this transverse oscillation decays over time. Interestingly, the different plasma threads within the jet's body exhibit decayless transverse oscillations, and these decayless oscillations are related to the main decaying transverse oscillation. The transverse velocity of these decayless oscillations ranges from 66 to 30 km/s, the amplitudes from 8.52 to 2.74 Mm, and periods from 811 to 406 s. In addition, the spectroscopic analysis reveals Si~{\sc iv} lines are forming in the optically thick conditions in high electron density regions (i.e., near the base of the blowout jet). Lastly, we mention that two weak C-class flares occurred during this event, and further, one CME also occurred, which propagated with the speed of $\sim$250 km/s.

The diversity of pulsar light-curves and radio polarisation properties originates in the structure of the magnetic field close to the stellar surface. For millisecond pulsars, this complexity is particularly puzzling. Fortunately, some means exist to uncover the magnetic field topology which indeed impacts the emission within the magnetosphere but also on the surface through its hot spot thermal radiation. We aim at deducing a plausible magnetic field geometry for the millisecond pulsar J0437$-$4715 by using combined information from the soft X-ray hot spot geometry deduced from NICER observations by pulse profile modelling and from radio and $\gamma$-ray pulse profile fitting. We also check the consistency between the geometry obtained and the radio polarisation data. Our $\gamma$-ray light-curve shapes rely on the striped wind model, whereas the radio polarisation fits rely on the rotating vector model. The magnetosphere structure is obtained from dipolar force-free magnetosphere simulations. We demonstrate that a slightly off-centred dipole augmented by a small scale dipole located on one polar cap explains simultaneously the shape of the hot spot and the radio and $\gamma$-ray data with a magnetic obliquity of $\alpha \approx (42\pm5) \degr$ and a line-of-sight inclination angle of $\zeta \approx (136 \pm5) \degr$. Our simple dipole model reproduces all the radio and $\gamma$-ray characteristics of PSR~J0437$-$4715, including its radio polarisation data. It shows that the radio emission could be produced in regions where the magnetic field is mainly of dipolar nature.

Sunspots are crucial for exploring how magnetic fields and plasma flows interact in the solar atmosphere, spanning from the stable photosphere to the shock-dominated chromosphere. To determine the thermal, magnetic, and kinematic properties of a sunspot across these layers and to investigate transient phenomena like umbral flashes, we analyzed high-resolution spectropolarimetric data from the CRISP instrument at the Swedish Solar Telescope. By applying the FIRTEZ inversion code, which incorporates non-local thermodynamic equilibrium (non-LTE) and 3D magneto-hydrostatic (MHS) equilibrium, to full Stokes measurements of multiple spectral lines (Mg I, Na I, Fe I, and Ca II), we successfully mapped the atmospheric parameters in a 3D domain. Our analysis reveals that the photospheric Evershed flow actually reverses into an inflow in the upper photosphere. In contrast, the surrounding moat flow persists as an outflow at similar heights, indicating that it is not a direct continuation of the Evershed flow. Furthermore, observations of an umbral flash event uncovered supersonic upflows (Mach numbers $\|M\|\geq 1.5$) and thermodynamic conditions characteristic of shock fronts. Ultimately, combining 3D MHS equilibrium and non-LTE effects across multiple spectral lines proves highly effective for simultaneously constraining parameters in both the photosphere and chromosphere. These findings provide clear evidence of shock dynamics in umbral flashes, supporting the theory that converging supersonic flows act as the primary driving mechanism while shifting optical depth iso-surfaces.

Physical models aimed to reproduce basic features of the solar sunspot cycle are typically based on the solar dynamo mechanism. Usually qualitative arguments are used to define parameters of the model, among which a challenging component is the nonlinear form of quenching of the alpha-effect governing regeneration of the magnetic field. We propose a novel approach, in which the functional form of the alpha-quenching is represented by a neural network model embedded into neural differential dynamo equations trained on observational data. For demonstration, we consider a low-mode dynamo model and find a wide set of alpha-quenching functions and corresponding dynamo numbers that provide an accurate fit to the average profile of the solar cycle data given by sunspot numbers. Within this set, we observe a strong relationship between the dynamo number and the shape of the alpha-quenching function indicating that additional magnetic field data or constraints are essential to unambiguously infer parameters of the dynamo model. In our opinion, the neural differential approach opens a new prospect for data-driven investigation of the closure problem in dynamo theory.

Aims. We aim to determine how ion neutral coupling and ambipolar diffusion affect the linear and the nonlinear growth of the RTinstability under astrophysically relevant conditions, and to identify the coupling regimes in which departures from the classical single fluid picture become significant. Methods. We perform high resolution two fluid numerical simulations using the MPI AMRVAC code, spanning a wide range of perturbation wavelengths, coupling strengths, from uncoupled to strongly coupled passing by intermediate or ambipolar diffusion dominated regimes, and magnetic field configurations. The linear theory is revisited using a physically consistent formulation with different ion neutral coupling strengths across the interface and validated against the simulations. We investigate the physics of the instability using morphology based diagnostics of the mixing layer to compare simulations at equivalent nonlinear stages, complemented by spectral, force, and energy budgets analyses. Results. In the linear regime, theoretical growth rates are recovered over a wide range of wavelengths, from the single fluid limit to intermediate bi fluid coupling. In the nonlinear regime, ambipolar diffusion modifies the classical quadratic growth and introduces a coupling dependent evolution. For multi wavelength perturbations, the nonlinear dynamics becomes strongly scale dependent: intermediate coupling enhances fragmentation in hydrodynamic configurations, while magnetised cases exhibit a non monotonic reorganisation of the interface, with the smoothest morphologies occurring at intermediate coupling. Spectral and energetic diagnostics indicate that these behaviours correlate with changes in the relative contributions of ion neutral drift and magnetic stresses during thenonlinear evolution

Extended radio sources present unique challenges for automated detection and classification in wide-field radio surveys. With current surveys such as the Evolutionary Map of the Universe (EMU), robust and scalable methods are essential to identify and catalogue these complex sources. We apply three automatic approaches to detect complex radio emission in EMU observations of the Galaxy And Mass Assembly (GAMA) 09 field (EMU-G09) in order to evaluate their relative strengths and limitations in preparation for large-scale application across future EMU data releases. These include DRAGNHunter, designed to detect likely DRAGNs (Double Radio sources associated with Active Galactic Nuclei) from a component catalogue; coarse-grained complexity, a metric designed to highlight regions of complex emission; and RG-CAT, a machine learning pipeline trained on radio sources identified in the EMU pilot survey. We find that together, the three methods recover nearly all extended sources in EMU-G09 but identify largely distinct, partially-overlapping subsets, with only 375 sources identified by all finders. This demonstrates that a combination of complementary techniques will be required to achieve a complete census of extended radio sources in future large-scale surveys.

Remnant black holes (BHs) of massive stars (``light seeds'') are a potential origin for supermassive black holes (SMBHs). We use magnetohydrodynamic simulations to study the formation and growth of light seeds in star-forming giant molecular clouds (GMCs) with masses $10^5$--$10^9\,M_\odot$, which evolve for $\sim 10$--$30\,\rm Myr$ and form compact star clusters, akin to high-redshift nuclear star clusters. In particular, the simulations resolve very massive stars (VMSs, 100--$300\,M_\odot$), including their radiative and mechanical feedback, and model feedback-regulated accretion onto remnant BHs. We find that, even in compact GMCs capable of forming deep potential wells, the gas reservoir is expelled by sustained stellar feedback and rapidly dispersed after supernova explosions. Remnant BH populations emerge $\sim 3\,\rm Myr$ after the starburst and concentrate at the cluster center (where $\rho_{\rm BH}\sim 10^4$--$10^6\,M_\odot\,{\rm pc}^{-3}$). With our fiducial sub-grid BH accretion/feedback model, in-situ BH accretion is inefficient for forming heavy seeds: some direct-collapse BHs briefly accrete at $\sim$\,(1--10)$\times$ the Eddington rate, but they reach only $\sim 400$--$500\,M_\odot$. A top-heavy initial mass function or natal kicks do not change this conclusion. Runaway accretion is only possible if the sub-grid BH model allows a high fraction of Bondi inflow to reach the BH, in which case a few seeds can grow to $\sim 10^6\,M_\odot$. We also discuss multiple-generation star formation that may be intrinsically correlated with remnant BH accretion.

As the first asteroid deflection test, NASA's successfully hit asteroid Dimorphos (secondary of the binary asteroid 65803 Didymos) with the DART kinetic impactor on September 26, 2022. To fully characterise the physical properties of the objects, and measure precisely the effects of this impact in the context of planetary defence, ESA launched the Hera mission on 7 October 2024, with scheduled arrival at Didymos in fall 2026. Among the core payload of the mission, the Asteroid Framing Cameras are two identical imaging systems that will support navigation and scientific activities, by acquiring images from various distances and observing geometries during the course of the mission. Built by ena-Optronik (Germany), the cameras match the requirements designed by the science team and will provide data that supports a wide range of investigations: hazard detection, system dynamics, shape reconstruction, surface morphology and mapping, and surface photometry. Each instrument is a panchromatic camera equipped with a 5.5 x 5.5 degree field of view, and an angular resolution of 93.7 micro-radians per pixel. The cameras shall provide the necessary data to address the mission requirements through a global mapping of the two components of the binary system at spatial scales of 2-3 m/pixel in the Early Characterisation Phase, 1-2 m/pixel in the Detailed Characterisation Phase, and 0.5-2 m/pixel in the Close Operation Phase. Dedicated flybys will bring the resolution down to < 10 cm/pixel on specific areas of interest on Dimorphos, such as the DART impact site and the JUVENTAS cubesat landing site. Here, we present the technical specifications of the camera, as well as the status of the calibration. We then summarise the planned operations in cruise and at the asteroids. Finally, we provide examples of the scientific investigations and products that will make use of the data returned by the cameras

The accretion of metal-poor gas sustains galactic star formation. In the Milky Way, this process is fueled by high-velocity clouds (HVCs), yet their fundamental properties have remained elusive in the absence of stellar tracers. Here we report a binary open cluster within HVC Complex H. With an age of 11.2 +- 0.6 Myr and a subsolar metallicity of 0.05(+0.05-0.02) Zsun, the clusters provide a direct stellar distance anchor to the cloud at 13.8 +- 0.6 kpc. Their proper motions indicate Complex H is on a prograde, south-to-north orbit through the outer Galactic disk. The resulting interaction produces a 'slow-fast-slow' velocity gradient, with the cloud's outer layers decelerating as they merge into the disk. Orbit integration suggests the clusters formed from an internal cloud-cloud collision. This triggering mechanism implies other HVCs could similarly produce high-velocity stars. The scarcity of previous stellar detections in HVCs is explained by the rapid escape of young stars (< 20 Myr), while CO non-detections may stem from weak emission due to low metallicity and gas dispersal. This work reveals that the circumgalactic medium can sustain star formation, offering a tangible laboratory to probe the physical conditions of accreting gas before it merges with the Galactic disk.

We introduce a new class of interacting dark sector models that couple the intrinsic entropy of dark matter to scalar field dark energy. Using the Lagrangian formulation for relativistic perfect fluids, we construct consistent covariant actions that incorporate algebraic and derivative entropy couplings. These interactions leave the expansion history unchanged, rendering the background cosmology indistinguishable from $\Lambda$CDM or uncoupled quintessence. At the level of cosmological perturbations, the entropy couplings generate scale-dependent modifications to the dark matter Euler equation, while the continuity equation remains unaltered at linear order. The resulting interactions correspond to a pure-momentum exchange within the dark sector. We show that intrinsic entropy perturbations can carry primordial scale dependence, and non-minimal couplings can lead to a scale-dependent suppression or enhancement of structure growth. Finally, we demonstrate that these models are generically compatible with current Cosmic Microwave Background observations, while inducing distinctive signatures in large-scale structure. The framework provides a theoretically well-motivated and observationally viable extension to the standard cosmological model, opening new directions to explore novel interactions in the dark sector.

Recently, a number of methods have emerged to detect galaxy clusters solely through their weak lensing signal. Using the recently-introduced wavelet multi-scale detection method, we focus here on the potential for the use of tomographic information of the source galaxies to increase the number of weak lensing detections. We apply the $z_{s,\mathrm{min}}$-cut technique, consisting of the combination of weak lensing peak detections emerging from lensing maps obtained using different source redshift bins, to mock data sets of progressively increasing sophistication. The source redshift distribution is chosen to be $Euclid$-like, with a maximum depth of $z_{s,\mathrm{max}}=3$, and overlapping tomographic redshift bins are constructed by progressively increasing the minimum source redshift $z_{s,\mathrm{min}}$. Considering all possible detection combinations from one to four tomographic bins, we find that a single source redshift bin, with $z_{s,\mathrm{min}}=0.4$, performs as well as the combination of multiple redshift bins. By running detections on synthetic clusters of varying complexity -- from isolated Navarro Frenk White haloes to haloes embedded in and formed within N-body cosmological simulations, and considering both true and photometric source redshifts -- we show that while large-scale structure contamination and photometric redshift errors reduce the potential gains of the tomographic approach, the dominant limitation is the accumulation of spurious detections across redshift bins, leading to decreased purity at a fixed detection threshold.

Plasma composition in flaring regions has been shown to have significant spatial and temporal variations, likely driven by dynamical processes that take place as a consequence of the sudden energy release at the reconnection site. The origins of these variations, as well as the effects they might, in turn, have on flare loops dynamics are not yet fully understood. In this work, we investigate the link between flare loop cooling times and plasma composition evolution in the loops formed during the M-class flare peaking at 13:56 UT on the 2022 April 2 using high cadence Hinode EIS spectroscopic observations. The analysis focuses on quantifying the cooling rate (using a series of emission lines covering a wide temperature range) and plasma composition evolution (using the Ca XIV 193.866 A/Ar XIV 194.401 A diagnostic) at the apex and footpoint of the flare loop arcade. Results show slower cooling and a FIP bias of 2.4 +/- 0.2 in the loop footpoint and faster cooling and a stronger FIP bias of 2.8 +/- 0.2 in the loop apex. The potential effects of plasma composition changes on the radiative cooling process of flare loops are also investigated by comparing observed loop cooling times to those predicted by simulations from the EBTEL 0D hydrodynamic model. The EBTEL simulations show that an higher FIP bias would lead to a faster radiative cooling rate and, therefore, shorter cooling times. This suggests that the variation in FIP bias observed in the two features could be responsible for the different cooling times observed.

The Lunar Electromagnetic Monitor in X-rays (LEM-X) is a proposed wide-field X-ray observatory designed for deployment on the Moon's surface. Its primary scientific goal is to enhance multi-messenger astrophysics by detecting, localizing, and monitoring high-energy transient phenomena and variable X-ray sources across the sky. Building on the heritage of the eXTP and LOFT mission proposals, LEM-X employs pairs of coded-aperture cameras equipped with large-area linear Silicon Drift Detectors (SDDs), offering excellent spectral resolution ($\leq$350 eV at 6 keV) over the 2-50 keV energy range. Each camera provides a field of view of ~1 steradian at 25% effective area and achieves a Point-Source Location Accuracy (PSLA) of 1 arcminute, with an on-axis sensitivity better than 5 mCrab in 50 ks and 700 mCrab in 1 s. In this paper we describe the experiment and focus on the detailed design and optimization of the LEM-X coded mask, analyzing its scientific performance, imaging capabilities, and thermo-mechanical properties. We describe the mask code generation, decoding algorithms, and the trade-offs involved in achieving the required angular resolution, sensitivity, and structural integrity. Imaging simulations and mechanical analyses confirm the effectiveness of the proposed design, demonstrating its suitability for high-precision, wide-field X-ray imaging devoted to multi-messenger astrophysics and transient events detection.

Estimating the meteoroid flux density at centimetre to metre sizes is notoriously difficult. Yet it is an important endeavour, as these sizes represent the transition between small meteoroids that pose a risk to spacecraft, and the Near-Earth Objects that are relevant for planetary defense. We present a novel automated methodology for debiasing meteor observations from multi-camera networks, applied to data from the Desert Fireball Network (DFN). Our approach utilizes the Hierarchical Equal Area isoLatitude Pixelisation (HEALPix) framework to partition the sky into equal-area pixels at 70 km altitude, enabling precise and convenient measurement of effective survey coverage and fireball counting across the network. We developed a comprehensive data processing pipeline that analyses millions of all-sky camera images to determine clear-sky conditions through automated star source detection and flux distribution analysis. As a case study, we apply this methodology to observations of the 2015 Southern Taurid meteor shower, during which there was significant fireball activity. Processing data from 33 cameras over a three-month period (October-December 2015), we calculate an effective observation coverage of $1.58 \times 10^{12}$ km$^2$.h and identified 54 Southern Taurid fireballs from 141 validated detections. Our results are consistent with the extrapolation of previous work done on the same meteor shower at smaller sizes, when we set a $\sim300$ kg.m$^{-3}$ mean meteoroid density, consistent with the cometary origin of the Taurid stream. The HEALPix-based approach successfully automates what was previously a labor-intensive manual process, providing a scalable solution for accurate flux measurements from distributed camera networks; it is directly applicable to other meteor surveys.

The recent release of high-precision cosmological data, particularly the small-scale cosmic microwave background (CMB) measurements from ACT and baryon acoustic oscillation (BAO) data from DESI, has opened a new landscape for probing the neutrino mass. In this work, we present updated constraints on the total neutrino mass, $\sum m_\nu$, and its hierarchy within the $\Lambda$CDM, $w$CDM, holographic dark energy (HDE), and $w_0w_a$CDM models, using the latest ACT DR6, DESI DR2, and DESY5 datasets. We find that the upper limits on $\sum m_\nu$ are critically governed by the evolutionary behavior of the dark energy equation of state. Specifically, models exhibiting early-time quintessence features (e.g., HDE) yield the most stringent constraints, whereas those allowing for early-time phantom behavior (e.g., $w_0w_a$CDM) result in significantly looser bounds. Despite these model-dependent variations, we observe a robust hierarchy dependence across all scenarios, where the inverted hierarchy consistently yields weaker constraints and the degenerate hierarchy consistently yields tightest constraints. Our analysis demonstrates that the improved small-scale CMB information from ACT, combined with high-precision BAO data, systematically tightens the limits on $\sum m_\nu$, providing a crucial benchmark for future neutrino mass measurement.

Pulsar timing of double neutron star (DNS) systems is one of the best methodologies to study the neutron star masses distribution. Here we report the discovery of a double neutron star system PSR J0641+0448 in the Five-hundred-meter Aperture Spherical radio Telescope (FAST) Galactic Plane Pulsar Snapshot (GPPS) survey. This pulsar has a 25.7 ms spin period and moves in a 3.73-days eccentric orbit with an eccentricity of 0.145. Using FAST observations, we obtained its phase-connected timing solution with periastron advance and Shapiro delay detected. Using $\chi^2$ analysis based on DDGR model, we constrain the pulsar mass to $1.319^{+0.021}_{-0.035}~M_\odot$, and the companion mass to $1.269^{+0.022}_{-0.016}~M_\odot$ with a 68.3\% confidence level. The low companion mass and mild orbital eccentricity is consistent with the correlation between neutron masses and orbital eccentricities.

The rapid growth of satellite constellations, particularly Starlink, is increasingly affecting ground-based astronomy. In this project, we developed a workflow to detect, identify, and measure the brightness of trails from artificial satellites and other orbiting objects in archival images from the Dark Energy Camera (DECam), available through the NOIRLab Data Archive. We filtered images with visible streaks, retrieved detector-level images, applied the Hough Transform (via satmetrics) to detect and align trails, and performed surface brightness photometry. We also used SatChecker to obtain likely identifications for each trail. Our sample of nine measured streaks includes Starlink satellites, a navigation satellite, a decommissioned science satellite, and a rocket body. Our results show that satellites and other orbiting objects are consistently detectable in DECam images, but their brightness varies significantly, reflecting design and operational differences across object types and models. While the methodology proved effective, detecting faint streaks was challenging, and short-lived glints remain an even harder problem for future work. This proof-of-concept establishes a foundation for larger statistical studies of satellite impacts on astronomical surveys. The code is available at this https URL

The recent evidence for a nanohertz gravitational wave background from Pulsar Timing Arrays highlights the urgent need to identify electromagnetic counterparts to supermassive binary black holes. Here, we perform global 3D general relativistic magnetohydrodynamic (GRMHD) simulations of a secondary black hole (mass ratio $q=0.1$) interacting with a Magnetically Arrested Disk around a primary black hole using a time-dependent superposed Kerr-Schild metric and post-processed general relativistic radiation transfer calculations based on thermal electron distribution function (eDF). We explore three orbital configurations: a vertical impact orbit, a coplanar embedded orbit, and a high-spin, eccentric, inclined scenario. Despite clear orbital periodicity and recurrent shock formation, the thermal synchrotron light curves frequently lack expected shock-induced flares. In vertical impacts, shock brightenings are typically sub-dominant to the stochastic MAD variability of the primary black hole, unless viewed at specific alignment phases. Conversely, coplanar orbits produce distinctive, rapid flares driven by gravitational self-lensing. We identify a frequency-dependent emission hierarchy: the primary black hole dominates sub-millimeter flux, while the secondary dominates near-infrared emission due to higher electron temperatures in thermal eDF. Finally, spin-orbit coupling drives Lense-Thirring precession, yielding twisted, wobbling jets reminiscent of OJ~287. Crucially, we demonstrate that intrinsic MAD turbulence may easily mask shock-induced flares at radio frequencies. We strongly advocate coordinated sub-millimeter and near-infrared monitoring to robustly isolate supermassive binary black hole self-lensing signatures.

We measure transverse proper motion velocities of LMC Wolf-Rayet (WR) stars using Gaia DR3 astrometry. The combined velocity distribution of WNh, O If*/WN, and WNL very massive stars ($>100\ M_\odot$; VMS) shows both slow, unejected objects ($v_\perp < 10$ $\rm km\ s^{-1}$) and stars dominated by fast, runaway velocities ($v_\perp > 24$ $\rm km\ s^{-1}$). This supports expectations that VMS ages are comparable to the dynamical ejection timescale ($\sim1.5$ Myr). These kinematics share similarities with those of lower-luminosity, classical WNh, O If*/WN, and WNL stars, as well as the SMC field OB stars, suggesting that dynamical ejections may also dominate these populations. In contrast, both single and binary WNE stars are ejected populations that show single-peaked velocity distributions, suggesting a different ejection mechanism(s). We speculate that single WNE stars might result from explosive mergers onto the shell-burning layer, thereby stripping the H envelope. Binary WC stars appear to be faster (median $v_\perp = 54$ $\rm km\ s^{-1}$) and have higher luminosities than singles (median $v_\perp = 38$ $\rm km\ s^{-1}$), suggesting that single WC stars are not descendants of the binaries. Thus, the binaries are probably stripped by mass transfer, while the WC singles likely originate from another process. The high velocities of binary WC stars are consistent with some predictions that lower mass clusters generate fast dynamical ejections. Single WC and WN3/O3 stars have ambiguous kinematics, but both show high $v_\perp$ (median $\sim 38$ $\rm km\ s^{-1}$), possibly linked to their lower masses.

The Peanut cluster (SRGe J023820.8+200556, SRGe CL0238.3+2005, $z_{spec}$ = 0.42) has recently emerged as a candidate for a rare, massive merger, potentially analogous to the Bullet cluster. We present the results of optical identification and spectroscopic redshift measurements for 31 galaxies in the Peanut cluster, including 26 new redshifts obtained with the 6-m telescope BTA (Big Telescope Alt-azimuthal) at SAO RAS between October 2024 and January 2025. The derived distribution of line-of-sight velocities reveals the possible presence of two subclusters with a line-of-sight velocity difference of ~2000 km/s. However, statistical tests and the Dressler-Schectman test show that the hypothesis that the observed velocity distribution can be described by a normal distribution for a single cluster cannot be ruled out, and the evidence for the existence of two gravitationally bound substructures remains ambiguous. Assuming a single cluster with the normal velocity distribution, the estimated galaxy velocity dispersion is $\sigma_{los} = 1455 \pm 83$ km/s, corresponding to the total cluster mass of $M_{200} = 2 \times 10^{15} M_\odot$ based on the mass-velocity dispersion scaling relation. In either scenario -- a single extremely massive cluster or an ongoing merger -- the Peanut cluster appears to be a very rare and peculiar object, comparable to such extreme systems as the Bullet cluster (1E 0657-56) or El Gordo (ACT-CL J0102-4915).

We aim to identify the physical mechanism responsible for the observed X-ray emission and polarization of Cyg X-1 in the soft state. We perform a detailed spectral analysis of X-ray data obtained by NICER, NuSTAR, and INTEGRAL during observations simultaneous with IXPE on 2023 June 20, supplemented at higher energies with archival CGRO data. We develop a new model, retBB, that describes thermal disk emission and its returning reflection, and apply it together with accurate models of Comptonization and relativistically broadened reflection. Using the resulting spectral solution, we compute the expected polarization signal and compare it with the IXPE measurements. Our spectropolarimetric modeling shows that the observed polarization is produced by Comptonization in a corona undergoing a semi-relativistic outflow with a velocity of ~0.3c. Our spectral solutions admit either low or high black-hole spin values, depending on the adopted model setup. However, the observed polarization strongly favors a low spin. At high spin, the polarization angle would inevitably rotate significantly across the energy band, which is not consistent with the observations. Apart from this rotation of the polarization angle, general relativistic effects do not play a significant role in producing the observed polarization. In particular, we find that, at most, returning disk radiation contributes only a minor amount.

The SPEARHEAD (specification, analysis, and re-calibration of high-energy particle data) project, funded by the European Union Horizon Europe programme, enhances the accuracy and usability of high-energy particle measurements. It investigates particle acceleration, release, and transport during solar eruptions by refining instrument response functions and cross-calibrating datasets. We present a comprehensive catalogue of greater than 100 MeV proton events identified from May 1996 to August 2024 using SOHO-ERNE penetrating-particle rates, together with associated solar phenomena derived from multi-instrument observations. The SEP events were detected through a systematic scan of ERNE-HED counter data and cross-calibrated with SOHO-EPHIN to derive peak fluxes and fluences. Each event was associated with its likely parent eruption using X-ray (XR) (GOES-XRS, RHESSI, SolO-STIX), radio (Wind-WAVES, STEREO-WAVES, ground-based observatories), and gamma-ray (Fermi-LAT) observations, CMEs (SOHO-LASCO), and ground-level enhancements (GLEs) (neutron monitors). Timings and physical properties were systematically compared to investigate the relationships between SEP onset, flare evolution, CME kinematics, and radio signatures. Statistical analyses reveal that most SEP releases closely follow flare and CME onsets, with moderate SEP-XR-CME correlations, and a strong SEP-GLE fluence link. These results indicate that high-energy SEP events are typically associated with strong solar activity signatures, with the observed intensities and timings strongly modulated by magnetic connectivity and coronal conditions. This catalogue provides the most extensive reference to date for high-energy SEP events over solar cycles 23-25, establishing a unified framework for future investigations of extreme particle acceleration.

Rocky exoplanets are particularly abundant around M-type stars. Their small radii and low luminosities provide favourable conditions for detecting transiting terrestrial planets and probing their atmospheric properties. We report the discovery and statistical validation of TOI-4616 b, an Earth-sized planet transiting a nearby mid-M dwarf observed by the Transiting Exoplanet Survey Satellite (TESS). We confirm the planetary nature of the signal and determine the system parameters by combining TESS photometry with ground-based multi-band transit observations, high-resolution imaging, and optical and near-infrared spectroscopy. The host star lies at a distance of 28.10 +(-) 0.07 pc and has a radius of 0.1889 +(-)0.0096 solar radii, a mass of 0.1881 +(-) 0.0094 solar masses, and an effective temperature of 3150 +(-) 75 K. TOI-4616 b has a radius of 1.22 Earth radii and an orbital period of 1.55 days. The planet receives an incident flux of approximately 40 times that of Earth, corresponding to an equilibrium temperature of about 525 K. This places TOI-4616 b in a regime intermediate between Earth-sized planets orbiting early M dwarfs and those around ultra-cool hosts. Statistical validation with the TRICERATOPS framework, supported by high-resolution imaging and chromatic transit constraints, yields a false-positive probability of 0.0135, below the recommended validation threshold of 0.015, confirming TOI-4616 b as a validated planet. Owing to its proximity to Earth, well-constrained stellar properties, and extensive multi-band follow-up, TOI-4616 b constitutes a valuable benchmark system for comparative studies of terrestrial planets around mid-M dwarfs and for future atmospheric investigations.

We report the discovery of the optical dipper and low-luminosity infrared stellar transient Gaia20ehk (hereafter, Gaia-GIC-1), which is currently undergoing high-amplitude variability due to transiting dusty material. In this work, we identify Gaia-GIC-1 as a likely young F-type star based on the spectral energy distribution before the onset of the high-amplitude optical variability. We detect a significant periodic modulation of 380.5 days in Gaia-G band before the onset of the infrared brightening, consistent with a $\sim$1.1 AU orbit assuming circular orbits and a 1.3 M$_{\odot}$ star. The system has remained in an infrared bright state for $>$4 years since the last near-infrared detection, confirmed by recent SPHEREx observations, while continuing to undergo large amplitude irregular optical dimming. We measure the dust temperature from the freshly generated debris to be $\sim$900 Kelvin based on available WISE photometry, and the dust clump size to have a minimum cross-sectional area of 0.13 AU$^{2}$, and the dust mass $4\times 10^{20}$ kg. Currently, optical follow-up spectroscopy has not revealed any prominent features in the system, likely due to its highly variable nature. We hypothesize that Gaia-GIC-1 represents debris recently formed in a planetary collision, which produced a clumpy dust cloud on a bound orbit, producing the observed dimming events. The ongoing collisional activity in this system presents a unique opportunity for understanding terrestrial planet formation.

The Shakura-Sunyaev $\alpha$-viscosity coefficient, defined as the ratio of total stress to total pressure, $\alpha= \mathbb{T}/p$, played an important role in the development of the accretion disc theory in the early 1970s. The origin of turbulence that causes the stress $\mathbb{T}$ was unknown at that time. Shakura and Sunyaev assumed $\alpha=$ const. Today we know that this was not quite realistic - the modern general relativistic magneto-hydrodynamic simulations (GRMHD) of black hole accretion discs revealed that $\alpha$ changes by about an order of magnitude within the disc, being smaller far away from the black hole and larger in the plunging region close in. It was found that the behaviour of $\alpha$ reflects some underlying, fundamental properties of the stress $\mathbb{T}$ itself. In particular, as argued by several authors, the stress must be zero at the black hole horizon. We notice that the stress calculated in GRMHD simulations by different authors, including us, has a maximum rather close to the location of the circular photon orbit. We propose a formula that accurately describes this universal behaviour of $\alpha$ in terms of the "gyration radius'', a physical characteristic of rotation well known in Newtonian dynamics and in the black hole case uniquely defined by the Kerr space-time geometry. Analytic and semi-analytic models of black hole accretion discs provide an invaluable insight into fundamental physics, and the GRMHD simulations do not aspire to replace them. Rather, simulations could help to improve analytic models by making them more realistic. For example, our $\alpha$-formula, deduced from the GRMHD simulations, may be handy in the construction of improved versions of thin and slim disc models.

Magnetic white dwarfs (MWDs) are key to understanding the origin and evolution of magnetic fields in compact stars. While large spectroscopic surveys such as SDSS have greatly expanded the known sample, the potential of LAMOST has not yet been fully explored. Our aim is to identify and characterize isolated MWDs in the LAMOST DR10 database. We cross-matched LAMOST DR10 spectra with white dwarf candidates from Gaia EDR3 and with recent SDSS-based catalogs of MWDs. Zeeman splitting in Balmer and helium absorption lines was used as the primary diagnostic to identify magnetic fields and to estimate their strengths. Reference objects from SDSS catalogs were used to test the detectability of MWDs in LAMOST low-resolution spectra. We identified 63 isolated MWDs in LAMOST DR10, of which 32 are new discoveries. Surface magnetic field strengths were measured from Zeeman splitting, covering a range from a few MG up to several tens of MG. For previously known SDSS MWDs, our LAMOST-based field measurements show mostly agreement with published values. This work demonstrates the capability of LAMOST low-resolution spectroscopy to identify and characterize isolated MWDs. The newly discovered objects expand the known population and provide valuable targets for future high-resolution spectroscopic and polarimetric follow-up studies. Our results highlight the potential of combining LAMOST with Gaia and other large surveys to build a more complete census of MWDs.

Recent work by Chae et al. (2026) reported a gravitational anomaly in 36 wide-binary pairs, finding a gravity boost factor of $\gamma \equiv G_{\rm eff}/G_{\rm N} \approx 1.60_{-0.14}^{+0.17}$ at low accelerations, consistent with predictions from Modified Newtonian Dynamics (MOND). We reanalyze the same dataset using a hierarchical Bayesian model that infers a global $\gamma$ across all systems while fitting three-dimensional orbital elements. Our model yields $\gamma = 1.12^{+0.27}_{-0.22}$, consistent with Newtonian gravity ($\gamma = 1$) at the $\sim0.4\sigma$ level. To identify the source of the discrepancy, we perform a test using an approach similar to Chae et al. (2026), replacing the semi-major axis with a geometric de-projection of the observed projected separation. This test yields $\gamma = 1.56^{+0.21}_{-0.18}$, closely matching the result of Chae et al. (2026). This suggests that the inferred value of $\gamma$ is sensitive to how the three-dimensional orbital separation is modeled, and including an independent semi-major axis parameter can account for velocity excesses that would otherwise be attributed to non-Newtonian gravity.

Oxygen is a well-studied biosignature. Studying potential abiotic pathways for O2 build-up in exoplanet atmospheres is essential for evaluating whether the detection of O2 would constitute a biosignature detection on other worlds. Previous modeling efforts in the literature demonstrated that detectable abiotic O2 and O3 can be produced through CO2 photolysis for rocky planets around M dwarf stars. Building on modeling approaches from previous studies, we use photochemical simulations to reassess the conditions under which O2 and O3 may accumulate through similar photochemical mechanisms. Using a Mars-like atmospheric composition and planetary parameters, we vary the hydrogen mole fraction to assess how changes in HOx chemistry can affect the resulting accumulation of abiotic O2 and O3. Across the range of hydrogen mole fractions explored, we obtain a maximum O2 abundance of ~2.7% for H = 0.0065 ppm, about an order of magnitude lower than reported in the literature. This reduction is consistent with the elevated water vapor abundance adopted in our simulations, which enhances HOx-driven recycling of CO and O and thereby suppresses the accumulation of O2 and O3. Our improved understanding of how this cycle results in atmospheric false positive biosignatures in crucial towards developing future exoplanet characterization strategies.

Observations of the Cosmic Microwave Background (CMB) B-mode polarisation provide a unique probe of inflationary physics. Extracting a reliable constraint on the tensor-to-scalar ratio $r$ nonetheless demands stringent suppression of diffuse Galactic foregrounds, whose residuals can bias the inferred signal. This work introduces and evaluates two extensions of the Needlet Internal Linear Combination (NILC) framework aimed at reducing foreground-induced biases on $r$. The first extension implements the deprojection of selected foreground moments directly within the component-separation step. The second performs a likelihood-level marginalisation over residual foreground power using a data-driven template. Using Simons Observatory Small Aperture Telescope (SO-SAT) - like simulations, we show that both methods effectively control residual contamination, yielding unbiased estimates of $r$ and a consistent reconstruction of the lensing B-mode amplitude. These results indicate that enhanced foreground-mitigation strategies will be useful for next-generation CMB polarisation analyses seeking a robust detection of primordial B-modes.

1FGL J0523.5$-$2529 is a Fermi selected redback millisecond pulsar candidate that exhibited luminous optical and X-ray flares in 2020-2021. We obtained a simultaneous X-ray and $U$-band observation with XMM-Newton in 2025, the first to cover the 16.5 hr orbit of 1FGL J0523.5$-$2529. The X-ray luminosity was in an intermediate state with a power-law photon spectral index of $\Gamma=1.53\pm0.02$. Frequent flares were superposed on a broad, single-peaked modulation, the latter characteristic of intrabinary shock models in which the shock front is wrapped around the pulsar. We speculate that density enhancements in the shocked companion wind cause flares, as well as variable optical recombination lines. The $U$-band light curve was dominated by ellipsoidal modulation of the nearly Roche lobe filling companion star, similar to that seen in ground-based optical photometry. We also used this effect in 10 years of ATLAS monitoring to improve the precision of the orbital period to 0.6881366(19) days. Considering that searches for radio pulsations from 1FGL J0523.5$-$2529 at all orbital phases have been unsuccessful, the shocked wind usually surrounds the pulsar.

We study the Stodolsky effect utilizing the most general form of neutrino interactions with electrons below the electroweak scale by considering all possible Lorentz invariant operators respecting SU(3)$\otimes$U(1) symmetry. We perform our calculation for both Dirac and Majorana neutrinos and find that in the most general setting, only the non-standard neutrino interactions and the tensor interaction terms provide a non-zero contribution, apart from the Standard Model contribution. We investigate the implications for the possible detection of the cosmic neutrino background (C$\nu$B) by analysing the energy shifts that are characteristic of the Stodolsky effect. We also discuss the implication of considerable asymmetry in the C$\nu$B on the present scenario.

Pulsars possess strong dipole magnetic fields that can source axion fields through the axion-photon interaction. Pulsars may therefore be surrounded by axion field configurations oscillating with the pulsar's rotational period. These axions could be detected by observing their effect on the polarization of the pular's emission. In this paper, we use time resolved observations of the optical polarization of the Crab pulsar to place bounds on the axion-photon coupling, demonstrating the potential of time resolved pulsar birefringence in the search for axions.

We investigate the phenomenology induced by generic PQ-breaking operators within the axion kinetic misalignment framework. We analyze their impact on the relic density of axion dark matter (DM), the PQ quality problem, axion-mediated fifth-force, as well as Big Bang Nucleosynthesis (BBN) and Cosmic Microwave Background (CMB) constraints. A nonzero initial axion velocity gives rise to brief periods of early matter domination and axion kinetic domination, leading to a nonstandard cosmological evolution. We compute the resulting gravitational wave (GW) signal from global cosmic strings and find that, because these nonstandard epochs are extremely short, the signal is highly suppressed and beyond the reach of existing experiments. Finally, we perform a parameter space scan, identify the regions and benchmark point that are consistent with all experimental constraints.

We formulate an extended holographic dark energy scenario based on a recently proposed two-parameter generalized entropic functional. Unlike constructions that phenomenologically impose modified entropy-area relations at the horizon level, the present framework is rooted in a microscopic entropy functional and the corresponding microstate counting. For bounded systems, the entropy acquires a generalized holographic scaling with two independent area contributions, recovering the Bekenstein-Hawking entropy in the appropriate limits. Implementing this entropy within the holographic principle, we derive a generalized dark energy density containing two distinct holographic sectors, naturally embedding standard holographic dark energy and $\Lambda$CDM as limiting cases. We analyze the cosmological evolution for both Hubble and future event horizon cutoffs and show that the model successfully reproduces the matter-to-dark-energy transition. The two entropic exponents enrich the dynamics, allowing for quintessence-like behavior or phantom regimes, while remaining compatible with the standard thermal history of the Universe.