Papers

We present a study that combines remote-sensing and in-situ observations of coronal mass ejections (CMEs) interacting with the nearby heliospheric current sheet (HCS). The sequence of eruptive events under study culminates in the largest directly observed flare of solar cycle 25 on 3 October 2024, producing a fast halo CME. Their source region can be linked to a so-called nested active region (or active longitude) that persisted over several solar rotations. Such long-lived regions reflect deep-seated magnetic structures that shape the global magnetic field configuration. By applying the drag-based CME propagation model, we connect the near-Sun observations from several CMEs during that activity period with in-situ measurements. While one of the CMEs propagated on the opposite side of the HCS from Earth, and therefore did not produce in-situ signatures near Earth, we detect, over the period October 5-10, 2024, a complex of HCS and CME structures propagating together with a shock ahead of them. The HCS seems to be locally replaced by the CME signatures, leading to a long-duration sector reversal of more than 48 hours. This event highlights the intrinsic connection between solar surface structures, the global magnetic field, and the evolution of complex eruptive events.

Identifying remnant radio-loud active galactic nuclei (AGNs) is challenging due to their diverse morphological and spectral characteristics. Using three-dimensional hydrodynamic simulations of 15 radio galaxies, we investigate how the spectral evolution of remnants depends on progenitor power, active lifetime, environment, and underlying dynamics. The simulations span low-density group and high-density cluster environments re-gridded from smooth-particle-hydrodynamic cosmological simulations. The resulting remnants exhibit a wide range of morphologies, from amorphous structures to double-lobed forms. We find that jet power correlates with the spectral slope. As the remnant lobes evolve, we find surface brightness depends strongly on environment: group remnants are systematically dimmer and more amorphous than cluster remnants, highlighting a potential observational bias against these low-surface-brightness sources. In our models, we estimate that the peak surface brightness of a low-redshift, 50 Myr-old remnant from a low-power progenitor in a 10^{13} M_sun group environment should be routinely detectable at the 3{\sigma} level with LOFAR, although 20-30% of the emission would remain undetectable within a reasonable integration time. We find young remnants exhibit low-frequency (150-1400 MHz) spectral indices that overlap with active sources, and follow a consistent and established spectral-evolution sequence: significant curvature ({\alpha}_{1400}^{6000} - {\alpha}_{150}^{1400} > 0.5) develops before an ultra-steep low-frequency index ({\alpha}_{150}^{1400} > 1.2). The results presented in this work are intended as a reference point for current and upcoming low-frequency studies of radio remnants.

Recent baryon acoustic oscillation (BAO) distance measurements, when combined with Cosmic Microwave Background (CMB) observations in the $\Lambda$CDM framework, lead to a preference for negative neutrino masses. We investigate whether this neutrino mass anomaly can be alleviated by a class of astrophysically motivated reionization histories. Using a frequentist analysis, we find that some reionization histories can move the best-fit value of $\sum m_\nu$ to a positive value and bring $\sum m_\nu\simeq0.06~{\rm eV}$ into the 95% confidence interval. To separate the effect of the total optical depth from that of the details of the reionization history, we compare a high-$\tau$ history with a two-step tanh-like reionization history of the same $\tau$. The resulting $\Delta\chi^2(\sum m_\nu)$ profiles are nearly identical. This indicates that the effect is mainly driven by the total optical depth, while the details of the reionization history play only a minor role.

This chapter reviews applications of machine learning (ML) to redshifted 21 cm cosmology, focusing on cosmic dawn, the Epoch of Reionization, and SKA-Low science. The redshifted 21 cm line directly probes diffuse neutral hydrogen, but the measured signal is not a simple astrophysical observable: density, ionization, heating, radiation backgrounds, foreground treatment, and instrumental response are coupled. The chapter first summarizes the physical ingredients needed in later sections, including the global signal, spatial fluctuations, morphology-sensitive statistics, and the 21 cm forest. It then discusses the main barriers to interpretation: bright foregrounds, radio-frequency interference, ionospheric and calibration effects, incomplete sampling, and the cost of forward modeling in high-dimensional parameter spaces. ML applications are organized by their role in the analysis chain. Observation-domain methods act on contaminated data products; theory-domain methods accelerate or compress forward modeling; and inference-domain methods connect observables to astrophysical and cosmological parameters.

