Papers for Tuesday, Jun 04 2024

The Zwicky Transient Facility (ZTF) was expected to detect more than one strong gravitationally-lensed supernova (glSN) per year, but only one event was identified in the first four years of the survey. This work investigates selection biases in the search strategy that could explain the discrepancy and revise discovery predictions. We present simulations of realistic lightcurves for lensed thermonuclear (glSNIa) and core-collapse supernova (glCCSN) explosions over a span of 5.33 years of the survey, utilizing the actual observation logs of ZTF. We find that the magnitude limit in spectroscopic screening significantly biases the selection towards highly magnified glSNe, for which the detection rates are consistent with the identification of a single object by ZTF. To reach the higher predicted rate of detections requires an optimization of the identification criteria for fainter objects. We find that around 1.36 (3.08) Type Ia SNe (CCSNe) are identifiable with the magnification method per year in ZTF, but when applying the magnitude cut of m < 19 mag, the detection rates decrease to 0.17 (0.32) per year. We compare our simulations with the previously found lensed Type Ia SNe, iPTF16geu and SN Zwicky, and conclude that considering the bias towards highly magnified events, the findings are within expectations in terms of detection rates and lensing properties of the systems. In addition, we provide a set of selection cuts based on simple observables to distinguish glSNe from regular, unlensed, supernovae to select potential candidates for spectroscopic and high-spatial resolution follow-up campaigns. We find optimal cuts in observed colours $g-r$, $g-i$, and $r-i$ as well as in the colour SALT2 fit parameter. The developed pipeline and the simulated lightcurves employed in this analysis can be found in the $LENSIT$ github repository.

With a particular focus on Scipy's minimize function the eclipse mapping method is thoroughly researched and implemented utilizing Python and essential libraries. Many optimization techniques are used, including Sequential Least Squares Programming (SLSQP), Nelder-Mead, and Conjugate Gradient (CG). However, for the purpose of examining photometric light curves these methods seek to solve the maximum entropy equation under a chi-squared constraint. Therefore, these techniques are first evaluated on two-dimensional Gaussian data without a chi-squared restriction, and then they are used to map the accretion disc and uncover the Gaussian structure of the Cataclysmic Variable KIC 201325107. Critical analysis is performed on the code structure to find possible faults and design problems. Additionally, the analysis shows how several factors impacting computing time and image quality are included including the variance in Gaussian weighting, disc image resolution, number of data points in the light curve, and degree of constraint.

Cyano-substituted polycyclic aromatic hydrocarbons (CN-PAHs) may contribute to the emission detected in the 7 - 9 $\mu$m (1430 - 1100 cm$^{-1}$) and 11 - 15 $\mu$m (900 - 670 cm$^{-1}$) regions of astronomical IR spectra. Anharmonic quantum chemical computations of 17 CN-PAH isomers for 4 small PAHs and Benzene reveal strong, broad absorption features across the entire 300 - 6200 cm$^{-1}$ (33 - 1.6 $\mu$m) frequency range. In particular, when a FWHM of 15 cm$^{-1}$ is applied, the composite CN-PAH spectrum is almost indistinguishable from the unsubstituted-PAH spectrum. At high resolution, however, the infrared absorption spectra reveal unique, identifiable features of CN-PAHs in the 700 - 950, 1100 - 1300, 2000 - 2500, and 3400 - 3600 cm$^{-1}$ ranges. The in-plane and out-of-plane CH bending vibrational frequencies of CN-PAHs are shifted when comparing isomers and to their unsubstituted counterparts, making their differentiation in mixed laboratory experiments possible. The overall aromatic CH stretch fundamental (2950 - 3200 cm$^{-1}$) and first overtone (5950 - 6200 cm$^{-1}$) regions are relatively unaffected by the cyano-substitution, with changes only to the breadth and intensity of the bands. Detailed spectroscopic data on the normal mode components of each state reported herein provide the means to directly assign future laboratory spectra and to guide direct IR observations of astronomical regions with, e.g., JWST.

Supernova (SN) 1987A is widely regarded as an excellent candidate for leveraging the capabilities of the freshly launched XRISM satellite. Recent researches indicate that the X-ray emission from SN 1987A will increasingly originate from its ejecta in the years to come. In a previous study, we thoroughly examined the proficiency of XRISM-Resolve in identifying signatures of shocked ejecta in SN 1987A, synthesizing the XRISM-Resolve spectrum based on a state-of-the-art magneto-hydrodynamic simulation. However, following the satellite's launch, a technical issue arose with the XRISM instrument's gate valve, which failed to open, thereby affecting observations with the Resolve spectrometer. Here, we update our analysis, reevaluating our diagnostic approach under the assumption that the gate valve remains closed. We find that, even with the reduced instrumental capabilities, it will be possible to pinpoint the ejecta contribution through the study of the line profiles in the XRISM-Resolve spectrum of SN 1987A.

We present an efficient technique for calculating the angular two-point correlation functions (or ''overlap reduction functions'') induced by gravitational waves in both the pulse arrival times of pulsars and in the angular deflections of distant sources. In the most general case, there are six auto- and cross-correlations for the pulse arrival times and the two components of the angular deflection. We provide results for spin-2 (i.e., general-relativistic) gravitational waves as well as the spin-1 modes that may arise in alternative-gravity theories. These calculations can be easily implemented for future analysis or study, and we provide code to do so.

Eclipsing binary systems are significant objects for astrophysics in that direct observations can determine the fundamental parameters of stars. In this study, we determined precisely the fundamental parameters of the binary component stars obtained by simultaneous analysis of radial velocities and the {\it TESS} light curve using the Wilson and Devinney code. Following the analysis, the masses and radii of the primary and secondary components were determined as $M_{1}= 1.58\pm 0.01M_\odot$, $M_{2}= 0.48\pm0.02M_\odot$, and $R_{1}=1.93\pm 0.01R_\odot$, $R_{2}= 1.14\pm 0.01 R_\odot$, respectively. Furthermore, the distance of IS CMa is calculated as $92.7\pm6.5$ pc. On the basis of the analysis of the mid-eclipse times, it was found that the variation in the orbital period is represented by an upward parabola. It has an increasing rate of $dP/dt$ = 1.09 $\times$ 10$^{-7}$ day yr$^{-1}$. Using PARSEC stellar evolutionary tracks and isochrones with solar metallicity were estimated the age of IS CMa as $1.3\pm0.1$ Gyr. Kinematic and Galactic orbital parameters of IS CMa were obtained from the astrometric and spectroscopic data of the system. The Galactic orbit analysis reveals that IS CMa formed inside the solar circle and it is a member of the young thin-disc population.

A long linear structure recently discovered could be the stellar wake produced by the passage of a runaway supermassive black hole (SMBH) or, alternatively, a bulgeless edge-on galaxy. We report on new very deep HST imaging that seems to be in tension with the SMBH runaway scenario but is consistent with the bulgeless edge-on galaxy scenario. The new observations were aimed at detecting two key features expected in the SMBH scenario, namely, the bow shock formed where the SMBH meets the surrounding medium, and a counter stellar wake created by another binary SMBH hypothesized as part of the ejection mechanism. Neither of these two features appears to be present in the new images, as would be expected in the edge-on galaxy scenario.

We study the bispectrum in Lagrangian perturbation theory. Extending past results for the power spectrum, we describe a method to efficiently compute the bispectrum in LPT, focusing on the Zeldovich approximation, in which contributions due to linear displacements are captured to all orders in a manifestly infrared (IR) safe way. We then isolate the effects of these linear displacements on oscillatory components of the power spectrum like baryon acoustic oscillations or inflationary primordial features and show that the Eulerian perturbation theory (EPT) prescription wherein their effects are resummed by a Gaussian damping of the oscillations arise as a saddle-point approximation of our calculation. These two methods of IR resummation are in excellent agreement at 1-loop in the bispectrum. At tree level, resummed EPT does less well to capture the nonlinear damping of the oscillations, and the LPT calculation does not require an artificial split of the power spectrum into smooth and oscillatory components, making the latter particularly useful for modeling exotic features. We finish by extending our analysis of IR resummation in LPT to N-point functions of arbitrary order.

The properties of the warm-hot intergalactic medium (WHIM) in the cosmic filaments are among the least quantified in modern astrophysics. eROSITA All Sky Survey (eRASS) provides us with a unique opportunity to study the X-ray emission of the WHIM. We applied both imaging and spectroscopic stacking techniques to the data of the first four eRASS scans to inspect the X-ray emissions from 7817 cosmic filaments identified from Sloan Digital Sky Survey optical galaxy samples. We obtained a 9-sigma significant detection of the total X-ray signal from filaments in the 0.3-1.2 keV band. We introduce a novel method to estimate the contamination fraction from unmasked X-ray halos, active galactic nuclei, and X-ray binaries that are associated with galaxies in filaments. We found an approximately 40% contamination fraction by these unmasked sources, suggesting that the remaining 60% of the signal could be from the WHIM and a 5.4-sigma detection significance of the WHIM. Moreover, we modeled the temperature and baryon density contrast of the detected WHIM by fitting the stacked spectrum and surface brightness profile. The best-fit temperature log(T/k)=6.84+/-0.07 obtained by using a single temperature model is marginally higher than the simulation results, which could be due to the fitting of a single temperature model on a multi-temperature spectrum. The best-fit baryon density contrast log\Delta_b=1.88+/-0.18 generally agrees with the X-ray emitting phases in the IllustrisTNG simulation, suggesting that the broadband X-ray emission traces the high temperature and density end of the entire WHIM population.

Recently, Biddle et al. (2024) claimed a non-detection of the protoplanet AB Aurigae b in Keck/NIRC2 Pa$\beta$ imaging. I reprocess these newly-public data and compare them to data from the extreme AO platform (SCExAO/CHARIS) used to discover AB Aur b. AB Aur b is decisively imaged with SCExAO/CHARIS at wavelengths covering Pa$\beta$. The Biddle et al. non detection of AB Aur b results from a far poorer image quality that is non competitive with SCExAO/CHARIS. Their contrast limits and thus constraints on accretion are overestimated due to an inaccurate AB Aur b source model. Consequentially, the revised Pa$\beta$ 2-$\sigma$ upper limit from these data is about three times higher than previously reported. Irrespective of image quality, single-band Pa$\beta$ imaging is ill suited to conclusively identifying accretion onto AB Aur b. Instead, high-resolution H$\alpha$ spectroscopy may provide accretion signatures. Aside from PDS 70, AB Aurigae remains the system with the strongest evidence for having a directly-imaged protoplanet.

Where and how flares efficiently accelerate charged particles remains an unresolved question. Recent studies revealed that a "magnetic bottle" structure, which forms near the bottom of a large-scale reconnection current sheet above the flare arcade, is an excellent candidate for confining and accelerating charged particles. However, further understanding its role requires linking the various observational signatures to the underlying coupled plasma and particle processes. Here we present the first study combining multi-wavelength observations with data-informed macroscopic magnetohydrodynamics and particle modeling in a realistic eruptive flare geometry. The presence of an above-the-looptop magnetic bottle structure is strongly supported by the observations, which feature not only a local minimum of magnetic field strength but also abruptly slowing down plasma downflows. It also coincides with a compact hard X-ray source and an extended microwave source that bestrides above the flare arcade. Spatially resolved spectral analysis suggests that nonthermal electrons are highly concentrated in this region. Our model returns synthetic emission signatures that are well-matched to the observations. The results suggest that the energetic electrons are strongly trapped in the magnetic bottle region due to turbulence, with only a small fraction managing to escape. The electrons are primarily accelerated by plasma compression and facilitated by a fast-mode termination shock via the Fermi mechanism. Our results provide concrete support for the magnetic bottle as the primary electron acceleration site in eruptive solar flares. They also offer new insights into understanding the previously reported small population of flare-accelerated electrons entering interplanetary space.

We introduce a pipeline that performs rapid image subtraction and source selection to detect transients, with a focus on identifying gravitational wave optical counterparts using the Dark Energy Camera (DECam). In this work, we present the pipeline steps from processing raw data to identification of astrophysical transients on individual exposures. We process DECam data and build difference images using the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines software, and we use flags and principal component analysis to select transients on a per-exposure basis, without associating the results from different exposures. Those candidates will be sent to brokers for further classification and alert distribution. We validate our pipeline using archival exposures that cover various types of objects, and the tested targets include a kilonova (GW170817), supernovae, stellar flares, variable stars (in a resolved galaxy or the Milky Way Bulge), and serendipitous objects. Overall, the data processing produces clean light curves that are comparable with published results, demonstrating the photometric quality of our pipeline. Real transients can be well selected by our pipeline when sufficiently bright (S/N $\gtrsim15$). This pipeline is intended to serve as a tool for the broader research community. Although this pipeline is designed for DECam, our method can be easily applied to other instruments and future LSST observations.

Exoplanets hosting M dwarfs are the best targets to characterize Earth-like or super-Earth planetary atmospheres with the James Webb Space Telescope (JWST). We determine detailed stellar parameters ($T_{\rm eff}$, log$g$, and $\xi$) and individual abundances of twelve elements for four cool M dwarfs hosting exoplanets TOI-1685, GJ 436, GJ 3470, and TOI-2445, scheduled for future observations by the JWST. The analysis utilizes high-resolution near-infrared spectra from the SDSS-IV APOGEE survey between 1.51-1.69$\mu$. Based on 1D-LTE plane-parallel models, we find that TOI-2445 is slightly metal-poor ([Fe/H] = -0.16$\pm$0.09), while TOI-1685, GJ 436 and GJ 3470 are more metal-rich ([Fe/H] = 0.06$\pm$0.18, 0.10$\pm$0.20 dex, 0.25$\pm$0.15). The derived C/O ratios for TOI-2445, TOI-1685, GJ 436, and GJ 3470 are 0.526$\pm$0.027, 0.558$\pm$0.097, 0.561$\pm$0.029, and 0.638$\pm$0.015, respectively. From results for 28 M dwarfs analyzed homogeneously from APOGEE spectra, we find exoplanet-hosting M dwarfs exhibit a C/O abundance ratio approximately 0.01 to 0.05 higher than those with non-detected exoplanets, at limits of a statistically significant offset. A linear regression of [Fe/H] \textit{vs.} C/O distribution reveals a noticeable difference in the angular coefficient between FGK dwarfs (0.27) and M dwarfs (0.13). Assuming our abundance ratios of Ca/Mg, Si/Mg, Al/Mg, and Fe/Mg, we determine a mass of 3.276$^{+0.448}_{-0.419}$$M_{\oplus}$ for TOI-2445 b, having density (6.793$^{+0.005}_{-0.099}$ this http URL$^{-3}$) and core mass fraction (0.329$_{-0.049}^{+0.028}$) very similar to Earth's. We also present an atlas of 113 well-defined spectral lines to analyze M dwarfs in the $H$-band and a comprehensive evaluation of uncertainties from variations in the atmospheric parameters, signal-to-noise, and pseudo-continuum.

