Papers for Tuesday, Mar 19 2024

Recent submillimeter observations have revealed signs of pc-scale molecular inflow and atomic outflow in the nearest Seyfert 2 galaxy, the Circinus galaxy. To verify the gas kinematics suggested by these observations, we performed molecular and atomic line transfer calculations based on a physics-based 3D radiation-hydrodynamic model, which has been compared with multi-wavelength observations in this paper series. The major axis position-velocity diagram (PVD) of CO(3-2) reproduces the observed faint emission at the systemic velocity, and our calculations confirm that this component originates from failed winds falling back to the disk plane. The minor-axis PVD of [CI]($^3P_1$-$^3P_0$), when created using only the gas with positive radial velocities, presents a sign of blue- and redshifted offset peaks similar to those in the observation, suggesting that the observed peaks indeed originate from the outflow, but that the model may lack outflows as strong as those in the Circinus galaxy. Similar to the observed HCN(3-2), the similar dense gas tracer HCO$^+$(3-2) can exhibit nuclear spectra with inverse P-Cygni profiles with $\sim$0.5 pc beams, but the line shape is azimuthally dependent. The corresponding continuum absorbers are inflowing clumps at 5-10 pc from the center. To detect significant absorption with a high probability, the inclination must be fairly edge-on ($\gtrsim$85$^\circ$), and the beam size must be small ($\lesssim$1 pc). These results suggest that HCN or HCO$^+$ and [CI] lines are effective for observing pc-scale inflows and outflows, respectively.

We present spot properties on 32 periodic young stellar objects in IC 5070. Long term, $\sim$5 yr, light curves in the $V$, $R$, and $I$-bands are obtained through the HOYS (Hunting Outbursting Young Stars) citizen science project. These are dissected into six months long slices, with 3 months oversampling, to measure 234 sets of amplitudes in all filters. We fit 180 of these with reliable spot solutions. Two thirds of spot solutions are cold spots, the lowest is 2150 K below the stellar temperature. One third are warm spots that are above the stellar temperature by less than $\sim$2000 K. Cold and warm spots have maximum surface coverage values of 40 percent, although only 16 percent of warm spots are above 20 percent surface coverage as opposed to 60 percent of the cold spots. Warm spots are most likely caused by a combination of plages and low density accretion columns, most common on objects without inner disc excess emission in $K-W2$. Five small hot spot solutions have $<3$ percent coverage and are 3000 - 5000 K above the stellar temperature. These are attributed to accretion, and four of them occur on the same object. The majority of our objects are likely to be accreting. However, we observe very few accretion hot spots as either the accretion is not stable on our timescale or the photometry is dominated by other features. We do not identify cyclical spot behaviour on the targets. We additionally identify and discuss a number of objects that have interesting amplitudes, phase changes, or spot properties.

The chemical composition of globular clusters (GCs) across the Local Group provides information on chemical abundance trends. The host galaxy, Sextans A, is a low-surface-brightness dwarf irregular galaxy located on the edge of the Local Group.We derive the dynamical mass of the GC together with the mass-to-light ratio and the abundances of the ${\alpha}$ (Mg, Ca, Ti, Si), Fe-peak (Fe, Cr, Mn, Sc, Ni), and heavy elements (Ba, Cu, Zn, Eu). Abundance ratios were determined from the analysis of an optical integrated-light spectrum of Sextans A GC-1, obtained with UVES on the VLT. We apply non-local thermodynamic equilibrium (NLTE) corrections to Mg, Ca, Ti, Fe, and Ni. The GC appears to be younger and more metal-poor than the majority of the GCs of the Milky Way, with an age of 8.6$\pm$2.7 Gyr and $\text{[Fe/H]}=-2.14\pm0.04$ dex. The calculated dynamical mass is $(5.18 \pm1.62) \times 10^5 M_{\odot}$, which results in an atypically high value of the mass-to-light ratio, 4.35$\pm$1.40 M$_{\odot}$/L$_{V \odot}$. Sextans A GC-1 has varying $\alpha$ elements - the Mg abundance is extremely low ($\text{[Mg/Fe]}=-0.79\pm0.29$), Ca and Ti are solar-scaled or mildly enhanced, and Si is enhanced. This makes the mean $\alpha$ abundance (excluding Mg) to be enhanced. The Fe-peak elements are consistent with scaled-solar or slightly enhanced abundances. Ba and Cu have sub-solar abundance ratios, while Zn and Eu are consistent with their upper limits being solar-scaled and enhanced. The composition of Sextans A GC-1 resembles the overall pattern and behaviour of GCs in the Local Group. The anomalous values are the mass-to-light ratio and the depleted abundance of Mg. There is no definite explanation for such an extreme abundance value. Variations in the initial mass function or the presence of an intermediate-mass black hole might explain the high mass-to-light ratio value.

Observations and high-resolution hydrodynamical simulations indicate that massive star clusters assemble hierarchically from sub-clusters with a universal power-law cluster mass function. We study the consequences of such assembly for the formation of intermediate-mass black holes (IMBHs) at low metallicities ($Z=0.01\;Z_\mathrm{\odot}$) with our updated N-body code BIFROST based on the hierarchical fourth-order forward integrator. BIFROST integrates few-body systems using secular and regularized techniques including post-Newtonian equations of motion up to order PN3.5 and gravitational-wave recoil kicks for BHs. Single stellar evolution is treated using the fast population synthesis code SEVN. We evolve three cluster assembly regions with $N_\mathrm{tot} = 1.70$--$2.35 \times 10^6$ stars following a realistic IMF in $\sim$1000 sub-clusters for $t=50$ Myr. IMBHs with masses up to $m_\bullet \sim 2200\:M_\mathrm{\odot}$ form rapidly mainly via the collapse of very massive stars (VMSs) assembled through repeated collisions of massive stars followed by growth through tidal disruption events and BH mergers. No IMBHs originate from the stars in the initially most massive clusters. We explain this by suppression of hard massive star binary formation at high velocity dispersions and the competition between core collapse and massive star life-times. Later the IMBHs form subsystems resulting in gravitational-wave BH-BH, IMBH-BH and IMBH-IMBH mergers with a $m_\bullet\sim1000\:M_\mathrm{\odot}$ gravitational-wave detection being the observable prediction. Our simulations indicate that the hierarchical formation of massive star clusters in metal poor environments naturally results in formation of potential seeds for supermassive black holes.

Even at modest amplification, the optical depth to gravitational lensing through the Galaxy is $<10^{-5}$. However, the large apparent isotropic-equivalent energy of GRB 221009A coupled with a path through low Galactic latitude suggests that the conditional probability that this particular GRB was lensed is greater than the very low a priori expectation. With the extreme brightness of the prompt emission, this Galactic lensing hypothesis can be constrained by autocorrelation analysis of Fermi photons on 0.1-1000 ms timescales. In relating lensing mass, magnification, and autocorrelation timescale, I show that a lensed-induced autocorrelation signature by stellar lenses falls below the minimal variability timescale (MVT) expected from a black hole central engine. However, lensing by Galactic dark matter MACHOs ($M_l > 10-1000\,M_\odot$) could be confirmed with this approach. Regardless, at a peak $\gamma$-ray photon rate of $>30$ ms$^{-1}$, GRB 221009A represents a prime opportunity to measure the smallest MVTs of GRBs.

We performed the optical spectroscopy of 16 classical Be stars in 11 open clusters older than 100 Myr. Ours is the first spectroscopic study of classical Be stars in open clusters older than 100 Myr. We found that the H alpha emission strength of most of the stars is less than 40 Angstrom, in agreement with previous studies. Our analysis further suggests that one of the stars, KW97 35 12, might be a weak H alpha emitter in nature, showing H alpha equivalent width of negative 0.5 Angstrom. Interestingly, we also found that the newly detected classical Be star LS III 47 37b might be a component of the possible visual binary system LS III 47 37, where the other companion is also a classical Be star. Hence, the present study indicates the possible detection of a binary Be system. Moreover, it is observed that all 16 stars exhibit a lesser number of emission lines compared to classical Be stars younger than 100 Myr. Furthermore, the spectral type distribution analysis of B type and classical Be stars for the selected clusters points out that the existence of CBe stars can depend on the spectral type distribution of B type stars present in these clusters.

Time-Delay Cosmography is a technique for masuring $H_0$ with strong gravitational lensing. It requires a correction for line of sight perturbations, and it is necessary to build tools to assess populations of these lines of sight efficiently. We aim demonstrate the techniques necessary to analyze line of sight effects at a population level, and investigate whether strong lenses fall in preferably overdense environments. We analyze a set of 25 galaxy-galaxy lens lines of sight in the Strong Lensing Legacy Survey sample using standard techniques, then perform a hierarchical analysis to constrain the population-level parameters. We introduce a new statistical model for these posteriors that may provide insight into the underlying physics of the system. We find the median value of $\kappa_{\rm{ext}}$ in the population model to be $0.033 \pm 0.010$. The median value of $\kappa_{\rm{ext}}$ for the individual lens posteriors is $0.008 \pm 0.015$. Both approaches demostrate that our systems are drawn from an overdense sample. The difference results form these two approaches show the importance of population models that do not multiply the effect of our priors.

