Papers for Thursday, May 15 2025

Sunspots are areas of strong magnetic fields driven by a convective dynamo. Rudolf Wolf devised his Sunspot Number (SN) series to describe its variation with time. Most other solar phenomena vary in concert with SN, in particular the microwave radiation. During the interval 2014-Jul to 2015-Dec, the variation of the SN with the microwave flux (e.g. F10.7) was anomalous, although microwave flux did not exhibit any anomaly with respect to the solar magnetic field measured by the HMI instrument on the SDO spacecraft. This points to a possible problem with the derivation of the SN using the 148,714 reports of sunspot observations received by the World Data Center Solar Influences (SIDC/SILSO at ROB) since 2011 (when the last public raw data was released). The unavailability of the raw data since then, contrary to avowed open data policy of ROB, prevents independent assessment and possible correction of the anomaly. We urge ROB to make the raw data available.

X-ray binaries (XRBs) display spectral state transitions that are accompanied by substantial changes in the hardness, luminosity, and structure of the accretion flow. We developed a GPU-accelerated cooling toolkit for general relativistic magnetohydrodynamic (GRMHD) simulations of accreting black holes that uses texture memory for fast retrieval of pre-computed values. The toolkit incorporates bremsstrahlung, synchrotron, inverse Compton radiation and Coulomb collision processes. We implemented our toolkit into a GRMHD code and used it to simulate a magnetically arrested disk in the context of the XRB low/hard state around a Kerr black hole. We explored the mass accretion rate in the $\sim (10^{-6}-0.3) \dot{M}_{\rm Edd}$ range, where $\dot{M}_{\rm Edd}$ is the Eddington accretion rate. Our simulations reveal that for low accretion rates ($\dot{M} \lesssim 0.01 \dot{M}_{\rm Edd}$), the flow settles into a geometrically thick, low-density, two-temperature hot accretion flow. At higher accretion rates, the flow turns into a cold single-temperature thin disk at $r_{\rm in} \gtrsim 50 r_g$. Inside, the disk breaks up into single-temperature thin filaments embedded into a two-temperature hot thick flow. Our GPU texture memory accelerated cooling prescription is $3-5$ times faster than the standard radiation M1 closure methods, and $\sim5$ times faster than storing the lookup table in global memory.

We explore the potential of the expanding ejecta method (EEM) as a cosmological probe, leveraging its ability to measure angular diameter distances to supernovae (SNe) with intensity interferometry. We propose three distinct applications of the EEM: (1) using Type IIP SNe as moderate-distance geometric anchors to calibrate Cepheids, replacing other local distance indicators; (2) directly calibrating Type Ia SNe, bypassing conventional calibration methods; (3) constructing a fully independent Hubble diagram with Type IIP (Type Ia) SNe, entirely decoupled from the traditional distance ladder. Incorporating realistic SN populations, we forecast a Hubble constant precision with next-generation intensity interferometers of $1.6\%$, $1.1\%$, and $9.3\% \,(3.6\%)$, respectively, for the three different proposed applications. Future intensity interferometry could yield improvements to $1.2\%$, $0.6\%$, and $1.5\%\,(0.4\%)$. The EEM thus offers a powerful geometric alternative for cosmic distance determination.

Supernovae, the explosive deaths of massive stars, create heavy elements and form black holes and neutron stars. These compact objects often receive a velocity at formation, a "kick" whose physical origin remains debated. We investigate kicks in Be X-ray binaries, containing a neutron star and a rapidly spinning companion. We identify two distinct populations: one with kicks below $10\,\rm{km}\,\rm{s}^{-1}$, much lower than theoretical predictions, and another with kicks around $100\,\rm{km}\,\rm{s}^{-1}$, that shows evidence for being aligned within 5 degrees of the progenitor's rotation axis. The distribution of progenitor masses for the two populations have medians around $2.3\,\rm{M}_\odot$ and $4.9\,\rm{M}_\odot$, corresponding to stars with birth masses of about $10\,\rm{M}_\odot$ and $15\,\rm{M}_\odot$. The second component matches the low-velocity mode observed in isolated pulsars. Combined with the known high-velocity component, which dominates isolated pulsars, this suggests three distinct kick modes. These results reveal previously unrecognized diversity in neutron-star formation.

For decades, studying quiescent galaxies beyond $z\sim1$ has been challenging due to the reliance on photometric spectral energy distributions (SEDs), which are highly susceptible to degeneracies between age, metallicity, dust, and star-formation history (SFH). Only recently has deep, rest-frame, optical spectroscopy made robust metallicity and age measurements possible, allowing us to empirically assess their effects on continuum shapes. To this end, we measure ages and metallicities of $\sim700$ massive ($10.2\lesssim\log(M_*/M_\odot)\lesssim11.8$), quiescent galaxies at $0.6\lesssim z\lesssim1.0$ from the Large Early Galaxy Astrophysics Census (LEGA-C) via continuum-normalized, absorption-line spectra, and compare with independent rest-frame $U-V$ and $V-J$ colours. Age increases along the quiescent sequence as both colours redden, consistent with stellar population synthesis (SPS) model predictions. Metallicity increases perpendicularly to the age trend, with higher metallicities at redder $U-V$ and bluer $V-J$ colours. Thus, age and metallicity behave differently in the $UVJ$ diagram. Moreover, this trend conflicts with SPS model expectations of increasing metallicity approximately along the quiescent sequence. Independent dynamical mass-to-light ratio trends also differ dramatically from SPS model predictions. These results demonstrate that relying on model fits to continuum shapes alone may lead to systematic biases in ages, metallicities, and stellar masses. The cause of these data-model disparities may stem from non-solar abundance patterns in quiescent galaxies or the treatment of evolved stellar phases in the models. Resolving these discrepancies is crucial, as photometric data remain central even with \textit{JWST}.

In this paper, we report on three of the largest (in terms of simulation domain size) and longest (in terms of duration) 3D general relativistic radiation magnetohydrodynamic simulations of super-critical accretion onto black holes. The simulations are all set for a rapidly rotating ($a_* = 0.9$), stellar-mass ($M_\mathrm{BH} = 6.62 M_\odot$) black hole. The simulations vary in their initial target mass accretion rates (assumed measured at large radius), with values sampled in the range $\dot{m}=\dot{M}/\dot{M}_\mathrm{Edd} = 1-10$. We find in practice, though, that all of our simulations settle close to a net accretion rate of $\dot{m}_\mathrm{net} = \dot{m}_\mathrm{in}-\dot{m}_\mathrm{out} \approx 1$ (over the radii where our simulations have reached equilibrium), even though the inward mass flux (measured at large radii) $\dot{m}_\mathrm{in}$ can exceed 1,000 in some cases. This is possible because the outflowing mass flux $\dot{m}_\mathrm{out}$ adjusts itself to very nearly cancel out $\dot{m}_\mathrm{in}$, so that at all radii $\dot{M}_\mathrm{net} \approx \dot{M}_\mathrm{Edd}$. In other words, these simulated discs obey the Eddington limit. We compare our results with the predictions of the slim disc (advection-dominated) and critical disc (wind/outflow-dominated) models, finding that they agree quite well with the critical disc model both qualitatively and quantitatively. We also speculate as to why our results appear to contradict most previous numerical studies of super-critical accretion.

Mass transfer in binary systems is the key process in the formation of various classes of objects, including merging binary black holes (BBHs) and neutron stars. Orbital evolution during mass transfer depends on how much mass is accreted and how much angular momentum is lost $-$ two of the main uncertainties in binary evolution. Here, we demonstrate that there is a fundamental limit to how close binary systems can get via stable mass transfer (SMT), that is robust against uncertainties in orbital evolution. Based on detailed evolutionary models of interacting systems with a BH accretor and a massive star companion, we show that the post-interaction orbit is always wider than $\sim10R_{\odot}$, even when extreme shrinkage due to L2 outflows is assumed. Systems evolving towards tighter orbits become dynamically unstable and result in stellar mergers. This separation limit has direct implications for the properties of BBH mergers: long delay times ($\gtrsim1 \rm Gyr$), and no high BH spins from the tidal spin-up of helium stars. At high metallicity, the SMT channel may be severely quenched due to Wolf-Rayet winds. The reason for the separation limit lies in the stellar structure, not in binary physics. If the orbit gets too narrow during mass transfer, a dynamical instability is triggered by a rapid expansion of the remaining donor envelope due to its near-flat entropy profile. The closest separations can be achieved from core-He burning ($\sim8-15R_{\odot}$) and Main Sequence donors ($\sim15-30R_{\odot}$), while Hertzsprung Gap donors lead to wider orbits ($\gtrsim30-50R_{\odot}$) and non-merging BBHs. These outcomes and mass transfer stability are determined by the entropy structures, which are governed by internal composition profiles. The formation of compact binaries is thus sensitive to chemical mixing in stars and may relate to the blue-supergiant problem.

Dwarf galaxies are uniquely sensitive to energetic feedback processes and are known to experience substantial mass and metal loss from their disk. Here, we investigate the circumgalactic medium (CGM) of 64 isolated dwarf galaxies ($6.0<$log(M$_*/M_{\odot}$)$<9.5$) at $z=0$ from the Marvel-ous Dwarfs and Marvelous Massive Dwarfs hydrodynamic simulations. Our galaxies produce column densities broadly consistent with current observations. We investigate these column densities in the context of mass and metal retention rates and the physical properties of the CGM. We find $48\pm11\%$ of all baryons within $R_{200c}$ reside in the CGM, with $\sim70\%$ of CGM mass existing in a warm gas phase, $10^{4.5}<T<10^{5.5}$ K that dominates beyond $r/R_{200c}\sim0.5$. Further, the warm and cool ($10^{4.0}<T<10^{4.5}$ K) gas phases each retain $5-10\%$ of metals formed by the dwarf galaxy. The significant fraction of mass and metals residing in the warm CGM phase provides an interpretation for the lack of observed low ion detections beyond $b/R_{200c}\sim0.5$ at $z\sim0$. We find a weak correlation between galaxy mass and total CGM metal retention despite the fraction of metals lost from the halo increasing from $\sim10\%$ to $>40\%$ towards lower masses. Our findings highlight the CGM (primarily its warm component) as a key reservoir of mass and metals for dwarf galaxies across stellar masses and underscore its importance in understanding the baryon cycle in the low-mass regime. Finally, we provide individual galaxy properties of our full sample and quantify the fraction of ultraviolet observable mass to support future observational programs, particularly those aimed at performing a metal budget around dwarf galaxies.

Debris disks are exoplanetary systems that contain planets, minor bodies (i.e., asteroids, Kuiper belt objects, comets, etc.), and micron-sized debris dust. Since water ice is the most common frozen volatile, it plays an essential role in the formation of planets and minor bodies. Although water ice has been commonly found in Kuiper belt objects and comets in the Solar System, no definitive evidence for water ice in debris disks has been obtained to date. Here, we report the discovery of water ice in the HD 181327 disk using the James Webb Space Telescope Near-Infrared Spectrograph. We detect the solid-state broad absorption feature of water ice at 3 $\mu$m and a distinct Fresnel peak feature at 3.1 $\mu$m, a characteristic of large water-ice particles. This implies the presence of a water-ice reservoir in the HD 181327 exoKuiper belt. Gradients of water-ice features at different stellocentric distances reveal a dynamic process of destroying and replenishing water ice in the disk, with estimated water-ice mass fractions ranging from 0.1% at ~85 au to 14% at ~113 au. It is highly plausible that the icy bodies that release water ice in HD 181327 could be the extra-solar counterparts of some of the Kuiper belt objects in our Solar System, supported by their spectral similarity.