A major goal of cosmology is to understand the nature of the field(s) which drove primordial Inflation. Through future observations, the statistics of large-scale structure will allow us to probe primordial non-Gaussianity of the curvature perturbation at the end of Inflation. We show how a new correlation statistic can significantly improve these constraints over conventional methods. Next-generation radio telescope arrays are under construction which will map the density field of neutral hydrogen to high redshifts. These telescopes can operate as an interferometer, able to probe small scales, or as a collection of single dishes, combining signals to map the large scales. We show how to fuse these these operating modes in order to measure the squeezed bispectrum with higher precision and greater economy. This leads to constraints on primordial non-Gaussianity that will improve on measurements by Planck, and out-perform other surveys such as Euclid. We forecast that $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}})\sim 3$, achieved by using a small subset, $\mathcal{O}(10^2-10^3)$, of the total number of accessible triangles. The proposed method can be used for any signal that peaks in squeezed configurations.

The coupling between the magnetic field and the gas during the collapsing phase of star-forming cores is strongly affected by the dust size distribution, which is expected to evolve. We aim to investigate the influence of key parameters on the evolution of the dust distribution as well as on the magnetic resistivities during the protostellar collapse. We perform collapsing single zone simulations with shark. The code computes the evolution of the dust distribution, accounting for different grain growth and destruction processes. It also computes the magnetic resistivities. We find that the dust distribution significantly evolves during the protostellar collapse, shaping the magnetic resistivities. The peak size of the distribution, the population of small grains and consequently the magnetic resistivities are controlled by both coagulation and fragmentation rates. Under standard assumptions, the small grains coagulate very early as they collide by ambipolar drift, yielding magnetic resistivities orders of magnitude away from the non-evolving dust case. In particular, the ambipolar resistivity \eta_AD is very high prior to nH=10^10 cm^-3, and as a consequence magnetic braking should be ineffective. In this case, large size protoplanetary disks should result, which is inconsistent with recent observations. To alleviate this tension, we identify mechanisms to reduce the ambipolar resistivity during the protostellar collapse. The most promising are namely: electrostatic repulsion and grain-grain erosion. The evolution of the magnetic resistivities during the protostellar collapse and consequently the shape of the magnetic field in the early life of the protoplanetary disk strongly depends on the possibility to repopulate the small grains or to prevent their early coagulation. Therefore, it is crucial to better constrain the collision outcomes and the dust grain elastic properties.

The H II region N88a in the Small Magellanic Cloud (SMC) is a spherical region of $1.5~\mathrm{pc}$ diameter, with high concentration of gas and dust, and at least four massive stars within it. Previous studies suggest that the four known sources may be insufficient to ionize the region and explain the nebular emission. In this contribution we analyze the ionizing photon production of the four known sources within the H II region. We compared the available photometry in the literature with the spectral energy distribution calculated for the ``Potsdam Wolf-Rayet'' (PoWR) models of massive stars for a wide range of dust extinction coefficient (Av). In particular, we selected models of OB-type atmospheres with SMC metallicity and compared the ionizing photon flux prediction with previous estimates based on the nebular emission of the H {\sc ii} region. We found that the Av values that best reproduce the photometry of each source vary from 0.82 to 3.84, increasing toward the dust band that runs through the nebula. In addition, all four sources are compatible with O-type stars with $\mathrm{T_{eff}} > 40 ~\mathrm{kK}$ and intrinsic ratios $\mathrm{F(LyC)/F(UV)}>1$. Lastly, the total ionizing photon rate predicted by the hottest models is $\mathrm{log(Q_H)} > 49.6 ~\mathrm{ph\ s}^{-1}$, which suggests that the stars could maintain the ionization state of the nebula.

Context. SN 2023ixf was discovered in Galaxy M101 in May 2023. Its proximity made it an extremely valuable opportunity for the scientific community to study the characteristics of the SN and its progenitor. A point source detected on archival images and hydrodynamical modelling of the bolometric light curve has been used to constrain the former star's properties. There is a significant variation in the published results regarding the initial mass of the progenitor. Nebular spectroscopy provides an independent tool to enhance our understanding of the supernova and its progenitor. Aims. We determine the SN progenitor mass by studying the first published nebular spectrum taken 259 days after the explosion. Methods. We analyze the nebular spectrum taken with GMOS at the Gemini North Telescope. Typical emission lines are identified, such as [O I], H{\alpha}, [Ca II], among others. Some species' line profiles show broad and narrow components, indicating two ejecta velocities and an asymmetric ejecta. We infer the progenitor mass of SN 2023ixf by comparing with synthetic spectra and by measuring the forbidden oxygen doublet flux. Results. Based on the flux ratio and the direct comparison with spectra models, the progenitor star of SN 2023ixf had a M_{ZAMS} between 12 and 15 M{_\odot}. In comparison with this, we find the use of the [O I] doublet flux to be less constraining of the progenitor mass. Our results agree with those from hydrodynamical modelling of the early light curve and pre-explosion image estimates pointing to a relatively low-mass progenitor.

A re-examination of high-resolution spectral monitoring of the W UMa-type binaries AW UMa and Epsilon CrA casts doubt on the widely utilized Lucy (1968a, 1968b) model of contact binaries. The detection of the very faint profile of the secondary component in AW UMa leads to a new spectroscopic determination of the mass ratio, q(sp) = 0.092 +/- 0.007, which is close to the previous, medium-resolution spectroscopic result of Pribulla & Rucinski (2008), q(sp) = 0.101 +/- 0.006, and remains substantially different from a cluster of generally accepted photometric results by several authors, concentrated around q(ph) = 0.080 +/- 0.005. The two approaches are independent, with the spectroscopic technique being more direct yet more demanding on telescope/spectrograph resources, while the photometric determinations are accessible to smaller telescopes but entirely dependent on the Lucy model. A survey of binaries with the best-determined values of the mass ratio shows a common tendency for q(ph) < q(sp). The tendency for systematically smaller values of q(ph) may result from the overfilling of the primary lobe and underfilling of the secondary lobe relative to the Roche model geometry, as predicted by the Stepien (2009) model; the tendency may be variable in time. Despite the observed moderate inter-systemic velocities, the photometric Lucy model may remain useful in providing approximate, though biased, results for the mass ratio. A complicating factor in detailed spectral analysis may be the occurrence of Enhanced Spectral-line Perturbations (ESP) projected over the secondary profiles, appearing in different numbers in the two studied binaries. The ESPs are tentatively identified within the Stepien model as collision fronts or fountains of hot, primary-component gas from the circumbinary, energy-carrying flow.

Regions of magnetic field with near-radial, Parker Spiral-like geometry known as quiescent regions have been observed in Parker Solar Probe data. These regions have very low $\delta B / \langle |B| \rangle$ compared to non-quiescent solar wind at the same heliocentric distances. Quiescent regions are observed to have lower solar wind bulk speeds, lower proton temperatures, and lower proton density, consistent with properties of the slow solar wind. Inside of 15 Rs, identified quiescent regions show distinct thermal properties, having higher proton temperature anisotropies and lower parallel plasma betas compared to switchback patches observed at the same heliocentric distances. When placed on $\mathcal{R}$ vs $\beta_{\parallel p}$ plots (where $\mathcal{R}$ is the proton temperature anisotropy), quiescent region solar wind is shown to be more stable to proton cyclotron and firehose instabilities than non-quiescent solar wind at the same heliocentric distances. It is shown that quiescent regions evolve similarly to the surrounding non-quiescent solar wind, but quiescent solar wind begins at a different location in the $\mathcal{R}$ vs $\beta_{\parallel p}$ parameter space, suggesting that these regions have separate origins than the more turbulent non-quiescent solar wind. Namely, enhanced temperature anisotropies and enhanced magnetic field strength may be consistent with magnetic field lines which have undergone less magnetic field expansion compared to non-quiescent wind at the same heliocentric distances.

WASP-107 b seems to be a poster child of the long-suspected high-eccentricity migration scenario. It is on a 5.7-day, polar orbit. The planet is Jupiter-like in radius but Neptune-like in mass with exceptionally low density. WASP-107 c is on a 1100-day, $e=0.28$ orbit with at least Saturn mass. Planet b may still have a residual eccentricity of $0.06\pm 0.04$: the ongoing tidal dissipation leads to the observed internally heated atmosphere and hydrodynamic atmospheric erosion. We present a population synthesis study coupling octopole Lidov-Kozai oscillations with various short-range forces, while simultaneously accounting for the radius inflation and tidal disruption of the planet. We find that a high-eccentricity migration scenario can successfully explain nearly all observed system properties. Our simulations further suggest that the initial location of WASP-107 b at the onset of migration is likely within the snowline ($<0.5\,{\rm AU}$). More distant initial orbits usually lead to tidal disruption or orbit crossing. WASP-107 b most likely lost no more than 20\% of its mass during the high-eccentricity migration, i.e. it did not form as a Jupiter-mass object. More vigorous tidally-induced mass loss leads to disruption of the planet during migration. We predict that the current-day mutual inclination between the planets b and c is substantial: at least 25-55$^\circ$ which may be tested with future Gaia astrometric observations. Knowing the current-day mutual inclination may further constrain the initial orbit of planet b. We suggest that the proposed high-eccentricity migration scenario of WASP-107 may be applicable to HAT-P-11, GJ-3470, HAT-P-18, and GJ-436 which have similar orbital architectures.

The next generation of spectroscopic surveys is expected to achieve an unprecedented level of accuracy in the measurement of cosmological parameters. To avoid confirmation bias and thereby improve the reliability of these results, blinding procedures become a standard practice in the cosmological analyses of such surveys. Blinding is especially crucial when the impact of observational systematics is important relative to the cosmological signal, and a detection of that signal would have significant implications. This is the case for local primordial non-gaussianity, as probed by the scale-dependent bias of the galaxy power spectrum at large scales that are heavily sensitive to the dependence of the target selection on the imaging quality, known as imaging systematics. We propose a blinding method for the scale-dependent bias signature of local primordial non-gaussianity at the density field level which consists in generating a set of weights for the data that replicate the scale-dependent bias. The applied blinding is predictable, and can be straightforwardly combined with other catalog-level blinding procedures that have been designed for the baryon acoustic oscillation and redshift space distortion signals. The procedure is validated through simulations that replicate data from the first year of observation of the Dark Energy Spectroscopic Instrument, but may find applications to other upcoming spectroscopic surveys.

Decameter hectometric (DH; 1-14 MHz) type-IV radio bursts are produced by flare-accelerated electrons trapped in post-flare loops or the moving magnetic structures associated with the CMEs. From a space weather perspective, it is important to systematically compile these bursts, explore their spectro-temporal characteristics, and study the associated CMEs. We present a comprehensive catalog of DH type-IV bursts observed by the Radio and Plasma Wave Investigation (WAVES) instruments onboard Wind and STEREO spacecraft, covering the period of white-light CME observations by the Large Angle and Spectrometric Coronagraph (LASCO) onboard the SOHO mission between November 1996 and May 2023. The catalog has 139 bursts, of which 73% are associated with a fast (>900 km/s) and wide (>60$^o$) CME, with a mean CME speed of 1301 km/s. All DH type-IV bursts are white-light CME-associated, with 78% of the events associated with halo CMEs. The CME source latitudes are within $\pm$45$^o$. 77 events had multi-vantage point observations from different spacecraft, letting us explore the impact of line of sight on the dynamic spectra. For 48 of the 77 events, there was good data from at least two spacecraft. We find that, unless occulted by nearby plasma structures, a type-IV burst is best viewed when observed within $\pm$60$^o$ line of sight. Also, the bursts with a duration above 120 min, have source longitudes within $\pm$60$^o$. Our inferences confirm the inherent directivity in the type-IV emission. Additionally, the catalog forms a sun-as-a-star DH type-IV burst database.

Considering a well motivated $f(R)$ modified-gravity model, in which an exponential function of the curvature is included, in this paper we implement a statistical data analysis to set constraints on the parameters of the model, taking into account an analytic approximate solution for the expansion rate, $H(z)$. From the Monte Carlo Markov Chain-based analysis of the expansion rate evolution, the standardized SN distance modulus and the redshift space distortion observational data, we find that the preferred value for the perturbative parameter, $b$, quantifying the deviation of the $f(R)$ model from $\Lambda$CDM, lives in a region which excludes $b = 0$ at $\gtrsim 4.5 \sigma$ C.L., and that the predicted current value of the Hubble parameter, $H_0$, locates in between the two observational results currently under scrutiny from Planck and SH0ES collaborations, indicating that the proposed model would alleviate the apparent tension. Under the implemented approximate solution, and with the constraints obtained for the parameters, the proposed $f(R)$ model successfully reproduce the observational data and the predicted evolution of interesting cosmological parameters resemble the results of $\Lambda$CDM, as expected, while an oscillatory behavior of the dark energy equation of state is observed, pointing to deviation from the concordance cosmological model. The results presented here reinforces the conclusion that the $f(R)$ modified-gravity model represents a viable alternative to describe the evolution of the Universe, evading the challenges faced by $\Lambda$CDM.