The Fornax galaxy cluster is the richest nearby (D ~ 20 Mpc) galaxy association in the southern sky. As such, it provides a wealth of oportunities to elucidate on the processes where environment holds a key role in transforming galaxies. Although it has been the focus of many studies, Fornax has never been explored with contiguous homogeneous wide-field imaging in 12 photometric narrow- and broad-bands like those provided by the Southern Photometric Local Universe Survey (S-PLUS). In this paper we present the S-PLUS Fornax Project (S+FP) that aims to comprehensively analyse the galaxy content of the Fornax cluster using S-PLUS. Our data set consists of 106 S-PLUS wide-field frames (FoV ~ 1.4 x 1.4 deg$^2$) observed in five SDSS-like ugriz broad-bands and seven narrow-bands covering specific spectroscopic features like [OII], CaII H+K, H$\delta$, G-band, Mg b triplet, H$\alpha$, and the CaII triplet. Based on S-PLUS specific automated photometry, aimed at correctly detecting Fornax galaxies and globular clusters in S-PLUS images, our dataset provides the community with catalogues containing homogeneous 12-band photometry for ~ 3 x 10$^6$ resolved and unresolved objects within a region extending over ~ 208 deg$^2$ (~ 5 Rvir in RA) around Fornax' central galaxy, NGC 1399. We further explore the EAGLE and IllustrisTNG cosmological simulations to identify 45 Fornax-like clusters and generate mock images on all 12 S-PLUS bands of these structures down to galaxies with M$\star \geq 10^8$ M$\odot$. The S+FP dataset we put forward in this first paper of a series will enable a variety of studies some of which are briefly presented.

Active galactic nuclei (AGNs) can power relativistic jets, which are called blazars when pointed close to our line of sight. Depending on the presence or absence of emission lines in their optical spectra, blazars are categorized into flat spectrum radio quasars (FSRQs) or BL Lacertae (BL Lacs) objects. According to the 'blazar sequence', as synchrotron peak frequency ($\nu^{sy}_{pk}$) shifts to higher energies, the synchrotron peak luminosity decreases. This means that BL Lacs as luminous as FSRQs, and with synchrotron peak frequencies $\nu^{sy}_{pk}>10^{15}$ Hz, should not exist. Detected as a high-synchrotron peak (HSP; $\nu^{sy}_{pk}>10^{15}$ Hz) BL Lac, 4FGL J1520.8-0348 shows high gamma-ray luminosity ($L_{\gamma}>10^{46}\,\rm erg~s^{-1}$), being at a high redshift of $z=$1.46. Since it is an outlier in the 'blazar sequence', the process of its jet acceleration and power may be different from bona fide BL Lacs. In this work, we constrain its spectral energy distribution (SED) by modeling the multi-wavelength data from infrared to $\gamma$-ray regime. Simultaneous X-ray data was obtained from X-ray Multi-Mirror Mission and Nuclear Spectroscopic Telescope Array to constrain the synchrotron emission and underlying electron distribution. On undertaking the SED modeling of the source, including the effect of extragalactic background light, we conclude that the source is more likely to be a 'blue FSRQ' or 'masquerading BL Lac' where the BL Lac is actually a FSRQ in disguise.

In the former studies, the time evolution information is missed in deducing the time-integrated polarizations of gamma-ray burst (GRB) prompt emission. Here, it is considered and the time-integrated polarizations is investigated through the summation of the time-resolved ones. The statistical properties of the distribution of the time-integrated polarization degree ($\Pi$) can be read from the $q-\Pi$ curve, where $q\equiv\theta_V/\theta_j$. $\theta_V$ and $\theta_j$ are the observational and jet half-opening angles, respectively. Hence, only the $q-\Pi$ curves are studied. In addition to a toroidal magnetic field in the radiation region, an aligned field is also discussed. We found the predicted time-integrated PD is around $(40-50)\%$ for High-energy Polarimetry Detector (HPD) on board POLAR-2 and is roughly $(30-40)\%$ for its Low-energy Polarimetry Detector (LPD). Therefore, $\Pi$ value detected by the HPD will be larger than that of the LPD in statistics and the result of the former estimations will underestimate the value of $\Pi$ in an ordered field. There are mainly two types of the $q-\Pi$ curve profiles, corresponding to two ordered magnetic field configurations.

In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $\Lambda$CDM cosmogony with $\Omega_m=0.3111$, $\Omega_{\Lambda}=0.6889$, and $\sigma_8=0.8102$. The simulation has a boxsize of $1~h^{-1} {\rm Gpc}$, and consists of $6144^3$ particles with mass resolution $3.723 \times 10^{8} h^{-1} M_\odot $. In order to take into account the effect of redshift evolution, we first use all 128 snapshots in the Jiutian simulation to generate a light-cone halo/subhalo catalog. Next, galaxy luminosities are assigned to the main and subhalo populations using the subhalo abundance matching (SHAM) method with the DESI $z$-band luminosity functions at different redshifts. Multi-band photometries, as well as images, are then assigned to each mock galaxy using a 3-dimensional parameter space nearest neighbor sampling of the DESI LS observational galaxies and groups. Finally, the CSST and DESI LS survey geometry and magnitude limit cuts are applied to generate the required MGRSs. As we have checked, this set of MGRSs can generally reproduce the observed galaxy luminosity/mass functions within 0.1 dex for galaxies with $L > 10^8 L_\odot$ (or $M_* > 10^{8.5} M_\odot$) and within 1-$\sigma$ level for galaxies with $L < 10^8L_\odot$ (or $M_* < 10^{8.5} M_\odot$). Together with the CSST slitless spectra and redshifts for our DESI LS seed galaxies that are under construction, we will set out to test various slitless observational selection effects in subsequent probes.

We present the results of our mosaic observations of a single Class 0 protostar IRAS 15398$-$3359 with Atacama Compact Array (ACA) in the CO $J=2\mbox{-}1$ line. The new observations covering a $\sim\!2'$ square region revealed elongated redshifted and blueshifted components, which are located at distances of $\sim\!30''\mbox{-}75''$ on the northern and southern sides of the protostar, respectively, in addition to the previously observed primary and secondary outflows. These elongated components exhibit Hubble-law like velocity structures, i.e., an increase of velocity with increasing distance from the protostar, suggesting that it is the third outflow associated with the protostar. Besides, a new redshifted component is detected at radii of $\sim\!40''\mbox{-}75''$ on the northwestern side of the protostar. This redshifted component also exhibits a Hubble-law like velocity profile, which could be the counterpart of the secondary outflow mostly detected at blueshifted velocities in a previous study. The three outflows are all misaligned by $\sim\!20\mbox{-}90^\circ$, and the dynamical timescale of the primary outflow is shorter than those of the other outflows approximately by an order of magnitude. These facts hint that the outflow launch direction has significantly changed with time. The outflow direction may change if the rotational axis and the magnetic field are misaligned, or if the dense core is turbulent. We favor the second scenario as the origin of the multiple outflows in IRAS 15398$-$3359 based on a comparison between the observational results and numerical simulations.

$\Lambda$CDM is challenged by observational tensions between late- and early-time cosmology, most dramatically in the Hubble constant $H_0$. $\Lambda$CDM hereby falls short as an effective infra-red (IR) limit of quantum cosmology. Consistency with general relativity in an effective dark energy $\Lambda=\alpha_p\Lambda_0\sim H^2/c^2$ obtains from an IR-coupling $\alpha_p\sim \hbar$ to the bare cosmological constant $\Lambda_0\sim 1/\hbar$, where $\hbar$ is the Planck constant. A path integral formulation with gauged total phase identifies $\Lambda$ with the trace of the Schouten tensor $J$. It predicts the scaling $H_0\simeq \sqrt{6/5}H_0^{\Lambda{\rm CDM}}$ inferred from the BAO, while preserving the age of the Universe. With no free parameters, this first principle model predicts a 9\% departure in $H_0$ between the Local Distance Ladder and the $\Lambda$CDM {\em Planck} measurements with no tension between late and early times. In galaxies, the same IR coupling predicts a sharp transition to anomalous rotation curves below the de Sitter acceleration $a_{dS}=cH$, where $c$ is the velocity of light. Excluded by galaxy models in $\Lambda$CDM, it points to ultra-light dark matter of mass $m_D\lesssim 3\times 10^{-21}$eV consistent with $m_D\gtrsim 10^{-22}$eV of wave-like $\psi$CDM.

Context: Molecular clouds are the prime locations of star formation. These clouds contain filamentary structures and cores which are crucial in the formation of young stars. Aims: In this work, we aim to quantify the physical properties of structural characteristics within the molecular cloud L1251 to better understand the initial conditions for star formation. Methods: We applied the getsf algorithm to identify cores and filaments within the molecular cloud L1251 using the Herschel multiband dust continuum image, enabling us to measure their respective physical properties. Additionally, we utilized an enhanced differential term algorithm to produce high-resolution temperature maps and column density maps with a resolution of ${13.5}''$. Results: We identified 122 cores in the region. Out of them, 23 are protostellar cores, 13 are robust prestellar cores, 32 are candidate prestellar cores (including 13 robust prestellar cores and 19 strictly candidate prestellar cores), and 67 are unbound starless cores. getsf also found 147 filament structures in the region. Statistical analysis of the physical properties (mass (M), temperature (T), size, and core brightness (hereafter, we are using the word luminosity (L)) for the core brightness) of obtained cores shows a negative correlation between core mass and temperature and a positive correlation between (M/L) and (M/T). Analysis of the filaments gives a median width of 0.14 pc and no correlation between width and length. Out of those 122 cores, 92 are present in filaments (75.4%) and the remaining were outside them. Out of the cores present in filaments, 57 (62%) cores are present in supercritical filaments ($M_{\rm line}>16M_{\odot }/{\rm pc}$).

