Papers for Wednesday, Apr 16 2025

Context: Fast Radio Bursts (FRBs) are enigmatic extragalactic bursts whose properties are still largely unknown, but based on their extremely small time duration, they are proposed to have a compact structure, making them candidates for wave-optical effects if gravitational lensed. If an FRB is lensed into multiple-images bursts at different times by a galaxy or cluster, a likely scenario is that only one image is detected, because the others fall outside the survey area and time frame. Aims: In this work we explore the FRB analog of quasar microlensing, namely the collective microlensing by stars in the lensing galaxy, now with wave optics included. The eikonal regime is applicable here. Methods. We study the voltage (rather than the intensity) in a simple simulation consisting of (a) microlensing stars, and (b) plasma scattering by a turbulent interstellar medium. Results: The auto-correlation of the voltage shows peaks (at order-microsecond separations) corresponding to wave-optical interference between lensed micro-images. The peaks are frequency dependent if plasma-scattering is significant. While qualitative and still in need of more realistic simulations, the results suggest that a strongly-lensed FRB could be identified from a single image. Conclusions: Microlensing could sniff out macro-lensed FRBs

The characterization and analysis of light curves are vital for understanding the physical and rotational properties of artificial space objects such as satellites, rocket stages, and space debris. This paper introduces the Light Curve Dataset Creator (LCDC), a Python-based toolkit designed to facilitate the preprocessing, analysis, and machine learning applications of light curve data. LCDC enables seamless integration with publicly available datasets, such as the newly introduced Mini Mega Tortora (MMT) database. Moreover, it offers data filtering, transformation, as well as feature extraction tooling. To demonstrate the toolkit's capabilities, we created the first standardized dataset for rocket body classification, RoBo6, which was used to train and evaluate several benchmark machine learning models, addressing the lack of reproducibility and comparability in recent studies. Furthermore, the toolkit enables advanced scientific analyses, such as surface characterization of the Atlas 2AS Centaur and the rotational dynamics of the Delta 4 rocket body, by streamlining data preprocessing, feature extraction, and visualization. These use cases highlight LCDC's potential to advance space debris characterization and promote sustainable space exploration. Additionally, they highlight the toolkit's ability to enable AI-focused research within the space debris community.

Strongly lensed supernovae are a promising new probe to obtain independent measurements of the Hubble constant (${H_0}$). In this work, we employ simulated gravitationally lensed Type Ia supernovae (glSNe Ia) to train our machine learning (ML) pipeline to constrain $H_0$. We simulate image time-series of glSNIa, as observed with the upcoming Nancy Grace Roman Space Telescope, that we employ for training an ensemble of five convolutional neural networks (CNNs). The outputs of this ensemble network are combined with a simulation-based inference (SBI) framework to quantify the uncertainties on the network predictions and infer full posteriors for the $H_0$ estimates. We illustrate that the combination of multiple glSN systems enhances constraint precision, providing a $4.4\%$ estimate of $H_0$ based on 100 simulated systems, which is in agreement with the ground truth. This research highlights the potential of leveraging the capabilities of ML with glSNe systems to obtain a pipeline capable of fast and automated $H_0$ measurements.

We use the precision measurements of the arrival time differences of the same fast radio burst (FRB) source along multiple sightlines to measure the primordial power spectrum and Non-Gaussianities. The anticipated experiment requires a sightline separation of 100 AU, achieved by sending three or more radio telescopes to the outer solar system. The Shapiro time delays, measured relatively between different telescopes, are sensitive to the gradient field of the gravitational potential between different sightlines. Since the arrival time difference is independent of when the transient signal is emitted from the source, every measurement of the detected FRB source can be correlated. With enough FRB sources discovered, we can map the gravitational potential across the sky. We further calculate the two-point and three-point correlation function of the arrival time difference between telescopes for different FRB sources in the sky. If $10^4$ FRBs were to be detected, our results suggest that this technique can test the inflationary scale-invariant power spectrum down to $\sim 10^3\,\rm Mpc^{-1}$ and primordial Non-Gaussianities at a level of $f_{\rm NL}\sim 1$.

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: "Cosmic Reionization On Computers" (CROC) and "Thesan". While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the radiation field in two simulations proceeds through the same evolutionary stages, but the timing of these stages is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

Radio-bright, jetted quasars at $z>5$ serve as unique laboratories for studying supermassive black hole activity in the early Universe. In this work, we present a sample of high-$z$ jetted quasars selected from the combination of the radio Rapid ASKAP Continuum Survey (RACS) with deep wide-area optical/near-infrared surveys. From this cross-match we selected 45 new high-$z$ radio quasar candidates with S$_{888MHz}>1$ mJy and mag$z<21.3$ over an area of 16000deg$^2$. Using spectroscopic observations, we confirmed the high-$z$ nature of 24 new quasars, 13 at $4.5<z<5$ and 11 at $z>5$. If we also consider similar, in terms of radio/optical fluxes and sky position, quasars at $z>5$ already reported in the literature, the overall $z>5$ RACS sample is composed by 33 powerful quasars, expected to be ~90% complete at mag$z<21.3$ and S$_{888MHz}>1$ mJy. Having rest-frame radio luminosities in the range $\nu L_{1.4GHz}=10^{41.5}-10^{44.4}$ erg s$^{-1}$, this sample contains the most extreme radio quasars currently known in the early Universe. We also present all X-ray and radio data currently available for the sample, including new, dedicated {\it Chandra}, uGMRT, MeerKAT and ATCA observations for a sub-set of the sources. from the modelling of their radio emission, either with a single power law or a broken power law, we found that these systems have a wide variety of spectral shapes with most quasars (22) having a flat radio emission (i.e., $-0.5<\alpha_{r}<0.5$). At the same time, the majority of the sources with X-ray coverage present a high-energy luminosity larger than the one expected from the X-ray corona only. Both the radio and X-ray properties of the high-$z$ RACS sample suggest that many of these sources have relativistic jets oriented close to our line of sight. (i.e., blazars) and can therefore be used to perform statistical studies on the entire jetted population at high redshift.

The baryonic features in the galaxy power spectrum offer tight, time-resolved constraints on the expansion history of the Universe but complicate the measurement of the broadband shape of the power spectrum, which also contains precious cosmological information. In this work we compare thirteen methods designed to separate the broadband and oscillating components and examine their performance. The systematic uncertainty between different de-wiggling procedures is at most $2$%, depending on the scale. The ShapeFit parameter compression aims to compute the slope $m$ of the power spectrum at large scales, sensitive to matter-radiation equality and the baryonic suppression. We show that the de-wiggling procedures impart large (50%) differences on the obtained slope values, but as long as the theory and data pipelines are set up consistently, this is of no concern for cosmological inference given the precision of existing and on-going surveys. However, it still motivates the search for more robust ways of extracting the slope. We show that post-processing the power spectrum ratio before taking the derivative makes the slope values far more robust. We further investigate eleven ways of extracting the slope and highlight the two most successful ones. We derive a systematic uncertainty on the slope $m$ of $\sigma_{m,\mathrm{syst}} = 0.023|m| + 0.001$ by studying the behavior of the slopes in different cosmologies within and beyond $\Lambda$CDM and the impact in cosmological inference. In cosmologies with a feature in the matter-power spectrum, such as in the early dark energy cosmologies, this systematic uncertainty estimate does not necessarily hold, and further investigation is required.

Direct collapse of pristine gas in early galaxies is a promissing pathway for forming supermassive black holes (BHs) powering active galactic nuclei (AGNs) at the epoch of reionization (EoR). This seeding mechanism requires suppression of molecular hydrogen (H$_2$) cooling during primordial star formation via intense far-ultraviolet radiation from nearby starburst galaxies clustered in overdense regions. However, non-detection of 21 cm signals from the EoR reported by the Hydrogen Epoch of Reionization Array (HERA) experiment suggests that such galaxies may also emit X-rays more efficiently than in the local universe, promoting H$_2$ production and thereby potentially quenching massive BH seed formation. In this study, we examine the thermal and chemical evolution of collapsing gas in dark matter halos using a semi-analytic model incorporating observationally calibrated X-ray intensities. We find that strong X-ray irradiation, as suggested by HERA, significantly suppresses direct collapse and leads most halos to experience H$_2$ cooling. Nevertheless, massive BH seeds with $M_\mathrm{BH} \gtrsim 10^4~M_\odot$ still form by $z\simeq 15$, particularly in regions with baryonic streaming motion, and their abundance reaches $\sim 10^{-4}~\mathrm{Mpc}^{-3}$ sufficient to explain the SMBHs identified by JWST spectroscopy at $3<z<6$. While the formation of highly overmassive BHs with masses comparable to their host galaxies is prohibited by X-ray ionization, our model predicts that BH-to-stellar mass ratios of $\simeq 0.01-0.1$ were already established at seeding.

In the coming age of precision neutrino physics, neutrinos from the Sun become robust probes of the conditions of the solar core. Here, we focus on $^8$B neutrinos, for which there are already high precision measurements by the Sudbury Neutrino Observatory and Super-Kamiokande. Using only basic physical principles and straightforward statistical tools, we calculate projected constraints on the temperature and density of the $^8$B neutrino production zone compared to a reference solar model. We outline how to better understand the astrophysics of the solar interior using forthcoming neutrino data and solar models.

Dust grains play a fundamental role in galaxies, influencing both their evolution and observability. As a result, incorporating dust physics into galaxy evolution simulations is essential. This is a challenging task due to the finite resolution of such simulations and the uncertainties on dust grain formation in stellar envelopes and their evolution in the Interstellar Medium (ISM). This report reviews some of the most commonly used techniques for modeling dust in galaxy evolution simulations, with a particular emphasis on developments from the past $\sim$10 years in both hydrodynamic and semi-analytic approaches. Key findings from these simulations are presented, discussed, and compared to the most recent available observations. These include the dust-to-gas vs metallicity relation, the abundance of dust within and outside galaxies, and predictions on the relative importance of various processes affecting dust. The analysis presented here highlights significant achievements as well as the limitations in our current modeling of dust evolution in galaxies.