Measuring fundamental stellar parameters is key to fully comprehending the evolution of stars. However, current theoretical models over-predict effective temperatures, and under-predict radii, compared to observations of K and M dwarfs (radius inflation problem). In this work, we developed a model independent method to infer precise radii of single FGK and M dwarfs using Gaia DR3 parallaxes and photometry, and we used it to study the radius inflation problem. We calibrated nine surface brightness-color relations for the three Gaia magnitudes and colors using a sample of stars with angular diameter measurements. We achieved an accuracy of 4% in our angular diameter estimations, which Gaia's parallaxes allow us to convert to a physical radii. We validated our method by comparing our radius measurements with literature samples and the Gaia DR3 catalog, which confirmed the accuracy of our method and revealed systematic offsets in the Gaia measurements. Moreover, we used a sample with measured Halpha equivalent width (HaEW), a magnetic activity indicator, to study the radius inflation problem. We demonstrated that active stars have larger radii than inactive stars, showing that radius inflation is correlated with magnetic activity. We found a correlation between the radius inflation of active stars and HaEW for the mass bin 0.5<M[Msun]<= 0.6, but we found no correlation for lower masses. This could be due to lack of precision in our radius estimation or a physical reason. Radius measurements with smaller uncertainties are necessary to distinguish between the two scenarios.

In this article we confront a class of $f(Q)$ gravity models with observational data of galaxy-galaxy lensing. Specifically, we consider the $f(Q)$ gravity models containing a small quadratic correction when compared with General Relativity (GR), and quantify this correction by a model parameter $\alpha$. To derive the observational constraints, we start by extracting the spherically symmetric solutions which correspond to the deviations from the Schwarzschild solution that depends on the model parameter in a two-fold way, i.e., a renormalized mass and a new term proportional to $r^{-2}$. Then, we calculate the effective lensing potential, the deflection angle, the shear component, and the effective Excess Surface Density (ESD) profile. After that, we employ the group catalog and shape catalog from the SDSS DR7 for the lens and source samples respectively. Moreover, we handle the off-center radius as a free parameter and constrain it using the MCMC. Concerning the deviation parameter from GR we derive $\alpha=1.202^{+0.277}_{-0.179}\times 10^{-6} {\rm Mpc}^{-2}$ at 1 $\sigma$ confidence level, and then compare the fitting efficiency with the standard $\Lambda$CDM paradigm by applying the AIC and BIC information criteria. Our results indicate that the $f(Q)$ corrections alongside off-center effects yield a scenario that is slightly favored.

An analysis of the Cosmological Constant $\Omega_\Lambda$ fitted to subsamples of the Pantheon+ Type Ia SN sample spanning 2$\pi$ sterradians for a grid of 432 pole positions covering the whole sky reveals two large scale asymmetries. One of them is closely aligned with the Galactic North-South direction and the other points approximately towards RA$\sim 217.5^\circ$ Dec$\sim -26.4^\circ$, $\sim$50.9 degrees from the CMB dipole Apex. The signal to noise ratio (S/N) of the multiple $\Omega_\Lambda$ measurements in these directions is $3.2\lesssim $ S/N $\lesssim 8.4$. The first asymmetry is puzzling, and would indicate a systematic effect related with the distribution of Pantheon+ SNe on the sky and, probably, how the correction for reddening in the Galaxy is calculated. The second one, which entails a 2.8-$\sigma$ tension between $\Omega_\Lambda$ measure in opposite directions, bears strong implications on its interpretation as Dark Energy: It is consistent with the prediction for tilted observers located in a Friedmann-Robertson-Walker universe who could measure an acceleration or a deceleration with a dipolar asymmetry, irrespective of what the universe as a whole is doing. In this case, $\Omega_\Lambda$ would not be a physical entity, a real Dark Energy, but an apparent effect associated with the relativistic frame of reference transformation.

Glitching pulsars are expected to be important sources of gravitational waves. In this paper, we explore six different models that propose the emission of transient continuous waves, lasting days to months, coincident with glitches. The maximal gravitational wave energy is calculated for each model, which is then used to determine whether associated gravitational waves could be detectable with LIGO-Virgo-KAGRA's O4 detectors. We provide an analytical approximation to calculate the signal-to-noise ratio which includes information about the source's sky position, improving on previous estimates that assume isotropic or sky and orientation averaged sensitivities. Applying the calculation to the entire glitching population, we find that certain models predict detectable signals in O4, whereas others do not. We also rank glitching pulsars in order of how significant a signal would be, based on archival data, and we find that for all models, the Vela pulsar (PSR J0835$-$4510) would provide the strongest signal. Moreover, PSR J0537$-$6910 is not expected to yield a detectable signal in O4, but will start becoming relevant for next generation detectors. Our analysis also extends to the entire pulsar population, regardless of whether they have glitched or not, and we provide a list of pulsars that would present a significant signal, if they were to glitch. Finally, we apply our analysis to the latest April 2024 Vela glitch and find that a signal should be detectable under certain models. The non-detection of a supposedly detectable signal would provide an efficiency factor that quantifies how much a model can contribute to gravitational wave emission, eventually leading to a differentiation of models and independent constraints on physical parameters.

Context. Ch-type asteroids are distinctive among other dark asteroids in that they exhibit deep negative polarization branches (NPBs). Nevertheless, the physical and compositional properties that cause their polarimetric distinctiveness are less investigated. Aims. We aim to investigate the polarimetric uniqueness of Ch-type asteroids by making databases of various observational quantities (i.e., spectroscopic and photometric properties as well as polarimetric ones) of dark asteroids.Methods. We conducted an intensive polarimetric survey of 52 dark asteroids (including 31 Ch-type asteroids) in the R$_\mathrm{C}$-band to increase the size of polarimetric samples. The observed data are compiled with previous polarimetric, spectroscopic, and photometric archival data to find their correlations. Results. We find remarkable correlations between these observed quantities, particularly the depth of NPBs and their spectroscopic features associated with the hydrated minerals. The amplitude of the opposition effect in photometric properties also shows correlations with polarimetric and spectral properties. However, these observed quantities do not show noticeable correlations with the geometric albedo, thermal inertia, and diameter of asteroids. Conclusions. Based on the observational evidence, we arrive at our conclusion that the submicrometer-sized structures (fibrous or flaky puff pastry-like structures in phyllosilicates) in the regolith particles could contribute to the distinctive NPBs of hydrated asteroids.

This work focuses on minimum-time low-thrust orbit transfers from a prescribed low Earth orbit to a specified low lunar orbit. The well-established indirect formulation of minimum-time orbit transfers is extended to a multibody dynamical framework, with initial and final orbits around two distinct primaries. To do this, different representations, useful for describing orbit dynamics, are introduced, i.e., modified equinoctial elements (MEE) and Cartesian coordinates (CC). Use of two sets of MEE, relative to either Earth or Moon, allows simple writing of the boundary conditions about the two celestial bodies, but requires the formulation of a multiple-arc trajectory optimization problem, including two legs: (a) geocentric leg and (b) selenocentric leg. In the numerical solution process, the transition between the two MEE representations uses CC, which play the role of convenient intermediate, matching variables. The multiple-arc formulation at hand leads to identifying a set of intermediate necessary conditions for optimality, at the transition between the two legs. This research proves that a closed-form solution to these intermediate conditions exists, leveraging implicit costate transformation. As a result, the parameter set for an indirect algorithm retains the reduced size of the typical set associated with a single-arc optimization problem. The indirect heuristic technique, based on the joint use of the necessary conditions and a heuristic algorithm (i.e., differential evolution in this study) is proposed as the numerical solution method, together with the definition of a layered fitness function, aimed at facilitating convergence. The minimum-time trajectory of interest is sought in a high-fidelity dynamical framework, with the use of planetary ephemeris and the inclusion of the simultaneous gravitational action of Sun, Earth, and Moon, along the entire transfer path.

Departing from standard slow-roll conditions is one way of putting the inflationary paradigm to test, and constraining the dynamics of the inflaton field with a constant-rate of roll of the inflaton field, a.k.a. the constant-roll scenario, is one way of exploring such deviation from the standard slow-roll dynamics. In this manuscript we explore such a possibility in a variant inflationary scenario, known as Warm Inflation. We construct and derive the conditions for having constant-roll WI models where inflation lasts at least for 60 $e-$folds, gracefully exits the constant-roll inflation phase, and maintains near thermal equilibrium of the system which is an essential feature of WI in the slow-roll regime. We show that while certain models of WI (the ones with dissipative coefficient as a function of temperature alone) can accommodate constant-roll dynamics, others (with dissipative coefficient as a function of temperature and the inflaton field both) fail to maintain thermal equilibrium once the constant-roll condition is imposed and hence cannot produce a constant-roll WI phase.

Diffusion tensor coefficients play a central role in describing cosmic-ray transport in various astrophysical environments permeated with magnetic fields, which are usually modeled as a fluctuating field on top of a mean field. In this article, a formal derivation of these coefficients is presented by means of the calculation of velocity decorrelation functions of particles. It relies mainly on expanding the 2-pt correlation function of the (fluctuating) magnetic field experienced by the particles between two successive times in the form of an infinite Dyson series and retaining a class of terms that converge to a physical solution. Subsequently, the velocity decorrelation functions, themselves expressed as Dyson series, are deduced from an iteration procedure that improves on the partial summation scheme. The results are shown to provide approximate solutions compared to those obtained by Monte-Carlo simulations as long as the Larmor radius of the particles is larger than at least one tenth of the largest scale of the turbulence.

The ~200 m/s impact of a single 400-kg Bjurböle L/LL ordinary chondrite meteorite onto sea ice resulted in the catastrophic disruption of the projectile. This resulted in a significant fraction of decimeter-sized fragments that exhibit power law cumulative size and mass distributions. This size range is underrepresented in impact experiments and asteroid boulder studies. The Bjurböle projectile fragments share similarities in shape (sphericity, and roughness at small and large scale) with asteroid boulders. However, the mean aspect ratio (3D measurement) and apparent aspect ratio (2D measurement) of Bjurböle fragment is 0.83 and 0.77, respectively, indicating that Bjurböle fragments are more equidimensional compared to both fragments produced in smaller scale impact experiments and asteroid boulders. These differences may be attributed either to the fragment source (projectile vs. target), to the high porosity and low strength of Bjurböle, to the lower impact velocity compared with typical asteroid collision velocities, or potentially to fragment erosion during sea sediment penetration or cleaning.

We present the results of a detailed study on the detectability of the High Frequency Feature (HFF) in core-collapse supernova (CCSN) gravitational wave (GW) signals. We applied Residual Neural Networks (ResNet50), one of the state-of-the-art deep learning architectures in computer vision, to perform a multi-class classification of image samples built from time-frequency Morlet wavelet scalograms of LIGO-Virgo noise plus CCSN GW signals. We consider three target labels for three consecutive and mutually exclusive intervals in which the (first-order) slope of the HFF can be located. We optimized, trained, and tested the ResNet50 model with phenomenological waveforms. Next, we tested the optimized ResNet50 model with GW signals from CCSN simulations. At galactic distances of $1$Kpc and $5$Kpc with H1 and L1 data and $1$Kpc with V1 data, we obtained highly accurate results (test accuracies from $0.8933$ to $0.9867$), which show the feasibility of our methodology. In the case of further distances, we observed declines in test performance until $0.8000$ with H1 and L1 data at $10$Kpc and until $0.5933$ with V1 data at $10$Kpc. Without assuming the continuity and/or discontinuity of the HFF slope values, our methodology is general enough to address, at an early stage, the characterization of the HFF.

In the star formation regions, the complex organic molecules that contain peptide bonds ($-$NH$-$C(=O)$-$) play a major role in the metabolic process because $-$NH$-$C(=O)$-$ is connected to amino acids (R-CHNH$_{2}$$-$COOH). Over the past few decades, many complex organic molecules containing peptide-like bonds have been detected in hot molecular cores, hot corinos, and cold molecular clouds, however, their prebiotic chemistry is poorly understood. We present the first detection of the rotational emission lines of formamide (NH$_{2}$CHO) and isocyanic acid (HNCO), which contain peptide-like bonds towards the chemically rich hot molecular core G358.93$-$0.03 MM1, using high-resolution and high-sensitivity Atacama Large Millimeter/submillimeter Array (ALMA) bands 6 and 7. We estimate that the column densities of NH$_{2}$CHO and HNCO towards G358.93$-$0.03 MM1 are (2.80$\pm$0.29)$\times$10$^{15}$ cm$^{-2}$ and (1.80$\pm$0.42)$\times$10$^{16}$ cm$^{-2}$ with excitation temperatures of 165$\pm$21 K and 170$\pm$32 K, respectively. The fractional abundances of NH$_{2}$CHO and HNCO towards G358.93$-$0.03 MM1 are (9.03$\pm$1.44)$\times$10$^{-10}$ and (5.80$\pm$2.09)$\times$10$^{-9}$. We compare the estimated abundances of NH$_{2}$CHO and HNCO with the existing three-phase warm-up chemical model abundance values and notice that the observed and modelled abundances are very close. We conclude that NH$_{2}$CHO is produced by the reaction of NH$_{2}$ and H$_{2}$CO in the gas phase towards G358.93$-$0.03 MM1. Likewise, HNCO is produced on the surface of grains by the reaction of NH and CO towards G358.93$-$0.03 MM1. We also found that NH$_{2}$CHO and HNCO are chemically linked towards G358.93$-$0.03 MM1.

In this study, we investigate the homologous collapse dynamics of a slowly rotating stellar nuclear core. We incorporate the virial theorem of hadron collisional relaxations to analyze the conversion of gravitational potential energy into kinetic energy of hadrons. Additionally, we consider the production of photons and electron-positron pairs resulting from hadron collisions. Our analysis reveals a gravo-thermal process wherein gravitational energy transforms into energy associated with photons and pairs. The presence of axial symmetric centrifugal potential reduces the gravitational potential energy gain. The presence of axial symmetric centrifugal potential mitigates the gain in gravitational potential energy. Consequently, we observe the formation of highly energetic and opaque photon-pair jets, whose size, peak spectra, energy and number densities are consistent with observed energetic features and empirical relationships of Gamma-Ray Burst (GRB) progenitors. Furthermore, we find that the angular momentum and mass ratio ($J/M$) of the collapsing core and binary coalescence play crucial roles in determining the jet angle and the differentiation between short and long bursts. Our findings contribute to a deeper understanding of the mechanisms governing GRBs' progenitor systems.