Geminga is the first pulsar around which a remarkable TeV gamma-ray halo extending over a few degrees was discovered by MILAGRO, HAWC and later by H.E.S.S., and by Fermi-LAT in the GeV band. More middle-aged pulsars have exhibited gamma-ray halos, and they are now recognized as an emerging class of Galactic gamma-ray sources. The emission appears in the late evolution stage of pulsars, and is most plausibly explained by inverse Compton scattering of CMB and interstellar photons by relativistic electrons and positrons escaping from the pulsar wind nebulae. These observations pose a number of theoretical challenges. Tackling these questions requires constraining the ambient magnetic field properties, which can be achieved through X-ray observations. If the gamma-ray halos originate from a distribution of highly energetic electrons, synchrotron losses in the ambient magnetic fields of the same particles are expected to produce a diffuse X-ray emission with a similar spatial extension. We present the most comprehensive X-ray study of the Geminga pulsar halo to date, utilising archival data from XMM-Newton and NuSTAR. Our X-ray analysis covers a broad bandwidth ($0.5\rm{-}79$ keV) and large field of view ($\sim 4^\circ$) for the first time. This is achieved by accurately measuring the background over the entire field of view, and taking into account both focused and stray-light X-ray photons with NuSTAR. We find no significant emission and set robust constraints on the X-ray halo flux. These are translated to stringent constraints on the ambient magnetic field strength and the diffusion coefficient by using a physical model considering particle injection, diffusion and cooling over the pulsar's lifetime, which is tuned by fitting multi-wavelength data. Our novel methodology for modelling and searching for synchrotron X-ray halos can be applied to other pulsar halo candidates.

Hot Jupiters, orbiting their host stars at extremely close distances, undergo tidal evolution, with some being engulfed by their stars due to angular momentum exchanges induced by tidal forces. However, achieving double synchronization can prolong their survival. Using the MESA stellar evolution code, combined with the magnetic braking model of Matt et al. (2015), we calculate 25,000 models with different metallicity and study how to attain the conditions that trigger the long-term double synchronization. Our results indicate that massive planets orbiting stars with lower convective turnover time are easier to achieve long-term double synchronization. The rotation angular velocity at the equilibrium point ($\Omega_{\mathrm{sta}}$) is almost equal to orbital angular velocity of planet ($\mathrm{n}$) for the majority of the main sequence lifetime if a system has undergone a long-term double synchronization, regardless of their state at this moment. We further compared our results with known parameters of giant planetary systems and found that those systems with larger planetary masses and lower convective turnover time seem to be less sensitive to changes in the tidal quality factor $Q'_{_*}$. We suggest that for systems that fall on the state of $\Omega_{\mathrm{sta}} \approx n$, regardless of their current state, the synchronization will persist for a long time if orbital synchronization occurs at any stage of their evolution. Our results can be applied to estimate whether a system has experienced long-term double synchronization in the past or may experience it in the future.

When there are many observations of an astronomical source - many images with different dithers, or many spectra taken at different barycentric velocities - it is standard practice to shift and stack the data, to (for example) make a high signal-to-noise average image or mean spectrum. Bound-saturating measurements are made by manipulating a likelihood function, where the data are treated as fixed, and model parameters are modified to fit the data. Traditional shifting and stacking of data can be converted into a model-fitting procedure, such that the data are not modified, and yet the output is the shift-adjusted mean. The key component of this conversion is a spectral model that is completely flexible but also a continuous function of wavelength (or position in the case of imaging) that can represent any signal being measured by the device after any reasonable translation (or rotation or field distortion). The benefits of a modeling approach are myriad: The sacred data never are modified. Noise maps, data gaps, and bad-data masks don't require interpolation. The output can take the form of an image or spectrum evaluated on a pixel grid, as is traditional. In addition to shifts, the model can account for line-spread or point-spread function variations, world-coordinate-system variations, and calibration or normalization variations. The noise in the output becomes uncorrelated across neighboring pixels as the shifts deliver good coverage in some sense. The only cost is a small increase in computational complexity over that of traditional methods. We demonstrate the method with a small data example and we provide open source sample code for re-use.

We report a simulation and quantification of the impact of the Starlink constellation on LSST in terms of the trail surface brightness using a BRDF-based satellite photometric model. A total of 11,908 satellites from the Gen1 and Gen2A constellations are used to focus on the interference to the initial phase of LSST operation. The all-sky simulation shows that approximately 69.33% of the visible satellites over station have an apparent brightness greater than 7 mag with a v1.5 satellite model. The impact of satellite streaks exhibit a non-monotonic relationship to the solar altitude, with the worst moments occurring around $-15^{\circ}$ solar altitude. The assessment based on simulated schedules indicates that no trails can reach the saturation-level magnitude, but 71.61% trails show a surface brightness brighter than the best-case crosstalk correctable limits, and this percentage increases as the dodging weight increases. Therefore, avoiding satellites in the scheduler algorithm is an effective mitigation method, but both the number of streaks and their brightness should be taken into account simultaneously.

Pulsars are known to manifest complex phenomena, such as nulling, sub-pulse drifting, and periodicity. Within the purview of this investigation, we have harnessed the wavelet analysis technique to scrutinize the multifaceted periodicities and sub-pulse drifting characteristics exhibited by PSR J1926-0652, discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Our analysis draws upon the rich dataset acquired from the FAST ultra-wide-bandwidth receiver (UWB), meticulously examining pulse attributes encompassing an entire pulse. It is notable that the pulse apex recurrently manifests approximately every 17.11 times the pulsar's period P, individual sub-pulses exhibit a drifting phenomenon, with a phase decrement of approximately $1.04^\circ$ over each P. Intriguingly, the central phase of each sub-pulse track gradually increments over temporal evolution. Furthermore, the relative offset distribution between successive sub-pulse tracks emanating from the trailing and leading components remains comparatively stable, with a central tendency of approximately $\sim$6.87$\pm$2.56 P. Most notably, derived from the outcomes of wavelet analysis, we ascertain a negative correlation of -0.98 between the periods of sub-pulses and their drifting rates, alongside the intrinsic period of sub-pulses identified at 28.14 seconds.

From a careful visual scrutiny of the radio structures of a well-defined sample of 2428 sources in the LoTSS DR2 survey made at 144 MHz with a 6" beam, we have selected a subset of 25 (i.e., 1%) sources showing highly unusual radio structures, not conforming to the prevalent radio morphological classification. Here we present and briefly discuss the basic properties of these rare morphological outliers and attempt to dissect their morphological peculiarities, based on multi-wavelength radio images and radio-optical overlays. Also, we underscore the need to accord due importance to such anomalous radio sources, considering the challenge they pose to the standard theoretical models and simulations of extragalactic double radio sources.

Precise measurements of black hole (BHs) masses are necessary to understand the coevolution of these sources and their host galaxies. Sometimes in the center of a galaxy there is not one, but two BHs. The BH duality of the quasar nucleus SDSS~J075217.84+193542.2 (herein SDSS~J0752) was recently proposed based on the observed strict periodicity of optical emission from the source. We tested this assumption using X-ray observations with Swift/XRT (2008--2010). We fitted the SDSS~J075217 spectrum using a Comptonization model and discovered soft X-ray variability in the 0.3--10~keV energy range. We pursued a scenario in which two supermassive BHs at the center of SDSS~J0752 form a pair; and the less massive (secondary) BH periodically crosses/punctures the disk around the more massive (primary) BH. We associate these periodic crossings with tidal disruptions of the disk and, as a consequence, with an increase in X-rays seen as a flare in SDSS~J0752. During such an X-ray flare event (2008--2010), we discovered a change in the source spectral states and the photon index saturation at the $\Gamma\sim3$ level with mass accretion rate $\dot M$. For BH mass scaling we used sources: OJ~287, M101 ULX--1 and HLX--1 ESO~243--49, as a reference ones, and found that M$_{SDSS}=9\times 10^7$ solar masses, assuming $d_{SDSS}= 500$ Mpc. Thus, we obtained a lower limit to a BH mass due the unknown inclination. In addition, we used the virial mass of the secondary BH based on $H_\alpha$-line measurements and we estimated the binary's inclination at SDSS~J0752, $i=80^{\circ}$, using a scaling technique.

The gravitational wave signals detected by the LIGO-Virgo-KAGRA collaboration can be explained by mergers of binary primordial black holes (PBHs) formed in the radiation dominated epoch. However, in early structures induced by the Poisson distribution of PBHs, a significant fraction of binaries are perturbed and avoid mergers. In addition, the internal dynamics of early halos lead to the formation of dense primordial black hole clusters within a few Hubble times from the moment of halo formation. In such clusters PBH binaries are effectively formed and their mergers potentially dominate in the modern era. We obtained that the PBH merger rate changes with redshift as $\mathcal{R} \propto (1 + z)^\beta$, where $\beta = 1.4 - 2.2$ reflects the influence of PBH clustering and depends on both $z$ and $f_{\rm PBH}$. The observed merger rate constraints the fraction of PBHs of tens solar masses in the composition of dark matter $f_{\rm PBH} \lesssim 0.001 - 0.1$ in dependence of the clustering efficiency.

We propose that gamma-ray burst pulses are produced when highly-relativistic jets sweep across an observer's line-of-sight. We hypothesize that axisymmetric jet profiles, coupled with special relativistic effects, produce the time-reversed properties of gamma-ray burst pulses. Curvature resulting from rapid jet expansion is responsible for much of the observed pulse asymmetry and hard-to-soft evolution. The relative obliqueness with which the jet crosses the line-of-sight explains the known GRB pulse morphological types. We explore two scenarios: one in which a rigid/semi-rigid jet moves laterally, and the other in which a ballistic jet sprays material from a laterally-moving nozzle. The ballistic jet model is favored based upon its consistency with standard emission mechanisms.