The gravitational capture of a small compact object by a supermassive black hole is one of the most intriguing sources of gravitational waves to be detected by space-borne observatories. Modeling gravitational waves is a challenging task. However, various approximations exist that enable us to claim detection and perform parameter extraction with high accuracy on synthetic data. The first numerical implementation of relativistic corrections for an astrophysical system was conducted in a post-Newtonian (PN) framework, and since then, most comparable programmes have followed a similar approach. Nevertheless, the PN approach has been developed for two particles in a vacuum. Here, we present the first results on the impact of cross-terms in the PN scheme for ``monochromatic'' sources (meaning that the peak frequency does not evolve over the observational time), which account for the interplay between different bodies. We present a conceptual and illustrative study of the cross-terms on the first PN term in a group of three bodies representing asymmetric binaries (X-MRIs, E-EMRIs, and X-MRIs with E-EMRIs). We find that these cross-terms can lead to a complete phase shift of the gravitational wave by a few times $2\pi$ over periods of time of less than one year for the close binaries. Ignoring the cross-terms in the relativistic semi-Keplerian system can significantly complicate parameter extraction for monochromatic sources.

We present the 4MOST IR AGN survey, the first large-scale optical spectroscopic survey characterizing mid-infrared (MIR) selected obscured active galactic nuclei (AGN). The survey targets $\approx 212,000$ obscured infrared (IR) AGN candidates over $\approx 10,000 \rm \: deg^2$ down to a magnitude limit of $r_{\rm AB}=22.8 \, \rm mag$ and will be $\approx 100 \times$ larger than any existing obscured IR AGN spectroscopic sample. We select the targets using a MIR colour criterion applied to the unWISE catalogue from the WISE all-sky survey, and then apply a $r-W2\geq 5.9 \rm \: mag$ cut; we demonstrate that this selection will mostly identify sources obscured by $N_{\rm H}>10^{22} \rm \: cm^{-2}$. The survey complements the 4MOST X-ray survey, which will follow up $\sim 1\rm M$ eROSITA-selected (typically unobscured) AGN. We perform simulations to predict the quality of the spectra that we will obtain and validate our MIR-optical colour-selection method using X-ray spectral constraints and UV-to-far IR spectral energy distribution (SED) modelling in four well-observed deep-sky fields. We find that: (1) $\approx 80-87\%$ of the WISE-selected targets are AGN down to $r_{\rm AB}=22.1-22.8 \: \rm mag$ of which $\approx 70\%$ are obscured by $N_{\rm H}>10^{22} \: \rm cm^{-2}$, and (2) $\approx 80\%$ of the 4MOST IR AGN sample will remain undetected by the deepest eROSITA observations due to extreme absorption. Our SED fitting results show that the 4MOST IR AGN survey will primarily identify obscured AGN and quasars ($\approx 55\%$ of the sample is expected to have $L_{\rm AGN,IR}>10^{45} \rm \: erg \: s^{-1}$) residing in massive galaxies ($M_{\star}\approx 10^{10}-10^{12} \rm \: M_{\odot}$) at $z\approx 0.5-3.5$ with $ \approx 33\%$ expected to be hosted by starburst galaxies.

We study the launching of magnetized jets from a resistive circumstellar disk within a binary system, employing a unique combination of 3D MHD jet launching simulations (PLUTO code) and post-processed 3D radiative transfer modeling (RADMC-3D code). Our findings reveal a well-defined jet originating from the inner region of the disk, extending to a larger disk area. While the model attains steady states for a single star, a binary system leads to the emergence of tidal effects such as the formation of ``spiral arms'' in the disk and inside the jet. Here we have consistently implemented a time-dependent Roche potential for the gravity of the binary. As a major step forward, we further present the first 3D radiation maps of the dust continuum for the disk-jet structure. In principle, this allows us to compare MHD simulation results to observed disk-outflow features. We, therefore, present convolved images of the dust continuum emission, employing exemplary point spread functions of the MIRI instrument (5~$\mu m$ band) and the ALMA array (320~$\mu m$ band). In these bands, we identify distinguishable features of the disk-jet structure, such as "spiral arms," which we have also seen in the MHD this http URL gas density increased by an order of magnitude, the disk become optically thick at 5~$\mu m$, but remains bright at 320~$\mu$m. At this wavelength, 320~$\mu$m, enhanced structural features in the disk and the base of the wind become more pronounced and are well resolved in the convolved image.

The dynamics of a binary neutron stars merger is governed by physics under the most extreme conditions, including strong spacetime curvature, ultra-high matter densities, luminous neutrino emission and the rapid amplification of the initial neutron star magnetic fields. Here we systematically explore how sensitive the magnetic field evolution is to the total mass of the merging binary, to the mass ratio of its components, the stellar spins and to the equation of state. For this purpose, we analyze 16 state-of-the-art GRMHD simulations that employ a subgrid-scale model to account for the unresolved small-scale turbulence. We find that strong and rapid amplification of the magnetic field to volume-averaged values of $\sim 10^{16}$~G in the high-density regions is a very robust outcome of a neutron star merger and this result is only marginally impacted by either mass, mass ratio, spin or equation of state.

The unusual light curve of the massive eclipsing binary W Ser was recently observed with high S/N and fast cadence by the NASA TESS mission. The TESS light curve records two eclipses and relatively fast variations outside of the eclipses. The eclipse timings verify the period increase of the binary, and the period derivative implies a mass transfer rate in excess of 10^{-5} solar masses per year.. The light curve shows a fading trend from just after an eclipse until the start of the next eclipse. The brightest flux source in the system is the accretion torus surrounding the mass gainer star, and we argue that these orbital-phase related fadings are the result of the injection of cooler gas from the mass donor entering the outskirts of the accretion torus. There are cyclic variations in the out-of-eclipse sections of the light curve that vary on a 2.8 day timescale. This equals the orbital period for gas in the outer regions of the accretion torus, so the photometric variations are probably the result of transitory, over-dense regions that form at the rim of the accretion torus.

Recent high-cadence flare campaigns by the Interface Region Imaging Spectrograph (IRIS) have offered new opportunities to study rapid processes characteristic of flare energy release, transport, and deposition. Here, we examine high-cadence chromospheric and transition region spectra acquired by IRIS during a C-class flare from 2022 September 25. Within the flare ribbon, the intensities of the Si IV 1402.77, C II 1334.53 and Mg II k 2796.35 lines peaked at different times, with the transition region Si IV typically peaking before the chromospheric Mg II line by 1 - 6 seconds. To understand the nature of these delays, we probed a grid of radiative hydrodynamic flare simulations heated by electron beams, thermal conduction-only, or Alfvén waves. Electron beam parameters were constrained by hard X-ray observations from the Gamma-ray Burst Monitor (GBM) onboard the Fermi spacecraft. Reproducing lightcurves where Si IV peaks precede those in Mg II proved to be a challenge, as only a subset of Fermi/GBM-constrained electron beam models were consistent with the observations. Lightcurves with relative timings consistent with the observations were found in simulations heated by either high-flux electron beams or by Alfvén waves, while the thermal conduction heating does not replicate the observed delays. Our analysis shows how delays between chromospheric and transition region emission pose tight constraints on flare models and properties of energy transport, highlighting the importance of obtaining very high-cadence datasets with IRIS and other observatories.

We conducted a $369\;{\rm ks}$ NuSTAR observation on the X-ray pulsar Centaurus X-3, which covered two consecutive orbital cycles of the source, including two eclipse durations. We investigated the orbital-phase spectral variability over the two orbital cycles. We divided the entire observation data into multiple segments, each covering an orbital interval of $\Delta\Phi=0.005$. The phenomenological spectral modeling applied to these orbital-phase-resolved spectra reveals that the photon index is the key parameter with the most variability and a strong correlation with the continuum flux. The photon index becomes softer during the high-flux phases and harder in the low-flux phases. The relation between the photon index and continuum flux remains consistent when investigating specific spin phases, suggesting that the spectral variability originates from extrinsic factors apart from the neutron star. Furthermore, the 3-5 keV pulse fraction also exhibits variability, being enhanced in the high-flux phases and suppressed in the low-flux phases, which indicates the presence of multiple emission components with different pulse fractions. These phenomenological analysis results enabled us to estimate the physical origin of the spectral variability. We successfully fitted the orbital-phase-resolved spectra with a physical model that assumes (1) stable emission from the neutron star, (2) attenuation by inhomogeneous, clumpy stellar wind, and (3) an additional non-pulsed emission component arising from thermal emission from the accretion disk. The thermal emission from the accretion disk can be described by a blackbody with a temperature of $kT\sim0.5\;{\rm keV}$ and a luminosity of $\sim10^{37}\;{\rm erg\;s^{-1}}$.

Enstatite chondrites (EC) are potential source material for the accretion of Mercury due to their reduced nature and enrichment in volatile elements. Understanding their melting properties is therefore important to better assess a scenario where Mercury formed from these chondrites. Here, we present experimental data on the partial melting of a modified EH4 Indarch EC, which was adjusted to have 18\% more metallic Si than SiO$_2$ in mass, yielding an oxygen fugacity of 3.7 below the iron--wüstite redox buffer and 12 wt\% Si in the metal. Experiments were performed from 0.5 to 5 GPa. Results indicate that the stability field of enstatite expands relative to olivine. This expansion is likely due to the presence of Ca--S and Mg--S complexes in the silicate melt, which enhance SiO$_2$ activity and promote enstatite crystallization. Additionally, sulfides show enrichment in Mg and Ca, up to 22 and 13 wt\% respectively, the main remaining cations being Fe, Cr, and Mn. These high Mg and Ca contents are observed at low temperatures and high silica content in the silicate melt, respectively. High-pressure melts (2 to 5 GPa, 160--400 km depth in Mercury) are Mg-rich, similar to those in Mercury's high-magnesium region (HMR), while low-pressure melts (0.5 to 1 GPa, 40--80 km depth) are Si-rich, comparable to the northern volcanic plains (NVP). Results suggest that a large fraction of Mercury's surface aligns compositionally with these melts, implying that Mercury's mantle could predominantly have a pyroxenitic composition. However, regions with differing compositions, such as aluminum-rich areas like the Caloris basin, suggest local variability in mantle geochemistry. Overall, our results show that if Mercury formed from materials similar to EC, batch melting of its primitive pyroxenite mantle would yield magmas with compositions resembling those of most rocks observed on the surface.

Models predict that more than half of all impacting meteoroids should be carbonaceous, reflecting the abundance of carbon-rich asteroids in the main belt and near-Earth space. Yet carbonaceous chondrites represent only about 4% of meteorites recovered worldwide. Here we analyse 7,982 meteoroid impacts and 540 potential meteorite falls from 19 global observation networks and demonstrate that intense thermal stress at low perihelion distances coupled with the filtering effect of Earth`s atmosphere explains this mismatch. Meteoroids repeatedly subjected to intense thermal cycling near the Sun fracture and weaken, removing the most friable objects even before atmospheric entry. Our data also show that tidally disrupted meteoroid streams produce especially fragile fragments that rarely survive to the ground. Consequently, compact, higher-strength, thermally cycled bodies dominate the meteorite record. These findings reconcile the predicted carbonaceous flux with its scarcity in collections, underscoring how orbital evolution and atmospheric filtering shape the materials that reach Earth`s surface.