The kinematics of the $\epsilon$ Cha young stellar association close to the Sun has been studied based on a list of candidate stars from the Dickson-Vandervelde work. The working sample consists of 26 stars with parallaxes, proper motions from the Gaia~DR3 catalog and radial velocities taken from literary sources. The orbits of the stars back to the past were constructed, and the moment when the association had a minimum spatial size was determined, as well as an analysis of the dependencies of the velocities $U,V,W$ on the coordinates $x,y,z$ was carried out. It is shown that the initial sample is divided into two parts with different kinematic properties. The first sample included 9 stars. Based on the construction of the orbits of these 9 stars, an age estimate of $t=4.9\pm0.8$ million years was obtained. An expansion coefficient in the $xz$ plane with the value $K_{xz}=135\pm19$ km/s/kpc was also found for them, on the basis of which another age estimate $t=7.2\pm1.0$ million years was obtained. The second sample included 17 stars. The construction of their orbits gave an estimate of age $t=0.2\pm0.3$ million years, and based on the gradient $\partial W/\partial z=707\pm248$ km/s/kpc, a second estimate of their age $t=1.4\pm0.5$ million years was obtained. This suggests that the $\epsilon$ Cha association either consists of two groupings of different ages, or a younger one arose as a result of a recent outbreak of star formation within a common star system. The question of the gravitational connection of the groupings has not been considered in the framework of this work.

Using a large sample of fast radio bursts (FRBs) from the first CHIME/FRB catalog, we apply the Lynden-Bell's c$^-$ method to study their energy function and formation rate evolutions with redshift. It is found with the non-parametric Kendell's $\tau$ statistics that the FRB energy strongly evolves with the cosmological redshift as $E(z)\propto(1 + z)^{5.23}$. After removing the redshift dependence, the local energy distribution can be described by a broken power-law form of $\Psi(E_{0})\propto E_{0}^{-0.38}$ for the low-energy segment and $\Psi(E_{0})\propto E_{0}^{-2.01}$ for the high-energy segment with a dividing line of $\sim2.1\times10^{40} \rm erg$. Interestingly, we find that the formation rate of CHIME FRBs also evolves with redshift as $\rho(z)\propto(1+z)^{-4.73\pm0.08}$. The local formation rate $\rho(0)$ of the CHIME FRBs is constrained to be about $ 1.25\times 10^4\rm{\,Gpc^{-3}yr^{-1}}$ that is comparable with some previous estimations. In addition, we notice the formation rate not only exceeds the star formation rate at the lower redshifts but also always declines with the increase of redshift, which does not match the star formation history at all. Consequently, we suggest that most FRBs could originate from the older stellar populations.

Understanding the origin of fast radio bursts (FRBs) has become the main science driver of recent dedicated FRB surveys. Between July 2019 and February 2022, we carried out ALERT, an FRB survey at 1370 MHz using the Apertif instrument installed at the Westerbork Synthesis Radio Telescope (WSRT). Here we report the detection of 18 new FRBs, and we study the properties of the entire 24 burst sample detected during the survey. For five bursts, we identify host galaxy candidates with >50% probability association. We observe an average linear polarisation fraction of $\sim$43% and an average circular polarisation fraction consistent with 0%. A third of the FRBs display multiple components. The sample next reveals a population of highly scattered bursts, which is most likely to have been produced in the immediate circumburst environment. Furthermore, two FRBs show evidence for high rotation measures, reaching |RM|>$10^3$ rad m$^{-2}$ in the source reference frames. Together, the scattering and rotation measures ALERT finds prove that a large fraction of FRBs are embedded in complex media such as star forming regions or supernova remnants. Through the discovery of the third most dispersed FRB so far, we show that one-off FRBs can emit at frequencies in excess of 6 GHz. Finally, we determine an FRB all-sky rate of $459^{+208}_{-155}$ sky$^{-1}$ day$^{-1}$ above a fluence limit of 4.1 Jy ms, and a fluence cumulative distribution with a power law index $\gamma=-1.23\pm0.06\pm0.2$, which is roughly consistent with the Euclidean Universe predictions. Through the high resolution in time, frequency, polarisation and localisation that ALERT featured, we were able to determine the morphological complexity, polarisation, local scattering and magnetic environment, and high-frequency luminosity of FRBs. We find all these strongly resemble those seen in young, energetic, highly magnetised neutron stars.

NGC 5128 (Cen A) is the nearest giant elliptical galaxy and one of the brightest extragalactic radio sources in the sky, boasting a prominent dust lane and jets emanating from its nuclear supermassive black hole. In this paper, we construct the star formation history (SFH) of two small fields in the halo of NGC 5128: a northeastern field (Field 1) at a projected distance of $\sim 18.8$ kpc from the center, and a southern field (Field 2) $\sim 9.9$ kpc from the center. Our method is based on identifying long period variable (LPV) stars that trace their sibling stellar population and hence historical star formation due to their high luminosity and strong variability; we identified 395 LPVs in Field 1 and 671 LPVs in Field 2. Even though the two fields are $\sim 28$ kpc apart on opposite sides from the center, they show similar SFHs. In Field 1, star formation rates (SFRs) increased significantly around $t\sim 800$ Myr and $t\sim 3.8$ Gyr; and in Field 2, SFRs increased considerably around $t\sim 800$ Myr, $t\sim 3.8$ Gyr, and $t\sim 6.3$ Gyr, where $t$ is look--back time. The increase in SFR $\sim 800$ Myr ago agrees with previous suggestions that the galaxy experienced a merger around that time. The SFH reconstructed from LPVs supports a scenario in which multiple episodes of nuclear activity lead to episodic jet-induced star formation. While there is no catalog of LPVs for the central part of NGC 5128, applying our method to the outer regions (for the first time in a galaxy outside the Local Group) has enabled us to put constraints on the complex evolution of this cornerstone galaxy.

We study the evolution of electric currents during the emergence of magnetic flux in the solar photosphere and the differences exhibited between solar active regions of different Hale complexity classes. A sample of 59 active regions was analyzed using a method based on image segmentation and error analysis to determine the total amount of non-neutralized electric current along their magnetic polarity inversion lines. The time series of the total unsigned non-neutralized electric current, $I_{NN,tot}$, exhibit intricate structure in the form of distinct peaks and valleys. This information is largely missing in the respective time series of the total unsigned vertical electric current $I_z$. Active regions with $\delta$- spots stand out, exhibiting 1.9 times higher flux emergence rate and 2.6 times higher $I_{NN,tot}$ increase. The median value of their peak $I_{NN,tot}$ is equal to $3.6\cdot10^12$ A, which is more than three times higher than that of the other regions of the sample. An automated detection algorithm was also developed to pinpoint the injection events of non-neutralized electric current. The injection rates and duration of these events were higher with increasing complexity of active regions, with regions containing $\delta$-spots exhibiting the strongest and longest events. These events do not necessarily coincide with increasing magnetic flux, although they exhibit moderate correlation. We conclude that net electric currents are injected during flux emergence, but are also shaped drastically by the incurred photospheric evolution, as active regions grow and evolve.

We report on the critical influence of small-scale flow structures (e.g., fronts, vortices, and waves) that immediately arise in hot-exoplanet atmosphere simulations initialized with a resting state. A hot, 1:1 spin-orbit synchronized Jupiter is used here as a clear example; but, the phenomenon is generic and important for any type of hot synchronized planet--gaseous, oceanic, or telluric. When the early-time structures are not captured in simulations (due to, e.g., poor resolution and/or too much dissipation), the flow behavior is markedly different at later times--in an observationally significant way; for example, the flow at large-scale is smoother and much less dynamic. This results in the temperature field, and its corresponding thermal flux, to be incorrectly predicted in numerical simulations, even when the quantities are spatially averaged.

The Cassiopée project aims to develop the key technologies that will be used to deploy very high-performance Adaptive Optics for future ELTs. The ultimate challenge is to detect earth-like planets and characterize the composition of their atmosphere. For this, imaging contrasts of the order of 109 are required, implying a leap forward in adaptive optics performance, with high density deformable mirrors (120x120 actuators), low-noise cameras and the control of the loop at few kHz. The project brings together 2 industrial partners: First Light Imaging and ALPAO, and 2 academic partners: ONERA and LAM, who will work together to develop a new camera for wavefront sensing, a new deformable mirror and their implementation in an adaptive optics loop. This paper will present the development of the fast large infrared e-APD camera which will be used in the wavefront sensor of the system. The camera will integrate the latest 512x512 Leonardo e-APD array and will benefit from the heritage of the first-light imaging's C-RED One camera. The most important challenges for the application are the autonomous operation, vibration control, background limitation, compactness, acquisition speed and latency.

The recent result from the DESI (Dark Energy Spectroscopic Instrument) data release 1 (DR1) combined with the cosmic microwave background and supernova data shows a preference for the dynamical dark energy over the cosmological constant [1]. For this analysis, the Chevallier-Polarski-Linder (CPL) parameterization of the dark energy equation of state has been used. Though CPL parametrization is the most popular parametrization of dark energy equation of state (EoS) recent studies have questioned its choice of prior and degeneracy. In this work, we propose a new parametrization of the dark energy equation of state motivated by its ability to cross the phantom barrier crossing. At a high redshift limit, this parametrization takes the form of CPL parametrization. We have used the Bayesian model selection criteria in which this parametrisation shows a moderate positive preference over the $\Lambda$CDM model. Contrary to the DESI2024 results our model suggests the current value of the dark energy equation of state to be phantom today with two phantom barrier crossings one at $z\simeq 0.4$ from deep phantom to quintessence which matches with the finding from the DESI2024 result and another phantom barrier crossing very recently at $z\simeq0.1$ from quintessence to phantom. The proposed model also weakens the Hubble tension significantly to $\approx 2 \sigma$ while comparing with the result from [2] and to $\approx 1.4 \sigma$ while comparing with [3].

By combining newly obtained deep GBT 21cm observations with optical spectroscopic data, we present an analysis of the gas content of BreakBRD galaxies, a population denoted by their blue star-forming centers and red quenched disks that do not appear to follow the typical inside-out evolution of spiral galaxies. We confirm previous results that the neutral atomic hydrogen (HI) gas fractions of BreakBRDs are on-average lower than those of typical galaxies on the star-forming sequence (SFS), and find that their \ion{H}{1} fractions are generally higher than Green Valley (GV) galaxies. HI depletion times for BreakBRDs are roughly an order of magnitude lower than those of SFS galaxies, in stark contrast with GV galaxies that typically have much longer depletion times than SFS galaxies. The nuclear gas-phase metallicities of BreakBRDs have a broader distribution than SFS galaxies and are skewed towards slightly above-average values. BreakBRDs are systematically offset from the Baryonic Tully-Fisher Relation towards lower baryonic mass at a given rotation velocity. They also have higher typical HI asymmetries than SFS galaxies, and of those galaxies with spatially resolved gas velocity fields from the SDSS-IV MaNGA survey, two-thirds are either highly distorted or completely misaligned relative to the stellar disk. Evidence supports a scenario where BreakBRDs are in an early phase of quenching, and there is mixed evidence that their behavior is related to past merger activity.

Complex organic molecules (COMs) have been widely observed in molecular clouds and protostellar environments. One of the formation mechanisms of COMs is radical reactions on the icy grain surface driven by UV irradiation. While many experiments have reported that various COMs can be synthesized under such ice conditions, the majority of the reaction processes are unclear. Complementary numerical simulations are necessary to unveil the synthetic process behind the formation of COMs. In this study, we develop a chemical reaction simulation using a Monte Carlo method. To explore the complex reaction network of COM synthesis, the model was designed to eliminate the need to prepare reaction pathways and to keep computational costs low. With this simulation, we investigate the chemical reactions occurring on icy dust surfaces during and after UV irradiation, assuming a protoplanetary disk environment. We aim to reveal the types of organic molecules produced in a disk and the formation mechanisms of COMs, in particular, amino acids and sugars. The results show that photodissociation and subsequent radical-radical reactions cause random rearrangement of the covalent bonds in the initial molecules composed of methanol, formaldehyde, ammonia, and water. Consequently, highly complex molecules such as amino acids and sugars were produced in a wide range of the initial conditions. We found that the final abundances of amino acids and sugars have extremely similar dependence on the atomic ratios of the initial molecules, which peak at C/H~0.1-0.3 and O/H~0.3-0.5, although the amino acids abundance is usually more than ten times higher than that of sugars. To understand this dependence, a semi-analytical formula was derived. Additionally, parameter surveys have suggested that the decomposition reactions of amino acids and sugars undergo a rapid transition within the threshold of a given parameter.

Clusters of galaxies are merging during the formation of large-scale structures in the Universe. Based on optical survey data, we identify a large sample of pre-mergers of galaxy clusters and merging subclusters in rich clusters. We find 39,382 partners within a velocity difference of 1500 km/s and a projected separation of 5r_{500} around 33,126 main clusters, where r_{500} is the radius of the main cluster. Based on the galaxy distribution inside rich clusters with more than 30 member galaxy candidates, we identify subclusters by modeling the smoothed optical distribution with a two-component profile, and a coupling factor is obtained for merging subclusters in 7845 clusters. In addition, we find 3446 post-collision mergers according to the deviations of brightest cluster galaxies from other member galaxies, most of which have been partially validated by using the Chandra and XMM-Newton X-ray images. Two new bullet-like clusters have been identified by using the optical and X-ray images. The large samples of merging clusters of galaxies presented here are important databases for studying the hierarchical structure formation, cluster evolution, and the physics of intergalactic medium.