Active Galactic Nuclei (AGN) are a key ingredient in galaxy evolution, possibly shaping galaxy growth through the generation of powerful outflows. Little is known regarding AGN-driven ionized outflows in moderate-luminosity AGN (logLbol[erg/s]<47) beyond cosmic noon (z>3). We present the first systematic analysis of the ionized outflow properties of a sample of X-ray-selected AGN (logLx[erg/s]>44) from the COSMOS-Legacy field at z~3.5 and with logLbol[erg/s]=45.2-46.7, by using JWST NIRSpec/IFU spectroscopic observations as part of the GA-NIFS program. We spectrally isolate and spatially resolve the ionized outflows, by performing a multi-component kinematic decomposition of the rest-frame optical emission lines. JWST/NIRSpec IFU data also revealed a wealth of close-by companions, of both non-AGN and AGN nature, and ionized gas streams likely tracing tidal structures and large-scale ionized gas nebulae, extending up to the circum-galactic medium. Ionized outflows are detected in all COS-AGN targets, which we compare with previous results from the literature up to z~3, opportunely (re-)computed for a coherent comparison. We normalize outflow energetics ($\dot{M}_{out}$, $\dot{E}_{out}$) to the outflow density to standardize the various assumptions that were made in the literature. Our choice is equal to assuming that each outflow has the same gas density. We find GA-NIFS AGN to show outflows consistent with literature results, within the large scatter shown by the collected measurements, suggesting no strong evolution with redshift in terms of total mass outflow rate, energy budget and outflow velocity for fixed bolometric luminosity. Moreover, we find no clear redshift evolution of the ratio of mass outflow rate and kinetic power over AGN bolometric luminosity beyond z>1. In general, our results indicate no significant evolution of the physics driving outflows beyond z~3. [abridged]

Early results from JWST suggest that epoch-of-reionization (EoR) galaxies produce copious ionizing photons, which, if they escape efficiently, could cause reionization to occur too early. We study this problem using JWST imaging and prism spectroscopy for 426 galaxies at 4.5 < z < 9.0. We fit these data simultaneously with stellar--population and nebular--emission models that include a parameter for the fraction of ionizing photons that escape the galaxy, f_esc. We find that the ionization production efficiency, $\xi_\mathrm{ion}$ = Q(H) / L(UV), increases with redshift and decreasing UV luminosity, but shows significant scatter, $\sigma( \log \xi_\mathrm{ion})$ ~= 0.3 dex. The inferred escape fractions averaged over the population are $\langle f_\mathrm{esc} \rangle$ = 3 +\pm 1% with no indication of evolution over 4.5 < z < 9.0. This implies that in our models nearly all of the ionizing photons need to be absorbed to account for the nebular emission. We compute the impact of our results on reionization, including the distributions for $\xi_\mathrm{ion}$ and f_esc, and the evolution and uncertainty of the UV luminosity function. Considering galaxies brighter than M(UV) < -16 mag, we would produce an IGM hydrogen-ionized fraction of x_e=0.5 at 5.5 < z < 7.1, possibly too late compared to constraints from from QSO sightlines. Including fainter galaxies, M(UV) < -14 mag, we obtain x_e = 0.5 at 6.7 < z < 9.6, fully consistent with QSO and CMB data. This implies that EoR galaxies produce plenty of ionizing photons, but these do not efficiently escape. This may be a result of high gas column densities combined with burstier star-formation histories, which limit the time massive stars are able to clear channels through the gas for ionizing photons to escape.

Bright-rimmed, cometary-shaped star-forming globules, associated with HII regions, are remnants of compressed molecular shells exposed to ultraviolet radiation from central OB-type stars. The interplay between dense molecular gas and ionizing radiation, analysed through gas kinematics, provides significant insights into the nature and dynamic evolution of these globules. This study presents the results of a kinematic analysis of the cometary globule, Lynds Bright Nebula (LBN) 437, focusing on the first rotational transition of 12CO and C18O molecular lines observed using the Taeduk Radio Astronomy Observatory (TRAO). The averaged 12CO spectrum shows a slightly skewed profile, suggesting the possibility of a contracting cloud. The molecular gas kinematics reveals signatures of infalling gas in the cometary head of LBN 437, indicating the initial stages of star formation. The mean infall velocity and mass infall rate towards the cometary head of LBN 437 are 0.25 km s^-1 and 5.08 * 10^-1 Msun yr^-1 respectively, which align well with the previous studies on intermediate or high-mass star formation.

We present JWST NIRCam observations of the emerging young star clusters (eYSCs) detected in the nearby spiral galaxy M83. The NIRcam mosaic encompasses the nuclear starburst, the bar, and the inner spiral arms. The eYSCs, detected in Pa$\alpha$ and Br$\alpha$ maps, have been largely missed in previous optical campaigns of young star clusters (YSCs). We distinguish between eYSCI, if they also have compact 3.3~$\mu$m PAH emission associated to them, and eYSCII, if they only appear as compact Pa$\alpha$ emitters. We find that the variations in the 3.3~$\mu$m PAH feature are consistent with an evolutionary sequence where eYSCI evolve into eYSCII and then optical YSCs. This sequence is clear in the F300M-F335M (tracing the excess in the \PAHlambda\ feature) and the F115W-F187N (tracing the excess in Pa$\alpha$) colors which become increasingly bluer as clusters emerge. The central starburst stands out as the region where the most massive eYSCs are currently forming in the galaxy. We estimate that only about 20~\% of the eYSCs will remain detectable as compact YSCs. Combining eYSCs and YSCs ($\leq$10 Myr) we recover an average clearing timescale of 6~Myr in which clusters transition from embedded to fully exposed. We see evidence of shorter emergence timescales ($\sim$5~Myr) for more massive ($>5\times10^3$ \msun) clusters, while star clusters of $\sim 10^3$ \msun\ about 7~Myr. We estimate that eYSCs remain associated to the \PAHlambda\ emission 3--4~Myr. Larger samples of eYSC and YSC populations will provide stronger statistics to further test environmental and cluster mass dependencies on the emergence timescale.

We show the variations of the CO J=2-1/1-0 line ratio (R21) across the barred spiral galaxy M83, using the 46 pc resolution data from ALMA. The R21 map clearly evidences the systematic large-scale variations as a function of galactic structures. Azimuthally, it starts from low R21<~0.7 in the interarm regions and becomes high ~>0.7 in the bar and spiral arms, suggesting that the density and/or kinetic temperature of molecular gas increase by about a factor of 2-3. This evolution is seen even in the parts of spiral arms without star formation, and R21 is often elevated even higher to ~0.8-1.0 when HII regions exist in the vicinity. Radially, R21 starts very high >~1.0 at the galactic center, remains low <~0.7 in the bar region, increases to >~0.7 around the bar end, and again decreases to <~0.7 in the rest of disk where the spiral arms dominate. The evolutionary sequence is synchronized with galactic rotation, and therefore, it is determined largely by the galactic structures and dynamics and is governed by the galactic rotation timescales. The R21 map also shows that the influence of stellar feedback is localized and limited. Massive, large, and non-star forming molecular structures have low R21, which also suggests that the bulk molecular gas in the disk is not regulated by stellar feedback, but more likely by galactic structures and dynamics. These results are consistent with suggestions by the earlier studies of the Milky Way and other barred spiral galaxies, and thus, are likely general among barred spiral galaxies in the local Universe.

Context. The crust composition of rocky exoplanets with a substantial atmosphere can not be observed directly. However, recent developments start to allow the observation and characterisation of their atmospheres. Aims. We aim to establish a link between the observable spectroscopic atmospheric features and the mineralogical crust composition of exoplanets. This allows to constrain the surface composition just by observing transit spectra. Methods. We use a diverse set of total element abundances inspired by various rock compositions, Earth, Venus, and CI chondrite as a basis for our bottom-to-top atmospheric model. We assume thermal and chemical equilibrium between the atmosphere and the planetary surface. Based on the atmospheric models in hydrostatic and chemical equilibrium with the inclusion of element depletion due to cloud formation theoretical transit spectra are calculated. Results. The atmospheric type classification allows constraints on the surface mineralogy especially with respect to sulphur compounds, iron oxides and iron hydroxides, feldspars, silicates and carbon species. Spectral features provide the possibility to differentiate the atmospheric types and thus allow some constraints on the surface composition.

Green Pea galaxies are a class of compact, low-mass, low-metallicity star-forming galaxies in the relatively local universe. They are believed to be analogues of high-redshift galaxies that re-ionised the universe and, indeed, the James Webb Space Telescope (JWST) is now uncovering such populations at record redshifts. Intriguingly, JWST finds evidence suggestive of active galactic nuclei (AGN) in many of these distant galaxies, including the elusive Little Red Dots, that broadly lack any detectable X-ray counterparts. Intuitively, one would expect to detect an AGN in their low-redshift analogues with X-rays, yet no study to date has conclusively identified an X-ray AGN within a Green Pea galaxy. Here we present the deepest X-ray campaign of a Green Pea galaxy performed to date, obtained with the goal of discerning the presence of a (potentially low-luminosity) AGN. The target $-$ SDSS J082247.66 +224144.0 (J0822+2241 hereafter) $-$ was previously found to display a comparable X-ray spectral shape to more local AGN ($\Gamma$ $\sim$ 2) and a high luminosity ($L_{2-10\,{\rm keV}}$ $\sim$ 10$^{42}$ erg s$^{-1}$). We show that over 6.2 years (rest-frame), the 2$-$10 keV luminosity of J0822+2241 is constant, whereas the soft 0.5$-$2 keV flux has decreased significantly by $\sim$60%. We discuss possible scenarios to explain the X-ray properties of J0822+2241, finding transient low-column density obscuration surrounding an AGN to be the only plausible scenario. J0822+2241 thus provides further evidence that low-luminosity AGN activity could have contributed to the epoch of reionisation, and that local analogues are useful to derive a complete multi-wavelength picture of black hole growth in high redshift low luminosity AGN.

Shock interaction in classical novae occurs when a fast outflow from the white dwarf > 1000 km s/s collides with a slower, cooler shell of gas released earlier in the outburst. The shocks radiate across the electromagnetic spectrum, from radio synchrotron to GeV gamma-rays. The hot shocked gas also emits >~ keV thermal X-rays, typically peaking weeks after the eruption, once the ejecta becomes transparent to photoelectric absorption. However, the observed hard X-ray luminosities are typically >4 orders of magnitude smaller than would be naively expected given the powerful shocks implied by the gamma-rays. We argue that a key missing piece to this puzzle is turbulence behind the shock, driven, e.g., by thin-shell and/or thermal instabilities. Turbulence efficiently mixes the hot X-ray emitting gas with cooler gas, sapping the hot gas of energy faster than it can directly radiate. Using analytic arguments motivated by numerical simulations, we show that energy losses due to turbulent mixing can easily balance shock heating, greatly reducing the volume of the hot gas and suppressing the X-ray luminosity. Equating the characteristic thickness of the X-ray emitting region to the minimum outer length scale of the turbulence capable of cooling the hot gas through mixing, we obtain X-ray luminosities consistent with nova observations if only ~1% of the shock's kinetic power goes into turbulent motions. A similar process may act to suppress thermal X-rays from other shock powered transients, such as interacting supernovae.