The observed spectra and light curves of the kilonova produced by the GW170817 binary neutron star merger provide complementary insights, but self-consistently modeling both the spectral- and time-domain has proven challenging. Here, we model the optical/infrared light curves of the GW170817 kilonova, using the properties and physical conditions of the ejecta as inferred from detailed modeling of its spectra. Using our software tool SPARK, we first infer the $r$-process abundance pattern of the kilonova ejecta from spectra obtained at 1.4, 2.4, 3.4, and 4.4 days post-merger. From these abundances, we compute time-dependent radioactive heating rates and the wavelength-, time-, and velocity-dependent opacities of the ejecta. We use these inferred heating rates and opacities to inform a kilonova light curve model, to self-consistently reproduce the observed early-time light curves and to infer a total ejecta mass of $M_{\mathrm{ej}} = {0.11}~M_{\odot}$, towards the higher end of that inferred from previous studies. The combination of a large ejecta mass from our light curve modeling and the presence of both red and blue ejecta from our spectral modeling suggests the existence of a highly magnetized hypermassive neutron star remnant that survives for $\sim$$0.01 - 0.5$ s and launches a blue wind, followed by fast, red neutron-rich winds launched from a magnetized accretion disk. By modeling both spectra and light curves together, we demonstrate how combining information from both the spectral and time domains can more robustly determine the physical origins of the ejected material.

Extracting additional information from old or incomplete fireball datasets remains a challenge. To address missing point-by-point observations, we introduce a method for estimating atmospheric flight parameters of meteoroids using metaheuristic optimization techniques. Using a fireball catalog from the European Fireball Network (EN), we reconstruct velocity profiles, meteoroid bulk densities, mass loss rates, and ablation and ballistic coefficients, based on the initial and terminal points' height, velocity, and mass with the purely dynamical $\alpha$-$\beta$ model. Additionally, the method's performance is compared to the Meteorite Observation and Recovery Project (MORP) derived fits, confirming the robustness of the computed parameters for objects with asteroidal compositions. Our findings show that $\alpha$-$\beta$ model yields parameters consistent with the photometric and dynamic mass estimates in the EN catalog for $P_E$ type I events. However, in the implementation proposed here, $\alpha$-$\beta$ model encounters limitations in accurately representing the final deceleration of more fragile high-velocity meteoroids. This is likely due to challenges in representing complex fragmentation processes by fitting only two points, even when initial and terminal residuals are minimal. The retrieved $\alpha$-$\beta$ distribution differs from the one derived from MORP data, likely due to the imposed mass constraints, which strongly influence the results, especially the bulk density. The results suggest that $P_E$ constraints reduce fitting accuracy (from 90\% to 44\%), while flexibility and freedom from assumptions improve $\alpha$-$\beta$ performance. The method yields 26\% of the events compatible with the catalog $P_E$ classification. Our approach is well-suited for interpreting historical or sparse datasets.

Recent observations expanded our understanding of galaxy formation and evolution, yet key challenges persist in the X-ray regime, crucial for studying Active Galactic Nuclei (AGN). These limitations drive the development of next-generation observatories such as ESA's NewAthena. Now in phase B (preliminary design), the mission requires extensive testing to ensure compliance with its scientific goals, particularly given the uncertainties surrounding high redshift AGN. This work leverages the IllustrisTNG cosmological simulation to build an X-ray AGN mock catalogue and assess the performance of NewAthena's WFI. We created a Super Massive Black Hole (SMBH) light cone, spanning 10 deg2, with corrections to account for the limited resolution of the simulation and X-ray properties derived in post-processing. The resulting catalogue reveals a 5* overabundance of faint AGN compared to current X-ray constraints, an inconsistency potentially resolved by invoking a higher Compton-thick (CTK) fraction and intrinsic X-ray weakness, as suggested by recent JWST findings. An end-to-end survey simulation using SIXTE predicts 250000 AGN detections, including 20,000 at z > 3 and 35 in the Epoch of Reionization (z > 6); notably, only AGN with LX > 43.5 erg/s are detectable at z > 6. The analysis also forecasts a significant population of detectable CTK AGN, even beyond z > 4. These findings suggest X-ray observations will, for the first time, probe a significant AGN population in the EoR, offering new insights into SMBH growth. They also provide key input for refining NewAthena's mission design and optimizing its survey strategy.

The energy spectra of cosmic rays (CRs) below tens of GeV are significantly modulated by solar activity within the heliosphere. To investigate the properties of Galactic CRs, it is important to determine the unmodulated local interstellar spectrum (LIS). Recent high-precision temporal measurements of CR energy spectra, released by the AMS-02 collaboration, provide a crucial observational foundation for this endeavor. In this study, we employ the widely used force-field approximation (FFA) model to analyze the AMS-02 data, and attempt to derive the LIS for CR protons and positrons. By applying a non-LIS method, we derive temporal variations of the relative solar modulation potential, $\Delta\phi$, for individual particle species. Our analysis demonstrates that the FFA provides sufficient accuracy in explaining the AMS-02 spectral measurements of all particles during the low solar activity period. Notably, the derived $\Delta\phi(t)$ for protons and positrons, as well as for electrons and antiprotons, exhibit excellent consistency, indicating that particles with the same charge sign can be effectively described within a unified FFA framework during the low solar activity period. Having established a well-constrained proton LIS and its associated modulation potential, we apply the common modulation behavior between positrons and protons to demodulate the AMS-02 positron measurements, and derive the positron LIS without necessitating prior knowledge of its characteristics. This LIS is useful for quantitative investigations into potential exotic origins of CR positrons.

The Fokker-Planck equation describing the transport of energetic particles interacting with turbulence is difficult to solve analytically. Numerical solutions are of course possible but they are not always useful for applications. In the past a subspace approximation was proposed which allows to compute important quantities such as the characteristic function as well as certain expectation values. This previous approach was applied to solve the one-dimensional Fokker-Planck equation which contains only a pitch-angle scattering term. In the current paper we extend this approach in order to solve the Fokker-Planck equation with a focusing term. We employ two- and three-dimensional subspace approximations to achieve a pure analytical description of particle transport. Additionally, we show that with higher dimensions, the subspace method can be used as a hybrid analytical-numerical method which produces an accurate approximation. Although the latter approach does not lead to analytical results, it is much faster compared to pure numerical solutions of the considered transport equation.

The origin and physics of X-ray intra-day variability (IDV) in blazars, which is a long-standing issue, is studied by modelling the broad-band X-ray spectrum, the light curves (LCs), and the Fourier time lags. We present the timing analysis of three archived XMM-Newton observations with a total exposure of $>80$ ks of PKS 2155-304, which is one of the brightest and most studied HBLs in the X-ray band. For each observation, we constructed averaged X-ray spectra in 0.5-10 keV band, as well as 100 s binned LCs in various sub-bands. We performed the Bayesian power spectral density (PSD) analysis and Fourier time-lag analyses of the variable LCs. The results are carefully modelled in the context of a multi-zone jet model. PSD analysis reveals that the X-ray variability can be characterised by red noise. The lag-frequency spectra measured in two observations show only the soft or negative lags, with the magnitude of the lags increasing as the frequency decreases. For another observation, the lag-frequency spectra are characterised by small positive or zero time lags at the lowest frequencies, which drops to negative values at higher frequencies. The magnitude of the soft lags ranges from $\sim5$ to $\sim40$ minutes, and increases with the energy difference of two compared LCs. The observed X-ray spectra and lag-frequency spectra can both be successfully described by our proposed two-zone model, with the physical parameters constrained in a fully acceptable space. Moreover, the LC profiles at different energy bands can be satisfactorily reproduced by only varying the injection rate of the energetic electrons. The IDV of PKS 2155-304 should be caused by the injection of energetic electrons, and accelerated by shocks formed in a weakly magnetised jet.

Cataclysmic variables (CVs), binary systems with a white dwarf accreting from a low-mass star, are significant Galactic X-ray sources. We present a systematic search for X-ray emitting CV candidates by cross-matching four X-ray catalogs (eROSITA, XMM-Newton, Swift, and ROSAT) with Gaia sources located in the bridge region between the main sequence and white dwarf cooling sequence in the Hertzsprung-Russell diagram. From 444 candidates (267 confirmed CVs and 177 new candidates), we detect orbital modulation in 56 sources using ZTF/TESS light curves. The eROSITA catalog contributes 51% of candidates, outperforming other surveys due to its wider sky coverage and higher sensitivity (~10^{-14} erg cm^{-2} s^{-1}). Our method demonstrates the efficiency of combining X-ray data with time-domain analysis for CV identification, with future eROSITA observations expected to expand the population of X-ray emitting CVs.

Extreme eclipsing binaries may harbor peculiar physical properties. In this work, we aim to identify a sample of such systems by selecting binaries with pronounced eclipsing light curves, characterized of large variability ($\Delta \mathrm{mag} > 0.3$ in ZTF $g$ band) and significant differences between primary and secondary eclipses (eclipse depth ratio $>$ 20 in ZTF $g$ band). We identified 23 candidates by combining the photometric data and the LAMOST spectroscopic survey. Spectroscopic analysis revealed that all of these systems are dominated by A-type stars in the optical band. Further investigation confirmed that all 23 candidates are Algol-type binaries, with 22 of them being newly discovered. Their orbital periods range from 2.57 to 19.21 days. These systems consist of low-luminosity, highly stripped subgiant donors and accreting A-type stars. The donor stars, with radii of $2.5-8.9~R_\odot$ and effective temperatures around 4000 K, have typical masses of $M_2 \sim 0.3~M_\odot$, indicating substantial mass loss through Roche-lobe overflow. The presence of ellipsoidal variability and H$\alpha$ emission provides strong evidence for ongoing mass transfer. By fitting the spectral energy distributions, spectra, and light curves, we found that most of the accretors have luminosities lower than expected from the mass-luminosity relation, aligning with the predicted faint phase for mass-gaining stars. Three objects of our sample exhibit pulsations with periods from 18 minutes to 8 hours, providing opportunities for asteroseismic studies. The low mass transfer rates and stability make the sample excellent systems for studying mass accretion, advancing our understanding of the Algol-type binary evolution.