By performing general relativistic radiative transfer calculations, we show the radio images of relativistic jets including highly magnetized regions inside jet funnels, based on steady, axisymmetric, and semi-analytic general relativistic magnetohydrodynamics models. It is found that multiple ring images appear at the photon frequency of 230 GHz for nearly pole-on observers, because of the strong light bending effect on photons generated at the separation surfaces which is the boundary between the inflow and outflow flows in the jet funnel. A bright teardrop-shaped component, which extends from the bright rings of the separation surface, also appears in the counter jet region. The diameter of the brightest outermost ring originated from the counter jet is $\sim 60 ~ \mu {\rm as}$, which is consistent with the ring-like images of M87 at 86 GHz observed with GMVA, ALMA and GLT, whose ring-diameter is $\sim 64^{+4}_{-8} ~\mu {\rm as}$. The thinner and smaller-diameter rings are exhibited when the black-hole spin magnitude is higher. These morphological features are expected to appear without being prominently affected by the detailed MHD-plasma parameters of our GRMHD jet model, since the location of the separation surfaces is mainly regulated by the black hole spin. Our GRMHD model and the emission features of the images in the horizon-scale, highly magnetized jet funnel may be tested by future observations, e.g., the next-generation Event Horizon Telescope and the Black Hole Explorer.

In the mixed dark matter scenarios consisting of primordial black holes (PBHs) and particle dark matter (DM), PBHs can accrete surrounding DM particles to form ultracompact minihalos (UCMHs or clothed PBHs) even at early epoch of the Universe. The distribution of DM particles in a UCMH follows a steeper density profile compared with classical DM halo. It is expected the DM annihilation rate is very large in UCMHs resulting in a contribution to, e.g., the extragalactic neutrino flux. In this work, we investigate the extragalactic neutrino flux from clothed PBHs due to DM annihilation, and then the muon flux for neutrino detection. Compared with the atmospheric neutrino flux, we derive the upper limits on the cosmological abundance of PBHs for 10 years exposure time of, e.g., the IceCube experiment. For the DM mass $m_{\chi}=100~(1000)$ GeV, the upper limits (2$\sigma$) on the fraction of DM in PBHs are $f_{\rm PBH}=1.2\times 10^{-3}~(8.9\times 10^{-5})$ for contained events and $f_{\rm PBH}=2.5\times 10^{-3}~(1.3\times 10^{-5})$ for upward events, respectively. Compared with other constrains, although the upper limits obtained by us are not the strongest, it is a different way to study the cosmological abundance of PBHs.

The local universe is highly inhomogeneous and anisotropic. We live in a relatively sparse region of the Laniakea supercluster at the edge of a large 80 Mpc-wide void. We study the effect of these inhomogeneities on the measured gravitational wave event rates. In particular, we estimate how the measured merger rate of compact binaries is biased by the local matter distribution. The effect of the inhomogeneities on the merger rate is suppressed by the low angular resolution of gravitational wave detectors coupled with their smoothly decreasing population-averaged sensitivity with distance. We estimate the effect on the compact binary coalescence event rate to be at most 6% depending of the chirp mass of the target binary system and the sensitivity and orientation of the detectors.

It is commonly believed that the dissipative properties of superdense matter play a negligible role in modeling gravitational waveforms from neutron star inspirals. This study aims to investigate whether this presumption holds true for the often neglected dissipative process associated with particle diffusion in superconducting neutron stars. As we demonstrate, diffusion effects can significantly impact the phase of the gravitational wave from the inspiral, manifesting at a magnitude of a few tens of milliradians at large orbit separations, equivalent to orbital frequencies of a few hertz. We also find that dissipation resulting from particle diffusion might increase the neutron star's temperature to approximately $10^7\rm K$ during the inspiral.

The evolution of Active Galactic Nuclei (AGN) is closely connected to their host galaxies and surroundings. Via feedback processes, AGN can counteract the cooling of the intracluster medium (ICM) and suppress star formation in their host galaxies. Radio observations at low frequencies provide a glimpse into the history of AGN activity. The Virgo cluster is a substantial reservoir of nearby galaxies and provides an ideal laboratory for the study of AGN as well as their feedback mechanisms. The aim of our work is to characterise the AGN population within the Virgo cluster down to low radio luminosities, constrain the AGN duty cycle and investigate environmental feedback in cluster member galaxies. We analyse 144 MHz and 1.3 GHz radio observations of early-type galaxies from the ACS Virgo Cluster Survey (ACSVCS) taken with LOFAR and MeerKAT. We detect 12 of these galaxies at 144 MHz, 5 of which show clearly extended radio emission. The radio luminosity shows a strong dependence on the stellar mass of the host galaxy, in agreement with previous results. As a notable outlier, the massive elliptical galaxy NGC 4365 ($M_* = 2.2 \times 10^{11} M_\odot$) is not detected as a compact source in the LOFAR observations. Instead, it is surrounded by diffuse, low-surface brightness emission, which hints towards a past phase of stronger nuclear activity. Furthermore, we find a cavity in NGC 4472 (= M 49) inflated by the wide-angle tail only visible in the LOFAR data, which implies that the cavity was created by a past outburst. The corresponding cavity power is of the same order of magnitude as the jet power in the present duty cycle of the AGN.

Located at the bottom of the main sequence, ultracool dwarf stars are widespread in the solar neighbourhood. Nevertheless, their extremely low luminosity has left their planetary population largely unexplored, and only one of them, TRAPPIST-1, has so far been found to host a transiting planetary system. In this context, we present the SPECULOOS project's detection of an Earth-sized planet in a 17 h orbit around an ultracool dwarf of M6.5 spectral type located 16.8 pc away. The planet's high irradiation (16 times that of Earth) combined with the infrared luminosity and Jupiter-like size of its host star make it one of the most promising rocky exoplanet targets for detailed emission spectroscopy characterization with JWST. Indeed, our sensitivity study shows that just ten secondary eclipse observations with the Mid-InfraRed Instrument/Low-Resolution Spectrometer on board JWST should provide strong constraints on its atmospheric composition and/or surface mineralogy.

The astrophysical origins of the majority of the IceCube neutrinos remain unknown. Effectively characterizing the spatial distribution of the neutrino samples and associating the events with astrophysical source catalogs can be challenging given the large atmospheric neutrino background and underlying non-Gaussian spatial features in the neutrino and source samples. In this paper, we investigate a framework for identifying and statistically evaluating the cross-correlations between IceCube data and an astrophysical source catalog based on the $k$-Nearest Neighbor Cumulative Distribution Functions ($k$NN-CDFs). We propose a maximum likelihood estimation procedure for inferring the true proportions of astrophysical neutrinos in the point-source data. We conduct a statistical power analysis of an associated likelihood ratio test with estimations of its sensitivity and discovery potential with synthetic neutrino data samples and a WISE-2MASS galaxy sample. We apply the method to IceCube's public ten-year point-source data and find no statistically significant evidence for spatial cross-correlations with the selected galaxy sample. We discuss possible extensions to the current method and explore the method's potential to identify the cross-correlation signals in data sets with different sample sizes.

The carbonaceous macromolecules imprinting in astronomical spectra the numerous absorptions called Diffuse Interstellar Bands (DIBs) are omnipresent in the Galaxy and beyond and represent a considerable reservoir of organic matter. However, their chemical formulae, formation and destruction sites remain open questions. Their spatial distribution and the local relation to other interstellar species is paramount to unravel their role in the lifecycle of organic matter. Volume density maps bring local instead of line-of-sight distributed information, and allow for new diagnostics. We present the first large-scale volume (3D) density map of a DIB carrier and compare it with an equivalent map of interstellar dust. The DIB carrier map is obtained through hierarchical inversion of about 202,000 measurements of the 8621 nm DIB obtained with the Gaia-RVS instrument. It covers about 4000 pc around the Sun in the Galactic plane. A dedicated interstellar dust map is built based on extinction towards the same target stars. At the 50 pc resolution of the maps, the 3D DIB distribution is found remarkably similar in shape to the 3D distribution of dust. On the other hand, the DIB-to-dust local density ratio increases in low-dust areas. It is also increasing away from the disk, however, the minimum ratio is found to be shifted above the Plane to Z = +50pc. Finally, the average ratio is also surprisingly found to increase away from the Galactic Center. We suggest that the three latter trends may be indications of a dominant contribution of material from the carbon-rich category of dying giant stars to the formation of the carriers. Our suggestion is based on recent catalogues of AGB stars and estimates of their mass fluxes of C-rich and O-rich ejecta.

We investigate the initiation of cosmic inflation, in full numerical relativity, from pre-inflationary scenarios with large tensor and vector fluctuations in the metric. These settings are characterized by having large values in the Weyl curvature tensor. In the matter sector, we consider a single scalar field with inhomogeneous field velocities, corresponding to a kination period. In the context of large-field inflation, it is shown that the onset of inflation continues to be robust to this type of initial conditions, and that during inflation the Universe successfully homogenizes and flattens any type of curvature, leading to a Friedmann-Lemaître-Robertson-Walker Universe after just a few tens of efolds of accelerating expansion.

We present significant improvements to our previous work on noise reduction in {\sl Herschel} observation maps by defining sparse filtering tools capable of handling, in a unified formalism, a significantly improved noise reduction as well as a deconvolution in order to reduce effects introduced by the limited instrumental response (beam). We implement greater flexibility by allowing a wider choice of parsimonious priors in the noise-reduction process. More precisely, we introduce a sparse filtering and deconvolution approach approach of type $l^2$-$l^p$, with $p > 0$ variable and apply it to a larger set of molecular clouds using {\sl Herschel} 250 $\mu $m data in order to demonstrate their wide range of application. In the {\sl Herschel} data, we are able to use this approach to highlight extremely fine filamentary structures and obtain singularity spectra that tend to show a significantly less $\log$-normal behavior and a filamentary nature in the less dense regions. We also use high-resolution adaptive magneto-hydrodynamic simulation data to assess the quality of deconvolution in such a simulated beaming framework.

We analyze more than 15 years of the \gr\ data, obtained with the Large Area Telescope (LAT) onboard {\it the Fermi Gamma-ray Space Telescope (Fermi)}, for the region of the young supernova remnant (SNR) RCW~103, since the nearby source 4FGL J1616.2$-$5054e, counterpart to HESS~J1616$-$518 and $\simeq$13\,arcmin away from the SNR, is determined to be extended in the more recent Fermi-LAT source catalog. Different templates for 4FGL J1616.2$-$5054e and RCW~103 are tested, and we find that a point source with a power-law (PL) spectrum at the southern limb of the SNR best describes the detected gamma-ray emission. The photon index of the PL emission is $\Gamma\simeq 2.31$, softer than the previously reported $\Gamma\simeq 2.0$ when the counterpart to HESS~J1616$-$518 was considered to be a point source (which likely caused mis-identification of extended emission at RCW~103). In order to produce the \gr\ emission in a hadronic scenario, we estimate that protons with an index$\sim$2.4 PL energy distribution are needed. These results fit with those from multi-wavelength observations that have indicated the remnant at the southern limb is interacting with a molecular cloud.

Recently, complex horizontal patterns of umbral oscillations have been reported, but their physical nature and origin are still not fully understood. Here we show that the two-dimensional patterns of umbral oscillations of slow waves are inherited from the subphotospheric fast-body modes. Using a simple analytic model, we successfully reproduced the temporal evolution of oscillation patterns with a finite number of fast-body modes. In this model, the radial apparent propagation of the pattern is associated with the appropriate combination of the amplitudes in radial modes. We also find that the oscillation patterns are dependent on the oscillation period. This result indicates that there is a cutoff radial mode, which is a unique characteristic of the model of fast-body modes. In principle, both internal and external sources can excite these fast-body modes and produce horizontal patterns of umbral oscillations.

The existence of intermediate-mass black holes (IMBHs) is crucial for understanding various astrophysical phenomena, yet their existence remains elusive, except for the LIGO-Virgo detection. We report the discovery of a high-velocity star J0731+3717, whose backward trajectory about 21 Myr ago intersects that of globular cluster M15 within the cluster tidal radius. Both its metallicity [Fe/H] and its alpha-to-iron abundance ratio [$\alpha$/Fe] are consistent with those of M15. Furthermore, its location falls right on the fiducial sequence of the cluster M15 on the color-absolute magnitude diagram, suggesting similar ages. These support that J0731+3717 is originally associated with M15 at a confidence level of $5.4\sigma$. We find that such a high-velocity star ($V_{\rm ej} = 548^{+6}_{-5}$ km s$^{-1}$) was most likely tidally ejected from as close as one astronomical unit to the center of M15, confirming an IMBH ($\ge 100 M_{\odot}$ with a credibility of 98%) as the exclusive nature of the central unseen mass proposed previously.

We estimate the progenitor and explosion properties of the nearby Type II SN 2023ixf using a synthetic model grid of Type II supernova light curves. By comparing the light curves of SN 2023ixf with the pre-existing grid of Type II supernovae containing about 228,000 models with different combinations of the progenitor and explosion properties, we obtain the chi2 value for every model and evaluate the properties of the models providing small values of chi2. We found that the light-curve models with the progenitor zero-age main-sequence mass of 10 Msun, the explosion energy of (2-3)e51 erg, the 56Ni mass of 0.04-0.06 Msun, the mass-loss rate of 1e-3 - 1e-2 Msun/yr with a wind velocity of 10 km/s, and the dense, confined circumstellar matter radius of (6-10)e14 cm match well to the observed light curves of SN 2023ixf. The photospheric velocity evolution of these models is also consistent with the observed velocity evolution. Although our parameter estimation is based on a pre-existing model grid and we do not perform any additional computations, the estimated parameters are consistent with those obtained by the detailed modeling of SN 2023ixf previously reported. This result shows that comparing the pre-existing model grid is a reasonable way to obtain a rough estimate for the properties of Type II supernovae. This simple way to estimate the properties of Type II supernovae will be essential in the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) era when thousands of Type II supernovae are expected to be discovered yearly.