We present 12CO(J=1-0) mapping observations over ~1/2 of the optical disk of 12 nearby galaxies from the Fundamental CO 1-0 Transition Survey of nearby galaxies (FACTS), using the ALMA Total Power array. Variations in the 12CO(J=2-1)/12CO(J=1-0) line ratio r21 are investigated. The luminosity-weighted r21 of the 11 sample galaxies ranges from 0.52 to 0.69 with an average of 0.61. We use position-velocity diagrams along the major axis and tilted ring models to separate the normal rotating galactic disk from kinematic outliers that deviate from pure circular rotation. We find that r21 is systematically higher in outliers compared to the disk. We compare r21 between SA, SAB and SB galaxies, and find no significant difference in the average r21 depending on the presence of galactic bars. We find, however, that the radial gradient in r21 is bimodal, where a group containing all SA galaxies prefer constant or very shallow r21 gradients out 40% of the optical radius, while another group containing all SB galaxies have a steep r21 gradient, decreasing by ~20% before 40% of the optical radius, which also corresponds to the radius of the stellar bar. After this radius, these galaxies become consistent with a constant or shallow trend in r21. The large scale trend in r21 can have implications for how we interpret observations made solely in the 12CO(J=2-1) line.

The detection of extensive air showers (EAS) induced by cosmic rays via radio signals has undergone significant advancements in the last two decades. Numerous ultra-high energy cosmic ray experiments routinely capture radio pulses in the MHz to GHz frequency range emitted by EAS. The Monte Carlo simulation of these radio pulses is crucial to enable an accurate reconstruction of the primary cosmic ray energy and to infer the composition of the primary particles. In this work, a comprehensive comparison of the predicted electric field in EAS simulated with CoREAS and ZHAireS was conducted to estimate the systematic uncertainties arising from the use of different simulation packages in the determination of two key shower observables namely, the electromagnetic energy of the EAS and the depth of maximum development ($X_{\rm max}$). For this comparison, input parameters and settings as similar as possible were used in both simulations, along with the same realistic atmospheric refractive index depending on altitude, which is crucial for the prediction of radio emission properties of EAS. In addition, simulated EAS with very similar values of depth of maximum development were selected. Good agreement was found between CoREAS and ZHAireS, with discrepancies in the dominant electric field components generally remaining below $10\%$ across the frequency range of a few MHz to hundreds of MHz, relevant for most radio detection experiments, translating into uncertainties in the determination of energy below $5\%$ and $\simeq 10\,\mathrm{g/cm^2}$ in $X_{\rm max}$.

Light echoes of stellar flares provide an intriguing option for exploring protoplanetary disks in young stellar systems. Previous work on light echoes of circumstellar disks made use of delta-function flares for modeling. We present a new model that incorporates echoes produced by extended, time-resolved flares. We then test this model on known disk-bearing stars with Kepler K2 data by estimating disk parameters from possible echo signals. We focus on two stars; the first appears to be a good candidate for use of this echo model, which predicts disk parameters that are consistent with known values. The second star turns out to be more problematic as a result of high brightness variability in its post-peak lightcurve. These two cases show both the promise and limitations of light echoes as a tool for exploring protoplanetary disks in the time domain

High-resolution spectroscopy (HRS) has grown into one of the main techniques to characterise the atmospheres of extrasolar planets. High spectral resolving power allows for the efficient removal of telluric and host-star contamination. Combined with the large collecting area of ground-based telescopes it enables detailed studies of atmospheric species, temperature structure, atmospheric loss, and global winds and circulation patterns. In this review, the wide range of HRS observation and data-analysis techniques are described and literature results discussed. Key findings include: * The highest irradiated planets show a rich spectrum of atomic and ionic species, just like stars. * Retrieval analyses of Hot Jupiters and directly imaged Super- Jupiters point to Solar metallicities and chemistry, but observed samples are still heterogeneous and incomplete. * There appears to be a clear dichotomy between Hot Jupiters with and without atmospheric inversions, depending on their equilibrium temperature. * Some highly irradiated planets exhibit enormous leading and/or trailing tails of helium gas, providing unique insights into planet evolution and atmospheric escape processes. * Minor isotopes of carbon and oxygen are now being detected in gas giant planets and brown dwarfs with the interesting potential to shed light on formation pathways. A list of potential pitfalls is provided for those new to the field, and synergies with JWST are discussed. HRS has a great future ahead with the advent of the extremely large telescopes, promising to bring temperate rocky exoplanets into view with their increase in HRS detection speed of up to three orders of magnitude.

We summarize the X-ray emission from young SNe. Having accumulated data on most observed X-ray SNe, we display the X-ray lightcurves of young SNe. We also explore the X-ray spectra of various SN types. The X-ray emission from Type Ib/c SNe is non-thermal. It is also likely that the emission from Type IIP SNe with low mass-loss rates (around 10$^{-7} \, Msun \,$ yr$^{-1}$) is non-thermal. As the mass-loss rate increases, thermal emission begins to dominate. Type IIn SNe have the highest X-ray luminosities, and are clearly thermal. We do not find evidence of non-thermal emission from Type IIb SNe. The aggregated data are used to obtain approximate mass-loss rates of the progenitor stars of these SNe. Type IIP's have progenitors with mass-loss rates $< 10^{-5}\, Msun \,$ yr$^{-1}$, while Type IIn progenitors generally have mass-loss rates $> 10^{-3}\, Msun $ yr$^{-1}$. However, we emphasize that the density of the ambient medium is the important parameter, and if it is due to a non-steady outflow solution, it can not be translated into a mass-loss rate.

Known for their large amplitude radial pulsations, classical Cepheids are critical standard candles in astrophysics. However, they also exhibit various pulsational irregularities and additional signals that provide deeper insights into their structure and evolution. These signals appear in spectroscopic observations as shape deformations of the spectral lines. Using semi-partial distance correlation periodograms, we analyse high-precision spectroscopic data from the VELOCE project for four stars: $\delta$ Cep, BG Cru, X Sgr, and Polaris. For $\delta$ Cep, our control star, only the main radial mode is detected, confirming its stability and suitability as a benchmark for the method. In BG Cru, a strong additional signal at $\sim 3.01$ d is identified, likely linked to line splitting. X Sgr exhibits dominant additional signals, notably one at $\sim 12.31$ d, also associated with significant line splitting. Polaris reveals multiple low-frequency signals, with the most prominent candidate at $\sim 59.86$ d, which might be linked to the star's rotation period. We explore the semi-partial distance correlation periodograms by incorporating CCFs and their variants, such as the median-subtracted CCFs, which improves the sensitivity to variations in line shape. In particular, the latter enables the faithful detection of primary and additional signals present in the 1D spectra of fainter stars and low-amplitude pulsators. The semi-partial distance correlation periodograms demonstrate their utility for isolating signals associated with line shape variations; although, the analyses are complicated by the presence of artefact subharmonics and a visible low-frequency power increase for Polaris and BG Cru. This study underscores the method's potential for finding new and unexpected signals as well as detailed analyses of Cepheid pulsations and opens new pathways for asteroseismic investigations.

Multiplicity among massive stars represents a major uncertainty in stellar evolution theory, because the extra physical processes that it introduces significantly impacts each star's structure. While multiplicity of O-type stars is fairly well constrained, for B-type stars it is not. B-type stars are more common and have longer lifetimes, thus providing an opportunity to characterize multiplicity at different ages. Moreover, young open clusters are advantageous for studying coeval and chemically homogeneous environments. Using a multi-epoch spectroscopic campaign with the HERMES spectrograph we determine multiplicity properties and rotation rates of 74 B-type stars in four Galactic open clusters: $h$ and $\chi$ Persei, NGC 457, NGC 581, and NGC 1960. We measure radial velocities with a cross-correlation method and determine tentative orbital solutions for 26 of the 28 identified binaries. We detect several Be stars, five of them being members of binary systems. We correct the observed binary fractions for observational biases and determine an average intrinsic binary fraction of 79$_{-16}^{+19}\%$. The consistency in binary fraction among the four clusters, which are between 15 and 30 Myr, suggests a reasonably homogenous binary fraction across this age range. We used TLUSTY atmosphere models to determine the projected rotational velocities, with an average value of 240 km s$^{-1}$ for both single and binary systems. Whereas, the Be stars show higher velocities between 200 and 360 km s$^{-1}$.

It is widely known that the Saha equation is not suitable for describing plasmas out of thermodynamic equilibrium. The primordial hydrogen recombination plasma is an example of this. In this work, we propose a theoretical modification to the standard Saha curve by introducing a non-Gaussian generalization of the Saha equation, motivated by Tsallis statistics. In particular, we explore the possibility that a time-dependent $q$-parameter may serve as an effective proxy for the evolving thermodynamic conditions during recombination, especially considering that hydrogen recombination occurs from excited states. Within this framework, the $q$-parameter could be interpreted as encoding departures from equilibrium and could play the role of effective time-dependent temperature. This indicates that the time evolution of the $q$-parameter could provide a phenomenological mechanism for incorporating non-equilibrium effects into the recombination history. Our findings suggest that the Tsallis parameterization provides an alternative path to fit the distribution of free electrons by using an effective temperature. Our approach enables us to reproduce hydrogen recombination results predicted by the Peebles kinetic approach, HYREC-2 code, or any other recombination model.

The present paper is devoted to studying the variability of the $\gamma$ Cas type star SAO 49725. Both the optical and X-ray spectra of the star are analyzed. Variability of line profiles in the spectra of SAO 49725 was discovered on short (70-223 minutes) scales. Based on the TESS photometric light curves of SAO 49725 the regular variations of light curves are detected with a period of 1.1989 days identified as the rotation period of the star. The pattern of photometric variability in SAO 49725, as observed by TESS, varies significantly across different epochs. The TESS components of the SAO 49725 light curves with periods of $\sim3-21\,$days may be instrumental.

The question of the nature of precursors of solar flares, as well as their relationship with subsequent flares, still has no unambiguous answer. This is due, in particular, to the lack of systematic statistical work, the relative incompleteness (in isolation from the context of the development of the entire active region) of studies of individual events and the ambiguity of the term 'precursor' itself. In this paper, we consider the dynamics of the NOAA 12230 active region (AR), in which a series of homologous flares (C5-C9) occurred on December 9, 2014 within 12 hours, with an average frequency of about 2 hours. This AR was characterized by a rapid increase in flare activity followed by a rapid decay, which can be considered as a good example for studying potential precursors of flare series. We investigated the evolution of AR 12230 over a relatively long period (several days) and its transition from the 'no-flare' state to the flare-active regime. For this purpose, we study the magnetic field dynamics using SDO/HMI vector magnetograms, UV-EUV images based on SDO/AIA data, and X-ray observations from GOES/XRS and RHESSI data. We conclude that the significant increase in the chromospheric radiation variations against the background of a small flux of soft X-ray and UV emission from the corona observed on December 8, 2014, together with the magnetic flux emergence, can be considered as a precursor to a series of flares. The results show the importance and prospects of applying new methods of synoptic observations of the Sun in the context of collecting statistics ("history") of the AR energy release in different ranges of the electromagnetic spectrum.