The detection of a $\simeq220$~PeV muon neutrino by the KM3NeT neutrino telescope offers an unprecedented opportunity to probe the Universe at extreme energies. We analyze the origin of this event under three scenarios, viz., a transient point source, a diffuse astrophysical emission, and line-of-sight interaction of ultrahigh-energy cosmic rays (UHECR; $E \gtrsim 0.1$~EeV). Our analysis includes the flux from both a KM3NeT-only fit and a joint fit, incorporating data from KM3NeT, IceCube, and Pierre Auger Observatory. If the neutrino event originates from transients, it requires a new population of transient that is energetic, gamma-ray dark, and more abundant than known ones. In the framework of diffuse astrophysical emission, we compare the required local UHECR energy injection rate at $\gtrsim4$ EeV, assuming a proton primary, with the rate derived from the flux measurements by Auger. This disfavors the KM3NeT-only fit at all redshifts, while the joint fit remains viable for $z\gtrsim 1$, based on redshift evolution models of known source populations. For cosmogenic origin from point sources, our results suggest that the luminosity obtained at redshifts $z \lesssim 1$ from the joint fit is compatible with the Eddington luminosity of supermassive black holes in active galactic nuclei.

Quasars serve as important cosmological probes and constructing accurate luminosity relations for them is essential for their use in cosmology. If the coefficients of quasar's luminosity relation vary with redshift, it could introduce biases into cosmological constraints derived from quasars. In this paper, we conduct a detailed analysis of the redshift variation in the X-ray luminosity and ultraviolet (UV) luminosity ($L_\mathrm{X}$-$L_\mathrm{UV}$) relations of quasars. For the standard $L_\mathrm{X}$-$L_\mathrm{UV}$ relation, we find that the relation coefficients exhibit a strong and linear correlation with redshift, which is not attributable to the selection effect. Additionally, we examine two three-dimensional, redshift-evolving $L_\mathrm{X}$-$L_\mathrm{UV}$ relations and find that the inclusion of a redshift-dependent term does not eliminate the impact of redshift evolution, as the relation coefficients continue to evolve with redshift. Finally, we construct a new $L_\mathrm{X}$-$L_\mathrm{UV}$ relation in which the redshift evolution of the relation coefficients is nearly eliminated. Calibrating the luminosity relations using Hubble parameter measurements, we demonstrate that quasars utilizing our new relation yield effective constraints on cosmological parameters that are consistent with results from Planck CMB data, unlike constraints derived from the standard relation.

The thickness of a galaxy's disk provides a valuable probe of its formation and evolution history. Observations of the Milky Way and local galaxies have revealed an ubiquitous disk structure with two distinctive components: an old thick disk and a relatively young thin disk. The formation of this dual-disk structure and the mechanisms that develop the thickness of the disk are still unclear. Whether the disk thickness inherit from the birth environment or is established through secular dynamical heating after formation is under debate. In this work we identify a relatively young ($\sim$6.6 billion years old) geometric thick disk in the Milky Way, with a scale height of $0.64$ kpc at the Solar Circle. This young thick component exhibits comparable thickness and flaring strength to the canonical old thick disk but is more radially extended and systematically younger. We also identify thin disk components that formed before and after this young thick disk. Detailed analysis of the solar vicinity structure suggests that the young thick disk marks the onset of a new phase of upside-down disk formation. These findings strongly discount the role of secular dynamical heating and support a turbulent, bursty birth environment as the primary mechanism behind thick disk formation. The existence of two thick disk components suggests that the Milky Way has undergone at least two episodes of turbulent and bursty star formation, likely triggered by galaxy mergers.

The Vainu Bappu Telescope (VBT) is a 2.34-m reflector, primarily supported on-axis field of view, offering high-resolution and low-to-medium resolution spectroscopic observations in its prime and Cassegrain configurations. This study presents the design and fabrication of a compact, lightweight, three-element wide-field corrector (WFC) utilizing three spherical lenses to cover a polychromatic wavelength range over a 30$'$ FoV at prime focus. The WFC design was optimized using ZEMAX, ensuring precision in aberrations, tolerances, and atmospheric dispersion. The fabricated lenses met stringent tolerances, with a $\pm$1 mm deviation in radius of curvature and 2 mm deviation in center thickness. A mechanical mount was developed to integrate all the WFC lenses, and wavefront error testing for the WFC system was performed using ZYGO interferometry, yielding a Wavefront Error of 0.05 $\lambda$. Laboratory performance tests were designed and conducted using a dedicated setup with achromatic lenses and 100 $\mu m$ fiber-coupled polychromatic light source showed a deviation of 0.1 pixel on-axis and 0.5 pixel at the extreme off-axis field compared to the ZEMAX design, demonstrating that the optical performance of WFC is with minimal aberrations across the entire FoV. The successful integration of the WFC at the VBT prime focus will increase the FoV, enabling the multi-fiber, multi-spectrograph setup in 30 arcmin field that will facilitate both OMR and Echelle spectrograph to be used on the same night along with the addition of new multi-object spectrograph and an integral field unit instrument. This will mark a significant upgrade for the VBT, broadening its research potential, and expanding its observational versatility.

X-ray light curves of gamma-ray burst (GRB) afterglows exhibit various features, with the shallow decay phase being particularly puzzling. While some studies report absence of the X-ray shallow decay for hyper-energetic GRBs, recently discovered GRB 240529A shows a clear shallow decay phase with an isotropic gamma-ray energy of \SI{2.2e54}{erg}, making it a highly unusual case compared to typical GRBs. In order to investigate the physical mechanism of the shallow decay, we perform the \textit{Fermi}-LAT analysis of GRB 240529A along with \textit{Swift}-XRT analysis. We find no jet break feature in the X-ray light curve and then give the lower bound of the collimation-corrected jet energy of $>10^{52}$~erg, which is close to the maximum rotational energy of a magnetar. Our LAT data analysis reveals evidence of GeV emission with a statistical significance of $4.5\sigma$ during the shallow decay phase, which can be interpreted as the first case for hyper-energetic GRBs with a typical shallow decay phase. The GeV to keV flux ratio is calculated to be $4.2\pm2.3$. Together with X-ray spectral index, this indicates an inverse Compton origin of the GeV emission. Multiwavelength modeling based on time-dependent simulations tested two promising models, the energy injection and wind models. Both models can explain the X-ray and gamma-ray data, while our modeling demonstrates that gamma-ray observations, along with future GeV--TeV observations by CTAO, will distinguish between them.

GJ1061 is a very nearby M star hosting three low-mass temperate planets detected from radial velocity variations. The close to 4:2:1 period commensurability of the planets, the available long-term monitoring of the system and new very high-precision radial velocity measurements from ESPRESSO enable the determination of masses from the planet-planet interaction. Using nested sampling, we derived parameter distributions for a co-planar configuration. The three planets (Mb =1.07 +- 0.11M_Earth, Pb =3.2073 +- 0.0003 d, Mc=1.76 +- 0.13M_Earth, Pc=6.6821 +- 0.0008 d, Md =1.55 +- 0.17M_Earth, Pd =13.066 +- 0.002 d) are potentially all rocky with equilibrium temperatures between 360 K and 240 K. This makes the GJ1061 system one of the prime targets for future ground or space based instruments suitable for a direct detection of the planetary atmospheres.

Dust is expected to form on a year timescale in core-collapse supernova (SN) ejecta. Its existence is revealed through an infrared brightening, an optical dimming, or a blue-red emission-line profile asymmetry. To investigate how the dust location and amount impact observations, we computed ultraviolet-to-optical spectra of interacting and standard, noninteracting Type II SNe using state-of-the-art models -- for simplicity we adopted 0.1micron silicate grains. These models account for the full ejecta and treat both radioactive decay and shock power that arises from interaction of the ejecta with circumstellar material. In a Type IIn SN such as 1998S at one year, approximately 3e-4Msun of dust within the dense shell reproduces the broad, asymmetric Halpha profile. It causes an optical dimming of ~2mag (which obscures any emission from the inner, metal-rich ejecta) but, paradoxically, a more modest dimming of the ultraviolet, which originates from the outer parts of the dense shell. In Type II SNe with late-time interaction, such as SN2017eaw, dust in the low-mass, fast outer ejecta dense shell tends to be optically thin, impacting little the optical spectrum for dust masses of order 1e-4Msun. In such SNe II with interaction, dust in the inner metal-rich ejecta has negligible effect on observed spectra in the ultraviolet and optical. In noninteracting SNe II, dust within the metal-rich ejecta preferentially quenches the [OI]6300,6364 and [CaII]7291,7323 metal lines, biasing the emission in favor of the H-rich material which generates the Halpha and FeII emission below 5500A. Our model with 5e-4Msun of dust below 2000km/s matches closely the optical spectrum of SN1987A at 714d. Modeling historical SNe requires treating both the ejecta material and the dust, as well as multiple power sources, although interaction power will generally dominate.

Spatial data fusion is a bottleneck when it meets the scale of 10 billion records. Cross-matching celestial catalogs is just one example of this. To challenge this, we present a framework that enables efficient cross-matching using Learned Index Structures. Our approach involves a data transformation method to map multi-dimensional data into easily learnable distributions, coupled with a novel search algorithm that leverages the advantages of model pairs, significantly enhancing the efficiency of nearest-neighbor search. In this study, we utilized celestial catalog data derived from astronomical surveys to construct the index and evaluated the speed of the cross-matching process. Using the HEALPix segmentation scheme, we built an independent model object for each tile and developed an end-to-end pipeline to construct a framework with semantic guarantees for record retrieval in query and range search. Our results show that the proposed method improves cross-matching speed by more than four times compared to KD-trees for a radius range between 1 milli-arcseconds and 100 arcseconds.

After more than three decades of investigation, the distribution of Ba stars in the e-log P diagram still defies our understanding. Recent smooth particle hydrodynamic simulations involving an asymptotic giant branch (AGB) primary have shown that a circumbinary disk (CB) can form around the binary and that the presence of dust in the wind of evolved low- and intermediate-mass stars can significantly affect the systemic angular momentum loss and mass accretion onto the companion through the wind Roche lobe overflow (WRLOF) phase. We used the binary evolution code BINSTAR, where we updated the modeling of the progenitors of Ba stars including a CB disk, the WRLOF, tidally enhanced wind mass loss, and non-conservative RLOF with their effects on the orbital evolution. In our approach, we considered that a CB disk forms when WRLOF is activated. The coupling between the CB disk and the binary follows the standard resonant interaction theory. We constructed grids of 2.0 + 1.0 $M_\odot{}$ and 1.2 + 0.8 $M_\odot{}$ binaries for initial orbital parameters that result in WRLOF, and evolved these systems until the end of the primary's AGB phase. WRLOF resulted in a significant shrinkage of the orbital separation during the AGB phase, leading to binaries with initial periods on the order of $\lesssim 12000$ d undergoing Roche lobe overflow (RLOF). The combination of WRLOF, eccentricity pumping from the CB disk, and/or tidally enhanced wind mass loss can lead to RLOF on eccentric orbits down to periods of $P_\mathrm{orb} \sim 3000$d. Non-conservative RLOF enabled a reduction of the period before circularization down to $\sim 2000$d, provided at least 50 percent of the transferred mass left the system. Our models still cannot account for the eccentricity distribution of Ba stars with periods shorter than $P_\mathrm{orb} \lesssim 2000$d, where a common envelope evolution appears unavoidable.