Due to advances in observation and imaging technologies, modern astronomical satellites generate large volumes of data. This necessitates efficient onboard data processing and high-speed data downlink. Reflecting this trend is the VERTECS 6U Astronomical Nanosatellite. Designed for the observation of Extragalactic Background Light (EBL), this mission is expected to generate a substantial amount of image data, particularly within the confines of CubeSat capabilities. This paper introduces the VERTECS Camera Control Board (CCB), an open-source payload interface board leveraging Commercial Off-The-Shelf (COTS) components, with a Raspberry Pi Compute Module 4 at its core. The VERTECS CCB hardware and software have been designed from the ground up to serve as the sole interface between the VERTECS bus system and astronomical imaging payload, while providing compute capability not usually seen in nanosatellites of this class. Responsible for mission data processing, it will facilitate high-speed data transfer from the imaging payload via gigabit Ethernet, while also providing a high-bitrate serial connection to the payload X-band transmitter for mission data downlink. Additional interfaces for secondary payloads are provided via USB-C and standard 15-pin camera connectors. The Raspberry Pi embedded within the VERTECS CCB operates on a standard Linux distribution, streamlining the software development process. Beyond addressing the current mission's payload control and data handling requirements, the CCB sets the stage for future missions with heightened data demands. Furthermore, it supports the adoption of machine learning and other compute-intensive applications in orbit. This paper delves into the development of the VERTECS CCB, offering insights into the design and validation of this next-generation payload interface, to ensure that it can survive the rigors of space flight.

In modeling a relativistic disk around a compact object, the self-gravity of the disk is often neglected while it needs to be incorporated for more accurate descriptions in several circumstances. Extending the Komatsu-Eriguchi-Hachisu self-consistent field method, we present numerical models of a rapidly rotating neutron star with a self-gravitating disk in stationary equilibrium. In particular, our approach allows us to obtain numerical solutions involving a massive disk with the rest mass $O(10^{-1})-O(10^0) M_\odot$ closely attached to a rotating neutron star. We also assess the impact of self-gravity on the internal structure of the disk and the neutron star. These axisymmetric, stationary solutions can be employed for simulations involving the neutron star-disk system in the context of high-energy transients and gravitational wave emissions.

The China Space Station Telescope (CSST) is a two-meter space telescope with multiple back-end instruments. The Fine Guidance Sensor (FGS) is an essential subsystem of the CSST Precision Image Stability System to ensure the required absolute pointing accuracy and line-of-sight stabilization. In this study, we construct the Main Guide Star Catalog for FGS. To accomplish this, we utilize the information about the FGS and object information from the Gaia Data Release 3. We provide an FGS instrument magnitude and exclude variables, binaries, and high proper motion stars from the catalog to ensure uniform FGS guidance capabilities. Subsequently, we generate a HEALPix index, which provides a hierarchical tessellation of the celestial sphere, and employ the Voronoi algorithm to achieve a homogeneous distribution of stars across the catalog. This distribution ensures adequate coverage and sampling of the sky. The performance of the CSST guide star catalog was assessed by simulating the field of view of the FGS according to the CSST mock survey strategy catalog. The analysis of the results indicates that this catalog provides adequate coverage and accuracy. The catalog's performance meets the FGS requirements, ensuring the functioning of the FGS and its guidance capabilities.

Recent studies have increasingly recognized the value of analyzing the post-reconstruction galaxy power spectrum for investigations into redshift space distortion (RSD) effects. In this paper, we present a novel theoretical model for the post-reconstruction galaxy power spectrum, designed for RSD analyses. In particular, we emphasize the importance of accounting for discrete effects arising from the reconstruction displacement vector, which have been overlooked. We specifically calculate these discrete effects, i.e., shot noise terms, within the framework of standard perturbation theory at the 1-loop level. In addition, we adopt a formulation that accounts for infrared (IR) effects to accurately model the non-linear damping of the Baryon Acoustic Oscillation (BAO) signal. Our model comprehensively integrates key physical phenomena relevant to the post-reconstruction galaxy power spectrum, such as gravitational non-linearities, RSD effects, bias effects, reconstruction effects, reconstruction-specific shot-noise effects, and the non-linear damping of the reconstructed BAO signal, thereby making it applicable to post-reconstruction RSD analyses.

To date, rings are found around a Centaur (10199) Chariklo, trans-neptunian objects (TNOs) (136108) Haumea, and (50000) Quaoar. These discoveries suggest that asteroidal ring systems may be common, particularly in the outer solar system. Since collisions are a ubiquitous and fundamental evolutionary process throughout the solar system, we conjecture that asteroidal ring systems must have experienced close encounters with small objects as part of their evolutionary process. Here, we investigate the response of ring systems when they experience gravitational disturbance by a close encounter with another small object, by calculating the change in eccentricity and the fraction of lost ring particles. We find that a perturber needs to be as massive as or more massive than the ringed object, and needs to pass in the immediate vicinity of the ring in order to cause significant disruption. The change in eccentricity expected for Chariklo's inner ring and Quaoar's outer ring agrees with the analytical expression derived from the impulse approximation, while that for Haumea's ring agrees with the analytical expression called "exponential regime". If we define a lifetime of a ring as the mean time to experience disruptive close encounters that can raise the eccentricity of ring particles greater than 0.1, the lifetime for ring systems around Chariklo, Haumea, and Quaoar are >10^4 Gyr. We conclude that ring systems around Chariklo, Haumea, and Quaoar are highly unlikely to suffer from close encounters with another small object even if those systems are as old as 4 Gyr. We also find that the lifetime is shorter for smaller ringed objects, and it still exceeds 4 Gyr for km-sized ringed objects in the outer solar system. Therefore, regardless of the size of ringed objects, asteroidal ring systems in the outer solar system are unlikely to suffer severe damage by close encounter with a small object.

Energy equipartition and the energy budget in the jet are import issues for the radiation mechanism of blazars. Early work predominantly concentrated on flat-spectrum radio quasars and a limited number of BL Lacertae objects (BL Lacs). In this paper, we compile 348 high-frequency peaked BL Lac objects (HBLs) based on the catalog of active galactic nuclei (4LAC-DR3) from Fermi-LAT, and employ \textit{JetSet} to fit the spectral energy distributions (SEDs) of these HBLs in the framework of the one-zone lepton model. We aim to determine whether the energy budget is reasonable and whether the energy equipartition is satisfied in HBLs. The results of the statistical analysis suggest that: (1) SEDs of HBLs can be reproduced well by using the one-zone lepton model; however it cannot achieve the energy equalization, and the relativistic electron energy density is far greater than the magnetic field energy density, $U_{e} \gtrsim100 U_{B}$; (2) the majority of the HBLs are located in the $t_{cool}$$<$$t_{dyn}$ region (where the horizontal coordinate represents the jet power of electrons, while the ordinate indicates the ratio between the dynamic time scale to the cooling timescale), and the jet kinetic power of HBLs is greater than the jet power of radiation; there is a very low radiation efficiency, we deduce that HBLs may have optically thin advection-dominated accretion flows; (3) the $\log\epsilon_{B}$ of HBLs is less than zero, which indicates that the jet kinetic power of HBLs is not affected by Poynting flux; (4) the relationships with $U_{e} >U_{Syn}\sim U_{B}$, $L_{e}\sim L_{p}>L_{B}\sim L_{rad}$ and $\log\epsilon_{e}>0.5$ are established. These relations indicate that most of the energy of HBLs is stored in the population of low-energy electrons.

Understanding the impact of the astrophysical environment on Type Ia supernova (SN Ia) properties is crucial to minimize systematic uncertainties in cosmological analyses based on this probe. We investigate the dependence of the SN Ia SALT2.4 light-curve stretch on the distance from their nearest galaxy cluster to study a potential effect of the intracluster medium (ICM) environment on SN Ia intrinsic properties. We use the largest SN Ia sample to date and cross-match it with existing X-ray, Sunyaev-Zel'dovich, and optical cluster catalogs in order to study the dependence between stretch and distance to the nearest detected cluster from each SN Ia. We model the underlying stretch distribution with a Gaussian mixture with relative amplitudes that depend on redshift and cluster-centric distance. We find a significant improvement of the fit quality of the stretch distribution if we include the distance-dependant term in the model with a variation of the Akaike information criterion $\rm{\Delta AIC} = -10.2$. Because of the known correlation between galaxy age and distance from cluster center, this supports previous evidence that the age of the stellar population is the underlying driver of the bimodial shape of the SN Ia stretch distribution. We further compute the evolution of the fraction of quenched galaxies as a function of distance with respect to cluster center from our best-fit model of the SNe Ia stretch distribution and compare it to previous results obtained from $H\alpha$ line measurements, optical broadband photometry, and simulations. We find our estimate to be compatible with these results. The results of this work indicate that SNe Ia searches at high redshift targeted towards clusters to maximize detection probability should be considered with caution as the stretch distribution of the detected sample would be strongly biased towards the old sub-population of SNe Ia.

Aims: The aim of the present study is to explore how to disentangle energy-dependent time delays due to a possible Lorentz invariance violation (LIV) at Planck scale from intrinsic delays expected in standard blazar flares. Methods: We first characterise intrinsic time delays in BL Lacs and Flat Spectrum Radio Quasars in standard one-zone time-dependent synchrotron self-Compton or external Compton models, during flares produced by particle acceleration and cooling processes. We simulate families of flares with both intrinsic and external LIV-induced energy-dependent delays. Discrimination between intrinsic and LIV delays is then investigated in two different ways. A technique based on Euclidean distance calculation between delays obtained in the synchrotron and in the inverse-Compton spectral bumps is used to assess their degree of correlation. A complementary study is performed using spectral hardness versus intensity diagrams in both energy ranges. Results: We show that the presence of non-negligible LIV effects, which essentially act only at very high energies (VHE), can drastically reduce the strong correlation expected between the X-ray and the VHE gamma-ray emission in leptonic scenarios. The LIV phenomenon can then be hinted at measuring the Euclidean distance $d_{E}$ from simultaneous X-ray and gamma-ray flare monitoring. Large values of minimal distance $d_{E,min}$ would directly indicate the influence of non-intrinsic time delays possibly due to LIV in SSC flares. LIV effects can also significantly modify the VHE hysteresis patterns in hardness-intensity diagrams and even change their direction of rotation as compared to the X-ray behaviour. Both observables could be used to discriminate between LIV and intrinsic delays, provided high quality flare observations are available.

Recent advances in numerical simulations of magnetically arrested accretion onto supermassive black holes have shed light on the formation and dynamics of magnetospheric current sheets near the black hole horizon. By considering the pair magnetization $\sigma_{\rm e}$ in the upstream region and the mass accretion rate $\dot{m}$ (in units of the Eddington mass accretion rate) as free parameters we estimate the strength of the magnetic field and develop analytical models, motivated by recent three-dimensional particle-in-cell simulations, to describe the populations of relativistic electron and positrons (pairs) in the reconnection region. Applying our model to M87*, we numerically compute the non-thermal photon spectra for various values of $\sigma_e$. We show that pairs that are accelerated up to the synchrotron radiation-limited energy while meandering across both sides of the current sheet, can produce MeV flares with luminosity of $\sim 10^{41}$~erg s$^{-1}$ -- independent of $\sigma_e$ -- for a black hole accreting at $\dot{m}=10^{-5}$. Pairs that are trapped in the transient current sheet can produce X-ray counterparts to the MeV flares, lasting about a day for current sheets with length of a few gravitational radii. We also show that the upstream plasma can be enriched due to photon-photon pair creation, and derive a new equilibrium magnetization of $\sigma_e \sim 10^3-10^4$ for $\dot{m}=10^{-6} - 10^{-5}$. Additionally, we explore the potential of magnetospheric current sheets to accelerate protons to ultra-high energies, finding that while acceleration to such energies is limited by various loss mechanisms, such as synchrotron and photopion losses from the non-thermal emission from pairs, maximal proton energies in the range of a few EeV are attainable in magnetospheric sheets forming around supermassive sub-Eddington accreting black holes.

The brightest Gamma-ray burst (GRB) ever, GRB 221009A, displays ultra-long GRB (ULGRB) characteristics, with a prompt emission duration exceeding 1000 s. To constrain the origin and central engine of this unique burst, we analyze its prompt and afterglow characteristics and compare them to the established set of similar GRBs. To achieve this, we statistically examine a nearly complete sample of Swift-detected GRBs with measured redshifts. Categorizing the sample to Bronze, Silver, and Gold by fitting a Gaussian function to the log-normal of T$_{90}$ duration distribution and considering three sub-samples respectively to 1, 2, and 3 times of the standard deviation to the mean value. GRB 221009A falls into the Gold sub-sample. Our analysis of prompt emission and afterglow characteristics aims to identify trends between the three burst groups. Notably, the Gold sub-sample (a higher likelihood of being ULGRB candidates) suggests a collapsar scenario with a hyper-accreting black hole as a potential central engine, while a few GRBs (GRB 060218, GRB 091024A, and GRB 100316D) in our Gold sub-sample favor a magnetar. Late-time near-IR (NIR) observations from 3.6m Devasthal Optical Telescope (DOT) rule out the presence of any bright supernova associated with GRB 221009A in the Gold sub-sample. To further constrain the physical properties of ULGRB progenitors, we employ the tool MESA to simulate the evolution of low-metallicity massive stars with different initial rotations. The outcomes suggest that rotating ($\Omega \geq 0.2\,\Omega_{\rm c}$) massive stars could potentially be the progenitors of ULGRBs within the considered parameters and initial inputs to MESA.