We present a comprehensive analysis of Vela X-1 using two new NuSTAR observations, placed in the context of four earlier datasets obtained between 2012 and 2020. The energy-resolved pulse profiles demonstrate a significant transformation from an asymmetric structure at low energies to distinct double peaks above 12 keV, whereas the pulse fraction escalates with photon energy but decreases with flux. Broadband spectra validate the Fe K alpha emission line and disclose both fundamental and harmonic cyclotron resonant scattering characteristics (CRSF). We observe no substantial link between CRSF energies and luminosity, contrary to previous findings; rather, the photon index and folding energy demonstrate distinct anti-correlations with flux, aligning with sub-critical accretion and increased Comptonization in the accretion column. Our results provide the first clear evidence that the harmonic CRSF in Vela X-1 does not follow the long-term decay previously claimed. The fundamental line energy also displays an irregular evolution, without a clear monotonic trend. Notably, the harmonic-to-fundamental energy ratio departs from the canonical value of two, suggesting that the line-forming regions are located at different heights within the accretion column. These results provide new constraints on the accretion geometry and magnetic field topology of Vela X-1, highlighting the importance of continued monitoring with current and future X-ray observatories.

Recent advances in astronomical observations have ushered in an era of remarkable discoveries. We now probe the Universe through multi-messenger signals, image the sky with unprecedented depth and resolution, and investigate individual sources using powerful large-aperture telescopes. Yet, a critical gap persists: the lack of wide-field, highly multiplexed spectroscopic capabilities needed to fully exploit the wealth of imaging data from current and upcoming surveys. In this review, we trace the historical development of large optical telescopes and spectroscopic surveys, assess the capabilities of ongoing and near-future facilities, and motivate the need for next-generation Stage-V spectroscopic experiments. As a representative example, we present the MUltiplexed Survey Telescope (MUST), the first Stage-V spectroscopic facility currently under construction. MUST is a 6.5-meter telescope designed to obtain optical spectra for over 20,000 targets simultaneously within a $\sim$5 deg$^2$ field, using a modular focal plane populated with 6.2-mm pitch fiber-positioning robots. Over an 8-year survey in the 2030s, MUST aims to build the most comprehensive 3D spectroscopic map of the Universe to date, measuring redshifts for over 100 million galaxies and quasars and opening new windows into cosmology, Galactic structure, and time-domain astrophysics.

Aims. We constrain the dynamical perturbations of the spiral arms in the Milky Way disk, based on the non-axisymmetric streaming motions of RGB stars revealed by APOGEE and Gaia. Methods. We develop a revised steady-state radial-velocity response model that incorporates both the \(V_{R,\sin}\) and the dynamically important \(V_{R,\cos}\) components for a two-armed logarithmic spiral potential. The model is validated using orbit integrations with AGAMA and Bayesian parameter recovery with dynesty, and is applied to the smoothed two-dimensional radial-velocity field of RGB stars while accounting for Lindblad and corotation resonances. Results. The revised model reproduces the phase and amplitude of the mock radial-velocity field to the \(\sim2\%\) level, substantially improving upon earlier \(V_{R,\sin}\)-only formulations. Applied to the observational data, it yields a robust pitch angle of \(p \simeq 10^\circ\) and a local surface density contrast of \(\xi \simeq 5\)–\(18\%\) at the solar radius. The radial scale length is less well-constrained (\(h_{R,1} \simeq 40\)–\(50\,\mathrm{kpc}\)) due to intrinsic parameter covariance. Resonance effects strongly shape the velocity field, thus affecting the fitting: the radial velocity becomes extremely large near the Lindblad resonances, whereas it vanishes close to the corotation resonance. Conclusions. Our results demonstrate that including both the \(V_{R,\sin}\) and \(V_{R,\cos}\) terms is essential for a physically consistent interpretation of stellar streaming motions induced by a spiral potential. The observed kinematics constrain the spiral pattern speed to \(\Omega_{p} \approx 10\)–\(20\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}\).

Intensity ratios of aromatic emission features are widely used to diagnose the size and ionization state of polycyclic aromatic hydrocarbons (PAHs) in astronomical environments. However, PAHs are known to typically carry aliphatic side chains, a structural feature that may compromise the reliability of traditional diagnostic methods. This study systematically investigates the effects of aliphatic components on the aromatic emission properties of PAHs. Based on theoretical data from the NASA Ames PAH IR Spectroscopic Database, we compare the emission behavior of purely aromatic PAHs with those containing aliphatic substituents, revealing that aliphatic functionalization may modify the intensity ratio of the 11.2 $\mu$m band relative to the 7.7 $\mu$m and 3.3 $\mu$m bands. This potentially leads to misidentification of their ionization state if molecular structural effects are neglected. Further analysis indicates that the impact of aliphatic components on diagnostic band ratios strongly depends on PAH size: small PAHs exhibit significant emission ratio shifts, deviating from traditional size/ionization trends, while larger PAHs are minimally affected. Despite these shifts, the classic $(I_{11.2/7.7})$ versus $(I_{11.2/3.3})$ diagnostic grid remains largely applicable to mixed aromatic-aliphatic PAHs, although some systematic calibration may be needed. Our findings emphasize the necessity for caution when interpreting PAH band ratios in aliphatic-rich environments, as variations in PAH molecular composition may distort inferences about physical conditions.