We introduce the Mirya-m1 Cosmic Ray Detector, the largest and only cosmic ray detector in Turkiye designed for space weather research. Mirya-m1, modeled and built after the Muon Impact Tracer and Observer (MITO) (Ayuso et al. 2021), is located at the Eastern Anatolia Observatory (DAG) site of the Turkiye National Observatories in Erzurum, Turkiye, at an altitude of 3099 meters. This elevation positions Mirya-m1 among the highest-altitude cosmic ray detectors globally. The detector consists of two stacked scintillator counters, each measuring 1x1 meters, separated by a vertical distance of 1.36 meters. Each scintillator is monitored by four H1411 Hamamatsu photomultiplier tubes, enabling precise detection and measurement of light by incident cosmic rays. In this study, we present the data collected throughout 2024, which includes the detection of two Forbush decrease events in March and May 2024. These significant detections demonstrate the capability of Mirya-m1 to contribute valuable data for space weather research, establishing its potential as a critical instrument for cosmic ray studies in the region.

There is a growing number of Earth's co-orbital bodies being discovered. At least five of them are known to be temporarily in quasi-satellite orbits. One of those, 469219 Kamo'oalewa, was identified as possibly having the same composition as the Moon. We explore the conditions necessary for lunar ejecta to evolve into Earth's co-orbital bodies, with particular attention to quasi-satellite orbits. We investigate the parameter space of ejection velocity and geographic launch location across the lunar surface. The study employs numerical simulations of the four-body problem (Sun-Earth-Moon-particle) with automated classification for co-orbital states. Particles are ejected from randomly distributed points covering the lunar surface with velocities from 1.0 to 2.6 times the Moon's escape velocity. Trajectories co-orbital to Earth are found to be common, with approximately 6.68% of particles evolving into Earth co-orbital motion and 1.92% exhibiting quasi-satellite behavior. We identify an optimal ejection velocity (1.2v$_{esc}$) for quasi-satellite production, yielding over 6% conversion efficiency. Successful ejections show a strong preference for the equatorial regions of the trailing hemisphere. Collisions with Earth or Moon occur for only 4% of the sample. Extended integrations reveal long-lived configurations, including tadpole orbits persisting for 10,000 years and horseshoe co-orbitals maintaining stability for 5,000 years. Our results strengthen the plausibility of lunar origin for Earth's co-orbital bodies, including quasi-satellites like Kamo'oalewa and 2024PT5. We identify both "prompt" and "delayed" co-orbital formation mechanisms, with a steady-state production regime that could explain the presence of lunar-derived objects in Earth's co-orbital regions despite infrequent major lunar impacts.

When several solar flares with comparable classes occur successively at the same location and exhibit similar morphological features, they are called homologous flares. During 2012 May 8-10, five M-class homologous circular-ribbon flares associated with no coronal mass ejection occurred in active region (AR) 11476. The formation process of these homologous confined flares, particularly the homologous aspect, is unclear and inconclusive. This paper is dedicated to studying how the energy for this series of flares was accumulated and whether there existed null points responsible for the flare energy release. Before and during the five flares, the sunspots with opposite polarities sheared against each other and also rotated individually. Before each flare, the magnetic fields at the polarity inversion line were highly sheared and there existed a magnetic flux rope overlain by arch-shaped loops. For the first four flares, we find magnetic null points in the fan-spine topology, situated at about 3.8 Mm, 5.7 Mm, 3.4 Mm, and 2.6 Mm above the photosphere, respectively. For the fifth flare, no null point is detected. However, in the (extreme-)ultraviolet images, the evolution behaviors of all the flares were almost identical. Therefore, we speculate that a null point responsible for the occurrence of the fifth flare may have existed. These results reveal that, for these homologous flares in AR 11476, the sunspot rotation and shearing motion play important roles in energy accumulation, the null point of the fan-spine topology is crucial for energy release through magnetic reconnection therein, and large-scale magnetic loops prevent the erupting material from escaping the Sun, thus forming the observed homologous confined major circular-ribbon flares. This study provides clear evidence for the drivers of successive, homologous flares as well as the nature of confined events.

This paper presents the spherically-averaged 21 cm power spectrum derived from Epoch of Reionization (EoR) observations conducted with the Murchison Widefield Array (MWA). The analysis uses EoR0-field data, centered at (RA$=0h$, DEC$=-27^{\circ}$), collected between 2013 and 2023. Building on the improved methodology described in Trott et al. 2020, we incorporate additional data quality control techniques introduced in Nunhokee et al. 2024. We report the lowest power level limits on the EoR power spectrum at redshifts $z=6.5$, $z=6.8$, and $z=7.0$. These power levels, measured in the East-West polarization, are $(30.2)^2$ mK$^2$ at $k=0.18\, h$ Mpc$^{-1}$, $(31.2)^2$ mK$^2$ at $k=0.18\, h$ Mpc$^{-1}$ and $(39.1)^2$ mK$^2$ at $k=0.21\, h$ Mpc$^{-1}$ respectively. The total integration time amounts to 268 hours. These results represent the deepest upper limits achieved by the MWA to date and provide the first evidence of heated intergalactic medium (IGM) at redshifts $z=6.5$ to $7.0$.

Tianyu telescope, an one-meter robotic optical survey instrument to be constructed in Lenghu, Qinghai, China, is designed for detecting transiting exoplanets, variable stars and transients. It requires a highly automated, optimally distributed, easily extendable, and highly flexible software to enable the data processing for the raw data at rates exceeding 500MB/s. In this work, we introduce the architecture of the Tianyu pipeline and use relative photometry as a case to demonstrate its high scalability and efficiency. This pipeline is tested on the data collected from Muguang observatory and Xinglong observatory. The pipeline demonstrates high scalability, with most processing stages increasing in throughput as the number of consumers grows. Compared to a single consumer, the median throughput of image calibration, alignment, and flux extraction increases by 41%, 257%, and 107% respectively when using 5 consumers, while image stacking exhibits limited scalability due to I/O constraints. In our tests, the pipeline was able to detect two transiting sources. Besides, the pipeline captures variability in the light curves of nine known and two previously unknown variable sources in the testing data. Meanwhile, the differential photometric precision of the light curves is near the theoretical limitation. These results indicate that this pipeline is suitable for detecting transiting exoplanets and variable stars. This work builds the fundation for further development of Tianyu software. Code of this work is available at this https URL.

Magnetic massive stars are stars of spectral types O, B and A that harbour $\sim$ kG strength (mostly dipolar) surface magnetic fields. Their non-thermal radio emission has been demonstrated to be an important magnetospheric probe, provided the emission is fully characterised. A necessary step for that is to build a statistically significant sample of radio-bright magnetic massive stars. In this paper, we present the `VAST project to study Magnetic Massive Stars' or VAST-MeMeS that aims to achieve that by taking advantage of survey data acquired with the Australian SKA Pathfinder telescope. VAST-MeMeS is defined under the `VAriable and Slow Transient' (VAST) survey, although it also uses data from other ASKAP surveys. We found radio detections from 48 magnetic massive stars, out of which, 14 do not have any prior radio detections. We also identified 9 `Main-sequence Radio Pulse Emitter' candidates based on variability and circular polarisation of flux densities. The expanded sample suggests a lower efficiency in the radio production than that reported in earlier work. In addition to significantly expanding the sample of radio-bright magnetic massive stars, the addition of flux density measurements at $\lesssim 1$ GHz revealed that the spectra of incoherent radio emission can extend to much lower frequencies than that assumed in the past. In the future, radio observations spanning wide frequency and rotational phase ranges should be conducted so as to reduce the uncertainties in the incoherent radio luminosities. The results from these campaigns, supplemented with precise estimations of stellar parameters, will allow us to fully understand particle acceleration and non-thermal radio production in large-scale stellar magnetospheres.

Evidence indicates that supermassive black holes exist at the centers of most galaxies. Their mass correlates with the galactic bulge mass, suggesting a co-evolution with their host galaxies, most likely through powerful winds. X-ray observations have detected highly ionized winds outflowing at sub-relativistic speeds from the accretion disks around supermassive black holes. However, the limited spectral resolution of current X-ray instruments has left the physical structure and location of the winds poorly understood, hindering accurate estimates of their kinetic power. Here, the first XRISM observation of the luminous quasar, PDS 456, is reported. The high-resolution spectrometer Resolve onboard XRISM enabled the discovery of five discrete velocity components outflowing at 20-30% of the speed of light. This demonstrates that the wind structure is highly inhomogeneous, which likely consists of up to a million clumps. The mass outflow rate is estimated to be 60-300 solar masses per year, with the wind kinetic power exceeding the Eddington luminosity limit. Compared to the galaxy-scale outflows, the kinetic power is more than 3 orders of magnitude larger, while the momentum flux is 10 times larger. These estimates disfavor both energy- and momentum-driven outflow models. It suggests that such wind activity occurs in less than 10% of the quasar phase and/or that its energy/momentum is not efficiently transferred to the galaxy-scale outflows due to the clumpiness of the wind and the interstellar medium.

The $S_8$ value inferred from the Subaru Hyper Suprime-Cam (HSC) Year 3 cosmic shear data, under the assumption of the flat $\Lambda$CDM model, is 2-3$\sigma$ lower than that inferred from observations of the early-time universe, such as cosmic microwave background (CMB) anisotropy data. Resolving the $S_8$ tension requires a scenario in which structure formation on small scales is suppressed in the late universe. As potential solutions, we consider extended models both within and beyond the $\Lambda$CDM model -- models that incorporate parameterized baryonic feedback effects, the effect of varying neutrino mass, and modified structure growth, each of which can lead to a suppression of structure growth at lower redshifts, with its own distinct scale- and redshift-dependencies. In particular, we consider phenomenological modified gravity models in which the suppression of structure growth is triggered at lower redshifts, as dark energy ($\Lambda$) begins to dominate the background expansion. We show that the modified growth factor models -- especially those featuring more rapid growth suppression at lower redshifts -- provide an improved fit to the combined datasets of the HSC-Y3 cosmic shear correlation functions, the Planck CMB, and the ACT DR6 CMB lensing, compared to the fiducial $\Lambda$CDM model and the models including the baryonic effects or the massive neutrino effect within the the $\Lambda$CDM framework.

High energy solar protons were observed by particle detectors aboard spacecraft in near-Earth orbit on May 11, 2024 and produced the 74th ground level enhancement (GLE74) event registered by ground-based neutron monitors. This study involves a detailed reconstruction of the neutron monitor response, along with the identification of the solar eruption responsible for the emission of the primary particles, utilizing both in situ and remote-sensing. Observations spanning proton energies from a few MeV to around 1.64 GeV, collected from the Solar and Heliospheric Observatory (SOHO), the Geostationary Operational Environmental Satellite (GOES), the Solar Terrestrial Relations Observatory (STEREO-A), and neutron monitors, were integrated with records of the associated solar soft X-ray flare, coronal mass ejection, and radio bursts, to identify the solar origin of the GLE74. Additionally, a time-shift analysis was conducted to link the detected particles to their solar sources. Finally, a comparison of GLE74 to previous ones is carried out. GLE74 reached a maximum particle rigidity of at least 2.4 GV and was associated with a series of type III, type II, and type IV radio bursts. The release time of the primary solar energetic particles (SEPs) with an energy of 500 MeV was estimated to be around 01:21 UT. A significant SEP flux was observed from the anti-Sun direction with a relatively broad angular distribution, rather than a narrow, beam-like pattern, particularly during the main phase at the particle peak flux. Comparisons with previous GLEs suggest that GLE74 was a typical event in terms of solar eruption dynamics.