We present integral field spectroscopy toward and around J1044+0353, a rapidly growing, low-metallicity galaxy which produces extreme [O III] line emission. A new map of the O32 flux ratio reveals a density-bounded ionization cone emerging from the starburst. The interaction of the hydrogen ionizing radiation, produced by the very young starburst, with a cavity previously carved out by a galactic outflow, whose apex lies well outside the starburst region, determines the pathway for global Lyman continuum (LyC) escape. In the region within a few hundred parsecs of the young starburst, we demonstrate that superbubble breakthrough and blowout contribute distinct components to the [O III] line profile, broad and very-broad emission line wings, respectively. We draw attention to the large [O III] luminosity of the broad component and argue that this emission comes from photoionized, superbubble shells rather than a galactic wind as is often assumed. The spatially resolved H eII 4686 nebula appears to be photoionized by young star clusters. Stellar wind emission from these stars is likely the source of line wings detected on the He II line profile. This broader He II component indicates slow stellar winds, consistent with an increase in stellar rotation (and a decrease in effective escape speed) at the metallicity of J1044+0353. At least in J1044+0353, the recent star formation history plays a critical role in generating a global pathway for LyC escape, and the anisotropic escape would likely be missed by direct observations of the LyC.

Observations have revealed a significant dark matter deficit in the ultra-diffuse galaxy NGC1052-DF2 (DF2). It is widely accepted that the formation of this unique galaxy can be attributed to the tidal stripping of its host galaxy, NGC1052. In this study, we simulate the evolution of a satellite system containing globular clusters (GCs) within an accreting host halo in the framework of self-interacting dark matter (SIDM). Our simulation results suggest that the heightened tidal stripping resulting from self-interaction can give rise to the transformation of a conventional dwarf galaxy into a galaxy resembling DF2 in all its observed properties. By comparing the simulation results with identical initial conditions in both the standard cold dark matter (CDM) and SIDM models, we find that the latter is more likely to replicate the properties of DF2. Furthermore, we demonstrate that a DF2 analog can also be produced on an orbit with a greater pericenter distance by increasing the strength of self-interaction. This suggests that the issue of extreme orbital parameters can be mitigated by implementing the SIDM model. The distributions of the GC population derived in our SIDM simulation are consistent with the observed characteristics of DF2. For comparison, we also explored the potential for achieving GC distributions in the context of CDM.

A proposed mechanism for solar chromospheric heating is that magnetohydrodynamic waves propagate upward along magnetic field lines and dissipate their energy in the chromosphere. In particular, compressible magneto-acoustic waves may contribute to the heating. Theoretically, the components below the cutoff frequency cannot propagate into the chromosphere; however, the cutoff frequency depends on the inclination of the magnetic field lines. In this study, using high temporal cadence spectral data of IRIS and Hinode SOT spectropolarimeter (SP) in plages, we investigated the dependence of the low-frequency waves on magnetic-field properties and quantitatively estimated the amount of energy dissipation in the chromosphere. The following results were obtained: (a) The amount of energy dissipated by the low-frequency component (3--6 mHz) increases with the field inclination below 40 degrees, whereas it is decreased as a function of the field inclination above 40 degrees. (b) The amount of the energy is enhanced toward $10^4 W/m^2$, which is the energy required for heating in the chromospheric plage regions, when the magnetic field is higher than 600 G and inclined more than 40 degree. (c) In the photosphere, the low-frequency component has much more power in the magnetic field inclined more and weaker than 400 G. The results suggest that the observed low-frequency components can bring the energy along the magnetic field lines and that only a specific range of the field inclination angles and field strength may allow the low-frequency component to bring the sufficient amount of the energy into the chromosphere.

We present an analysis of the galaxy merger rate in the redshift range $4.0<z<9.0$ (i.e. about 1.5 to 0.5 Gyr after the Big Bang) based on visually identified galaxy mergers from morphological parameter analysis. Our dataset is based on high-resolution NIRCam JWST data (F150W and F2000W broad-band filters) in the low-to-moderate magnification ($\mu<2$) regions of the Abell 2744 cluster field. From a parent set of 675 galaxies $(M_{UV}\in[-26.6,-17.9])$, we identify 64 merger candidates from the Gini, $M_{20}$ and Asymmetry morphological parameters, leading to a merger fraction $f_m=0.11\pm0.04$. There is no evidence of redshift evolution of $f_m$ even at the highest redshift considered, thus extending well into the epoch of reionization the constant trend seen previously at $z\lesssim 6$. Furthermore, we investigate any potential redshift dependent differences in the specific star formation rates between mergers and non-mergers. Our analysis reveals no significant correlation in this regard, with deviations in the studied redshift range typically falling within $0.25$ dex (logarithmic scale) that can be attributed to sample variance and measurement errors. Finally, we also demonstrate that the classification of a merging system is robust with respect to the observed (and equivalently rest-frame) wavelength of the high-quality JWST broad-band images used. This preliminary study highlights the potential for progress in quantifying galaxy assembly through mergers during the epoch of reionization, with significant sample size growth expected from upcoming large JWST infrared imaging datasets.

Post-asymptotic giant branch stars (post-AGB) in binary systems, with typical orbital periods between ~100 to ~1000 days, result from a poorly understood interaction that terminates their precursory AGB phase. The majority of these binaries display a photospheric anomaly called 'chemical depletion', thought to arise from an interaction between the circumbinary disc and the post-AGB star, leading to the reaccretion of pure gas onto the star, devoid of refractory elements due to dust formation. In this paper, we focus on a subset of chemically peculiar binary post-AGBs in the Galaxy and the Magellanic Clouds (MCs) whose high-resolution optical spectroscopic study revealed a carbon and s-process enrichment, contrary to the commonly observed photospheric chemical depletion. Using spectral energy distribution (SED) fitting and period-luminosity-colour (PLC) relation methods, we determine the luminosity of the targets ($2700-8300L_{\odot}$), which enables confirmation of their evolutionary phase and estimation of initial masses (as a function of metallicity) ($1-2.5M_{\odot}$). In conjunction with predictions from dedicated ATON stellar evolutionary models, our results indicate a predominant intrinsic enrichment of carbon and s-process elements in our binary post-AGB targets. We qualitatively rule out extrinsic enrichment and inherited s-process enrichment from the host galaxy as plausible explanations for the observed overabundances. Our chemically peculiar subset of intrinsic carbon and s-process enriched binary post-AGBs also hints at potential variation in the efficiency of chemical depletion between stars with C-rich and O-rich circumbinary disc chemistries. However, critical observational studies of circumbinary disc chemistry are necessary to address gaps in our current understanding of disc-binary interactions inducing chemical depletion in binary post-AGB systems.

The inflated radii observed in hundreds of Hot Jupiters represent a long-standing open issue. The observed correlation between radii and irradiation strength, and the occasional extreme cases, nearly double the size of Jupiter, remain without a comprehensive quantitative explanation. In this investigation, we delve into this issue within the framework of Ohmic dissipation, one of the most promising mechanisms for explaining the radius anomaly. Using the evolutionary code MESA, we simulate the evolution of irradiated giant planets, spanning the range 1 to 8 Jupiter masses, incorporating an internal source of Ohmic dissipation located beneath the radiative-convective boundary. Our modeling is based on physical parameters, and accounts for the approximated conductivity and the evolution of the magnetic fields, utilizing widely-used scaling laws. We compute the radius evolution across a spectrum of masses and equilibrium temperatures, considering varying amounts of Ohmic dissipation, calculated with the internal conductivity profile and an effective parametrization of the currents, based on the typical radius of curvature of the field lines. Our analysis reveals that this internal Ohmic dissipation can broadly reproduce the range of observed radii using values of radius of curvature up to about one order of magnitude lower than what we estimate from the Juno measurements of the Jovian magnetosphere and from MHD dynamo simulations presented herein. The observed trend with equilibrium temperature can be explained if the highly-irradiated planets have more intense and more small-scale magnetic fields. This suggests the possibility of an interplay between atmospherically induced currents and the interior, via turbulence, in agreement with recent box simulations of turbulent MHD in atmospheric columns.

The Mg II resonance doublet at 2796 {\AA} and 2803 {\AA} is an increasingly important tool to study cold, $T \sim 10^{4}\,$K, gas -- an observational driven development requiring theoretical support. We develop a new Monte Carlo radiative transfer code to systematically study the joined Mg II and Ly$\alpha$ escape through homogeneous and `clumpy' multiphase gas with dust in arbitrary 3D geometries. Our main findings are: (i) The Mg II spectrum differs from Ly$\alpha$ due to the large difference in column densities, even though the atomic physics of the two lines are similar. (ii) the Mg II escape fraction is generally higher than that of Ly$\alpha$ because of lower dust optical depths and path lengths -- but large variations due to differences in dust models and the clumpiness of the cold medium exist. (iii) Clumpy media possess a `critical covering factor' above which Mg II radiative transfer matches a homogeneous medium. The critical covering factors for Mg II and Ly$\alpha$ differ, allowing constraints on the cold gas structure. (iv) The Mg II doublet ratio $R_{\rm MgII}$ varies for strong outflows/inflows ($\gtrsim 700 \mathrm{km\,s}^{-1}$), in particular, $R_{\rm MgII}<1$ being an unambiguous tracer for powerful galactic winds. (v) Scattering of stellar continuum photons can decrease $R_{\rm MgII}$ from two to one, allowing constraints on the scattering medium. Notably, we introduce a novel probe of the cold gas column density -- the halo doublet ratio -- which we show to be a powerful indicator of ionizing photon escape. We discuss our results in the context of interpreting and modeling observations as well as their implications for other resonant doublets.