Wind is a key mechanism for supermassive black hole (SMBH) feedback to their host galaxies. In tidal disruption events (TDEs), black holes spend most of their time accreting at highly sub-Eddington rates, implying that feedback from persistent sub-Eddington winds could be significant. We investigate the effects of black hole mass, viscosity parameter and stellar debris temperature on the properties of winds from hot accretion flows in TDEs. We find that more massive black holes yield a higher accreted fraction and launch faster winds, while the debris temperature has a negligible influence on the accretion flow. For $\alpha=0.1$, the mildly-relativistic unbound winds ($\sim 0.1c$) are launched predominantly from the outside of the accretion flows along the equatorial plane, with a kinetic energy of $\sim10^{-4}L_\mathrm{Edd}$. In contrast, convective bound outflows dominate for $\alpha=0.01$, which differs from the true winds typically seen in active galactic nuclei and X-ray binaries. Potential applications for explaining delayed radio brightening in TDEs at $\sim10^3$ days and for searching for intermediate-mass black holes through radio and X-ray surveys are also discussed.

Modified Newtonian Dynamics (MOND) is a paradigm that can do away with dark matter at galaxy scales, but displays a residual missing mass discrepancy in galaxy clusters. Prompted by the updated JWST-based gravitational lens model of the Bullet Cluster, I confirm here that this cluster exhibits the same residual missing mass discrepancy as other clusters of similar mass in the MOND context. Moreover, this missing mass should be mostly collisionless, since it is centred on the galaxies of the Bullet Cluster.

Since the early 1990s, there has been significant interest in the physics of dusty plasmas, which has now become a new discipline in plasma science. Dusty plasma exhibits new and unusual behaviour, and provides a possibility for modified or entirely new collective modes of oscillations, instabilities as well as coherent nonlinear structures.\\ First, a review of the important recurring terms – The Cosmic Waves (CRs), the Alfven Waves (AWs), and the associated charged dust grains is presented. Starting from the basic composition of the CRs to their scattering mechanism, along with the different modes of scattering, is presented, along with the modes of confinement and a precise definition of each term. The paper also includes some useful diagrams and brief notes from the references. \\ Gravitation plays a significant role in the collapse of matter and the formation of cosmological structures. Unlike a static universe, this paper investigates the Jeans instability in a radiation-pressure-dominated expanding universe using the Einstein-de Sitter model for Euclidean geometry with zero curvature ($\kappa=0$). The fluid model for an expanding universe is constructed, and by taking small perturbations, the perturbed fluid equations are obtained. The dispersion relation of gravitational instability is derived using plane-wave solutions. In the static case (Newtonian cosmology), the classical Jeans instability criterion is revisited and modified in an expanding universe. The critical Jeans wave number of perturbations to excite Jeans instability depends upon the time-dependent expansion factor $S(t)$. It is found that short-wavelength perturbations are expected during the inflationary period of big bang cosmology, which are responsible for gravitational collapse and the formation of galaxies.

Tidal fragmentation at perihelion and nontidal fragmentation elsewhere cause the orbital distribution of Kreutz sungrazers of all sizes to be extremely complicated and highly nonuniform. Among the features are (largely fortuitous) clusters of bright (naked-eye) objects and clumps of dwarf objects (often closely genetically related, as their detection primarily by the SOHO coronagraphs suggests) on the one hand; and both spectacular and less brilliant sibling sungrazers, whose perihelion times are scattered over centuries, on the other hand. Investigation of four fragment nuclei of the Great September Comet of 1882, the products of a perihelion breakup of the comet's original nucleus, showed that their orbital periods followed a distinct pattern, which likewise applied to other tidally split sungrazers and was characterized by a specific value of the second difference of parameter u_frg of neighboring fragments' centers of mass. The algorithm has a potential for the prognostication of bright Kreutz sungrazers over the rest of the 21st century and beyond. However, because of its as yet unverified empirical character, the utmost caution should be exercised when applying the procedure.