Giant radio galaxies (GRG) are radio galaxies with physical sizes of their radio emission larger than 0.7 Mpc. Recently, the sample of GRGs has become large enough to study the extreme end of the GRG size distribution. We examine the properties of GRGs with largest linear sizes larger than 3 Mpc in order to shed light on the nature and origin of GRGs. We select, corroborate, and revise if necessary, the largest GRGs from literature. We add to these the GRGs we identified in our own search, combined with optical surveys and catalogues of spectroscopic and photometric redshifts in order to find their projected linear radio size. We study their radio power--size relation, the asymmetry in the lobes, their association with clusters of galaxies, as well as their bending angles. We present an unprecedented sample of 143 GRGs larger than 3 Mpc, of which 69 were newly found by us. The sample includes six GRGs with projected linear sizes clearly exceeding 5 Mpc and reaching up to 6.6 Mpc. We find that GRGs larger than 3 Mpc are distributed in redshift and radio luminosity indistinguishable from those of smaller GRGs. At most a single one of the GRGs larger than 3 Mpc can be classified as a clear Fanaroff-Riley (FR) type I source, and only 6 per cent deviate from a clear FR II radio morphology. One quarter of our GRGs show very diffuse lobes typical for remnant radio galaxies, and only 59 per cent show indications of hotspots in at least one lobe, with 38 per cent featuring a hotspot in both lobes. As in the case of smaller radio galaxies, the shorter lobe is most often also the brighter one. We find tentative evidence that the bending angle decreases with size of the GRG, but no trend with redshift is detected. The bending angles of GRGs > 3 Mpc in known clusters are larger than for those GRGs not associated with clusters.

We present the first results of a pilot 'TASmanian Search for Inclined Exoplanets' (TASSIE) program. This includes observations and analysis of five short-period exoplanet candidates using data from TESS and the Harlingten 50 cm telescope at the Greenhill Observatory. We describe the instrumentation, data reduction process and target selection strategy for the program. We utilise archival multi-band photometry and new mid-resolution spectra to determine stellar parameters for five TESS Objects of Interest (TOIs). We then perform a statistical validation to rule out false positives, before moving on to a joint transit analysis of the remaining systems. We find that TOI3070, TOI3124 and TOI4266 are likely non-planetary signals, which we attribute to either short-period binary stars on grazing orbits or stellar spots. For TOI3097, we find a hot sub-Jovian to Jovian size planet ($R_{3097Ab}$ = 0.89 $\pm$ 0.04 $R_{J}$, $P_{3097Ab}$ = 1.368386 $\pm$ 0.000006 days) orbiting the primary K dwarf star in a wide binary system. This system shows indications of low metallicity ([Fe/H] $\approx$ -1), making it an unlikely host for a giant planet. For TOI3163, we find a Jovian-size companion on a circular orbit around a late F dwarf star, with $R_{3163b}$ = 1.42 $\pm$ 0.05 $R_{J}$ and $P_{3163b}$ = 3.074966 $\pm$ 0.000022 days. In future, we aim to validate further southern giant planet candidates with a particular focus on those residing in the sub-Jovian desert/savanna.

We present a JWST and ALMA detailed study of the ISM properties of high-redshift galaxies. Our JWST/NIRSpec IFU spectroscopy targeting three galaxies at $z=6-7$ detects key rest-frame optical emission lines, allowing us to derive [OII]$\lambda\lambda$3726,3729-based electron densities of $n_\mathrm{e,optical}\sim1000$ cm$^{-3}$ on average and [OIII]$\lambda$4363-based metallicities of $\mathrm{12+log(O/H)}=8.0-8.2$ in two galaxies. New ALMA Band 9/10 observations detect the [OIII]52$\mu$m line in one galaxy but do not in the others, resulting in FIR-based densities of $n_\mathrm{e,FIR}\lesssim500$ cm$^{-3}$ from the [OIII]52$\mu$m/[OIII]88$\mu$m ratios, systematically lower than the optical [OII]-based measurements. These low FIR-based densities are comparable to those at both $z\sim0$ and $z>6$ in the literature, including JADES-GS-z14-0 at $z=14.18$, suggesting little evolution up to $z\sim14$, in contrast to the increasing trend of optical-based densities with redshift. By conducting a JWST and ALMA joint analysis using emission lines detected with both telescopes, we find that the observed FIR [OIII]52,88$\mu$m luminosities are too high to be explained by the optical-based densities at which they would be significantly collisionally de-excited. Instead, a 2-zone model with distinct high- and low-density regions is required to reproduce all observed lines, indicating that FIR [OIII] emission arises predominantly from low-density gas, while optical [OIII] and [OII] lines trace both regions. We further demonstrate that the direct-$T_\mathrm{e}$ method can sometimes significantly underestimate metallicities up to 0.8 dex due to the presence of the low-density gas not fully traced by optical lines alone, highlighting the importance of combining optical and FIR lines to accurately determine gas-phase metallicities in the early universe.

We investigate the effects of dark matter (DM) on neutron star (NS) properties using the neutron decay anomaly model within the relativistic mean-field (RMF) framework. Three nucleonic models (HCD0-HCD2) are developed, satisfying astrophysical constraints such as the maximum NS mass ($\geq 2 M_\odot$), the NICER mass-radius limits, and the tidal deformability constraint from the GW170817 event. The equation of states of the NS admixed with DM (DMANS) are calculated by incorporating the self-interactions between them. The macroscopic properties, such as mass, radius, and tidal deformability of the NSs, are obtained for HCD models along with five others by varying self-interaction strength. By combining NS observations with scattering cross-section constraints from galaxy clusters, we explore model-dependent trends in the DM self-interaction parameter space. While the quantitative bounds may vary with hadronic model choice, our analysis offers insights into the interplay between DM interactions and NS observables within the RMF framework.

The nature of dark matter is still mysterious despite various astronomical evidence. As a possible candidate, self-interacting dark matter (SIDM) can potentially resolve some issues appearing in cold dark matter paradigm. Here we investigate how SIDM around supermassive black holes (SMBH) in galaxy centers may form a density spike and imprint in the spectrum shape of stochastic gravitational-wave background from SMBH binaries (SMBHBs). Employing a refined dynamical friction formula and consistently evolving the orbital dynamics, we demonstrate that current pulsar timing arrays (PTAs) data is sensitive to the cross section of SIDM with $\sigma(v)/m_\chi\lesssim0.66\,\mathrm{cm}^2/\mathrm{g}$, comparable to other astrophysical probes. We also highlight the importance of including the eccentricity of SMBHBs in the parameter inference, which would affect the results significantly. Our findings reveal the promising potential of PTAs observations in probing the nature of dark matter.

Angular momentum transport by magnetic fields is important for formation and evolution of protoplanetary disks. The effects of magnetic fields are suppressed due to non-ideal magnetohydrodynamic (MHD) effects such as ambipolar diffusion and Ohmic dissipation, which depend on the degree of ionization. Cosmic rays (CRs) are the primary source of ionization in star-forming clouds, and their distribution is nonuniform as it is affected by gas density and magnetic fields. Therefore, CRs, magnetic fields, and gas interact with each other. In this work, we develop a new fully implicit cosmic ray transport module in Athena++ and perform three-dimensional simulations of disk formation from collapse of molecular cloud cores. Since CRs are strongly attenuated in the dense gas at the disk scale, distribution of magnetic fields is considerably altered compared to conventional models assuming a uniform ionization rate. While the total magnetic fluxes accreted onto the disks remain similar as the gas outside the disks remain sufficiently ionized and well coupled, the magnetic fields in the disks are less twisted due to the stronger non-ideal MHD effects. As a consequence, magnetic angular momentum transport is strongly suppressed at the disk scale, resulting in more gravitationally unstable disks with more prominent spiral arms. Our simulations demonstrate influence of non-uniform ionization resulting from CR transport and attenuation on the disk formation and evolution.

An automatic procedure to perform sub-clustering on large samples is presented. At each iteration, the most diverse cluster is sub-clustered, and the global diversity of the new classification is compared to the previous one. The process stops if no improvement is found. The key to our procedure is the use of a quantitative measure of diversity, called the Leinster-Cobbold index, that takes into account the similarity between clusters. While this procedure has been successfully applied on a large sample of spectra of galaxies, we illustrate its efficiency with two examples in this paper.

When a star is torn apart by the tidal forces of a supermassive black hole (a so-called TDE) a transient accretion episode is initiated and a hot, often X-ray bright, accretion disk is formed. Like any accretion flow this disk is turbulent, and therefore the emission from its surface will vary stochastically. As the disk has a finite mass supply (i.e., at most the initial mass of the disrupted star) the disk will also undergo long-timescale evolution, as this material is lost into the black hole. In this paper we combine theoretical models for this long time evolution of the disk with models for the stochastic variability of turbulent accretion flows which are correlated on short (orbital) timescales. This new framework allows us to demonstrate that (i) dimming events should be more prevalent than brightening events in long term TDE X-ray light curves (i.e., their log-luminosity distribution should be asymmetric), (ii) TDE X-ray light curves should follow a near- (but formally sub-)linear correlation between their root mean square variability and the mean flux, (iii) the fractional variability observed on short timescales across an X-ray observing band should increase with observing energy, and (iv) TDEs offer a unique probe of the physics of disk turbulence, owing to their clean spectra and natural evolutionary timescales. We confirm predictions (i) and (ii) with an analysis of the long timescale variability of two observed TDEs, and show strong support for prediction (iii) using the intra-observation variability of the same two sources.

We present a fully analytical propagator for the orbits of lunar artificial satellites in a lunar gravity and third-body model sufficiently precise for a wide range of practical applications. The gravity model includes the twelve most important lunar gravity harmonics as well as the Earth's quadrupole tidal terms with a precise representation of the Earth's lunicentric ephemeris, and it gives an accuracy comparable to the way more extended semi-analytical propagator SELENA [6] for satellite orbits at altitudes from 300 to 3000 km. Extra terms of a more complete gravity model are straightforward to include using the formulas of the presently discussed analytical theory. The theory is based on deriving an approximate analytical solution of the secular part of the equations of motion using a Hamiltonian normal form in closed form. In total, we have two types of element transformations: from osculating to mean elements (as in [6]), and from mean to proper elements. The solution of the problem in proper elements is trivial, and, through the inverses of the above transformations, it allows to recover the position and velocity of a satellite analytically at any time t given initial conditions of the osculating elements at time $t_0$ without any intermediate numerical propagation. The propagator model is valid in time spans of several decades, and for every initial condition leading to no-fall on the Moon's surface, except for identified thin zones around a set of secular resonances corresponding to commensurabilities between the satellite's secular frequencies and the secular frequencies of the lunicentric Earth's orbit. Open software python and symbolic routines implementing our propagator are provided in the repository [14]. Precision tests with respect to fully numerical orbital propagation in Cartesian coordinates are reported.