It is believed that the core of a neutron star can be host to various novel phases of matter, from nucleon superfluid phase to exotic high baryon density QCD phases. Different observational signals for such phase transitions have been discussed in the literature. Here, we point out a unique phenomenon associated with phase transition to a superfluid phase, which may be the nucleon superfluid phase or a phase like the CFL phase, allowing for superfluid vortices. In any superfluid phase transition, a random network of vortices forms via the so-called Kibble-Zurek mechanism, which eventually mostly decays away, finally leaving primarily vortices arising from the initial angular momentum of the core. This transient, random vortex network can have a non-zero net angular momentum for the superfluid component, which will generally be oriented in an arbitrary direction. This is in contrast to the final vortices, which arise from initial rotation and hence have the initial angular momentum of the neutron star. The angular momentum of the random vortex network is balanced by an equal and opposite angular momentum in the normal fluid due to the conservation of angular momentum, thereby imparting an arbitrarily oriented angular momentum component to the outer shell of the neutron star. This will affect the pulse timing and pulse profile of a pulsar. These changes in the pulses will decay away in a characteristic manner as the random vortex network decays, obeying specific scaling laws leading to universal features for the detection of superfluid transitions occurring in a pulsar core.

Direct imaging of exoplanets requires to separate the background noise from the exoplanet signals. Statistical methods have been recently proposed to avoid subtracting any signal of interest as opposed to initial self-subtracting methods based on Angular Differential Imaging (ADI). However, unless conservative thresholds are chosen to claim for a detection, such approaches tend to produce a list of candidates that include many false positives. Choosing high, conservative, thresholds leads to miss the faintest planets. We extend a statistical framework with a logistic regression to filter the list of candidates. Features with physical/optical meaning (in two wavelengths) are used, leading to a very fast and pragmatic approach. The overall method requires a simple edge detection (image processing) and clustering algorithm to work with sub-images. To estimate its efficiency, we apply our approach to targets observed with the ESO/SPHERE high contrast imager, that were previously used as tests for blind surveys. Experimental results with injected signals show that either the number of false detections is considerably reduced or faint exoplanets that would otherwise not be detected can be sometimes found. Typically, on the blind tests performed, we are now able to detect around 50% more of the injected planets with an SNR below 5, and with a very low number of additional candidates.

Aims. In this work, we intend to examine how environment influences the merger fraction, from the low density field environment to higher density groups and clusters. We also aim to study how the properties of a group or cluster, as well as the position of a galaxy in the group or cluster, influences the merger fraction. Methods. We identified galaxy groups and clusters in the North Ecliptic Pole using a friends-of-friends algorithm and the local density. Once identified, we determined the central galaxies, group radii, velocity dispersions, and group masses of these groups and clusters. Merging systems were identified with a neural network as well as visually. With these, we examined how the merger fraction changes as the local density changes for all galaxies as well as how the merger fraction changes as the properties of the groups or clusters change. Results. We find that the merger fraction increases as local density increases and decreases as the velocity dispersion increases, as is often found in literature. A decrease in merger fraction as the group mass increases is also found. We also find groups with larger radii have higher merger fractions. The number of galaxies in a group does not influence the merger fraction. Conclusions. The decrease in merger fraction as group mass increases is a result of the link between group mass and velocity dispersion. Hence, this decrease of merger fraction with increasing mass is a result of the decrease of merger fraction with velocity dispersion. The increasing relation between group radii and merger fraction may be a result of larger groups having smaller velocity dispersion at a larger distance from the centre or larger groups hosting smaller, infalling groups with more mergers. However, we do not find evidence of smaller groups having higher merger fractions.

In this paper, we revisit the data for the galaxy cluster ZwCl 0024.0+1652 provided by the GLACE survey and study the mass--metallicity function and its relationship with the environment. Here we describe an alternative way to reduce the data from OSIRIS tunable filters. This method gives us better uncertainties in the fluxes of the emission lines and the derived quantities. We present an updated catalogue of cluster galaxies with emission in H$\alpha$ and [N\,{\sc{ii}}] $\lambda\lambda$6548,6583. We also discuss the biases of these new fluxes and describe the way in which we calculated the mass--metallicity relationship and its uncertainties. We generated a new catalogue of 84 emission-line galaxies with reliable fluxes in [N\,{\sc{ii}}] and H$\alpha$ lines from a list of 174 galaxies. We find a relationship between the clustercentric radius and the density of galaxies. We derived the mass--metallicity relationship for ZwCl 0024.0+1652 and compared it with clusters and field galaxies from the literature. We find a difference in the mass--metallicity relationship when compared to more massive clusters, with the latter showing on average higher values of abundance. This could be an effect of the quenching of the star formation, which seems to be more prevalent in low-mass galaxies in more massive clusters. We find little to no difference between ZwCl 0024.0+1652 galaxies and field galaxies located at the same redshift.

Advancements in instrumentation have revealed a multitude of small-scale EUV events in the solar atmosphere. Our aim is to employ high-resolution magnetograms to gain a detailed understanding of the magnetic origin of such phenomena. We have used coordinated observations from SST, IRIS, and SDO to analyze an ephemeral magnetic flux emergence episode and the following chain of small-scale energetic events. These unique observations clearly link these phenomena together. The high-resolution (0."057/pixel) magnetograms obtained with SST/CRISP allows us to reliably measure the magnetic field at the photosphere and detect the emerging bipole that causes the subsequent eruptive atmospheric events. Notably, this small-scale emergence episode remains indiscernible in the lower resolution SDO/HMI magnetograms (0."5/pixel). We report the appearance of a dark bubble in Ca II K related to the emerging bipole, a sign of the canonical expanding magnetic dome predicted in flux emergence simulations. Evidences of reconnection are also found: first through an Ellerman bomb, and later by the launch of a surge next to a UV burst. The UV burst exhibits a weak EUV counterpart in the coronal SDO/AIA channels. By calculating DEM, its plasma is shown to reach a temperature beyond 1 MK and have densities between the upper chromosphere and transition region. Our study showcases the importance of high-resolution magnetograms to unveil the mechanisms triggering phenomena such as EBs, UV bursts, and surges. This could hold implications for small-scale events akin to those recently reported in EUV using Solar Orbiter. The finding of temperatures beyond 1 MK in the UV burst plasma strongly suggests that we are examining analogous features. Therefore, we signal caution regarding drawing conclusions from full-disk magnetograms that lack the necessary resolution to reveal their true magnetic origin.

Atomic Hydrogen-21 cm transition (HI) is an excellent tracer to study and understand the properties of the atomic gas in the Galaxy. Using the Westerbork Synthesis Radio Telescope (WSRT), we observed 12 quasar sightlines to detect galactic HI in absorption. We achieve an optical depth RMS of $\sim 1-2 \times 10^{-3}$, essential to detect the Warm Neutral Medium (WNM). We detect HI absorption in all our sightlines except along 1006+349, for which we set a strict upper limit on the spin temperature as $\langle T_s \rangle > 570$ K. We find around 50\% of our sightlines have $\langle T_s \rangle > 500$ K, indicating a WNM dominance. Further, we calculate an upper limit of the CNM fraction along our sightlines and find a median CNM fraction of $\sim 0.12$. With our observations, we reconfirm the existence of a threshold column density of $\sim 2 \times 10^{20} \ cm^{-2}$ to form CNM in the ISM. Using a two-temperature model of the HI disk, we explore the distribution of spin temperature in the Galactic ISM. We find that a simple fixed axisymmetric two-temperature model could not produce either the observed column density or the integral optical depth. This indicates the existence of a more complex distribution of spin temperatures in the Galaxy.

New interior models of Jupiter and Saturn suggest that both planets have "fuzzy cores". These cores should be viewed as central regions that are enriched with heavy elements but are not distinct from the rest of the deep interior. These cores may contain large amounts of hydrogen and helium though small pure-heavy element cores may also exist. New measurements along with advanced planetary modeling have revolutionized the way we think about the interiors of giant planets and provide important constraints for planet formation and evolution theories. These developments are also relevant for the characterization of giant exoplanets.

Narrow absorption features in nearby supernova (SN) spectra are a powerful diagnostic of the slow-moving material in the line of sight: they are extensively used to infer dust extinction from the host galaxies, and they can also serve in the detection of circumstellar material originating from the SN progenitor and present in the vicinity of the explosion. Despite their wide use, very few studies have examined the biases of the methods to characterize narrow lines, and not many statistical analyses exist. This is the first paper of a series in which we present a statistical analysis of narrow lines of SN spectra of various resolutions. We develop a robust automated methodology to measure the equivalent width (EW) and velocity of narrow absorption lines from intervening material in the line of sight of SNe, including Na I D , Ca II H&K, K i and diffuse interstellar bands (DIBs). We carefully study systematic biases in heterogeneous spectra from the literature by simulating different signal-to-noise, spectral resolution, slit size and orientation and present the real capabilities and limitations of using low- and mid-resolution spectra to study these lines. In particular, we find that the measurement of the equivalent width of the narrow lines in low-resolution spectra is highly affected by the evolving broad P-Cygni profiles of the SN ejecta, both for core-collapse and type Ia SNe, inducing a conspicuous apparent evolution. We present thus an easy way to detect and exclude those cases to obtain more robust and reliable measurements. Finally, after considering all possible effects, we analyse the temporal evolution of the narrow features in a large sample of nearby SNe to detect any possible variation in their EWs over time. We find no time evolution of the narrow line features in our large sample for all SN types