The leading-trailing asymmetry in lunar crater distribution provides a critical record of inner solar system dynamics, yet the long-standing discrepancy between the observed higher asymmetry and lower theoretical predictions indicates a gap in our understanding of the impactor population. This paper hypothesizes that lunar impact ejecta, which can enter Earth-like orbits and return, constitute a previously unaccounted-for component. Through numerical simulations, we find that ~25% of escaped ejecta will re-impact the Earth-Moon system within 3 Myr, with about 1.2% striking the Moon. Crucially, these lunar impacts exhibit an extreme leading-trailing asymmetry with a ratio of 5.9. Our results indicate that lunar ejecta, if comprising ~15% of total impactors, can fully explain the observed asymmetry, leading to their recognition as active agents shaping the lunar impact record. This work provides new constraints for understanding the impact environment of the Earth-Moon system, with direct relevance to the interpretation of lunar geology, the transport of lunar material to Earth, and ongoing space exploration missions.

We present RAYTHEIA, a high-performance reverse ray-tracing algorithm designed to efficiently solve three-dimensional direction-dependent equations in astronomical simulations. The algorithm uses a dual-grid framework in which the native simulation mesh – serving as the source grid for ray emission – and an adaptive mesh refinement (AMR) Cartesian contribution grid are constructed for efficient ray-walking and contribution accumulation. The core of the algorithm integrates a leaf-only linear-octree data structure to reduce memory overhead, the digital differential analyzer (DDA) traversal method to efficiently determine the ray-walking path, Morton Code indexing to fast leaf cell lookup during traversal, and the slab method to analytically compute the path length. Furthermore, RAYTHEIA employs a hybrid (MPI/OpenMP) distributed parallel framework with a chunk-to-chunk communication strategy, achieving exceptional, near-ideal linear speed-up ratio and delivering high-end performance. We integrate RAYTHEIA with the 3D-PDR code to solve the complex chemistry and radiation transfer in photodissociation regions (PDRs). This allowed the modelling of three-dimensional PDR chemistry in a turbulent, star-forming cloud at an unprecedented resolution of $512^3$ grid cells. The algorithm demonstrates accuracy and convergence even at low angular resolutions. We further showcase the capabilities of RAYTHEIA by producing high-resolution synthetic emission maps of key diagnostic lines of a star-forming region capturing physical effects such as [O I] $63\mu$m self-absorption, measuring the [C I]-bright but CO-dark molecular gas, and deriving a CO-to-H$_2$ conversion factor in agreement with observations.

Resonant populations of trans-Neptunian objects serve as crucial dynamical archives for unraveling the early migratory history of the Solar System. A quantitative assessment of the capture efficiency into various mean motion resonances (MMRs) during migration is essential for understanding the origins of these populations, constraining migration parameters, and reconstructing of the primordial planetesimal disk. Using numerical simulations, this study systematically investigates the capture capability of exterior MMRs during Neptune's outward migration in a planar model. For a specific p:q MMR, the small bodies can be captured only when their eccentricities surpass a certain threshold, which increases with faster migration rates, greater distances of MMRs, and higher resonance orders. On the other hand, as long as a particle's eccentricity is suitable, its capture efficiency shows little dependence on the migration rate; instead, it mainly depends on the p value and heliocentric distance, decaying exponentially as either parameter increases. Based on our simulation results, we derive for the first time a simple empirical expression to calculate eccentricity threshold and the capture efficiency. This research provides a systematic quantitative framework for understanding capture into Neptunian MMRs during migration. Future integrations of more comprehensive observational data will facilitate a more precise reconstruction of the Solar System's early dynamical evolution.

Nonlinear Force-free Field (NLFFF) models are widely used to investigate coronal magnetic field structure in solar active regions, but methods to validate them remain limited. Here, we use Gaussian separation, recently applied to solar vector magnetogram data, to assess the accuracy of NLFFF models constructed with two methods: optimization and the current-field iteration (CFIT) implementation of the Grad-Rubin method. Gaussian separation partitions the photospheric vector magnetic field into three components associated with currents flowing below, above, and passing through the photosphere, respectively. Comparing the photospheric field components due to coronal currents in an NLFFF model with those in the original vector magnetogram data provides a check on the accuracy of the model's coronal currents. We consider NLFFF models constructed for the active region AR 11429. The photospheric signatures of coronal currents in both the models and the vector magnetogram data indicate currents flowing above and parallel to central, sheared polarity inversion lines (PILs), consistent with other recent studies. We find that while both models reproduce the coronal current signatures along the upper section of the main PIL, the CFIT model significantly alters the signature of a flux rope along the lower section of the PIL, including shifting its positive-polarity footpoint. These differences arise from modifications to the vector magnetogram boundary data when solving the NLFFF equations, and from the assumptions underlying the models. We propose Gaussian separation as a useful tool to validate coronal magnetic field models, in addition to existing methods.