Gravitational tidal interactions drive long-term rotational and orbital evolution in planetary systems, in multiple (particularly close binary) star systems and in planetary moon systems. Dissipation of tidal flows in Earth's oceans is primarily responsible for producing gradual expansion of the Moon's orbit at a few centimetres per year as the Earth's day lengthens by a few milliseconds per century. Similar processes occur in many astrophysical systems. For example, tidal dissipation inside (slowly rotating) stars hosting short-period planets can cause the orbits of these planets to decay, potentially leading to planetary destruction; tidal dissipation inside stars in close stellar binary systems -- and inside short-period planets such as hot Jupiters in planetary systems -- can cause initially eccentric orbits to become circular. To model these processes, explain many current observational results, and make predictions for future observations, we require a detailed theoretical understanding of tidal flows and the mechanisms by which -- and how efficiently -- they are dissipated inside stars and planets. This article will introduce our current understanding of tidal flows and dissipation inside stars (and to a lesser extent giant planets), as well as highlight some unsolved problems.

The Imaging X-ray Polarimetry Explorer is an X-ray observatory measuring the X-ray polarization in the 2-8 keV energy range. Highly sensitive to the system's geometry, X-ray polarization is a unique method to probe the structure of X-ray binaries. The Imaging X-ray Polarimetry Explorer observed the High-Mass X-ray Binary Cygnus X-1 and the Low-Mass X-ray Binary Swift J1727.8-1613 in different accretion states: in the hard state and in the soft state. The X-ray polarimetry analysis of both sources shows a linear polarization degree increasing with energy, with higher values in the hard state than in the soft state. However, the linear polarization angle stays similar in both states and is aligned with the radio jet within $5^\circ$. Furthermore, the Low-Mass X-ray Binary Swift J1727.8-1613 has a lower optical intrinsic polarization and a lower X-ray polarization degree for a softer spectrum. The similarities observed in this analysis between the X-ray polarization results of different types of X-ray Binaries show that the innermost accretion processes are independent of the companion star's type.

Pan and Daphnis are embedded in Saturn's rings and opening a gap with satellite wakes at the gap edges. Furthermore, in the case of Daphnis, pronounced vertical wall structures casting shadows on the rings are also observed in the satellite wakes. In this paper, we perform a global 3D N-body simulation with non-zero $e_{\rm s}$ or non-zero $i_{\rm s}$ of the satellite orbit to investigate how they affect the gap edge structures. We found that the effect of satellite eccentricity is important both in the satellite wakes and the vertical walls at the gap edges. The non-sinusoidal sawtooth-like satellite wakes and azimuthally more localized vertical walls observed by Cassini are simultaneously reproduced in the detailed structures and spatial scales. Both of them periodically vary due to the satellite excursions between the apocenter and the pericenter. The ring particles in outer (inner) rings that undergo closest encounters with the satellite near the apocenter (pericenter) are excited the most highly. Because the excited eccentricities of the ring particles are converted to the inclinations through physical collisions, the conversion is the most active for the particles that acquire the highest eccentricities, resulting in the azimuthally more localized vertical wall structures. The predicted height of the tallest vertical walls is $\sim 0.2$ times the satellite Hill radius in the case of the satellite eccentricity comparable to Daphnis when adopting Hill scaling, which is twice as much as the height obtained in the case of the circular satellite orbit and is quantitatively more consistent with the Cassini observation. These results show that the observed vertical walls are actually formed by the satellite wakes followed by their conversion to the vertical motions through inter-particle collisions, rather than by the out-of-plane perturbation from the satellite in an inclined orbit.

Radio-emitting active galactic nuclei (AGNs) are common in elliptical galaxies and AGN feedback is one of the possible mechanisms for regulating star formation in massive galaxies. It is unclear if all passive galaxy populations host radio AGNs and if AGN feedback is a plausible mechanism for truncating or regulating star formation in these galaxies. To determine if radio AGNs are common in passive spiral galaxies, we have measured the radio emission of 38 low-redshift passive spiral galaxies using RACS-low at 887.5 MHz and VLASS at 3 GHz. We selected a subset of 2MRS galaxies with negligible WISE 12 $\mu$m emission from warm dust, and spiral morphologies from HyperLeda, RC3, 2MRS and manual inspection. In contrast to comparable early-type galaxies, our sample has no significant radio detections, with radio flux densities below 1 mJy, implying that radio AGNs are rare or non-existent in passive spirals. Using the combined radio images and assuming radio luminosity is proportional to $K$-band luminosity, we find ${\rm log}~L_\nu \lesssim 9.01-0.4~M_K$. This falls below the radio luminosities of passive elliptical galaxies, implying radio luminosity in passive galaxies is correlated with host galaxy morphology and kinematics.

In this work, we present a comprehensive and systematic study of the statistical complexity, originally introduced by López-Ruiz, Mancini, and Calbet [Phys. Lett. A 209, 321-326 (1995)], across a broad range of compact star models. We explore how complexity correlates not only with macroscopic observables such as mass and radius, but also with the microscopic characteristics of the underlying equation of state. By incorporating both realistic equations of state and analytical solutions to Einstein's field equations, we demonstrate that gravitational mass plays a dominant role in determining the behavior of complexity. Furthermore, we show that strong phase transitions within the stellar interior, such as those hypothesized in hybrid stars, can manifest as distinct features in the complexity profile, offering a potential informational signature of such transitions. This work offers new insights into the link between information theory and compact object physics, highlighting complexity's potential as a diagnostic tool in astrophysics.

Quasi-periodic eruptions (QPEs) are intense repeating soft X-ray bursts with recurrence times about a few hours to a few weeks from galactic nuclei. More and more analyses show that (at least a fraction of) QPEs are the result of collisions between a stellar mass object (SMO, a stellar mass black hole or a main sequence star) and an accretion disk around a supermassive black hole (SMBH) in galactic nuclei. Previous studies have shown that the SMO trajectory can be reconstructed from QPE timing data, consequently the SMBH mass can be robustly measured from tracing a single SMO. In this Letter, we construct a comprehensive Bayesian framework for implementing the QPE timing method, explore the optimal QPE observation strategy for measuring SMBH masses, and forecast the measurement precision expected in the era of multi-target X-ray telescope, Chasing All Transients Constellation Hunters (CATCH). As a result, we find the QPE timing method is of great potential in precise measurement of SMBH masses (and spins), especially in the lower mass end ($\lesssim 10^7 M_\odot$) where QPEs prevail and relevant dynamical timescales are reasonably short to be measured.

The search for long-term variability of compact components of radio sources B0821+394 and B1812+412 over an interval of 10 years was carried out. The LPA LPI radio telescope with an operating frequency of 111 MHz was used for observations. According to our estimates, the characteristic time of variability for both sources is 1.5-2.5 years. It is shown that the observed variability is not related to intrinsic variations in the radiation flux, but is due to refractive scintillation on inhomogeneities of the interstellar medium. From the obtained upper estimates of the apparent angular dimensions of the sources, it follows that the main contribution to the scattering of radio emission is made by turbulent plasma concentrated in sufficiently thin screens, the distance to which does not exceed 300-400 pc.

The disk around PDS 70 hosts two directly imaged protoplanets in a gap. Previous VLT/SPHERE and recent James Webb Space Telescope/NIRCam observations have hinted at the presence of a third compact source in the same gap at ~13 au, interior to the orbit of PDS 70 b. We reduce seven published and one unpublished VLT/SPHERE datasets in YJH and K bands, as well as an archival VLT/NaCo dataset in L' band, and an archival VLT/SINFONI dataset in H+K band. We combine angular-, spectral- and reference star differential imaging to search for protoplanet candidates. We recover the compact source in all epochs, consistent with the JWST detection, moving on an arc that can be fit by Keplerian motion of a protoplanet which could be in a resonance with PDS 70 b & c. We find that the spectral slope is overall consistent with the unresolved star and inner disk emission at 0.95--1.65$\mu$m, which suggests a dust scattering dominated spectrum. An excess beyond 2.3$\mu$m could be thermal emission from either a protoplanet or heated circumplanetary dust, variability, or inner disk contamination, and requires confirmation. While we currently cannot rule out a moving inner disk feature or a dust clump associated with an unseen planet, the data supports the hypothesis of a third protoplanet in this remarkable system.

Here we use the samples of spiral and elliptical galaxies, in order to investigate theoretically some of their properties and to test the empirical relations, in the light of modified gravities. We show that the baryonic Tully-Fisher relation can be described in the light of $f(R)$ gravity, without introducing the dark matter. Also, it is possible to explain the features of fundamental plane of elliptical galaxies without the dark matter hypothesis.

The evolution of the universe together with the galaxies is one of the fundamental issues that we humans are most interested in. Both the observations of tidal streams from SDSS and the theory of $\Lambda$CDM support the hierarchical merging theory. The study of high redshift celestial bodies contributes to a more in-depth study of cosmology. The LAMOST low resolution search catalog DR11 v1.0 has released 11,939,296 spectra, including 11,581,542 stars, 275,302 galaxies, and 82,452 quasars, and so on. The data of 28,780 stellar population synthesis of galaxies and some high redshift quasars are used to do a preliminary statistical research. We selected the data with small errors for analysis and obtained some basic statistical conclusions. Older galaxies have relatively larger stellar velocity dispersions. The larger the metallicity, the greater the stellar velocity dispersion. These statistical results are reasonable and consistent with previous work. Because the stellar velocity dispersion is driven by the total mass of a galaxy at the first order and more massive galaxies have older ages and greater metallicities. The spectra of high redshift quasars show clear Gunn-Peterson trough and Lyman-$\alpha$ forest. The identified emission lines and high redshift celestial spectra released by LAMOST can be used for cosmological research.

High-energy physics often motivates multi-field inflationary scenarios where stochastic effects play a crucial role. Peculiar to multi-field models, the noise-induced centrifugal force results in a longer duration of inflation depending on the number of fields, even when the stochastic noises themselves are small. We show that, in such small-noise regimes, the number of fields generically discriminates whether inflation successfully terminates or lasts forever. Our results indicate that inflation with an extremely large number of fields may fail to realise our observable Universe.