In the disk-corona model, the X-ray spectrum of a stellar-mass black hole in an X-ray binary is characterized by three components: a thermal component from a thin and cold accretion disk, a Comptonized component from a hot corona, and a reflection component produced by illumination of the cold disk by the hot corona. In this paper, we present ZIJI, which is a model to calculate the X-ray spectrum of black hole X-ray binaries in the disk-corona model. The accretion disk is described by the Novikov-Thorne model. The reflection spectrum is produced by the direct radiation from the corona as well as by the returning radiation of the thermal and reflection components and is calculated considering the actual spectrum illuminating the disk. If we turn the corona off, the reflection spectrum is completely generated by the returning radiation of the thermal component, as it may happen for some sources in soft spectral states. The user decides the radial density profile of the accretion disk and the ionization parameter is calculated self-consistently at any radial coordinate of the disk from the illuminating X-ray flux and the local electron density. We show the predictions of ZIJI in different regimes. We discuss current limitations of our model as well as the next steps to improve it.

Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. Whereas the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius and $\gamma$-ray emission is believed to be produced in the striped wind, outside the light-cylinder, their non-thermal X-ray production site remains unknown. The aim of this letter is to localize the non-thermal X-ray emission region based on multi-wavelength pulse profile fitting for PSR J2229+6114, a particularly good candidate due to its high X-ray brightness. Based on the geometry deduced from the joint radio and $\gamma$-ray pulse profiles, we fix the magnetic axis inclination angle and the line of sight inclination angle but we leave the region of X-ray emission unlocalised, somewhere between the surface and the light-cylinder. We localize this region and its extension by fitting the X-ray pulse profile as observed by the NICER, NuSTAR and RXTE telescopes in the ranges 2-7 keV, 3-10 keV and 9.4-22.4 keV, respectively. We constrain the non-thermal X-ray emission to arise from altitudes between $0.2\,r_L$ and $0.55\,r_L$ where $r_L$ is the light cylinder radius. The magnetic obliquity is approximately $\alpha \approx 45°-50°$ and the line of sight inclination angle $\zeta \approx 32°-48°$. This letter is among the first works to tightly constrain the location of the non-thermal X-ray emission from pulsars. We plan to apply this procedure to several other good candidates to confirm this new result.

Observations find that some fast radio bursts (FRBs) have extremely narrow-band spectra, i.e., $\Delta\nu/\nu_0 \ll 1$. We show that when the angular size of the emission region is larger than the Doppler beaming angle, the observed spectral width ($\Delta\nu/\nu_0$) exceeds 0.58 due to the {\it high latitude effects} for a source outside the magnetosphere, even when the spectrum in the source's comoving frame is monochromatic. The angular size of the source for magnetospheric models of FRBs can be smaller than the Doppler beaming angle, in which case this geometric effect does not influence the observed bandwidth. We discuss various propagation effects to determine if any could transform a broad-spectrum radio pulse into a narrow-spectrum signal at the observer's location. We find that plasma lensing and scintillation can result in a narrow bandwidth in the observed spectrum. However, the likelihood of these phenomena being responsible for the narrow observed spectra with $\Delta\nu/\nu_0 < 0.58$ in the fairly large observed sample of FRBs is exceedingly small.

A powerful tool to probe the gas content at high redshift is the [C II] 158 $\mu$m sub-millimeter emission line, which, due to its low excitation potential and luminous emission, is considered a possible direct tracer of star forming gas. In this work we investigate the origin, evolution and environmental dependencies of [C II] 158 $\mu$m emission line, as well as its expected correlation with stellar mass and star formation activity of the high-redshift galaxies observed by JWST. We use a set of state-of-the-art cold-gas hydrodynamic simulations (ColdSIM) with fully coupled time-dependent atomic and molecular non-equilibrium chemistry and self-consistent [C II] emission from metal enriched gas. We accurately track the evolution of H I, H II and $H_2$ in a cosmological context and predict both global and galaxy-based [C II] properties. For the first time, we predict the cosmic mass density evolution of [C II] and find that it is in good agreement with new measurements at redshift z = 6 from high-resolution optical quasar spectroscopy. We find a correlation between [C II] luminosity, $L_{[C II]}$, and stellar mass, consistent with results from ALMA high-redshift large programs. We predict a redshift evolution in the relation between $L_{[C II]}$ and the star formation rate, SFR, and provide a fit to relate $L_{[C II]}$ to SFR which can be adopted as a more accurate alternative to the currently used linear relation. Our findings provide physical grounds to interpret high-redshift detections in contemporary and future observations, such as the ones performed by ALMA and JWST, and to advance our knowledge on structure formation at early times.

Moderate resolution spectra of four globular clusters in the dwarf spheroidal galaxy IKN obtained with the 6-m telescope of the Special Astrophysical Observatory have been used to determine the radial velocities, ages, and metallicities of the clusters, and also to derive the first approximate estimates of the abundances of Mg, Ca, and C. Cross-correlation with radial-velocity standards, fitting of the observed spectra with model spectra, diagnostic diagrams based on the Lick absorption indices, and comparison of the spectra and absorption indices with those of Galactic globular clusters are applied. The integrated-light (IL) spectrum of the two bright clusters IKN4 and IKN5, which are close to the center of the galaxy in projection on the celestial sphere, yields the heliocentric radial velocity $V_h= 38\pm30$ km/s, age $T=12.6\pm2$ Gyr, metallicity $[Fe/H]=-2.1\pm0.2$ dex and abundance of $\alpha$-process elements $[\alpha/Fe] \sim 0.5$ dex. The IL spectrum of the two weaker clusters IKN1 and IKN3, which are far from the center of the galaxy, yields the radial velocity $V_h=-39 \pm 50$ km/s. Despite of the low signal to noise ratio in the summary spectrum of IKN1 and IKN3, one can conclude from the comparison of the results of different used methods, that IKN1 and IKN3 have seemingly the same age and metallicity as IKN4 and IKN5. According to the measured Lick indices H$_{\delta_{\rm F}}$ and H$_{\beta}$, the studied globular clusters in IKN have blue horizontal branches.

Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observations of the \textsuperscript{12}CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionised-gas emission lines from SITELLE at the Canada-France-Hawaii Telescope. Although constrained by the resolution of the ALMA observations ($\approx37$~pc), we identify only $52$ GMCs with radii ranging from $30$ to $140$~pc, a low mean molecular gas mass surface density $\langle\Sigma_{\rm gas}\rangle\approx125$~M$_\odot$~pc$^{-2}$ and a high mean virial parameter $\langle\alpha_{\rm obs,vir}\rangle\approx5.3$. We measure spatially-resolved molecular gas depletion times ($\tau_{\rm dep}\equiv1/{\rm SFE}$) with a spatial resolution of $\approx100$~pc within a galactocentric distance of $1.5$~kpc. The global depletion time is $\approx2.0$~Gyr but $\tau_{\rm dep}$ increases toward the galaxy centre, with a maximum $\tau_{\rm dep,max}\approx5.2$~Gyr. However, no pure \ion{H}{II} region is identified in NGC~524 using ionised-gas emission-line ratio diagnostics, so the $\tau_{\rm dep}$ inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC~524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC~524.

We measure the radius-velocity phase-space edge profile for Abell 1063 using galaxy redshifts from arXiv:1409.3507 and arXiv:2109.03305. Combined with a cosmological model and after accounting for interlopers and sampling effects, we infer the escape velocity profile. Using the Poisson equation, we then directly constrain the gravitational potential profile and find excellent agreement between three different density models. For the NFW profile, we find log$_{10}$(M$_{200},{\rm crit}$)= $15.40^{+0.06}_{-0.12}$M$_{\odot}$, consistent to within $1\sigma$ of six recently published lensing masses. We argue that this consistency is due to the fact that the escape technique shares no common systematics with lensing other than radial binning. These masses are 2-4$\sigma$ lower than estimates using X-ray data, in addition to earlier velocity dispersion estimates. We measure the 1D velocity dispersion within r$_{200}$ to be $\sigma_{v} = 1477^{+87}_{-99}$ km/s, which combined with our escape velocity mass, brings the dispersion for A1063 in-line with hydrodynamic cosmological simulations for the first time.

Mergers of compact objects (binary neutron stars, BNS, or neutron star-black hole, NSBH) with a substantial mass ratio ($q>1.5$) are expected to produce a mildly relativistic ejecta within $\sim20^\circ$ from the orbital plane. We present a semi-analytic approach to calculate the expected synchrotron emission observed from various viewing angles, along with the corresponding radio maps, that are produced by a collisionless shock driven by such ejecta into the interstellar medium. This method is supported by 2D numerical calculations of the full relativistic hydrodynamics with an agreement in the observed emission of $\sim30\%$. We consider a toroid ejecta with an opening angle of $15^\circ\leq\theta_ \text{open}\leq30^\circ$ and broken power-law mass distribution, $M(>\gamma\beta)\propto(\gamma\beta)^{-s}$ with $s=s_{\rm KN}$ at $\gamma\beta<\gamma_0\beta_0$ and $s=s_{\rm ft}$ at $\gamma\beta>\gamma_0\beta_0$ (where $\gamma$ is the Lorentz factor). While the peak flux is dimmer by a factor of $\sim$2-3, and the peak time remains roughly the same (within $20\%$), for various viewing angles compared to isotropic equivalent ejecta ($\theta_\text{open}=90^\circ$) considered in preceding papers, the radio maps are significantly different from the spherical case. The semi-analytic method can provide information on the ejecta geometry and viewing angle from future radio map observations and, consequently, constrain the ejection mechanism. Due to the larger ejected mass in NSBH merger, this late non-thermal signal can be observed to distances of $\lesssim 200$Mpc for typical parameter values.

Aims: We address two issues for the adoption of convex optimization in radio interferometric imaging. First, a method for a fine resolution setup is proposed which scales naturally in terms of memory usage and reconstruction speed. Second, a new tool to localize a region of uncertainty is developed, paving the way for quantitative imaging in radio interferometry. Methods: The classical $\ell_1$ penalty is used to turn the inverse problem into a sparsity-promoting optimization. For efficient implementation, the so-called Frank-Wolfe algorithm is used together with a \textit{polyatomic} refinement. The algorithm naturally produces sparse images at each iteration, leveraged to reduce memory and computational requirements. In that regard, PolyCLEAN reproduces the numerical behavior of CLEAN while guaranteeing that it solves the minimization problem of interest. Additionally, we introduce the dual certificate image, which appears as a numerical byproduct of the Frank-Wolfe algorithm. This image is proposed as a tool for uncertainty quantification on the location of the recovered sources. Results: PolyCLEAN demonstrates good scalability performance, in particular for fine-resolution grids. On simulations, the Python-based implementation is competitive with the fast numerically-optimized CLEAN solver. This acceleration does not affect image reconstruction quality: PolyCLEAN images are consistent with CLEAN-obtained ones for both point sources and diffuse emission recovery. We also highlight PolyCLEAN reconstruction capabilities on observed radio measurements. Conclusions: PolyCLEAN can be considered as an alternative to CLEAN in the radio interferometric imaging pipeline, as it enables the use of Bayesian priors without impacting the scalability and numerical performance of the imaging method.

Asteroseismology has been shown to be, together with stellar modelling, an invaluable tool in constraining properties of novel physics. In this work, we study for the first time the influence of axionic production in the evolution of a late main-sequence star, comparing computational models with observational data in order to constrain the axion-photon $g_{a\gamma}$ coupling parameter. We first perform a high-precision calibration of a stellar model to our target star, in order to obtain a benchmark for our other diagnostics. We then apply a two-stage test, first using global quantities and then resorting to precision seismic ratios. We find that seismology allows us to place an independent upper bound of $g_{a\gamma} \leq 0.98\times 10^{-10} $ GeV$^{-1}$ at a $68\%$ confidence level (CL), in the same order of magnitude as both the most recent constraints from the observation of globular clusters and previous bounds obtained through stellar modelling, but more stringent than most current direct axion detections. We also suggest a more conservative limit of $g_{a\gamma} \leq 1.38\times 10^{-10} $ GeV$^{-1}$ at a $95\%$ CL. Moreover, this new diagnostic method can be applied to stellar data that will be obtained in future asteroseismic projects.

In order to predict the black hole mass distributions at high redshift, we need to understand whether very massive single stars ($M>40$ M$_\odot$) at low metallicity $Z$ lose their hydrogen-rich envelopes, like their metal-rich counterparts, or whether a binary companion is required to achieve this. To test this, we undertake a deep spectroscopic search for binary companions of the seven apparently single Wolf-Rayet (WR) stars in the Small Magellanic Cloud (SMC; $Z \simeq 1/5 Z_\odot$). For each of them, we acquired six high-quality VLT-UVES spectra spread over 1.5 years. By using the narrow N V lines in these spectra, we monitor radial velocity (RV) variations to search for binary motion. We find low RV variations between 6 and 23 km/s for the seven WR stars, with a median standard deviation of $5$ km/s. Our Monte Carlo simulations imply probabilities below ~5% for any of our target WR stars to have a binary companion more massive than ~5 M$_\odot$ at orbital periods of less than a year. We estimate that the probability that all our target WR stars have companions with orbital periods shorter than 10 yr is below ~10$^{-5}$, and argue that the observed modest RV variations may originate from intrinsic atmosphere or wind variability. Our findings imply that metal-poor massive stars born with $M \gtrsim 40$ M$_\odot$ can lose most of their hydrogen-rich envelopes via stellar winds or eruptive mass loss, which strongly constrains their initial mass - black hole mass relation. We also identify two of our seven target stars (SMC AB1 and SMC AB11) as runaway stars with a peculiar radial velocity of ~80 km/s. Moreover, with all five previously detected WR binaries in the SMC exhibiting orbital periods of below 20 d, a puzzling absence of intermediate-to-long-period WR binaries has emerged, with strong implications for the outcome of massive binary interaction at low metallicity.

Understanding Type Ia supernovae (SNe~Ia) and the empirical standardisation relations that make them excellent distance indicators is vital to improving cosmological constraints. SN~Ia ``siblings", i.e. two or more SNe~Ia in the same host or parent galaxy offer a unique way to infer the standardisation relations and their diversity across the population. We analyse a sample of 25 SN~Ia pairs, observed homogeneously by the Zwicky Transient Factory (ZTF) to infer the SNe~Ia light curve width-luminosity and colour-luminosity parameters $\alpha$ and $\beta$. Using the pairwise constraints from siblings, allowing for a diversity in the standardisation relations, we find $\alpha = 0.218 \pm 0.055 $ and $\beta = 3.084 \pm 0.312$, respectively, with a dispersion in $\alpha$ and $\beta$ of $\leq 0.195$ and $\leq 0.923$, respectively, at 95$\%$ C.L. While the median dispersion is large, the values within $\sim 1 \sigma$ are consistent with no dispersion. Hence, fitting for a single global standardisation relation, we find $\alpha = 0.228 \pm 0.029 $ and $\beta = 3.160 \pm 0.191$. We find a very small intrinsic scatter of the siblings sample $\sigma_{\rm int} \leq 0.10$ at 95\% C.L. compared to $\sigma_{\rm int} = 0.22 \pm 0.04$ when computing the scatter using the Hubble residuals without comparing them as siblings. Splitting the sample based on host galaxy stellar mass, we find that SNe~Ia in both subsamples have consistent $\alpha$ and $\beta$. The $\beta$ value is consistent with the value for the cosmological sample. However, we find a higher $\alpha$ by $\sim 2.5 - 3.5 \sigma$. The high $\alpha$ is driven by low $x_1$ pairs, potentially suggesting that the slow and fast declining SN~Ia have different slopes of the width-luminosity relation. We can confirm or refute this with increased statistics from near future time-domain surveys. (abridged)