We explore the density profile, shape, and virial mass of the Milky Way's dark matter halo using K giants (KG) from LAMOST and SDSS/SEGUE, as well as blue horizontal branch (BHB) stars from SDSS. Incorporating Gaia DR3 proper motions, we first investigate the velocity ellipsoid distribution within the $(R, |z|)$ space. The ellipsoids projected onto the $(v_R, v_z)$ plane exhibit near-spherical alignment. We then probe the underlying dark matter distribution using the axisymmetric Jeans equations with multi-Gaussian expansion (MGE) and the spherically aligned Jeans anisotropic modelling (JAM${\rm sph}$), allowing for different flattened dark matter density models. For each model, we apply two fitting approaches: fitting the KGs and BHBs separately or fit them simultaneously as two dynamical tracers in one gravitational potential. We find consistent results on the dark matter density profiles, $r_{200}$, and $M_{200}$ within a 1-$\sigma$ confidence region for models constrained by KGs, BHBs, and both. We find the strongest consistency between KGs and BHBs in constraining dark matter profiles for models incorporating radially varying halo flattening ($q(r_{\rm gc})$), which suggests the Milky Way's dark matter halo shape evolves with Galactocentric distance ($r_{\rm gc}$). Specifically, the halo flattening parameter $q_h$ decreases within $r_{\rm gc} < 20$ kpc and increases for $r_{\rm gc} > 20$ kpc. In this model, $M_{\rm tot} (< 60~{\rm kpc}) = 0.533^{+0.061}_{-0.054} \times 10^{12}$ $M_{\odot}$, $r_{200}$ is $188\pm15$ kpc, with $M_{200}$ estimated at $0.820^{+0.210}_{-0.186} \times 10^{12} M_{\odot}$.

The well-known discrepancy in galactic rotation curves refers to the mismatch between observed rotational velocities and the velocities predicted by baryonic matter. In this study we investigate a potential pattern in the discrepancy, which may point to an underlying pattern in dark-halo distributions. By looking at rotational-velocity curves from an alternate perspective, the angular-velocity curves, it appears that the observed angular velocities and their corresponding baryonic predictions differ by a constant shift. That is, the discrepancy may be reduced to a constant angular-velocity term, independent of the radius. We test the generality of the suggested property by analysing 143 high-quality rotation curves. The property appears significant as it performs equally well (or better) than well-established dark-halo models. Compared to a Burkert profile, it is preferred in 60% of the cases, while relative to a Navarro Frenk White profile (NFW), it is superior in 73% of the cases. Next, by including the new phenomenological property within the dynamical equations, we find an explicit expression for the dark-halo profile. The new single-parameter profile is characterised by a remarkable property: it is intrinsically related to the baryonic distribution. Thus, information regarding the cuspy or cored nature of a particular dark halo, according to this profile, is encoded (and explicitly determined) by the respective baryonic behaviour.

Theoretical frameworks for reflection and emission spectroscopy of exoplanet surfaces are becoming increasingly important for the characterization of rocky exoplanets, especially with the rapid growth of the detected exoplanet population and observational capabilities. The Hapke theory of reflectance and emittance spectroscopy has been widely adopted in the exoplanet community, yet a key physical effect - the opposition surge enhancement at small phase angles - remains largely neglected. This phenomenon, driven by shadow hiding and coherent backscattering, introduces a significant brightening that depends on wavelength, particle size, and surface morphology. In this paper, I propose an alternative formulation for opposition surge modeling, ensuring a smooth-to-sharp transition at small phase angles, dictated by wavelength-dependent scattering properties. I evaluate the impact of opposition surge on phase curves and surface spectra, comparing a family of models with increasing simplifications, ranging from a full wavelength-dependent opposition effect to its complete omission. My results indicate that neglecting opposition effects can introduce systematic deviations in retrieved albedos, spectral features, and phase curves, with errors reaching up to 20%-30% in certain spectral bands. Upcoming JWST observations will probe phase angles below ~10° for rocky exoplanets around M dwarfs; thus, accounting for opposition effects is crucial for accurate surface characterization. Proper treatment of this effect will lead to improved retrievals of surface albedo, mineralogical composition, and roughness properties. This study establishes a physically consistent framework for exoplanet phase-curve modeling and provides a foundation for future retrieval algorithms aimed at interpreting exoplanet surfaces.

The frequencies of gravity mode oscillations are determined by the chemical, thermal, and structural properties of stellar interiors, facilitating the study of internal mixing mechanisms in stars. We investigate the impact of discontinuities in the chemical composition induced by the formation of an adiabatic semiconvection region during the core helium (He)-burning phase of evolution of hot subdwarf B-type (sdB) stars. We scrutinize the asteroseismic attributes of the evolutionary stages and assess the core He-burning phase by evaluating the parameter linked to the average interval between the deep trapped modes in both sdB evolutionary models and the observations of KIC 10001893. We perform evolutionary and asteroseismic analyses of sdB stars using MESA and GYRE to examine the properties of the semiconvection region. We address the challenges of relying solely on average interval between oscillation mode periods with consecutive radial orders to identify the core He-burning stage. To enhance identification, we propose a new parameter representing the average interval between deep trapped modes during some of the stages of sdB evolutionary models. Additionally, our results show that integrating convective penetration with convective premixing improves our models and yields comparable outcomes without the need for additional model parameters. Our results can advance the development of detailed evolutionary models for sdB stars by refining internal mixing schemes, enhancing the accuracy of pulsation predictions, and improving alignment with observational data.

We report a novel mechanism where two families of primordial black holes (PBHs) may form at nearly the same comoving scales but at two different epochs. It is realized in two-stage inflation where a non-inflationary stage is sandwiched by the two inflationary stages. In this case, smaller PBHs form when the comoving scale of interest re-enters the horizon during the break period, and larger PBHs form when the scale re-enters the horizon after inflation. This mechanism may realize both reheating of the universe through the evaporation of ultralight PBHs formed during the break stage and the dark matter by those formed after inflation. We show that this scenario may give rise to a distinctive signature in the stochastic gravitational wave background that can be tested by the near-future gravitational wave observatories such as LISA and DECIGO. Our work thus provides a unified observational window into the physics of inflation, reheating, and dark matter.

In the study of stellar flares, traditional method of calculating total energy emitted in the continuum assumes the emission originating from a narrow chromospheric condensation region with a constant temperature of 10000 K and variable flare area. However, based on multicolor data from 7 new flares observed in Białków and Shumen observatory and 8 flares from (Howard et al., 2020) observed on 10 main-sequence stars (spectral types M5.5V to K5V) we show that flare areas had a relative change in the range of 10% - 61% (for more than half of the flares this value did not exceed 30%) throughout the events except for the impulsive phase, and had values starting from 50 +/- 30 ppm to 300 +/- 150 ppm for our new flares and from 380 +/- 200 ppm to 7600 +/- 3000 ppm from (Howard et al., 2020) data, while their temperature increased on average by the factor 2.5. The peak flare temperatures for our seven observed flares ranged from 5700 +/- 450 K to 17500 +/- 10050 K. Five of these flares had their temperatures estimated using the Johnson-Kron-Cousins B filter alongside TESS data, one flare was analyzed using the SLOAN g' and r' bandpasses, and another was evaluated using both the SLOAN g' and r' bandpasses and TESS data. Using flare temperature and area from our data and data from Howard et al., 2020, along with the physical parameters of stars where the flares occurred, we developed a semi-empirical grid that correlates a star's effective temperature and flare amplitude in TESS data with the flare's peak temperature. This allows interpolation of a flare's peak temperature based on the star's effective temperature (ranging from 2700 K to 4600 K) and flare amplitude from TESS observations. Applying this grid to 42257 flares from TESS survey, we estimated peak flare temperatures between 5700 K and 38300 K, with most flares showing peak black-body temperatures around 11100 +/- 2400 K.

Accretion disks around black holes host extreme conditions where general relativity and magnetohydrodynamics dominate. These disks exhibit two distinct dynamical regimes -- Standard and Normal Evolution (SANE) and Magnetically Arrested Disk (MAD). In the MAD regime, these systems exhibit magnetic fields up to $10^8$ G and variability on gravitational timescales $t_g \sim 10^{-4}$ s for stellar-mass black holes. While classical magnetohydrodynamics has been extensively applied, quantum effects in these high-energy environments remain unexplored. Here, we employ quantum field theory in background gauge fields (QFTBGF) to demonstrate that the dynamic magnetic fields of MADs drive significant pair production via the Schwinger mechanism. The resulting pairs emit non-thermal (synchrotron) radiation with a peak frequency tunable across $ \sim 1 - 3000$ MHz, depending on the magnetic field strength (peaking at higher frequencies for stronger fields). For $ B \sim 10^8 $ G, our model predicts a peak spectral flux density of $ \sim 1 - 100$ mJy, detectable with next-generation radio telescopes (e.g., SKA, ngVLA). This work provides a direct and observable signatures of quantum effects in black hole accretion disks.

The Lyman-$\alpha$ forest opacity fluctuations observed from high redshift quasar spectra have been proven to be extremely successful in order to probe the late phase of the reionization epoch. For ideal modeling of these opacity fluctuations, one of the main challenges is to satisfy the extremely high dynamic range requirements of the simulation box, resolving the Lyman alpha forest while probing the large cosmological scales. In this study, we adopt an efficient approach to model Lyman-$\alpha$ opacity fluctuations in coarse simulation volume, utilizing the semi-numerical reionization model SCRIPT (including inhomogeneous recombination and radiative feedback) integrated with a realistic photoionization background fluctuation generating model. Our model crucially incorporates ionization and temperature fluctuations, which are consistent with the reionization model. After calibrating our method with respect to high-resolution full hydrodynamic simulation, Nyx, we compare the models with available observational data at the redshift range, $z=5.0-6.1$. With a fiducial reionization model (reionization end at $z=5.8$), we demonstrate that the observed scatter in the effective optical depth can be matched reasonably well by tuning the free parameters of our model i.e. the effective ionizing photon mean free path and mean photoionization rate. We further pursue an MCMC-based parameter space exploration utilizing the available data to put constraints on the above free parameters. Our estimation prefers a slightly higher photoionization rate and slightly lower mean free path than the previous studies which is also a consequence of temperature fluctuations. This study holds significant promise for efficiently extracting important physical information about the Epoch of Reionization, utilizing the wealth of available and upcoming observational data.

We present general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers, exploring the conditions to launch relativistic outflows compatible with gamma-ray burst (GRB) jets. Employing an unprecedentedly low numerical density floor decreasing as the sixth power of radial distance, we obtain for the first time an evolution where such a floor has no effect on any outflow produced during and after merger. In models forming an accreting black hole (BH) 25 ms after merger, incipient jets are launched with terminal Lorentz factors and Poynting-flux luminosities compatible with GRBs. The rather high jet velocities reached in a few tens of ms, when compared to previous simulations adopting a much higher and uniform density floor, show the importance of removing any influence of the latter as the jet evolves toward increasingly large time and spatial scales. In absence of collapse, we find the emergence of much heavier and slower polar outflows, incapable of producing a GRB. We thus favor a BH origin for GRB jets.

Context. The sizes of many asteroids, especially slowly rotating, low-amplitude targets, remain poorly constrained due to selection effects. These biases limit the availability of high-quality data, leaving size estimates reliant on spherical shape assumptions. Such approximations introduce significant uncertainties propagating, e.g. into density determinations or thermophysical and compositional studies, affecting our understanding of asteroid properties. Aims. This work targets poorly studied main-belt asteroids, most of which lacked shape models. Using only high-quality dense light curves, thermal IR observations (incl. WISE), and stellar occultations, we aimed to produce reliable shape models and scale them via two independent techniques, allowing size comparison. We conducted two campaigns to obtain dense photometric light curves and to acquire multi-chord stellar occultations. Methods. Shape and spin models were reconstructed using lightcurve inversion. Sizes were determined by (1) thermophysical modeling with the Convex Inversion Thermophysical Model (CITPM), optimizing spin and shape models to visible lightcurve and IR data, and (2) scaling shape models with stellar occultations. Results. We obtained precise sizes and shape models for 15 asteroids. CITPM- and occultation-derived sizes agree within 5% in most cases, demonstrating the modeling's reliability. Larger discrepancies usually stem from incomplete occultation chord coverage. The study also gives insights into surface properties incl. albedo, roughness and thermal inertia. Conclusions. Using high-quality data and an advanced TPM integrating thermal and visible data with shape adjustment enabled precise size estimates comparable to those from multi-chord stellar occultations. We resolved major inconsistencies in previous size estimates, providing solid input for future studies on asteroid densities and surfaces.