Context. The High-Mass X-ray Binary (HMXB) system GX 301$-$2 is a persistent source with a well-known variable cyclotron line centered at 35 keV. Recently, a second cyclotron line at 50 keV has been reported with a presumably different behavior than the 35 keV line. Aims. We investigate the presence of the newly discovered cyclotron line in the phase-averaged and phase-resolved spectra at higher luminosities than before. We further aim to determine the pulse-phase variability of both lines. Methods. We analyze a NuSTAR observation of GX 301$-$2 covering the pre-periastron flare, where the source luminosity reached its peak of ${\sim} 4 \times 10^{37}\,\mathrm{erg}\,\mathrm{s}^{-1}$ in the 5-50 keV range. We analyze the phase-averaged spectra in the NuSTAR energy range from 3.5-79 keV for both the complete observation and three time segments of it. We further analyze the phase-resolved spectra and the pulse-phase variability of continuum and cyclotron line parameters. Results. We confirm that the description of the phase-averaged spectrum requires a second absorption feature at $51.5^{+1.1}_{-1.0}$ keV besides the established line at 35 keV. The statistical significance of this feature in the phase-averaged spectrum is $>99.999\%$. We further find that the 50 keV cyclotron line is present in three of eight phase bins. Conclusions. Based on the results of our analysis, we confirm that the detected absorption feature is very likely to be a cyclotron line. We discuss a variety of physical scenarios which could explain the proposed anharmonicity, but also outline circumstances under which the lines are harmonically related. We further present the cyclotron line history of GX 301$-$2 and evaluate concordance among each other. We also discuss an alternative spectral model including cyclotron line emission wings.

We present the new ARTEMIS emulator suite of high resolution (baryon mass of $2.23 \times 10^{4}$ $h^{-1}$M$_{\odot}$) zoom-in simulations of Milky Way mass systems. Here, three haloes from the original ARTEMIS sample have been rerun multiple times, systematically varying parameters for the stellar feedback model, the density threshold for star formation, the reionisation redshift and the assumed warm dark matter (WDM) particle mass (assuming a thermal relic). From these simulations emulators are trained for a wide range of statistics that allow for fast predictions at combinations of parameters not originally sampled, running in $\sim 1$ms (a factor of $\sim 10^{11}$ faster than the simulations). In this paper we explore the dependence of the central haloes' stellar mass on the varied parameters, finding the stellar feedback parameters to be the most important. When constraining the parameters to match the present-day stellar mass halo mass relation inferred from abundance matching we find that there is a strong degeneracy in the stellar feedback parameters, corresponding to a freedom in formation time of the stellar component for a fixed halo assembly history. We additionally explore the dependence of the satellite stellar mass function, where it is found that variations in stellar feedback, the reionisation redshift and the WDM mass all have a significant effect. The presented emulators are a powerful tool which allows for fundamentally new ways of analysing and interpreting cosmological hydrodynamic simulations. Crucially, allowing their free (subgrid) parameters to be varied and marginalised, leading to more robust constraints and predictions.

Context. Identifying minerals on asteroid surfaces is difficult as space weathering modifies the minerals infrared spectra. This shouldbe better understood for proper interpretation. Aims. We simulated the space weathering effects on a meteorite and recorded the alterations of the crystalline structure, such as the change in peak positions and full width at half maximum values. Methods. We used proton irradiation to simulate the effects of solar wind on a sample of NWA 10580 CO3 chondrite meteorites. After irradiation in three gradually increased steps with 1 keV ion energy, we used infrared microscopic reflectance and diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) to identify and understand the consequences of irradiation. Results. We find negative peak shifts after the first and second irradiations at pyroxene and feldspar minerals, similarly to the literature, and this shift was attributed to Mg loss. However, after the third irradiation a positive change in values in wavenumber emerged for silicates, which could come from the distortion of SiO4 tetrahedra, resembling shock deformation. The full width at half maximum values of major bands show a positive (increasing) trend after irradiations in the case of feldspars, using IR reflection measurements. Comparing DRIFTS and reflection infrared data, the peak positions of major mineral bands were at similar wavenumbers, but differences can be observed in minor bands. Conclusions. We measured the spectral changes of meteorite minerals after high doses of proton irradiation for several minerals. We show the first of these measurements for feldspars; previous works only presented pyroxene, olivine, and phyllosilicates.