We simulate star formation and star cluster assembly inside a molecular cloud with parameters we derive directly from observations of the Aquila Rift. We model the evolution of stars and gas together while resolving close encounters between stars, the formation of new stars, and stellar feedback to follow cluster formation up to the expulsion of the surrounding gas. We find that star formation takes place in clumps spaced unevenly along Serpens South and that these clumps accrete surrounding gas to grow and form new stars. Gas flows along the filament promote the merger of these clumps into a star cluster inside the Serpens South filament. The imprints of these mergers are seen in the dynamics of the Serpens South cluster in the form of velocity space anisotropies, cluster rotation, and cluster expansion. Before gas is removed from the simulation, the Serpens South cluster merges with the nearby cluster W40 non-monolithically resulting in a fractal cluster at the end of the simulation. The dynamics inherited from the mergers throughout the simulation are still seen in the final bound stellar system after the gas has been removed. We compare these results with recent observations of Milky Way clusters to comment on their formation histories. We also study how our results change when lowering the mass resolution of our simulation and removing observations of dense gas tracers from our initial condition setup. Each of the three simulations result in different final cluster configurations pointing towards the importance of gas in cluster assembly.

Recent JWST observations suggest that galaxies at z > 9 may be dominated by low-opacity SNe-produced dust before efficient ISM grain growth is established. This transition in dust origin and opacity could explain both the prevalence of galaxies with extremely low dust attenuation and the excess of UV-bright galaxies relative to most pre-JWST predictions. We investigate whether this transition, combined with evolving star-formation efficiency, can reproduce these observed properties. We develop a physically motivated attenuation framework combining (i) extinction laws for reverse-shock-processed SNe dust, (ii) metallicity- and dust-to-metal-dependent opacity scalings, and (iii) porous radiative-transfer geometries allowing partial UV-photon leakage. Unlike outflow-driven scenarios requiring large-scale gas evacuation, our approach preserves gas reservoirs while reducing effective UV opacity through dust composition and geometry. We introduce extinction-based, gas-based, and hybrid attenuation prescriptions linking SNe-dominated and ISM grain-growth dust regimes. We find that the observed A_FUV-M_star relation at z > 9 is best reproduced for an intrinsic FUV dust opacity kappa_UV(dust)=10^3 - 10^4 cm2/g, characteristic of low-opacity SNe dust, naturally producing very low attenuation even in gas-rich galaxies. This regime reproduces galaxies with extremely low dust attenuation (GELDAs), which dominate observed samples at z > 9. Applied to intrinsic UV luminosity function models, our SNe-dominated and hybrid prescriptions mainly suppress the brightest galaxies, bringing predictions into agreement with JWST measurements without requiring extreme star-formation efficiencies or dust-free interstellar media. Our results suggest that the UV-bright galaxy excess at z > 9 reflects a transition in dust origin and opacity during the earliest phases of galaxy evolution.

We investigate the morphology and structural properties of 246 massive compact quiescent galaxies (MCGs; $\log M_{\star} \sim 10$-$11$, $\sigma_{\mathrm{e}} \sim 150$-$350\,$km\,s$^{-1}$, $R_{\mathrm{e}} \sim 0.7$-$2.5\,$kpc) at $z \sim 0$, selected as outliers in the stellar mass-velocity dispersion and velocity dispersion-size relations, using $g$-, $r$-, and $i$-band Hyper Suprime-Cam images. We compare them to a control sample of average-sized quiescent galaxies (CSGs) matched in stellar mass, star formation rate, redshift, and $g-i$ color. Both samples are dominated by S0 galaxies, comprising $93\%$ of MCGs and $71\%$ of CSGs, while ellipticals account for $4\%$ and $11\%$, respectively. The fraction of interacting or morphologically disturbed systems is low in both samples ($13\%$ for MCGs and $16\%$ for CSGs). Multi-component decompositions of the $g$- and $r$-band images show that $75\%$ of MCGs require a three-component model (bulge, disk, and envelope), while $21\%$ are best fit by two components and $4\%$ by a single Sérsic profile. Two-component MCGs are preferentially low-inclination systems, suggesting that the three-component fraction represents a lower limit. In contrast, only $7\%$ of CSGs exhibit a comparable three-component structure. Bars are present in $29\%$ of CSGs but are absent in MCGs. For three-component systems, MCGs and CSGs have similar bulge ($R_\mathrm{e}=0.39$ vs.\ $0.45$\,kpc) and envelope ($R_\mathrm{e}=6.4$ vs.\ $5.8$\,kpc) sizes, while MCG disks are significantly more compact ($R_\mathrm{e}=1.9$ vs.\ $3.3$\,kpc). The envelope component shows a broad ellipticity distribution ($\epsilon_\mathrm{Envelope} \sim 0.0$-$0.6$), which we interpret as corresponding to either a stellar halo or a thick disk.