Recent Type Ia supernova (SN Ia) simulations featuring a double detonation scenario have managed to reproduce the overall trend of the Phillips relation reasonably well. However, most, if not all, multidimensional simulations struggle to reproduce the scatter of observed SNe around this relation, exceeding it substantially. In this study, we investigate whether the excessive scatter around the Phillips relation can be caused by an off-center ignition of the carbon-oxygen (CO) core in the double detonation scenario and if this can help constrain possible SN Ia explosion channels. We simulated the detonation of three different initial CO white dwarfs of $0.9$, $1.0$, and $1.1\,M_\odot$, artificially ignited at systematically offset locations using the Arepo code. After nucleosynthetic postprocessing, we generated synthetic observables using the Artis code and compared these results against observational data and models of other works. We find that our simulations produce synthetic observables well within the range of the observed data in terms of viewing angle scatter. The majority of the viewing angle variability seems to be caused by line blanketing in the blue wavelengths of intermediate-mass elements and lighter iron-group elements, which are asymmetrically distributed in the outer layers of the ashes. Our results suggest that although the off-center ignition of the CO introduces substantial line of sight effects, it is not responsible for the excessive viewing angle scatter observed in other models. Instead, this effect seems to be caused by the detonation ashes from the rather massive helium (He) shells in current state-of-the-art models. Further reducing the He-shell masses of double detonation progenitors may be able to alleviate this issue and yield observables that reproduce the Phillips relation.

This article questions the common assumption of cold dark matter (DM) by exploring the possibility of a non-zero equation of state (EoS) without relying on any parametric approach. In standard cosmological analyses, DM is typically modeled as pressureless dust with $w_{\rm DM} = 0$, an assumption that aligns with large-scale structure formation, supports the empirical success of the $\Lambda$CDM model, and simplifies cosmological modeling. However, there is no fundamental reason to exclude a non-zero $w_{\rm DM}$ from the cosmological framework. In this work, we explore this possibility through non-parametric and model-independent reconstructions based on Gaussian Process Regression. The reconstructions use Hubble parameter measurements from Cosmic Chronometers (CC), the Pantheon+ sample of Type Ia supernovae, and Baryon Acoustic Oscillation (BAO) data from DESI DR1 and DR2. Our findings suggest that a dynamical EoS for DM, although statistically mild, cannot be conclusively ruled out. Notably, we observe a tendency toward a negative $w_{\rm DM}$ at the present epoch, a result that is particularly intriguing and may have meaningful implications for modern cosmology.

A mechanism is proposed for synchronizing core-collapse supernova with a recent loss of a red supergiant (RSG) envelope in the common envelope regime. A perequisite for the synchronization is a moderate RSG expansion during final decade. This scenario is based on the phenomenon of preSN~II dense shell formed at the final stage of 10-20 yr as a result of powerfull mass loss. The energy deposition into the RSG envelope that powers the enormous mass loss rate is able to expand the RSG. The moderate expansion is sufficient for the close secondary component to plunge into the common envelope with a subsequent explosion of stripped helium core. Superluminous SN~2006gy is suggested to be the outcome of this scenario.

The Yarkovsky effect on real asteroids is complicated to calculate either by analytical or numerical methods, since they are generally irregular in shape. We propose an index to properly characterise the shape of any asteroid, through which the Yarkovsky effect can be easily calculated without the heavy computations of surface temperatures. By analysing the energy absorbed and then emitted by a surface element, we find that the effective working power produced by the radiation recoil force on this surface element and its contribution to the Yarkovsky effect are both proportional to the double projected area of the surface element. The normalized total projected area over the asteroid's surface is defined as the shape index ($S_1$). We model the Yarkovsky effects of different asteroids using multiphysics software COMSOL, and take the rate of semi-major axis drift ($d a/d t$) obtained in these numerical simulations as the measurement of the strength of Yarkovsky effect. A linear relationship between $d a/d t$ and $S_1$ is confirmed. The shape index is then improved by taking the shadowing effect into account. A much better linear relationship is found between $d a/d t$ and the improved index $S_2$. This linear relationship is obeyed very well in a wide range of thermal parameter values. The influences of scattering and self-heating effects on the linear relationship are found ignorable. Using the shape index and the linear relation obtained in this paper, the rate of semi-major axis migration due to the Yarkovsky effect can be calculated accurately. Compared with the full numerical modeling of surface temperature and then the thermal radiation on an irregularly shaped asteroid, it is very easy to compute the shape index, which brings great convenience to the estimation of Yarkovsky effect.

The surface of Europa experiences a competition between thermally-induced crystallization and radiation-induced amorphization processes, leading to changes of its crystalline structure. The non-linear crystallization and temperature-dependent amorphization rate, incorporating ions, electrons and UV doses, are integrated into our multiphysics surface model (MSM) LunaIcy, enabling simulations of these coupled processes on icy moons. Thirty simulations spanning 100 000 years, covering the full ranges of albedo and latitude values on Europa, explore the competition between crystallization and irradiation. This is the first modeling of depth-dependent crystallinity profiles on icy moons. The results of our simulations are coherent with existing spectroscopic studies of Europa, both methods showing a primarily amorphous phase at the surface, followed by a crystalline phase after the first millimeter depth. Our method provides quantitative insights into how various parameters found on Europa can influence the subsurface crystallinity profiles. Interpolating upon our simulations, we have generated a crystallinity map of Europa showing, within the top millimeter, highly crystalline ice near the equator, amorphous ice at the poles, and a mix of the two at mid-latitudes. Regions/depths with balanced competition between crystallization and amorphization rates are of high interest due to their periodic fluctuations in crystalline fraction. Our interpolated map reveals periodic variations, with seasonal amplitudes reaching up to 35% of crystalline fraction. These variations could be detected through spectroscopy, and we propose a plan to observe them in forthcoming missions.

We present the one-dimensional Lyman-alpha forest power spectrum measurement derived from the data release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI). The measurement of the Lyman-alpha forest power spectrum along the line of sight from high-redshift quasar spectra provides information on the shape of the linear matter power spectrum, neutrino masses, and the properties of dark matter. In this work, we use a Fast Fourier Transform (FFT)-based estimator, which is validated on synthetic data in a companion paper. Compared to the FFT measurement performed on the DESI early data release, we improve the noise characterization with a cross-exposure estimator and test the robustness of our measurement using various data splits. We also refine the estimation of the uncertainties and now present an estimator for the covariance matrix of the measurement. Furthermore, we compare our results to previous high-resolution and eBOSS measurements. In another companion paper, we present the same DR1 measurement using the Quadratic Maximum Likelihood Estimator (QMLE). These two measurements are consistent with each other and constitute the most precise one-dimensional power spectrum measurement to date, while being in good agreement with results from the DESI early data release.

Recurrent novae undergo thermonuclear-powered eruptions separated by less than 100 years, enabled by subgiant or red giant donors transferring hydrogen-rich matter at very high rates onto their massive white dwarf companions. The most-rapidly moving parts of envelopes ejected in successive recurrent nova events are predicted to overtake and collide with the slowest ejecta of the previous eruption, leading to the buildup of vast (~ 10 - 100 parsec) super-remnants surrounding all recurrent novae; but only three examples are currently known. We report deep narrowband imaging and spectroscopy which has revealed a ~ 70-parsec-diameter shell surrounding the frequently recurring nova RS Ophiuchi. We estimate the super-remnant mass to be ~ 20 - 200 solar masses, expanding at a few tens of km/s, with an age of order 50-100 kyr. Its extremely low surface brightness and large angular size help explain the hitherto surprising absence of nova super-remnants. Our results support the prediction that ALL recurrent novae are surrounded by similar extended structures.

In the first part of this work, we provide a curated overview of the theoretical framework necessary for incorporating dephasing due to environmental effects (EE) in gravitational wave (GW) templates. We focus in particular on the relationship between orbital perturbations in the time-domain and the resulting dephasing in both time and frequency domain, elucidating and resolving some inconsistencies present in the literature. We discuss how commonly studied binary environments often result in several sources of dephasing that affect the GW signal at the same time. This work synthesizes insights from two decades of literature, offering a unified conceptual narrative alongside a curated reference of key formulas, illustrative examples and methodological prescriptions. It can serve both as a reference for researchers in the field as well as a modern introduction for those who wish to enter it. In the second part, we derive novel aspects of dephasing for eccentric GW sources and lay the foundations for consistently treating the full problem. Our results highlight the unique potential of modelling and searching for EE in eccentric binary sources of GWs.

Pulsar timing at low frequencies offers a powerful tool for studying the interstellar medium. Additionally, pulsar observations in the ecliptic enables us to study the effects of the solar wind which becomes much more prominent at low radio frequencies. The Irish station of the LOw Frequency ARray (I-LOFAR) is a sensitive low-frequency radio telescope, capable of delivering high-precision data for pulsar studies. We present a comprehensive dataset of times-of-arrival, timing solutions and dispersion measure (DM) time series for seven ecliptic pulsars observed over two-to-three years with I-LOFAR. The primary objectives are to investigate time-dependent dispersion effects and provide high-precision timing data for pulsar timing experiments. We measure DM variations through pulsar timing and analysed these across different ecliptic latitudes to assess the impact of the solar wind on each pulsar. We model the intrinsic pulse-profile variability as a function of frequency. The high-precision DM time series for all seven pulsars exhibit clear variations dependent on their ecliptic latitudes, revealing the impact of the solar wind. Some pulsars show significant changes in their pulse widths across the frequency band, while others remain stable. We examine and quantify the pulse-nulling present in PSR J0826+2637, we report evidence for DM chromaticity in PSR J1645-0317, and we describe how PSR J2145-0750's DM precision is such that it could resolve the ionospheric DM contribution. This makes it a target of interest for telescopes in areas of the globe where the ionospheric electron density is higher, e.g. the Murchison Radio Observatory in Australia. This data release underscores the potential of I-LOFAR, or any standalone international LOFAR station, for advancing low-frequency pulsar studies, particularly in analyses of dispersion in the interstellar medium, the solar wind and the ionosphere.

Debris discs provide valuable insights into the formation and evolution of exoplanetary systems. Their structures are commonly attributed to planetary perturbations, serving as probes of as-yet-undetected planets. However, most studies of planet-debris disc interactions ignore the disc's gravity, treating it as a collection of massless planetesimals. Here, using an analytical model, we investigate how the vertical structure of a back-reacting debris disc responds to secular perturbations from an inner, inclined planet. Considering the disc's axisymmetric potential, we identify two dynamical regimes: planet-dominated and disc-dominated, which may coexist, separated by a secular-inclination resonance. In the planet-dominated regime ($M_d/m_p\ll1$), we recover the classical result: a transient warp propagates outward until the disc settles into a box-like structure centered around the planetary orbit's initial inclination $I_p(0)$, with a distance-independent aspect ratio $\mathcal{H}(R)\approx I_p(0)$. In contrast, in the disc-dominated regime ($M_d/m_p\gtrsim1$), the disc exhibits dynamical rigidity, remaining thin and misaligned, with significantly suppressed inclinations and a sharply declining aspect ratio, $\mathcal{H}(R)\propto I_p(0)R^{-7/2}$. In the intermediate regime ($M_d/m_p\lesssim1$), the system exhibits a secular-inclination resonance, leading to long-lived, warp-like structures and a bimodal inclination distribution, containing both dynamically hot and cold populations. We provide analytic formulae describing these effects as a function of system parameters. We also find that the vertical density profile is intrinsically non-Gaussian and recommend fitting observations with non-zero slopes of $\mathcal{H}(R)$. Our results may be used to infer planetary parameters and debris disc masses based on observed warps and scale heights, as demonstrated for HD110058 and $\beta$ Pic.