Validating modeling choices through simulated analyses and quantifying the impact of different systematic effects will form a major computational bottleneck in the preparation for 3$\times$2 analysis with Stage-IV surveys such as Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). We can significantly reduce the computational requirements by using machine learning based emulators, which allow us to run fast inference while maintaining the full realism of the data analysis pipeline. In this paper, we use such an emulator to run simulated 3$\times$2 (cosmic shear, galaxy-galaxy lensing, and galaxy clustering) analyses for mock LSST-Y1/Y3/Y6/Y10 surveys and study the impact of various systematic effects (galaxy bias, intrinsic alignment, baryonic physics, shear calibration and photo-$z$ uncertainties). Closely following the DESC Science Requirement Document (with several updates) our main findings are: {\it a)} The largest contribution to the `systematic error budget' of LSST 3$\times$2 analysis comes from galaxy bias uncertainties, while the contribution of baryonic and shear calibration uncertainties are significantly less important. {\it b)} Tighter constraints on intrinsic alignment and photo-$z$ parameters can improve cosmological constraints noticeably, which illustrates synergies of LSST and spectroscopic surveys. {\it c)} The scale cuts adopted in the DESC SRD may be too conservative and pushing to smaller scales can increase cosmological information significantly. {\it d)} We investigate the impact of photo-$z$ outliers on 3$\times$2 pt analysis and find that we need to determine the outlier fraction to within $5-10\%$ accuracy to ensure robust cosmological analysis. We caution that these findings depend on analysis choices (parameterizations, priors, scale cuts) and can change for different settings.

'Debris disks' are collections of small bodies around stars, such as the Asteroid Belt and Kuiper Belt in our Solar System. These disks are composed of objects smaller than planets, including asteroids, comets, dust, and dwarf planets. We detect debris disks around a significant fraction of stars, and these disks appear to be common components of planetary systems. Extrasolar debris disks have a broad range of locations, shapes and features. This chapter provides an introduction to debris disks around main-sequence stars. It summarises our understanding of the field, and covers a wide range of concepts from observations and theory. It describes how we detect extrasolar debris disks, what we see, and what these observations tell us. It also describes how debris disks evolve, and how they interact with planets. The chapter concludes by discussing several unsolved questions in debris-disk science.

We present a detailed analysis of the X-ray luminosity (Lx) dependence of the cyclotron absorption line energy (Ecyc) for the X-ray binary pulsar 1A 0535+262 during its 2020 giant outburst based on pulse-to-pulse analysis. By applying this technique to high cadence observations of Insight-HXMT, we reveal the most comprehensive Ecyc-Lx correlation across a broad luminosity range of ~(0.03-1.3)*10^38 erg/s. Apart from the positive and negative correlations between cyclotron line energy and luminosity at Lx~(1-3)*10^37 erg/s and ~(7-13)*10^37 erg/s, which are expected from the typical subcritical and supercritical accretion regimes, respectively, a plateau in the correlation is also detected at ~(3-7)*10^37 erg/s^-1. Moreover, at the lowest luminosity level (Lx<10^37 erg/s), the positive Ecyc-Lx correlation seems to be broken, and the pulse profile also occurs a significant transition. These discoveries provide the first complete view on the correlation between luminosity and the centriod energy of the cyclotron line, and therefore are relevant for understanding how accretion onto magnetized neutron stars depends on luminosity.

Recently LHAASO Collaboration gives precise measurements of cosmic rays (CR) all particle energy spectrum and mean logarithmic mass $\left \langle \ln A \right \rangle$ from 0.3 PeV to 30 PeV. Combining the CR measurements by AMS-02 and DAMPE in space and that by LHAASO and Auger on the ground we construct a model to recover all these measurements from tens of GeV to tens of EeV. We find the LHAASO measurement is crucial in the model construction by connecting the Galactic component to the extragalactic component. The precise measurements of CR spectra for individual species by AMS-02 and DAMPE together with the newest LHAASO results clearly indicates three Galactic CR components, that is, a soft low energy background, a hard high energy component, and a local source contribution. However, the LHAASO data show that above $\sim 10^{16}$ eV a nonnegligible extragalactic component must be included. Combining the Auger results and the LHAASO results we figure out the extragalactic CRs which need at least two components at lower and higher energies. Thanks to the precise measurements by LHAASO the constraints on the model parameters are quite stringent. The spectra features and mass measurements in all energy range are all well reproduced in the model.

In this paper we study the role of radiative cooling in a two-fluid model consisting of coupled neutrals and charged particles. We first analyze the linearized two-fluid equations where we include radiative losses in the energy equation for the charged particles. In a 1D geometry for parallel propagation and in the limiting cases of weak and strong coupling, it can be shown analytically that the instability conditions for the thermal mode and the sound waves, the isobaric and isentropic criteria, respectively, remain unchanged with respect to one-fluid radiative plasmas. For the parameters considered in this paper, representative for the solar corona, the radiative cooling produces growth of the thermal mode and damping of the sound waves. When neutrals are included and are sufficiently coupled to the charges, the thermal mode growth rate and the wave damping both reduce by the same factor, which depends on the ionization fraction only. For a heating function which is constant in time, we find that the growth of the thermal mode and the damping of the sound waves are slightly larger. The numerical calculation of the eigenvalues of the general system of equations in a 3D geometry confirm the analytic results. We then run 2D fully nonlinear simulations which give consistent results: a higher ionization fraction or lower coupling will increase the growth rate. The magnetic field contribution is negligible in the linear phase. Ionization-recombination effects might play an important role because the radiative cooling produces a large range of temperatures in the system. In the numerical simulation, after the first condensation phase, when the minimum temperature is reached, the fraction of neutrals increases four orders of magnitude because of the recombination.

We assess the consistency of cosmological models that alter the size of the sound horizon at last scattering to resolve the Hubble tension with data from ACT + Planck CMB lensing, Big Bang Nucleosynthesis, and supernova data from Pantheon or Pantheon+. We use early dark energy (EDE) as an example model but conclude that the results apply to other similar models. We constrain $\Lambda$CDM and EDE with these data finding that while they can constrain $\Lambda$CDM very tightly, EDE opens up the parameter space significantly and allows $H_0 > 72$ km s$^{-1}$ Mpc$^{-1}$. We combine these data with measurements from ACT + Planck TT650TEEE CMB primary anisotropy and galaxy baryon acoustic oscillations, and find that overall, EDE fits these data better than $\Lambda$CDM at $\approx 2\sigma$. However, the fit to specifically the sound-horizon-independent measurements is worse for EDE than $\Lambda$CDM. We assess this increase in $\chi^2$ coming from the sound-horizon-independent measurements and find that the best-fit model is still consistent with a random statistical fluctuation even with $H_0$ values around $72$ km s$^{-1}$ Mpc$^{-1}$. We conclude that these specific sound-horizon-independent data cannot rule out the possibility of a miscalibration of the size of the sound horizon, but leave open the possibility that other current or future sound-horizon-independent data sets could rule out such a miscalibration.

Long gamma-ray burst (GRB) prompt emission shows a correlation between the intrinsic peak energy, $E_{\mathrm{p,i}}$, of the time-average $\nu F_{\nu}$ spectrum and the isotropic-equivalent peak gamma-ray luminosity, $L_{{\rm p,iso}}$, as well as the total released energy, $E_{\rm iso}$. The same correlation is found within individual bursts, when time-resolved $E_{\rm p,i}$ and $L_{\rm iso}$ are considered. These correlations are characterised by an intrinsic dispersion, whose origin is still unknown. Discovering the origin of the correlation and of its dispersion would shed light on the still poorly understood prompt emission and would propel GRBs to powerful standard candles. We studied the dispersion of both isotropic-equivalent and collimation-corrected time-resolved correlations. We also investigated whether the intrinsic dispersion computed within individual GRBs is different from that obtained including different bursts into a unique sample. We then searched for correlations between key features, like Lorentz factor and jet opening angle, and intrinsic dispersion, when the latter is treated as one of the characterising We performed a time-resolved spectral analysis of 20 long Type-II or collapsar-candidate GRBs detected by the Fermi Gamma-ray Burst Monitor with known redshift and estimates of jet opening angle and/or Lorentz factor. The collimation-corrected correlation appears to be no less dispersed than the isotropic-equivalent one. Also, individual GRBs are significantly less dispersed than the whole sample. We excluded (at $4.2 \sigma$ confidence level) the difference in samples' sizes as the possible reason, thus confirming that individual GRBs are {\em intrinsically} less dispersed than the whole sample. No correlation was found between intrinsic dispersion and other key properties for the few GRBs with available information.

The measurement of the flux of muons produced in cosmic ray air showers is essential for the study of primary cosmic rays. Such measurements are important in extensive air shower detectors to assess the energy spectrum and the chemical composition of the cosmic ray flux, complementary to the information provided by fluorescence detectors. Detailed simulations of the cosmic ray air showers are carried out, using codes such as CORSIKA, to estimate the muon flux at sea level. These simulations are based on the choice of hadronic interaction models, for which improvements have been implemented in the post-LHC era. In this work, a deficit in simulations that use state-of-the-art QCD models with respect to the measurement deep underwater with the KM3NeT neutrino detectors is reported. The KM3NeT/ARCA and KM3NeT/ORCA neutrino telescopes are sensitive to TeV muons originating mostly from primary cosmic rays with energies around 10 TeV. The predictions of state-of-the-art QCD models show that the deficit with respect to the data is constant in zenith angle; no dependency on the water overburden is observed. The observed deficit at a depth of several kilometres is compatible with the deficit seen in the comparison of the simulations and measurements at sea level.

We demonstrate the effectiveness of one of the many multi-tracer analyses enabled by Optimal Transport (OT) reconstruction. Leveraging a semi-discrete OT algorithm, we determine the displacements between initial and observed positions of biased tracers and the remaining matter field. With only redshift-space distorted final positions of biased tracers and a simple premise for the remaining mass distribution as input, OT solves the displacement field. This extracted field, assuming asymptotically uniform density and a gradient flow displacement, enables reconstruction of the initial overdensity fluctuation field. We show that the divergence of the OT displacement field is a good proxy of the linear density field, even though the method never assumes the linear theory growth. Additionally, this divergence field can be combined with the reconstructed protohalos to provide a higher signal-to-noise measurement of the BAO standard ruler than was possible with either measurement individually.

G352.7$-$0.1 is a mixed-morphology (MM) supernova remnant (SNR) with multiple radio arcs and has a disputed supernova origin. We conducted a spatially resolved spectroscopic study of the remnant with XMM-Newton X-ray data to investigate its explosion mechanism and explain its morphology. The global X-ray spectra of the SNR can be adequately reproduced using a metal-rich thermal plasma model with a temperature of $\sim 2$ keV and ionization timescale of $\sim 3\times 10^{10}~{\rm cm^{-3}~s}$. Through a comparison with various supernova nucleosynthesis models, we found that observed metal properties from Mg to Fe can be better described using core-collapse supernova models, while thermonuclear models fail to explain the observed high Mg/Si ratio. The best-fit supernova model suggests a $\sim 13$ $M_\odot$ progenitor star, consistent with previous estimates using the wind bubble size. We also discussed the possible mechanisms that may lead to SNR G352.7$-$0.1 being an MMSNR. By dividing the SNR into several regions, we found that the temperature and abundance do not significantly vary with regions, except for a decreased temperature and abundance in a region interacting with molecular clouds. The brightest X-ray emission of the SNR spatially matches with the inner radio structure, suggesting that the centrally filled X-ray morphology results from a projection effect.

Crystallization of dense matter in neutron star crusts and white dwarf cores may be similar to epitaxial crystal growth in terrestrial laboratories. However in stellar crystals, the spacing between horizontal planes has to gradually increase with the outward movement of the crystallization front, tracing decrease of the electron density. This process produces Coulomb crystals with stretched rather than cubic elementary cells. We extend the analysis of the elastic and breaking properties of such crystals to the face-centered (fc) lattice. Shear deformations orthogonal to the stretch direction have been studied for 22 crystallographic shear planes. A common property for all these planes is a reduction and eventual nulling of the breaking shear strain with deviation from the unstretched configuration. The effective shear moduli for deformations orthogonal to the stretch direction have been calculated. It is possible that the epitaxial crystallization in compact stars results in a formation of large-scale crystallites or, at least, in growth of the whole crystallization front perpendicular to particular crystallographic planes. For fc structure growth orthogonal to the $\{111\}$ planes, we expect that, at any density, $\sim 5\%$ ($\sim 0.5\%$) of crystallite height is occupied by layers one (two) orders of magnitude weaker than the bulk of the crystallite. This may be important for realistic modeling of crustquakes on neutron stars.

Strong gravitational lens system catalogues are typically used to constrain a combination of cosmological and empirical lens mass model parameters, even though the simplest singular isothermal sphere (SIS) models yield a $\chi^2$ per degree of freedom $\simeq 2$. To date, this problem has been alleviated by introducing additional empirical parameters in the extended power law (EPL) models and constraints from high resolution imagery. The EPL parameters are taken to vary from lens to lens, rather than defining universal lens profiles. We investigate these lens models using Bayesian methods through a novel alternative that treats spatial curvature via the non-FLRW Timescape cosmology. We apply Markov Chain Monte Carlo methods using the catalogue of 161 lens systems of Chen et al (arXiv:1809.09845) to simulate large mock catalogues for: (i) the standard $\Lambda$CDM model with zero spatial curvature; and (ii) the Timescape model. Furthermore, this methodology can be applied to any cosmological model. In agreement with previous results we find that in combination with SIS parameters, models with zero FLRW spatial curvature fit better as the free parameter approaches an unphysical empty universe, $\Omega_{\rm M0}\to0$. By contrast, the Timescape cosmology is found to prefer parameter values in which its cosmological parameter, the present void fraction, is driven to $f_{\rm v0}\to0.73$ matches, close to values found to best fit independent cosmological data sets: supernovae Ia distances and cosmic microwave background. This conclusion holds for a large range of seed values $f_{\rm v0}\to0.73\in\{0.1,0.9\}$, and for Timescape fits to both Timescape and FLRW mocks. Regardless of cosmology, unphysical estimates of the distance ratios given from power-law lens models result in poor goodness of fit. Nonetheless, the results are consistent with non-FLRW spatial curvature evolution.

We have identified nearly a hundred close triply eclipsing hierarchical triple star systems from data taken with the space telescope TESS. These systems are noteworthy in that we can potentially determine their dynamical and astrophysical parameters with a high precision. In the present paper, we report the comprehensive study of seven new compact triply eclipsing triple star systems taken from this larger sample: TICs 133771812, 176713425, 185615681, 287756035, 321978218, 323486857, and 650024463. Most of the data for this study come from TESS observations, but two of them have Gaia measurements of their outer orbits, and we obtained supplemental radial velocity (RV) measurements for three of the systems. The eclipse timing variation curves extracted from the TESS data, the photometric light curves, the RV points, and the spectral energy distribution (SED) are combined in a complex photodynamical analysis to yield the stellar and orbital parameters of all seven systems. Four of the systems are quite compact with outer periods in the range of 41-56 days. All of the systems are substantially flat, with mutual inclination angles of < ~2 degrees. Including the systems reported in this work, we have now studied in considerable detail some 30 triply eclipsing triples with TESS, and are accumulating a meaningful census of these systems.

Observations of Type Ia supernovae (SNe\,Ia) reveal diversity, even within assumed subcategories. Here, the composition of the peculiar iPTF16abc (SN\,2016bln) is derived by modeling a time series of optical spectra. iPTF16abc's early spectra combine traits of SNe 1999aa and 1991T known for weak \SiII\ $\lambda$ 6355 and prominent \FeIII\ features. However, it differs with weak early \FeIII\ lines, and persistent \CII\ lines post-peak. It also exhibits a weak \CaII\ H\&K feature aligning it with SN\,1991T, an observation supported by their bolometric light curves. The early attenuation of \FeIII\ results from abundance effect. The weakening of the \SiII\ $\lambda$ 6355 line, stems from silicon depletion in the outer shells, a characteristic shared by both SNe 1999aa and 1991T, indicating a common explosion mechanism that terminates nuclear burning at around 12000 \kms\, unseen in normal events. Beneath a thin layer of intermediate mass elements (IMEs) with a total mass of 0.18 \Msun, extends a \Nifs\ rich shell totaling 0.76 \Msun\ and generating a bolometric luminosity as high as ${L_{\mathrm{peak}}}=1.60 \pm 0.1 \times$ $10^{43}$ ergs s$^{-1}$. Inner layers, typical of SNe\,Ia, hold neutron-rich elements, (\Feff\ and \Nife), totaling 0.20 M${\odot}$. Stable iron, exceeding solar abundance, and carbon, coexist in the outermost layers, challenging existing explosion models. The presence of carbon down to $v\approx$ 9000\,\kms, totalling $\sim$ 0.01 \Msun\, unprecedented in this class, links iPTF16abc to SN\,2003fg-like events. The retention of 91T-like traits in iPTF16abc underscores its importance in understanding the diversity of SNe\,Ia.

Even after dark matter chemically freezes out in the early universe, electromagnetic cascades from dark matter annihilation can still perturb the background photon spectrum when the universe temperature cools down to 0.5~keV. We revisit the CMB spectrum distortions caused by $s$-wave dark matter annihilation under the updated Planck data and the future CMB sensitivity, concluding that $s$-wave annihilation cannot create observable distortions under forecast sensitivities of the (Super-)PIXIE missions. We further detail the case of $p$-wave dark matter annihilation, demonstrating the observability of the primordial $\mu$-distortion. Taking current constraints from primordial light elements, structure formations, cosmic electron-positron rays, and gamma rays, we find that the $\mu$-distortion reaching the observational limit as large as $\mu\simeq 3\times 10^{-8}$ can only be realized by a dark matter mass at 10--50~MeV and kinetic decoupling at keV temperature. The upper bound of the $p$-wave annihilation cross section can be strengthened by an order of magnitude if the $\mu$-distortion is not detected.

We study the physics of photon rings in a wide range of axisymmetric black holes admitting a separable Hamilton-Jacobi equation for the geodesics. Utilizing the Killing-Yano tensor, we derive the Penrose limit of the black holes, which describes the physics near the photon ring. The obtained plane wave geometry is directly linked to the frequency matrix of the massless wave equation, as well as the instabilities and Lyapunov exponents of the null geodesics. Consequently, the Lyapunov exponents and frequencies of the photon geodesics, along with the quasinormal modes, can be all extracted from a Hamiltonian in the Penrose limit plane wave metric. Additionally, we explore potential bounds on the Lyapunov exponent, the orbital and precession frequencies, in connection with the corresponding inverted harmonic oscillators and we discuss the possibility of photon rings serving as holographic horizons in a holographic duality framework for astrophysical black holes. Our formalism is applicable to spacetimes encompassing various types of black holes, including stationary ones like Kerr, Kerr-Newman, as well as static black holes such as Schwarzschild, Reissner-Nordstr\"om, among others.

These lecture notes provide a pedagogical introduction to some aspects of the inflationary cosmology, including the background scalar field dynamics, generation of primordial seed perturbations via quantum fluctuations during inflation, and the process of reheating after inflation in the single-field inflationary paradigm.

Shear strength and cohesion of granular materials are important geotechnical properties that play a crucial role in the stability and behavior of lunar and Martian regolith, as well as their terrestrial analog materials. To characterize and predict shear strength and cohesion for future space missions, it is also important to understand the effects of particle size distribution and density on these fundamental geotechnical properties. Generalized equations have been established using empirical data from direct shear measurements of lunar and Martian regolith simulants to quantify the effects of particle size distribution and density on cohesion and shear strength. Preliminary results are also presented highlighting the effects of atmospheric absorbed water on shear strength and cohesion when conducting experiments in atmospheric conditions on Earth. The results of this study show that cohesion increases exponentially with bulk density, while the exponential growth constant is also dependent on particle size distribution.

We perform the numerical simulation of primordial black hole formation from a nonspherical profile of the initial curvature perturbation $\zeta$. We consider the background expanding universe filled with the perfect fluid with the linear equation of state $p=w\rho$ ($w=1/3$ or $1/5$), where $p$ and $\rho$ are the pressure and the energy density, respectively. The initial condition is set in a way such that the principal directions of the second derivatives of $\zeta$ and $\triangle \zeta$ at the central peak are misaligned, where $\triangle$ is the Laplacian. In this setting, since the linearized density is proportional to $\triangle \zeta$, the inertia tensor and deformation tensor $\partial_i\partial_j \zeta$ are misaligned. Thus tidal torque may act and the spin of a resultant primordial black hole would be non-zero in general, although it is estimated to be very small from previous perturbative analyses.As a result, we do not find a finite value of the spin within our numerical precision, giving support for the negligibly small value of the black hole spin for $1/5\lesssim w \lesssim 1/3$. More specifically, our results suggest that the dimensionless PBH spin $s$ is typically so small that $s\ll0.1$ for $w\gtrsim0.2$.

In view to scrutinise the idea that nonlocal modifications of \GR~could dynamically address the dark energy problem, we investigate the evolution of the Universe at infrared scales as an \IDG~model of the Ricci scalar, without introducing the cosmological constant $\Lambda$ or any scalar field. The accelerated expansion of the late Universe is shown to be compatible with the emergence of nonlocal gravitational effects at sufficiently low energies. A technique for circumventing the mathematical complexity of the nonlocal cosmological equations is explained and, after drawing a connection with the Starobinsky gravity, verifiable predictions are considered, like a possible decreasing in the strength of the effective gravitational constant.

We argue that the axionic domain-wall with a QCD bias may be incompatible with the NANOGrav 15-year data on a stochastic gravitational wave (GW) background, when the domain wall network collapses in the hot-QCD induced local CP-odd domain. This is due to the drastic suppression of the QCD bias set by the QCD topological susceptibility in the presence of the CP-odd domain with nonzero $\theta$ parameter of order one which the QCD sphaleron could generate. We quantify the effect on the GW signals by working on a low-energy effective model of Nambu-Jona-Lasinio type in the mean field approximation. We find that only at $\theta=\pi$, the QCD bias tends to get significantly large enough due to the criticality of the thermal CP restoration, which would, however, give too big signal strengths to be consistent with the NANOGrav 15-year data and would also be subject to the strength of the phase transition at the criticality.

A simple and natural extension of the standard $\Lambda$CDM model is to allow relic neutrinos to have non-zero degeneracy. We confront this $\Lambda$CDM$\xi$ model, $\Lambda$CDM with neutrino mass $M_\nu$ and degeneracy $\xi_3$ as additional parameters, with the \textit{Planck} TT, lowT, plik--lensing, BAO, and DES datasets, and we observe a strong preference (Bayes factor $\log_{10}B=1.9$) for it over the standard $\Lambda$CDM model. Both the $H_0$ and $S_8$ tensions are resolved to within 1$\sigma$ with the same set of neutrino parameters, along with 3$\sigma$ evidence for nonzero neutrino mass ($M_\nu=0.58^{+0.17}_{-0.13}\ \mathrm{eV}$) and degeneracy ($\xi_3=1.27^{+0.42}_{-0.22}$). Furthermore, our analysis favors the scalar index $n_s$ to be slightly larger than 1, compatible with some hybrid inflation models, as well as a significantly larger optical depth $\tau$ than the standard Planck value, indicating an earlier onset of reionization.

In this work, we analyze the observational properties of static, spherically symmetric boson stars with fourth and sixth-order self-interactions, using the Julia-based general-relativistic radiative transfer code Skylight. We assume the boson stars are surrounded by an optically thick, geometrically thin accretion disk. We use the Novikov-Thorne model to compute the energy flux, introducing a physically based accretion model around these boson star configurations. Additionally, we calculate the relativistic broadening of emission lines, incorporating a lamppost corona model with full relativistic effects for the first time around a boson star. Our results show distinct observational features between quartic-potential boson stars and Schwarzschild black holes, owing to the presence of stable circular orbits at all radii around the former. On the other hand, compact solitonic boson stars, which possess an innermost stable circular orbit, have observational features closely similar to black holes. This similarity emphasizes their potential as black-hole mimickers. However, the compact boson stars, lacking an event horizon, have complex light-ring structures that produce potentially observable differences from black holes with future generations of experiments.

Comparing different dark matter (DM) models, we explore the DM influence on black hole (BH) accretion disk physics, considering corotating and counterrotating thick accretion tori orbiting a central spinning BH. Our results identify accretion onto a central BH as a good indicator of DM presence, signaling possible DM tracers in accretion physics. We analyze accretion around a spinning BH immersed in perfect-fluid dark matter, cold dark matter and scalar field dark matter. Our investigation addresses observational evidence of distinctive DM effects on toroidal accretion disks and protojet configurations, proving that BH accretion tori immersed in DM can present characteristics, such as interdisk cusp or double tori, which have usually been considered as tracers for superspinars and naked singularity attractors. Therefore, in this context DM influence on the BH geometry could manifest as superspinar mimickers. DM also affects the central spinning attractor energetics associated with accretion physics, and its influence on accretion disks can be searched for in a variation of the central BH energetics as an increase of the mass accretion rates.