We investigate the thermal evolution of isolated neutron stars containing hyperon–mixed kaon–condensed matter, focusing on the role of proton superconductivity. The equation of state utilized for cooling calculation is based upon the minimal relativistic mean–field framework supplemented by chiral SU(3) dynamics for kaon condensation with an additional component on the three-baryon force, which ensures stiffness at high densities enough to meet astrophysical constraints on neutron-star masses and radii. We show that the nucleonic direct Urca processes operate at relatively low stellar masses ($M \gtrsim 1.3\,M_\odot$), erasing any observable signature of strangeness in the absence of superfluidity. However, if the proton $^1{\rm S}_0$ superconductivity works, because of suppression of fast neutrino cooling processes, the cooling scenario could become relevant with the strangeness, depending on the density regions of the pairing gap. In particular, if the proton superconductivity is so strong in high-density regions ($T_{c,p}\sim10^{10}~{\rm K}$), the nucleon and hyperon direct Urca processes shut down, which makes the kaon-induced Urca processes dominant in massive neutron stars. This scenario is in good agreement with several cold isolated neutron stars identified recently recently, such as the Vela Jr., PSR J0205+6449, and PSR B2334+61. Hence, we suggest that strong proton superconductivity can render kaon condensation observationally visible through cold neutron-star observations, providing a potential signature of strangeness in dense matter.