I share fond memories of my former PhD advisor Alexei Starobinsky with whom I was closely associated for nearly 45 years. I reflect upon my early years in Moscow when I worked with him on my thesis, and touch upon the seminal work on inflation which he did during that period. Alexei visited India often and actively interacted with Indian scientists and students on issues relating to inflation, large scale structure and dark energy. This extensive collaboration, which lasted several decades, resulted in the publication of over a dozen important papers, several PhD's, and the development of the Statefinder and Om diagnostics, which I briefly discuss.

In this work, we systematically study gravitational wave (GW) production during both the inflationary and post-inflationary epochs. While inflationary GWs can be readily derived from tensor perturbations during inflation, post-inflationary GWs arise from a variety of processes during reheating and require detailed treatment for quantitative analysis. We consider four distinct production channels: $(i)$ pure inflaton annihilation, $(ii)$ graviton bremsstrahlung from inflaton decay, $(iii)$ radiation-catalyzed inflaton-graviton conversion, and $(iv)$ scattering among fully thermalized radiation particles. For each channel, we solve the corresponding Boltzmann equation to obtain the GW spectrum and derive a simple yet accurate analytical expression for it. By employing a consistent treatment of all production channels, our analysis yields for the first time the full spectrum of GWs produced during the inflationary and post-inflationary epochs. We find that, while inflationary GWs dominate at low frequencies, post-inflationary processes generally produce high-frequency GWs with considerably high energy densities that may significantly exceed that of inflationary GWs.

The late time response of vacuum black holes in General Relativity is notoriously governed by power law tails arising from the wave scattering off the curved spacetime geometry far from the black hole. While it is known that such tails are universal to a certain extent, a precise characterization of their key ingredients is missing. Here we provide an analytical proof that the tail fall-off is universal for any effective potential asymptotically decaying as $1/r^2$, while the power law decay is different if the potential decays as $1/r^\alpha$ with $1<\alpha<2$. This result extends and revises some previous work and is in agreement with numerical analyses. Our proof is based on an analytical evaluation of the branch cut contribution to the Green function, and includes charged black holes, different kinds of perturbations, Teukolsky equation for the Kerr metric, exotic compact objects, extensions of General Relativity, and environmental effects. In the latter case, our results indicate that tails are largely insensitive to a wide range of physically motivated matter distributions around black holes, including the Navarro Frenk White profile commonly used to model dark matter.

We analyse properties of general stationary and axisymmetric spacetimes, with a particular focus on circularity -- an accidental symmetry enjoyed by the Kerr metric, and therefore widely assumed when searching for rotating black hole solutions in alternative theories of gravity as well as when constructing models of Kerr mimickers. Within a gauge specified by seven (or six) free functions, the local existence of which we prove, we solve the differential circularity conditions and translate them into algebraic relations among the metric components. This result opens the way to investigating the consequences of circularity breaking in a controlled manner. In particular, we construct two simple analytical examples of non-circular deformations of the Kerr spacetime. The first one is "minimal", since the horizon and the ergosphere are identical to their Kerr counterparts, except for the fact that the horizon is not Killing and its surface gravity is therefore not constant. The second is "not so minimal", as the horizon's profile can be chosen arbitrarily and the difference between the horizon and the so-called rotosurface can be appreciated. Our findings thus pave the way for further research into the phenomenology of non-circular stationary and axisymmetric spacetimes.

The characterization of exoplanetary atmospheres through spectral analysis is a complex challenge. The NeurIPS 2024 Ariel Data Challenge, in collaboration with the European Space Agency's (ESA) Ariel mission, provided an opportunity to explore machine learning techniques for extracting atmospheric compositions from simulated spectral data. In this work, we focus on a data-centric business approach, prioritizing generalization over competition-specific optimization. We briefly outline multiple experimental axes, including feature extraction, signal transformation, and heteroskedastic uncertainty modeling. Our experiments demonstrate that uncertainty estimation plays a crucial role in the Gaussian Log-Likelihood (GLL) score, impacting performance by several percentage points. Despite improving the GLL score by 11%, our results highlight the inherent limitations of tabular modeling and feature engineering for this task, as well as the constraints of a business-driven approach within a Kaggle-style competition framework. Our findings emphasize the trade-offs between model simplicity, interpretability, and generalization in astrophysical data analysis.

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

In warmer climates, cities face the prospect of higher air temperatures in summer compared to historical averages due to the urban heat island effect. An approach intended to address this problem is the application of "cool pavement" treatments (CPT) to city streets to make them more reflective to sunlight. Raising the albedo of roadways for this purpose may also have the effect of increasing the amount of street light that is reflected into the night sky. The simplest hypothesis explaining the relationship between CPT application and upward radiance is that CPT applied to road surfaces in areas where street lighting is dominant should increase the upward radiance of neighborhoods where the treatments are applied. A simple model predicted radiance increases of 2-6% immediately after CPT application. To test the hypothesis and model predictions, we looked for radiance changes coinciding with the application of CPT in residential neighborhoods of Phoenix, U.S., since 2020. We obtained time series radiance measurements from Visible Infrared Imaging Radiometer Suite Day-Night Band (VIIRS-DNB) data from Radiance Light Trends for Phoenix neighborhoods receiving CPT and nearby "control" neighborhoods not receiving CPT. At the 95% confidence level, we found that any increases in nighttime radiances from treated neighborhoods did not exceed about 14%. As a consequence, we cannot rule out either the expected radiance increases from our model or the possibility that CPT application yielded no change in radiance. We therefore cannot draw robust conclusions about the potential influence of CPT deployment on skyglow given the limitations of the DNB as a data source.

Boson stars, hypothetical astrophysical objects bound by the self-gravity of a scalar field, have been widely studied as a type of exotic compact object that is horizonless and provides a testing ground for physics beyond the Standard Model. In particular, many previous works have demonstrated methods for distinguishing compact boson stars from black holes in general relativity through gravitational wave observations. However, the formation scenario of compact boson stars within the age of the universe remains unclear. In this paper, we explore a possible scenario for the formation of compact boson stars. The model we consider requires two coupled scalar fields: a complex scalar field that forms a boson star and a spatially homogeneous background field, as formation of a compact boson star cannot be achieved in a single filed model. Using the adiabatic approximation, we show that non-relativistic boson clouds can evolve into compact boson stars through the cosmological time-evolution of the background field. In our model the background field evolves to increase the effective mass of the scalar field, and as a result compact boson stars can form within the cosmological timescale, if the variation of the background field is as large as the Planck scale. However, further investigation is required because the required initial states are not the configurations that can be described by the well-studied Schrödinger-Poisson system.

Interplanetary coronal mass ejections (ICMEs) are major drivers of space weather disturbances, posing risks to both technological infrastructure and human activities. Automatic detection of ICMEs in solar wind in situ data is essential for early warning systems. While several methods have been proposed to identify these structures in time series data, robust real-time detection remains a significant challenge. In this work, we present ARCANE - the first framework explicitly designed for early ICME detection in streaming solar wind data under realistic operational constraints, enabling event identification without requiring observation of the full structure. Our approach evaluates the strengths and limitations of detection models by comparing a machine learning-based method to a threshold-based baseline. The ResUNet++ model, previously validated on science data, significantly outperforms the baseline, particularly in detecting high-impact events, while retaining solid performance on lower-impact cases. Notably, we find that using real-time solar wind (RTSW) data instead of high-resolution science data leads to only minimal performance degradation. Despite the challenges of operational settings, our detection pipeline achieves an F1 score of 0.53, with an average detection delay of 21.5% of the event's duration while only seeing a minimal amount of data. As more data becomes available, the performance increases significantly. These results mark a substantial step forward in automated space weather monitoring and lay the groundwork for enhanced real-time forecasting capabilities.

Important scientific discoveries should be backed by high statistical significance. In the 2030s, multiple space-based gravitational wave detectors are expected to operate. While many works aim to achieve quick and reliable detection and parameter estimation of millihertz gravitational wave sources, dedicated studies are lacking to assess the significance of space-based detectors. In this work, we propose a framework to assess the statistical significance of massive black hole binaries (MBHBs) detections with space-based gravitational wave detectors. We apply this algorithm to simulated data with Gaussian stationary noise and the complex LDC-2a dataset to measure the false alarm rate and significance of MBHB signals. We also analyze factors affecting the significance of MBHBs and design a method to mitigate multi-source confusion interference. In Gaussian noise conditions, MBHBs with a signal-to-noise ratio of about 7 can achieve $3 \sigma$ significance, and those with a signal-to-noise ratio of about 8 achieve $4 \sigma$. Our analysis demonstrates that all MBHB signals in the LDC-2a dataset have a significance exceeding $4.62 \sigma$.

We present a unified Bayesian framework to jointly constrain the Hubble constant $H_0$ and the post-Newtonian parameter $\gamma$, a key probe of deviations from general relativity, using the population characteristics of strongly lensed gravitational wave (GW) events from binary black hole mergers. Unlike traditional methods that rely on electromagnetic counterparts or GW waveform modeling, our approach exploits the time-delay distribution and the total number of lensed events, achievable with third-generation detectors such as the Einstein Telescope. Assuming a flat $\Lambda$CDM cosmology, we demonstrate that this method can achieve precision levels of $0.4\% - 0.7\%$ for$ H_0$ and $0.5\% - 3.3\%$ for $\gamma$ at $68\%$ credibility, significantly outperforming existing joint constraints. These results underscore the power of lensed GW population statistics as a robust and efficient probe of both cosmic expansion and the nature of gravity.

Accurate and efficient modeling of the Laser Interferometer Space Antenna (LISA) response is crucial for gravitational-wave (GW) data analysis. A key computational challenge lies in evaluating time-delay interferometry (TDI) variables, which require projecting GW polarizations onto the LISA arms at different retarded times. Without approximations, the full LISA response is computationally expensive, and traditional approaches, such as the long-wavelength approximation, accelerate the response calculation at the cost of reducing accuracy at high frequencies. In this work, we introduce a novel hybrid time-domain response for LISA that balances computational efficiency and accuracy across the binary's evolution. Our method implements a fast low-frequency approximation during the early inspiral$\unicode{x2013}$where most binaries spend most of the time in the sensitive frequency band of LISA$\unicode{x2013}$while reserving the computationally intensive full-response calculations for the late inspiral, merger, and ringdown phases. The low-frequency approximation (LFA) is based on Taylor expanding the response quantities around a chosen evaluation time such that time delays correspond to central finite differences. Our hybrid approach supports CPU and GPU implementations, TDI generations 1.5 and 2.0, and flexible time-delay complexity, and has the potential to accelerate parts of the global fit and reduce energy consumption. We also test our LFA and hybrid responses on eccentric binaries, and we perform parameter estimation for a "golden" binary. Additionally, we assess the efficacy of our low-frequency response for "deep alerts" by performing inspiral-only Bayesian inference.