We present a MeerKAT survey of portions of the Milky Way bulge. The survey covers 172.8 square degrees in two contiguous mosaics above and below the Galactic Center as well as 32 single pointing fields at higher longitudes. The resolution of the images is $\sim$8\asec\ at a frequency of 1333 MHz with a typical Stokes I RMS of 20 $\mu$Jy Beam$^{-1}$. Most of the emission seen is from background extragalactic sources but many compact Galactic objects are identifiable by their polarization properties. Apparent polarized emission resulting from fine scale Faraday rotation in the ISM is widespread in this region of the Galaxy. The survey is used to search for background Giant Radio Galaxies, $>$700 kpc in size, identifying 17 such objects. Data products include FITS images of Stokes I, Q, U and V as well as a Faraday analysis and lists of compact total intensity and polarized sources.

Recent findings on retrograde co-orbital mean-motion resonances in the Earth-Moon system, highlight the potential use of spacecraft in retrograde resonances. Based on these discoveries, this study investigates retrograde co-orbital resonances within the Earth-Moon system, focusing on both optimal and sub-optimal orbital transfers to such configurations. The paper provides a comprehensive analysis of retrograde co-orbital resonances, optimization techniques to evaluate and enhance the performance of bi-impulsive transfers to these configurations. The results reveal the feasibility of low-cost transfers, which could support a range of future missions, including space exploration and satellite deployment. Combining advanced optimization processes, we obtained solutions for orbital transfers for different arrival points in retrograde co-orbitals improving mission efficiency and offering a cost-effective approach to interplanetary exploration.

Hot- and cold-start planet formation models predict differing luminosities for the young, bright planets that direct imaging surveys are most sensitive to. However, precise mass estimates are required to distinguish between these models observationally. The presence of two directly imaged planets, PDS 70 b and c, in the PDS 70 protoplanetary disk provides us a unique opportunity for dynamical mass measurements, since the masses for these planets are currently poorly constrained. Fitting orbital parameters to new astrometry of these planets, taken with VLTI/GRAVITY in the $K$~band, we find $2\sigma$ dynamical upper mass limits of 4.9 $M_{\rm Jup}$ for b and 13.6 $M_{\rm Jup}$ for c. Adding astrometry from the newly proposed planet candidate PDS 70 d into our model, we determine $2\sigma$ dynamical upper mass limits of 5.3 $M_{\rm Jup}$, 7.5 $M_{\rm Jup}$ and 2.2 $M_{\rm Jup}$ for b, c, and the candidate d respectively. However, $N$-body analysis of the orbits fit in this case suggest that the inclusion of $d$ makes the system unstable. Using the upper mass limits for b and c we rule out the coldest-start formation models for both planets, calculating minimum post-formation entropies ($S_i$) of 9.5 $k_{\rm B}/{\rm baryon}$ and 8.4 $k_{\rm B}/{\rm baryon}$ respectively. This places PDS 70 b and c on the growing list of directly-imaged planets inconsistent with cold-start formation.

We observed three recurrent blowout jets in an active regio with Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). Using Helioseismic Magnetic Imager (HMI) data. We found that the magnetic flux of an emerging negative pole increases steadily before declining just as the jets erupt. Certain physical quantities, like the total unsigned vertical current, align with the periodicity of the jets. The differential affine velocity of the vector magnetograms reveals strong shear around the negative pole. The Doppler velocity map, calculated from the H$\alpha$ spectra observed by the Chinese H$\alpha$ Solar Explorer (CHASE), shows upflows with large initial velocity before it can be observed by AIA. The magnetic field derived from the nonlinear force-free field (NLFFF) model suggests a topology akin to fan-spine structure, consistent with AIA images. We calculated the evolution of volumetric helicity ratio using the NLFFF model and found its phase aligns with the jet flux in AIA 171 Å. These results suggest that recurrent jets may be triggered by the accumulation and release of energy and helicity, driven by emergence, shearing and cancellation of photospheric magnetic field.

Over the past decade, there has been a significant increase in the reporting of extensive and luminous star-forming regions associated with explosive outflows. Nevertheless, there is still a lack of understanding of the possible physical mechanisms that produce such energetic and isotropic events. Then, we propose a gravitational interaction as a likely mechanism that could trigger explosive outflows in dense star-forming regions. This could constrain the physical conditions that generate an explosive outflow produced by the close dynamical encounter of a runaway star with a clump cluster in dynamical equilibrium. Then, we have produced a set of $N$-body simulations that account for the collision of a 10M$_\odot$ stellar object with a cluster of particles with a mass that ranges from 0 to 50 \msun. We propose a parameter to describe the interaction, the evaporation parameter, that represents the fraction of stars that become unbound. The main result is that, when the cluster mass is less than, or up to a few times the stellar mass, the collision will produce an explosive outflow, ejecting a significant fraction of the cluster members with velocities larger than the impact velocity. All of our models produce an explosive outflow, with different characteristics, which increases the probability that a close encounter could be responsible for producing the observed flows.

We present a comprehensive analysis of 2RXP J130159.6-635806, a persistent low-luminosity Be/X-ray pulsar, focusing on its transition to a spin equilibrium state and the discovery of a bimodal luminosity distribution revealing possibly a new accretion regime. Using data from NuSTAR, Swift, XMM-Newton, and Chandra observatories, we investigate changes in the pulsar's timing and spectral properties. After more than 20 years of continuous spin-up, the pulsar's spin period stabilized, marking the onset of spin equilibrium. This transition was accompanied by the emergence of a previously unobserved accretion regime at $L_{\rm bol} = (2.0_{-1.0}^{+2.3})\times 10^{34}$ erg s$^{-1}$, an order of magnitude lower than its earlier quiescent state. After that, the source occasionally switched between these regimes, remaining in each state for extended periods, with the transition time from a luminosity of $10^{35}$ erg s$^{-1}$ to $10^{34}$ erg s$^{-1}$ taking less than 2.3 day. The analysis of the spectral data collected during this new low-luminosity state revealed a two-hump shape which is different from the cutoff power-law spectra observed at higher luminosities. The discovery of pulsations in this state, together with the hard spectral shape, demonstrates ongoing accretion. We estimate the magnetic field strength to be $\sim 10^{13}$ G based on indirect methods. Additionally, we report a hint of a previously undetected $\sim$90-day orbital period in the system.

Asteroseismic modelling is crucial for upcoming missions like PLATO, CubeSpec, and Roman. Despite significant progress, discrepancies between observations and theoretical predictions introduce biases in stellar characterisation at the precision required by PLATO. Current models typically ignore magnetic activity, assuming its effects are hidden within surface effects. However, recent studies have shown significant impacts of magnetic activity on the Sun's asteroseismic characterisation using forward modelling. Using GOLF and BiSON observations of two full solar activity cycles, we quantified the impact of magnetic activity on solar mean density and acoustic radius inversions. Observations were segmented into yearly overlapping snapshots, each offset by 91.25 days. Inversions were performed for each snapshot to determine mean density and acoustic radius, tracking their temporal evolution and estimating systematic uncertainty due to magnetic activity. We observed a clear imprint of the magnetic activity cycle on solar mean density and acoustic radius through helioseismic inversions, consistent across GOLF and BiSON datasets. This imprint is the largest source of systematic uncertainty in solar asteroseismic characterisation. Including low radial-order modes mitigates these effects more significantly than previously measured for other stellar variables. We recommend asteroseismic values for solar mean density (1.4104 \pm 0.0051 g/cm3) and acoustic radius (3722.0 \pm 4.1 s), averaged over two activity cycles. These values account for major systematic errors, achieving high precision (0.36% for mean density and 0.11% for acoustic radius). These results are promising for high-precision characterisation of Sun-like stars, a better-constrained mean density being able to enhance the precision of stellar radius estimate, which is crucial for exoplanetary system characterisation.

We explore the cosmological implications of generalized entropic models within the framework of Gravity-Thermodynamics (GT) approaches. These models, characterized by three or four additional free parameters, are designed to capture deviations from the standard Bekenstein-Hawking entropy and can reproduce well-known entropic formulations, including Tsallis, Rényi, Sharma-Mittal, Barrow, Kaniadakis, and Loop Quantum Gravity entropies in various analytical limits. We implement the corresponding cosmological models using a fully numerical GT approach to constrain the model parameters and to study the evolution of the dark energy equation of state as a function of the scale factor. Our Bayesian analysis, which incorporates the Pantheon+ and DESy5 supernovae data alongside the recently released DESI-DR2/DR1 Baryon Acoustic Oscillation (BAO) measurements, shows that the data favor the standard Bekenstein-Hawking entropy, leading to a $\Lambda$CDM-like late-time behavior. In this context, the three-parameter ($\mathcal{S}_3$) entropic model appears to be sufficient to capture the observed dark energy phenomenology. Furthermore, a direct comparison of the Bayesian evidence indicates that the three-parameter model is preferred over the four-parameter ($\mathcal{S}_4$) variant by a factor of $\Delta\log\mathcal{B} \sim -6$, while the GT approach as a whole is significantly disfavored relative to the $\Lambda$CDM model with at least $\Delta\log\mathcal{B} \sim -8$ ($\mathcal{S}_3$) to $\Delta\log\mathcal{B} \sim -13$ ($\mathcal{S}_4$), when using the DESy5 and DESI-DR2 datasets.

We investigated the impact of a spatially varying matter potential $\lambda$, coming from neutrino-electron forward scattering, on the emergence of fast neutrino flavor conversion (FFC) triggered by the presence of zero crossings in the angular distribution of the neutrino electron lepton number (ELN). We find that FFC can be significantly affected as the spatial variation rate of $\lambda$ increases, and strong spatial variations can completely stabilize initially unstable systems. Using stability analysis based solely on initial conditions, we identified for the first time a critical variation rate above which no FFC occurs even if the flavor instability exists. By analyzing several representative matter profiles based on an 18 $M_{\odot}$ SN model, we show that spatially inhomogeneous $\lambda$ can suppress the occurrence of FFC associated with shallow ELN zero crossings in most of the SN's radial region, especially during the accretion phase. Our finding highlights the need to consider the impact of matter inhomogeneity in the development of improved SN models that aim to include the effect of neutrino flavor conversions.