We study in detail a scenario in which the inflaton scalar field couples to both a visible sector (VS) and a hidden sector (HS). The VS is assumed to contain the Standard Model (SM), while the HS contains a dark matter (DM) candidate. We are in particular interested in a scenario in which the inflaton decays dominantly into the HS degrees of freedom. The DM candidate is taken to be a dark Dirac fermion $\chi$, coupled to a massive dark photon $\gamma'$, a popular model for a HS also known as Dark QED. The inflaton decays into particles of both sectors generate an initial asymmetry between the SM and HS fermion abundances, which we model as being proportional to the ratio of effective Yukawa couplings, $y$ and $y'$. We pay particular attention to the process of thermalisation of the HS, with temperature $T'$, as a function of $y'$ and $\alpha'$, the HS fine structure constant. We investigate the several possible ways of producing the observed DM relic abundance, and their interplay with the reheating of the HS and the transfer of energy between the HS and the VS. Key results, beyond the systematic character of our analysis, include: a new mechanism for DM production, which occurs when DM particles annihilate while still being produced by the inflaton decay; a study of the temperature ratio $\xi = T'/T$ and its relation with the initial energy asymmetry between the HS and VS, as parameterized by $\xi_i = \sqrt{y'/y}$; a reassessment of the domain of viable DM candidates, taking into account the constraints set by unitarity and the thermalisation of the HS, accounting for the LPM effect; and, in cases where the HS does not reach thermal equilibrium, an analysis of how non-thermal DM production fits within the domain of thermal DM candidates.

Neutron star (NS) low-mass X-ray binary (LMXB) systems depict a variety of X-ray spectral and timing features, which can be useful to probe the accretion-ejection mechanism in the strong gravity regime. Here, we study the relatively unexplored and faint NS LMXB 1A 1246-588, which is also an ultra-compact X-ray binary (UCXB) with a white dwarf donor. We investigate its temporal and spectral behavior using pointed NICER and AstroSat observations, supported by long-term MAXI/GSC monitoring. The MAXI light curve shows modest, recurrent outburst-like enhancements, providing the long-term flux context for interpreting the pointed observations. During the AstroSat observations in 2017, the source exhibits an absorbed 0.4-20 keV flux of $(1.18 \pm 0.02)$ x $10^{-10}$ $erg$ $cm^{-2}$ $s^{-1}$, while during the NICER observations in 2019, it spans an absorbed 0.5-10 keV flux range of $(0.7-3.7)$ x $10^{-10}$ $erg$ $cm^{-2}$ $s^{-1}$ and traces an atoll-like pattern in the hardness-intensity diagram. Broadband spectral modeling shows that the emission is well described by a soft blackbody and a hard Comptonized component, with no statistically required multicolor disk contribution. The blackbody temperature increases from 0.28 to 0.39 keV, with an emitting radius consistent within 6.9-13.8 km, while the Comptonization photon index varies from 1.8 to 2.3. We find that the observed spectral-state evolution is driven by a redistribution of accretion power between thermal emission from the NS boundary layer and Comptonized emission, consistent with atoll-type behavior. These results provide the first quantitative, multi-epoch view of accretion-state evolution in 1A 1246-588, revealing systematic changes in the thermal boundary-layer emission and the Comptonizing region in this UCXB system.

Eclipse shadow bands are transient intensity patterns that can appear on the ground near solar eclipse totality. This study presents a reusable image-processing framework for analyzing shadow-band video recordings collected with consumer-grade cameras. The framework quantifies band prominence, band orientation, band separation, band velocity and power spectral density from video recordings. Applied to two eclipse datasets, the method detected statistically significant shadow-band activity during eclipse windows that align with the scintillation theory for shadow bands. The findings also indicate that eclipse shadow bands can consist of simultaneous superimposed band modes with different orientations and velocities. This study demonstrates that consumer-grade cameras can support quantitative analysis of shadow bands and may support future observational and atmospheric studies

We investigate the optical appearance of a Kerr-Bertotti-Robinson (Kerr-BR) black hole illuminated by a geometrically and optically thin accretion disk. Instead of using a phenomenological power-law emissivity, we adopt a magnetically driven synchrotron emissivity proxy coupled to the local electromagnetic environment. With a backward ray-tracing framework, we examine the effects of the spin $a$, magnetic parameter $B$, and observer inclination $\theta_O$ on the ray-classification maps, redshift distributions, and specific-intensity images. We show that the ISCO position is modified by both $a$ and $B$, and that rapidly rotating prograde configurations can develop an additional model-dependent inner cutoff when the magnetically dominated approximation underlying the emissivity prescription ceases to be applicable. High-resolution one-dimensional intensity profiles further separate the direct image, the $n=1$ lensing-ring contribution, and the higher-order $n\geq 2$ photon-ring subimages, while quantifying the Doppler-induced brightness asymmetry. Retrograde disks exhibit a wider emission-depleted central region because of the outwardly shifted ISCO, making the higher-order lensed components more clearly distinguishable from the direct emission. These results show that the disk inner boundary and the magnetic-field-dependent emissivity can substantially influence the observable appearance of Kerr-BR black holes.

We study the postinflationary dynamics of an Einstein-Cartan-Holst gravity-motivated inflationary scenario, known as Einstein-Cartan pseudoscalaron inflation, coupled to a type-I seesaw extension of the Standard Model with three heavy right-handed Majorana neutrinos. In particular, we show that nonthermal leptogenesis emerges as a necessary and self-consistent mechanism for generating the observed baryon asymmetry of the Universe, mainly because of the universal coupling of the inflaton to the additional heavy Majorana fermions. The resulting framework provides theoretical predictions that are fully compatible with the latest cosmological constraints from the Cosmic Microwave Background, Baryon Acoustic Oscillations, and Big Bang Nucleosynthesis, as well as with neutrino oscillation experiments, for a wide range of the fundamental Barbero-Immirzi model parameter $\gamma$, which controls the inflationary and postinflationary phases. In particular, for $\gamma \sim -1/100$ and a lightest Majorana-neutrino mass of order $10^{13}$ GeV, we find a scalar spectral index $n_s \sim 0.970$, a tensor-to-scalar ratio $r \sim 0.004$, for a number of e-folds before the end of inflation $N_e \lesssim 60$, and a baryon-to-entropy ratio $n_B/s \sim 8.7 \times 10^{-11}$.

The origin of the highly ionized iron emission (Fe XXV at $6.7\,\mathrm{keV}$) characterizing the Galactic ridge X-ray emission (GRXE) remains a fundamental puzzle in high-energy astrophysics. Although the GRXE continuum is largely resolved into discrete populations of cataclysmic variables and coronally active stars, these sources exhibit Fe XXV equivalent widths significantly lower than that of the total GRXE, leaving the intense iron line emission unexplained. In this work, we cross-correlated the XMM-Newton survey of the inner Galactic disk with Gaia DR3 astrometry to identify and characterize hard X-ray sources ($>2\,\mathrm{keV}$) with reliable stellar counterparts. We selected 107 X-ray sources located within the red giant branch of the color-magnitude diagram, many of which are verified long-period variables. These sources exhibit high X-ray luminosities ($L_{\mathrm{X}} \approx 10^{31}$–$10^{33}\,\mathrm{erg~s^{-1}}$), significantly exceeding the typical coronal saturation levels of single giants. Their X-ray spectra are notably harder than those of quiescent stellar coronae, with plasma temperatures reaching up to $kT \approx 6\,\mathrm{keV}$ and a prominent emission feature at $\sim 6.7\,\mathrm{keV}$. The combination of high $L_{\mathrm{X}}$, hard spectra, and intense Fe XXV emission identifies this population as accretion-powered binaries associated with late-type giants. Our analysis demonstrates that this population contributes $\sim 20\%$ of the total GRXE continuum and $\sim 40\%$ of its iron line emission, providing a key component to resolving the Galactic X-ray background puzzle.

In the last year, several pieces of evidence have pointed to a possible deviation from the standard cosmological model, $\Lambda$CDM. The recent work by the Dark Energy Survey (DES) collaboration reports a preference in the ballpark of $3\sigma$ in favor of dynamical dark energy against the standard cosmological model. For that, it used its final analyses of Baryonic Acoustic Oscillations and type Ia Supernovae, both sensitive to the expansion history of the Universe, in combination with the Cosmic Microwave Background (CMB) from Planck. This adds to the growing debate about the nature of dark energy. Published as a Perspective in Nature Astronomy in August 2025.

Spectropolarimetric interpretation of chromospheric lines requires solving the radiative transfer problem under non-local thermodynamic equilibrium (non-LTE) conditions. This means computing atomic-level populations self-consistently with the radiation field. While traditional inversion codes employ 1.5D approximations, they neglect horizontal radiative transfer, which can be significant near magnetic structures and in the chromosphere. We present a method to solve 3D atomic-level populations using Graph Neural Networks (GNNs), extending prior 1.5D work to the full 3D domain. By discretizing the solar atmosphere as a directed graph, in which nodes encode physical properties and edges encode spatial distances, an Encode-Process-Decode GNN propagates information to efficiently capture radiative coupling. The network is trained on a Bifrost simulation using Ca II populations from Multi3D as ground truth. The trained GNN accurately predicts populations of the five-level Ca II atom plus continuum. Correlations exceed 0.99 in the photosphere and mid-chromosphere; errors in the upper chromosphere remain unbiased. Inference is $\sim 10^6$ times faster than traditional iterative solvers. Spectral synthesis of the Ca II 8542 Å\ line yields intensity profiles with $< 2 \%$ mean residuals relative to the full 3D solution. This framework bypasses the computational bottleneck of iterative solvers while preserving essential non-LTE physics, including horizontal transfer, paving the way toward routine 3D non-LTE inversions.

Assuming that the compact source at the Galactic center, Sagittarius A*, is a massive boson star, we fit the near-infrared flare astrometry data. We consider 12 discrete boson star configurations and model the flare as a hotspot on a circular equatorial orbit. The analysis is performed in a Bayesian framework using nested sampling, yielding the marginal posterior distributions of all parameters as well as the Bayesian evidence for each model. For comparison, the same procedure is applied to a Schwarzschild black hole. The resulting Bayesian evidence values differ only marginally between the boson star and black hole cases, and the well-determined mass of Sgr~A* (${\sim}4.296\times 10^6\,M_\odot$) falls within the 68% highest density interval in every configuration. We conclude that, under current near-infrared astrometric constraints and within the considered parameter ranges, a massive boson star and a Schwarzschild black hole remain statistically indistinguishable as the compact object at the Galactic center.

In this Paper we consider twisted accretion disks in supermassive binary black hole by analytical and numerical means. It is assumed that the disk orbiting around the more massive rotating component and that the disk rings are inclined with respect to the orbital plane. We use orbital parameters of the binary often employed in the precessing massive (PM) model of the well-known blazar OJ 287. Unlike our previous investigation of a similar problem, here we consider disks with both small and relatively large relative thicknesses $\delta=h/r$, where $h$ is the disk's height at a typical radius $r$, as well as a range of values of the viscosity parameter, $\alpha$, including the cases when $\alpha \lesssim \delta$. Similar to our previous results, we find that the twisted disk relaxes to a quasi-stationary state in the frame precessing with the Lense-Thirring frequency of the orbit. However, its shape is qualitatively different from that corresponding to the case of $\delta=10^{-3}$ and $\alpha=0.1$ considered in our previous work. In a disk with $\delta=10^{-3}$ but $\alpha \le 2\cdot 10^{-2}$, we find the new effect of generation of a twisting spiral wave near the resonance between a forcing frequency associated with the presence of the secondary and the Lense-Thirring frequency of a particular disk ring defined in the precessing frame. We propose an analytic theory of it, which is in a good agreement with our numerical results. This effect leads to multiple crossings of the orbit with the disk per one orbital period, which contradicts the PM model. When $\delta \gtrsim 0.1$, a typical disk's inclination within the orbit of the binary is smaller than that of the orbit which results in only two crossings of the orbit with the disk per one orbital period. We suggest that the additional heating of the disk gas by the secondary-disk collisions may result in $\delta \sim 0.1$.