We present a study of the radio variability of bright, $S_{1.4}\geq100$ mJy, high-redshift quasars at $z\geq3$ on timescales up to 30-40 years. The study involved simultaneous RATAN-600 measurements at frequencies of 2.3, 4.7, 8.2, 11.2, and 22.3 GHz in 2017-2020. In addition, data from the literature were used. We have found that the variability index, $V_S$, which quantifies the normalized difference between the maximum and minimum flux density while accounting for measurement uncertainties, ranges from 0.02 to 0.96 for the quasars. Approximately half of the objects in the sample exhibit a variability index within the range of 0.25 to 0.50, comparable to that observed in blazars at lower redshifts. The distribution of $V_S$ at 22.3 GHz is significantly different from that at 2.3-11.2 GHz, which may be attributed to the fact that a compact AGN core dominates at the source's rest frame frequencies greater than 45 GHz, leading to higher variability indices obtained at 22.3 GHz (the $V_S$ distribution peaks around 0.4) compared to the lower frequencies (the $V_S$ distribution at 2.3 and 4.7 GHz peaks around 0.1-0.2). Several source groups with distinctive variability characteristics were found using cluster analysis of quasars. We propose 7 new candidates for gigahertz peaked-spectrum (GPS) sources and 5 new megahertz peaked-spectrum (MPS) sources based on their spectrum shape and variability features. Only 6 out of 23 sources previously reported as GPS demonstrate a low variability level typical of classical GPS sources ($V_{S} < 0.25$) at 4.7-22.3 GHz. When excluding the highly variable peaked-spectrum blazars, we expect no more than 20% of the sources in the sample to be GPS candidates and no more than 10% to be MPS candidates.

In this paper, we investigate scalar-induced gravitational waves (GWs) generated in the post-inflationary universe to infer new limits on warm inflation. We specifically examine the evolution of primordial GWs produced by scalar perturbations produced during the radiation-dominated epoch. For this purpose, we assume a weak regime of warm inflation under the slow-roll approximation, with a dissipation coefficient linearly dependent on the temperature of the radiation bath. We then derive analytical expressions for the curvature power spectrum and the scalar index, in the cases of chaotic and exponential potentials of the inflationary field. Subsequently, we compare the theoretical predictions regarding the relic energy density of GWs with the stochastic GW background signal recently detected by the NANOGrav collaboration through the use of pulsar timing array measurements. In so doing, we obtain numerical constraints on the free parameters of the inflationary models under study. Finally, we conduct a model selection analysis through the Bayesian inference method to measure the statistical performance of the different theoretical scenarios.

We derive new constraints on the CIB monopole spectrum from reprocessed COBE-DIRBE sky maps with lower instrumental and astrophysical contamination than the legacy DIRBE maps. These maps have been generated through a global Bayesian analysis framework that simultaneously fits cosmological, astrophysical, and instrumental parameters, described in a series of papers referred to as Cosmoglobe DR2. We have applied this method to the DIRBE Calibrated Individual Observations, complemented by selected HFI and FIRAS sky maps to break key astrophysical degeneracies, as well as WISE and Gaia compact object catalogs. In this paper, we focus on CIB constraints that result from this work. We report detections of an isotropic signal in five out of the ten DIRBE bands (1.25, 2.2, 3.5, 140, and 240 $\mu$m). For the 2.2 $\mu$m channel, we find an amplitude of $13\pm3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 36% lower than reported from the official DIRBE release. For the 240 $\mu$m channel, we find $10\pm3\,\mathrm{nW\,m^{-2}\,sr^{-1}}$, 26% lower than the official DIRBE release. We interpret these lower values as resulting from improved zodiacal light and Galactic foreground modeling. For the bands between 4.9 and 100 $\mu$m, the presence of excess radiation in solar-centric coordinates reported in a companion paper precludes the definition of lower limits. However, the analysis still provides well-defined upper limits. For the 12 $\mu$m channel, we find an upper 95% confidence limit of 103 $\mathrm{nW\,m^{-2}\,sr^{-1}}$, a factor of four lower than the corresponding legacy result of 468 $\mathrm{nW\,m^{-2}\,sr^{-1}}$. The results presented in this paper redefines the state-of-the-art CIB monopole constraints from COBE-DIRBE, and it provides a real-world illustration of the power of global end-to-end analysis of multiple complementary data sets which is the foundational idea of the Cosmoglobe project.

The TEMPO (Transiting Exosatellites, Moons, and Planets in Orion) Survey is a proposed 30-day observational campaign using the Nancy Grace Roman Space Telescope. By providing deep, high-resolution, short-cadence infrared photometry of a dynamic star-forming region, TEMPO will investigate the demographics of exosatellites orbiting free-floating planets and brown dwarfs -- a largely unexplored discovery space. Here, we present the simulated detection yields of three populations: extrasolar moon analogs orbiting free-floating planets, exosatellites orbiting brown dwarfs, and exoplanets orbiting young stars. Additionally, we outline a comprehensive range of anticipated scientific outcomes accompanying such a survey. These science drivers include: obtaining observational constraints to test prevailing theories of moon, planet, and star formation; directly detecting widely separated exoplanets orbiting young stars; investigating the variability of young stars and brown dwarfs; constraining the low-mass end of the stellar initial mass function; constructing the distribution of dust in the Orion Nebula and mapping evolution in the near-infrared extinction law; mapping emission features that trace the shocked gas in the region; constructing a dynamical map of Orion members using proper motions; and searching for extragalactic sources and transients via deep extragalactic observations reaching a limiting magnitude of $m_{AB}=29.7$\,mag (F146 filter).

Mesospheric CO2 clouds are one of two types of carbon dioxide clouds known on Mars. We present observations of mesospheric CO2 clouds made by Atmospheric Chemistry Suit (ACS) onboard the ESA-Roscosmos ExoMars Trace Gas Orbiter (TGO). We analyzed 1663 solar occultation sessions of Thermal InfraRed (TIRVIM) and Middle InfraRed (MIR) channels of ACS covering more than two Martian years that contain spectra of 2.7 {\mu}m carbon dioxide ice absorption band. That allowed us to unambiguously discriminate carbon dioxide ice aerosols from mineral dust and water ice aerosols, not relying on the information of atmospheric thermal conditions. CO2 clouds were detected in eleven solar occultation observations at altitudes from 39 km to 90 km. In five cases, there were two or three layers of CO2 clouds that were vertically separated by 5-15 km gaps. Effective radius of CO2 aerosol particles is in the range of 0.1-2.2 {\mu}m. Spectra produced by the smallest particles indicate a need for a better resolved CO2 ice refractive index. Nadir optical depth of CO2 clouds is in the range 5*10^{-4}-4*10^{-2} at both 2.7 {\mu}m and 0.8 {\mu}m. Asymmetrical diurnal distribution of detections observed by ACS is potentially due to local time variations of temperature induced by thermal tides. Two out of five cases of carbon dioxide cloud detections made by the TIRVIM instrument reveal the simultaneous presence of CO2 ice and H2O ice aerosols. Temperature profiles measured by the Near InfraRed (NIR) channel of ACS are used to calculate CO2 saturation ratio S at locations of carbon dioxide clouds. Supersaturation S > 1 is detected in only 5 out of 19 cases of CO2 cloud layers; extremely low values of S < 0.1 are found in 9 out of 19 cases.

We present a method for modelling the cluster-galaxy correlation function in redshift-space, down to ~ Mpc scales. The method builds upon the so-called Galaxy Infall Kinematics (GIK) model, a parametric model for the pairwise velocities of galaxies with respect to nearby galaxy clusters. We fit the parameters of the GIK model to a suite of simulations run with different cosmologies, and use Gaussian Processes to emulate how the GIK parameters depend upon cosmology. This emulator can then be combined with knowledge of the real-space clustering of clusters and galaxies, to predict the cluster-galaxy correlation function in redshift space. Fitting this model to an observed correlation function enables the extraction of cosmological parameter constraints, and we present forecasts for a DESI-like survey. We also perform tests of the robustness of our constraints from fitting to mock data extracted from N-body simulations, finding that fitting to scales < 3 Mpc/h leads to a biased inference on cosmology, due to model misspecification on these scales. Finally, we discuss what steps will need to be taken in order to apply our method to real data.

Recent studies of wide binary stars based on Gaia DR3 suggest that the relative orbital velocities of objects with separations s > 3'000 astronomical units are statistically larger than the standard Newtonian predictions. Obviously there is no Dark Matter halo arround binary stars that could be invoked to explain these high velocities. However, we explore the properties of two-body systems in the framework of scale invariant vacuum theory, focusing on the case of objects with extreme separations. In this regime, the additional acceleration term present in the modified Newton equation with scale invariance becomes important, and may even dominate the dynamical evolution at very low gravities. Comparisons with Gaia DR3 observations of wide binaries are performed and suggest that binaries with separations s > 3'000 astronomical units have experienced such an evolution for a few Gyr, accounting well for the observed velocity excesses.

This paper presents maps of the J=2-1 transition of CO toward the Draco Nebula Intermediate Velocity Cloud (IVC). The maps cover 8500 square arcmin with a velocity resolution of 0.33 km~s$^{-1}$ and angular resolution of 38", or 0.11 pc at the cloud distance of 600 pc. The mapped area includes all the emission detected by the {\it Herschel} satellite with 250 $\mu$m intensity >5 MJy/sr. Previously published observations of the far-IR emission and the 21 cm line of HI are used to derive the column density distribution of H$_2$ and the abundance ratio CO/H$_2$, as well as the distribution of the molecular fraction of hydrogen, which approaches 90\% over much of the brighter parts of the nebula. The CO emission is highly clumpy and closely resembles the structures seen in far-IR images. The kinematics of the CO show supersonic motions between clumps but near-thermal to trans-sonic motions within clumps, consistent with model predictions that the scale length for dissipation of supersonic turbulence should be $\sim0.1$ pc, mediated by kinematic viscosity and/or ambipolar diffusion. Different parts of the nebula show evidence for a spread of molecular formation timescales of a few 10$^5$ years, comparable to the dynamical timescale of the infalling gas. The IVC will likely merge with the Galactic interstellar medium in $\sim 10^7$ years, and the densest clumps may form an unbound cluster of low-mass stars.

The gradient technique is a promising tool with theoretical foundations based on the fundamental properties of MHD turbulence and turbulent reconnection. Its various incarnations use spectroscopic, synchrotron, and intensity data to trace the magnetic field and measure the media magnetization in terms of Alfven Mach number. We provide an analytical theory of gradient measurements and quantify the effects of averaging gradients along the line of sight and over the plane of the sky. We derive analytical expressions that relate the properties of gradient distribution with the Alfven Mach number $M_A$. We show that these measurements can be combined with measures of sonic Mach number or line broadening to obtain the magnetic field strength. The corresponding technique has advantages to Davis-Chandrasekhar-Fermi way of obtaining the magnetic field strength.

The Sun's history is still a subject of interest to modern astrophysics. Observationally constrained CME rates of young solar analogues are still lacking, as those require dedicated monitoring. We present medium resolution optical spectroscopic monitoring of a small sample of bright and prominent solar analogues over a period of three years using the 0.5m telescope at observatory Lustbühel Graz (OLG) of the University of Graz, Austria. The aim is the detection of flares and CMEs from those spectra. In more than 1700 hours of spectroscopic monitoring we found signatures of four flares and one filament eruption on EK Dra which has been reported in previous literature, but we complementarily extended the data to cover the latter phase. The other stars did not reveal detectable signatures of activity. For these non-detections we derive upper limits of occurrence rates of very massive CMEs, which are detectable with our observational setup, ranging from 0.1 to 2.2 per day , but these may be even smaller than the given rates considering observational biases. Furthermore, we investigate the detectability of flares/CMEs in OLG spectra by utilizing solar 2D H{\alpha} spectra from MEES solar observatory. We find that solar-sized events are not detectable within our observations. By scaling up the size of the solar event, we show that with a fractional active region area of 18% in residual spectra and 72% in equivalent width time series derived from the same residuals that solar events are detectable if they had hypothetically occurred on HN Peg.

The bright GRB 210610B was discovered simultaneously by Fermi and Swift missions at redshift 1.13. We utilized broadband Fermi-GBM observations to perform a detailed prompt emission spectral analysis and to understand the radiation physics of the burst. Our analysis displayed that the low energy spectral index ($\alpha_{\rm pt}$) exceeds boundaries expected from the typical synchrotron emission spectrum (-1.5,-0.67), suggesting additional emission signature. We added an additional thermal model with the typical Band or CPL function and found that CPL + BB function is better fitting to the data, suggesting a hybrid jet composition for the burst. Further, we found that the beaming corrected energy (E$_{\rm \gamma, \theta_{j}}$ = 1.06 $\times$ 10$^{51}$ erg) of the burst is less than the total energy budget of the magnetar. Additionally, the X-ray afterglow light curve of this burst exhibits achromatic plateaus, adding another layer of complexity to the explosion's behavior. Interestingly, we noted that the X-ray energy release during the plateau phase (E$_{\rm X,iso}$ = 1.94 $\times$ 10$^{51}$ erg) is also less than the total energy budget of the magnetar. Our results indicate the possibility that a magnetar could be the central engine for this burst.