The detectors of the JWST Mid-Infrared Instrument (MIRI) Medium Resolution Spectrometer (MRS) form low-finesse resonating cavities that cause periodic count rate modulations (fringes) with peak amplitudes of up to 15% for sources external to MIRI. To detect weak features on a strong continuum and reliably measure line fluxes and line-flux ratios, fringe correction is crucial. This paper describes the first of two steps implemented in the JWST Science Calibration Pipeline, which is the division by a static fringe flat that removes the bulk of the fringes for extended sources. Fringe flats were derived by fitting a numerical model to observations of spatially extended sources. The model includes fringes that originate from two resonating cavities in the detector substrate (a third fringe component that originates from the dichroic filters is not included). The model, numerical implementation, and resulting fringe flats are described, and the efficiency of the calibration was evaluated for sources of various spatial extents on the detector. Flight fringe flats are obtained from observations of the planetary nebula NGC 7027. The two fringe components are well recovered and fitted by the model. The derived parameters are used to build a fringe flat for each MRS spectral band, except for 1A and 1B due to the low signal-to-noise ratio of NGC 7027 in these bands. When applied to extended sources, fringe amplitudes are reduced to the sub-percent level on individual spaxels. For point sources, they are reduced to amplitudes of 1 to 5% considering individual spaxels and a single dither position, and decrease to the 1 to 2% level after two-dimensional residual fringe correction. The fringe flats derived from this work are the reference files currently in use by the JWST Science Calibration Pipeline. They provide an efficient calibration for extended sources, and are less efficient for point sources.

I examine the morphology of the core-collapse supernova (CCSN) remnant (SNR) G0.9+0.1 and reveal a point-symmetrical morphology that implies shaping by three or more pairs of jets, as expected in the jittering jets explosion mechanism (JJEM). The large northwest protrusion, the ear (or lobe), has two bright rims. I compare this ear with its rims to an ear with three rims of a jet-shaped planetary nebula and jets from an active galactic nucleus that shaped several rims on one side. Based on this similarity, I argue that two jets or more shaped the northwest ear of SNR G0.9+0.1 and its two rims. I identified the bright region south of the main shell of SNR G0.9+0.1 as a jet-shaped blowout formed by a jet that broke out from the main SNR shell. I base this on the similarity of the blowout of SNR G0.9+0.1 with that of SNR G309.2-00.6, argued in the past to be shaped by jets. I identify four symmetry axes along different directions that compose the point-symmetric morphology of SNR G0.9+0.1. I show that the morphological features of holes, granular texture, and random filaments exist in CCSNe and planetary nebulae and are unlikely to result from some unique processes in CCSNe. These structures result from similar instabilities in the JJEM and the neutrino-driven explosion mechanism and, unlike a point-symmetric morphology, cannot determine the explosion mechanism. Identifying SNR G0.9+0.1 as a new point-symmetric CCSN strengthens the JJEM as the primary explosion mechanism of CCSNe.

The isotopic composition of meteorites linked to S-complex asteroids has been used to suggest that these asteroids originated in the terrestrial planet's region, i.e., within 1.5 au, and later got implanted into the main asteroid belt (MAB). Dynamical models of planet formation support this view. Yet, it remains to be demonstrated whether the currently observed size-frequency distribution (SFD) of S-complex bodies in the MAB can be reproduced via this implantation process. Here we studied the evolution of the SFD of planetesimals during the accretion of terrestrial planets with the code LIPAD self-consistently accounting for growth and fragmentation of planetesimals. In our simulations we vary the initial surface density of planetesimals, the gaseous disk lifetime, and the power slope of the initial planetesimals' SFD. We compared the final SFDs of leftover planetesimals in the terrestrial planet region with the SFD of observed S-complex MAB objects (D $>$ 100km). We found that the SFDs of our planetesimal populations and that of S-complex MAB objects show very similar cumulative power index (i.e., q $\approx$ 3.15 in N($>$D)$~\propto$ D$^{-q}$) for slopes in the diameter range 100 km $<$ D $<$ 400 km by the end of our simulations. Our results support the hypothesis of S-complex MAB implantation from the terrestrial planet forming region, assuming implantation is size-independent, and implies that implantation efficiency is smaller than $\mathcal{O}$(10$^{\rm -2}$--10$^{\rm -4}$) to avoid over-implantation of (4) Vesta-sized objects or larger.

We explore SN 2023pel, the most recent event associated with gamma-ray bursts (GRBs), specifically GRB 230812B. SN 2023pel has a high luminosity and low expansion velocities compared to other GRB-SNe. These properties seem difficult to reconcile with a single nickel power source. We searched for models that can explain the properties of this event. We calculated a grid of hydrodynamic models based on pre-SN structures derived from evolutionary calculations. We compared our models with observations of SN~2023pel and selected our preferred model using statistical analysis, taking both light curves and expansion velocities into account. This allowed us to derive a set of physical properties for SN~2023pel. Our models suggest that the most probable scenario involves a millisecond magnetar as the primary power source, supplemented by energy from radioactive decay. Our preferred model has a spin period of P = 3.2 ms, a magnetic field of B = 28 x 10^14 G, an explosion energy of 2.3 foe, a nickel mass of M_Ni= 0.24 solar masses, and an ejected mass of 3.4 solar masses. Alternatively, we find that a purely nickel-powered model also provides a good match with the observations, though M_Ni > 0.8 solar masses are always required. However, the combination of such high values of M_Ni and low M_ej is difficult to reconcile, indicating that this scenario is less probable. We have also identified a specific region within the peak luminosity-velocity plane where an additional energy source beyond nickel may be necessary to power SNe with characteristics similar to SN~2023pel. Our study indicates that an additional energy source beyond radioactive decay is essential to explain the high brightness and relatively low expansion velocities of SN 2023pel. A magnetar-powered model, similar to the models proposed for the very luminous GRB-SN 2011kl, aligns well with these characteristics.

2MASS J04381486+2611399 (or J0438) is one of the few young brown dwarfs (BD) with a highly inclined ($i\!\sim\!70^\circ$) disk. Here we report results from JWST-MIRI MRS, HST-ACS and ALMA Band 7 observations. Despite its late spectral type (M7.25), the spectrum of J0438 resembles those of inner disks around earlier-type stars (K1-M5, T Tauri stars), with a volatile reservoir lacking hydrocarbons (except for acetylene, C$_2$H$_2$) and dominated by water. Other identified species are H$_2$, CO$_2$, HCN, [Ar$^{+}$], and [Ne$^{+}$]. The dominance of water over hydrocarbons is driven by multiple factors such as disk dynamics, young disk age, low accretion rate and possible inner disk clearing. J0438 appears highly dynamic, showing a seesaw-like variability and extended emission in H$_2 \,\,\, S$(1), $S$(3), $S$(5), [Ne$^{+}$] and CO ($J=3-2$). Interestingly, the CO emission reaches up to 400 au from the brown dwarf, suggesting ongoing infalling/outflowing activity impacting the disk chemistry. These observations underscore the combined power of JWST, HST and ALMA in characterizing the chemical diversity and dynamics of brown dwarf disks.

Infrared observations of the inner disks around very low-mass stars (VLMS, $<$0.3$\,M_{\odot}$) have revealed a carbon-rich gas composition in the terrestrial planet-forming regions. Contrary to the typically water-rich T Tauri disk spectra, only two disks around VLMS have been observed to be water-rich among more than ten VLMS disks observed so far with JWST/MIRI. In this letter, we systematically search for the presence of water and other oxygen-bearing molecules in the JWST/MIRI spectra of ten VLMS disks from the MIRI mid-INfrared Disk Survey (MINDS). In addition to the two previously reported detections of water emission in this VLMS sample, we detect water emission in the spectra of three other sources and tentatively in one source, and we provide strong evidence for water emission in the remaining disks in the MINDS sample, most of which have bright emission from carbon-bearing molecules. We show that the $\rm C_2H_2$ emission is much stronger than that of water for sources with low luminosities, and the hydrocarbons outshine the water emission in such conditions. We propose that the appearance of water-rich vs. hydrocarbon-rich spectra is related to the location of the water reservoir in the disk relative to the main hydrocarbon reservoir. Our findings indicate that the terrestrial planet forming regions in VLMS disks have high carbon-to-oxygen ratios (C/O$>$1), but can still harbor ample water similar to those in the T Tauri disks.

We present a new framework for modeling gravitational wave diffraction near fold caustics using the Uniform Approximation (UA), focusing on binary mass lenses-axially asymmetric systems with complex caustic structures. Full-wave methods based on the Kirchhoff integral become impractical in this regime due to highly oscillatory integrands. The UA provides a robust and accurate description of the wave field near folds, resolving the breakdown of Geometrical Optics at caustics and improving upon Transitional Asymptotics-based on Airy function approximations-which lack global validity. Central to our approach is the concept of the caustic width, $d_c$, a characteristic length scale defining the region where diffraction significantly alters wave propagation. We find that $d_c$ scales universally with the gravitational wavelength as ~ $ \lambda^{2/3}$ and inversely with the redshifted lens mass as ~ $ M_{Lz}^{-2/3}$. The wave amplification near the fold grows as ~ $ d_c^{-1/4}$, substantially enhancing the signal and potentially playing a key role in the detection of gravitational waves lensed near caustics. Notably, for lens masses below the galactic scale, the caustic width for gravitational waves is not negligible compared to the Einstein radius-as it is in electromagnetic lensing-making the UA essential for accurately capturing wave effects.

The scattering of extremely energetic cosmic rays with both cosmic microwave background and extragalactic background light, can produce $\mathcal{O}(10^{18} \,{\rm eV})$ neutrinos, known as cosmogenic neutrinos. These neutrinos are the only messengers from the extreme cosmic accelerators that can reveal the origin of the most energetic cosmic rays. Consequently, much effort is being devoted to achieving their detection. In particular, the GRAND project aims to observe the $\nu_\tau$ and $\bar \nu_\tau$ components of the cosmogenic neutrino flux in the near future using radio antennas. In this work, we investigate how the detection of cosmogenic neutrinos by GRAND can be used to probe beyond the Standard Model physics. We identify three well-motivated scenarios which induce distinct features in the cosmogenic neutrino spectrum at Earth: neutrino self-interactions mediated by a light scalar ($\nu$SI), pseudo-Dirac neutrinos (PD$\nu$) and neutrinos scattering on ultra-light Dark Matter ($\nu$DM). We show these scenarios can be tested by GRAND, using 10 years of cosmogenic neutrino data, in a region of parameter space complementary to current experiments. For the $\nu$SI model,, we find that GRAND can constrain the coupling to $\nu_\tau$ in the range [$10^{-2},10^{-1}$] for a scalar with mass in the range 0.1 to 1 GeV. For PD$\nu$, we find that GRAND is sensitive to sterile-active mass squared splitting in the range [$10^{-15},10^{-13}$] ${\rm eV}^2$. Finally, for the $\nu$DM model, assuming a heavy mediator, GRAND can do substantially better than the current limits from other available data. These results rely on the fact that the actual cosmogenic flux is around the corner, not far from the current IceCube limit.

Precise estimation of the Tully-Fisher relation is compromised by statistical biases and uncertain inclination corrections. To account for selection effects (Malmquist bias) while avoiding individual inclination corrections, I introduce a Bayesian method based on likelihood functions that incorporate Sine-distributed scatter of rotation velocities, Gaussian scatter from intrinsic dispersion and measurement error, and the observational selection function. However, tests of unidirectional models on simulated datasets reveal an additional bias arising from neglect of the Gaussian scatter in the independent variable. This additional bias is identified as a generalized Eddington bias, which distorts the data distribution independently of Malmuqist bias. I introduce two extensions to the Bayesian method that successfully mitigate the Eddington bias: (1) analytical bias corrections of the dependent variable prior to likelihood computation, and (2) a bidirectional dual-scatter model that includes the Gaussian scatter of the independent variable in the likelihood function. By rigorously accounting for Malmquist and Eddington biases in a latent-inclination regression analysis, this work establishes a framework for unbiased distance estimates from standardizable candles, critical for improving determinations of the Hubble constant.

The freeze-in mechanism describes the out-of-equilibrium production of dark matter (DM) particles via feeble couplings or non-renormalisable interactions with large suppression scales. In the latter case, predictions suffer from a strong sensitivity to the initial conditions of the universe, such as the details of reheating. In this work, we investigate how this sensitivity is altered in the presence of a cosmological first-order phase transition. We show that freeze-in via non-renormalisable interactions is not always dominated by the highest temperatures of the Standard Model (SM) thermal bath, but instead may be governed by the period immediately after the phase transition, during which the decaying scalar field transfers its energy density to the SM radiation. We refer to this alternative production regime as DM $\textit{phase-in}$. Using numerical and approximate analytical solutions of the relevant Boltzmann equations, we determine the conditions that under which phase-in or conventional freeze-in production dominates the final DM abundance in terms of the type of interaction between the DM and SM particles, the amount of supercooling before and the evolution of the scalar field after the phase transition. In the phase-in regime, the DM abundance is correlated with the peak frequency of the gravitational wave signal associated with the phase transition, opening up new observational possibilities.

This report, summarising work achieved in the context of the LHC Dark Matter Working Group, investigates the phenomenology of $t$-channel dark matter models, spanning minimal setups with a single dark matter candidate and mediator to more complex constructions closer to UV-complete models. For each considered class of models, we examine collider, cosmological and astrophysical implications. In addition, we explore scenarios with either promptly decaying or long-lived particles, as well as featuring diverse dark matter production mechanisms in the early universe. By providing a unified analysis framework, numerical tools and guidelines, this work aims to support future experimental and theoretical efforts in exploring $t$-channel dark matter models at colliders and in cosmology.

Proposed half a century ago, the quantum chromodynamics (QCD) axion explains the lack of charge and parity violation in the strong interactions and is a compelling candidate for cold dark matter. The last decade has seen the rapid improvement in the sensitivity and range of axion experiments, as well as developments in theory regarding consequences of axion dark matter. We review here the astrophysical searches and theoretical progress regarding the QCD axion. We then give a historical overview of axion searches, review the current status and future prospects of dark matter axion searches, and then discuss proposed dark matter axion techniques currently in development.

Below about 2.3 $\mu$m, the nighttime emission of the Earth's atmosphere is dominated by non-thermal radiation from the mesosphere and thermosphere. As this airglow can even outshine scattered moonlight in the near-infrared regime, the understanding of the Earth's night-sky brightness requires good knowledge of the complex airglow emission spectrum and its variability. As airglow modelling is very challenging, the comprehensive characterisation of airglow emission requires large data sets of empirical data. For fixed locations, this can be best achieved by archived spectra of large astronomical telescopes with a wide wavelength coverage, high spectral resolving power, and good temporal sampling. Using 10 years of data from the X-shooter echelle spectrograph in the wavelength range from 0.3 to 2.5 $\mu$m and additional data from the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Cerro Paranal in Chile, we have succeeded to build a comprehensive spectroscopic airglow model for this low-latitude site under consideration of theoretical data from the HITRAN database for molecules and from different sources for atoms. The Paranal Airglow Line And Continuum Emission (PALACE) model comprises 9 chemical species, 26,541 emission lines, and 3 unresolved continuum components. Moreover, there are climatologies of relative intensity, solar cycle effect, and residual variability with respect to local time and day of year for 23 variability classes. Spectra can be calculated with a stand-alone code for different conditions, also including optional atmospheric absorption and scattering. In comparison to the observed X-shooter spectra, PALACE shows convincing agreement and is significantly better than the previous, widely used airglow model for Cerro Paranal.

We briefly review the recent developments in magnetohydrodynamics, which in particular deal with the evolution of magnetic fields in turbulent plasmas. We especially emphasize (i) the necessity of renormalizing equations of motion in turbulence where velocity and magnetic fields become Hölder singular; (ii) the breakdown of Laplacian determinism (spontaneous stochasticity) for turbulent magnetic fields; and (iii) the possibility of eliminating the notion of magnetic field lines, using instead magnetic path lines as trajectories of Alfvenic wave-packets. These methodologies are then exemplified with their application to the problem of magnetic reconnection -- rapid change in magnetic field pattern that accelerates plasma -- a ubiquitous phenomenon in astrophysics and laboratory plasmas. The necessity of smoothing out rough velocity and magnetic fields on a finite scale L implies that magnetohydrodynamic equations should be regarded as effective field theories with running parameters depending upon the scale L.

We derive the analytical expression of the coordinate time $t$ in terms of the eccentric anomaly $u$ at the second post-Newtonian order in General Relativity for a compact binary system moving on eccentric orbits. The parametrization of $t$ with $u$ permits to reduce at the minimum the presence of discontinuous trigonometric functions. This is helpful as they must be properly connected via accumulation functions to finally have a smooth coordinate time $t(u)$. Another difficulty relies on the presence of an infinite sum, about which we derive a compact form. This effort reveals to be extremely useful for application purposes. Indeed, we need to truncate the aforementioned sum to a certain finite threshold, which strongly depends on the selected parameter values and the accuracy error we would like to achieve. Thanks to our work, this analysis can be easily carried out.

We consider the $f\left( Q \right) $-theory for the description of dark energy with a non-trivial connection defined in the non-coincident gauge. The resulting field equations form a two-scalar-field, quintom-like gravitational model. For the power-law model $f\left( Q \right) \simeq Q^{\frac{n}{n-1}}$, we construct an analytic expression for the dynamical evolution of dark energy, which depends on the parameter $n$. We constrain this dark energy model using the the baryon acoustic oscillations from DESI DR2, and gamma-ray bursts. The cosmological data provides $n \simeq0.33 $. The $f\left( Q \right) $-model challenges the $\Lambda$CDM by providing a smaller value for $\chi_{\min}^{2}$, while there is weak evidence in favor of the $f\left( Q \right) $-theory.

The solar wind is a medium characterized by strong turbulence and significant field fluctuations on various scales. Recent observations have revealed that magnetic turbulence exhibits a self-similar behavior. Similarly, high-resolution measurements of the proton density have shown comparable characteristics, prompting several studies into the multifractal properties of these density fluctuations. In this work, we show that low-resolution observations of the solar wind proton density over time, recorded by various spacecraft at Lagrange point L1, also exhibit non-linear and multifractal structures. The novelty of our study lies in the fact that this is the first systematic analysis of solar wind proton density using low-resolution (hourly) data collected by multiple spacecraft at the L1 Lagrange point over a span of 17 years. Furthermore, we interpret our results within the framework of non-extensive statistical mechanics, which appears to be consistent with the observed nonlinear behavior. Based on the data, we successfully validate the q-triplet predicted by non-extensive statistical theory. To the best of our knowledge, this represents the most rigorous and systematic validation to date of the q-triplet in the solar wind.

A large primordial density perturbation of the Hubble scale will gravitationally collapse, generating an outgoing sound shell whether or not a primordial black hole (PBH) is formed. In this Letter, we report a new source of the stochastic gravitational wave background induced by the collision of sound shells in the early Universe. The peak frequency and amplitude in the GW spectrum depend on the Hubble horizon and the abundance of sound shells. Abundant density perturbations would lead to GW backgrounds potentially detectable for future pulsar timing arrays and ground-based/space-borne detectors. For those perturbations that collapse into PBHs, future null detection of the corresponding high-frequency GW background could put new observational constraints on those PBHs that have already evaporated.

The Laser Interferometer Space Antenna (LISA) will detect gravitational waves from the population of merging massive black holes binaries (MBHBs) throughout the Universe. The LISA data stream will feature many superposed signals from different astrophysical sources, requiring a global fit procedure. Most of the MBHB signals will be loud enough to be detected days or even weeks before the merger; and for those sources LISA will be able to predict the time of the merger well in advance of the coalescence, as well as an approximate position in the sky. In this paper, we present a fast detection and signal reconstruction scheme for massive black hole binaries in the LISA observation band. We propose: (i) a detection scheme for MBHB mergers allowing a first subtraction of these signals for the purpose of a global fit, and (ii) an efficient early detection scheme providing a time-of-merger estimate for a pre-merger signal, that will allow to trigger a protection period, placing LISA in ``do not disturb'' mode and enabling more detailed analysis that will facilitate multi-messenger observations. We highlight the effect of confusion of several overlapping in time MBHB signals in the pre-merger detection.

Domain walls are the simplest type of topological defects formed at cosmological phase transitions, and one of the most constrained. Their studies typically assume a quartic double well potential, but this model is not fully representative of the range of known or plausible particle physics models. Here we study the cosmological evolution of domain walls in two other classes of potentials. The Sine-Gordon potential allows several types of walls, interpolating between different pairs of minima (which demands specific numerical algorithms to separately measure the relevant properties of each type). The Christ-Lee potential parametrically interpolates between sextic and quartic behavior. We use multiple sets of simulations in two and three spatial dimensions, for various cosmological epochs and under various choices of initial conditions, to discuss the scaling properties of these networks. In the Sine-Gordon case, we identify and quantify deviations from the usual scaling behavior. In the Christ-Lee case, we discuss conditions under which walls form (or not), and quantify how these outcomes depend on parameters such as the energy difference between the false and true vacua and the expansion rate of the Universe. Various biased initial conditions are also addressed in appendices. Finally, we briefly comment on the possible cosmological implications of our results.