Papers for Thursday, Feb 23 2023

We investigate the dependence of elemental abundances on physical constants, and the implications this has for the distribution of complex life for various proposed habitability criteria. We consider three main sources of abundance variation: differing supernova rates, alpha burning in massive stars, and isotopic stability, and how each affects the metal-to-rock ratio and the abundances of carbon, oxygen, nitrogen, phosphorus, sulfur, silicon, magnesium, and iron. Our analysis leads to several predictions for which habitability criteria are correct by determining which ones make our observations of the physical constants, as well as a few other observed features of our universe, most likely. Our results indicate that carbon-rich or carbon-poor planets are uninhabitable, slightly magnesium-rich planets are habitable, and life does not depend on nitrogen abundance too sensitively. We also find suggestive but inconclusive evidence that metal-rich planets and phosphorus-poor planets are habitable. These predictions can then be checked by probing regions of our universe that closely resemble normal environments in other universes. If any of these predictions are found to be wrong, the multiverse scenario would predict that the majority of observers are born in universes differing substantially from ours, and so can be ruled out, to varying degrees of statistical significance.

We investigate the nucleosynthesis and kilonova properties of binary neutron-star (NS) merger models which lead to intermediate remnant lifetimes of ~0.1-1seconds until black-hole (BH) formation and describe all components of material ejected during the dynamical merger phase, NS-remnant evolution, and final viscous disintegration of the BH torus after gravitational collapse. To this end we employ a combination of hydrodynamics, nucleosynthesis, and radiative-transfer tools to achieve a consistent end-to-end modeling of the system and its observables. We adopt a novel version of the Shakura-Sunyaev scheme allowing to vary the approximate turbulent viscosity inside the NS remnant independently of the surrounding disk. We find that asymmetric progenitors lead to shorter remnant lifetimes and enhanced ejecta masses, although the viscosity affects the absolute values of these characteristics. The integrated production of lanthanides and heavier elements in such binary systems is sub-solar, suggesting that the considered scenarios contribute in a sub-dominant fashion to r-process enrichment. One reason is that BH-tori formed after delayed collapse exhibit less neutron-rich conditions than typically found, and often assumed in previous BH-torus models, for early BH formation. The outflows in our models feature strong anisotropy as a result of the lanthanide-poor polar neutrino-driven wind pushing aside lanthanide-rich dynamical ejecta. Considering the complexity of the models, the estimated kilonova light curves show promising agreement with AT2017gfo after times of several days, while the remaining inconsistencies at early times could possibly be overcome in binary configurations with a more dominant neutrino-driven wind relative to the dynamical ejecta.

Astrophysical transients with rapid development on sub-hour timescales are intrinsically rare. Due to their short durations, events like stellar superflares, optical flashes from gamma-ray bursts, and shock breakouts from young supernovae are difficult to identify on timescales that enable spectroscopic followup. This paper presents the Evryscope Fast Transient Engine (EFTE), a new data reduction pipeline designed to provide low-latency transient alerts from the Evryscopes, a North-South pair of ultra-wide-field telescopes with an instantaneous footprint covering 38% of the entire sky, and tools for building long-term light curves from Evryscope data. EFTE leverages the optical stability of the Evryscopes by using a simple direct image subtraction routine suited to continuously monitoring the transient sky at minute cadence. Candidates are produced within the base Evryscope two-minute cadence for 98.5% of images, and internally filtered using VetNet, a convolutional neural network real-bogus classifier. EFTE provides an extensible, robust architecture for transient surveys probing similar timescales, and serves as the software testbed for the real-time analysis pipelines and public data distribution systems for the Argus Array, a next generation all-sky observatory with a data rate 62x higher than Evryscope.

The imminent launch of XRISM will usher in an era of high-resolution X-ray spectroscopy. For active galactic nuclei (AGN) this is an exciting epoch that is full of massive potential for uncovering the ins and outs of supermassive black hole accretion. In this work, we review AGN research topics that are certain to advance in the coming years with XRISM and prognosticate the possibilities with Athena and Arcus. Specifically, our discussion focuses on: (i) the relatively slow moving ionised winds known as warm absorbers and obscurers; (ii) the iron emitting from different regions of the inner and outer disc, broad line region, and torus; and (iii) the ultrafast outflows that may be the key to understanding AGN feedback.

The relationship linking a galaxy cluster's total mass with the concentration of its mass profile and its redshift is a fundamental prediction of the Cold Dark Matter (CDM) paradigm of cosmic structure formation. However, confronting those predictions with observations is complicated by the fact that simulated clusters are not representative of observed samples where detailed mass profile constraints are possible. In this work, we calculate the Symmetry-Peakiness-Alignment (SPA) morphology metrics for maps of X-ray emissivity from THE THREE HUNDRED project hydrodynamical simulations of galaxy clusters at four redshifts, and thereby select a sample of morphologically relaxed, simulated clusters, using observational criteria. These clusters have on average earlier formation times than the full sample, confirming that they are both morphologically and dynamically more relaxed than typical. We constrain the concentration-mass-redshift relation of both the relaxed and complete sample of simulated clusters, assuming power-law dependences on mass ($\kappa_m$) and $1+z$ ($\kappa_\zeta$), finding $\kappa_m = -0.12 \pm 0.07$ and $\kappa_\zeta = -0.27 \pm 0.19$ for the relaxed subsample. From an equivalently selected sample of massive, relaxed clusters observed with ${\it Chandra}$, we find $\kappa_m = -0.12 \pm 0.08$ and $\kappa_\zeta = -0.48 \pm 0.19$, in good agreement with the simulation predictions. The simulated and observed samples also agree well on the average concentration at a pivot mass and redshift providing further validation of the $\Lambda$CDM paradigm in the properties of the largest gravitationally collapsed structures observed. This also represents the first clear detection of decreasing concentration with redshift, a longstanding prediction of simulations, in data.

We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-bolometric luminosity of log L$_{\text{max}} = 45.7$ erg s$^{-1}$ from slowly fading ztf-$\it{g}$ and ztf-$r$ light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ($g > 23.3$ mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra from Lick and Keck Observatories exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si III] $\lambda$1892, C III] $\lambda$1909, and C II] $\lambda$2325. Typical nebular lines in AGN spectra from ions such as [O II] and [O III] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity, and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a $\approx 14 M_{\odot}$ star and a supermassive black hole of mass $M_{\text{BH}} \sim$ $10^{8} M_{\odot}$. Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy.

We select and characterise a sample of massive (log(M$_{*}/$M$_{\odot})>10.6$) quiescent galaxies (QGs) at $3<z<5$ in the latest COSMOS2020 catalogue. QGs are selected using a new rest-frame colour selection method, based on their probability of belonging to the quiescent group defined by a Gaussian Mixture Model (GMM) trained on rest-frame colours ($NUV-U, U-V, V-J$) of similarly massive galaxies at $2<z<3$. We calculate the quiescent probability threshold above which a galaxy is classified as quiescent using simulated galaxies from the SHARK semi-analytical model. We find that at $z\geq3$ in SHARK, the GMM/$NUVU-VJ$ method out-performs classical rest-frame $UVJ$ selection and is a viable alternative. We select galaxies as quiescent based on their probability in COSMOS2020 at $3<z<5$, and compare the selected sample to both $UVJ$ and $NUVrJ$ selected samples. We find that although the new selection matches $UVJ$ and $NUVrJ$ in number, the overlap between colour selections is only $\sim50-80\%$, implying that rest-frame colour commonly used at lower redshifts selections cannot be equivalently used at $z>3$. We compute median rest-frame SEDs for our sample and find the median quiescent galaxy at $3<z<5$ has a strong Balmer/4000 Angstrom break, and residual $NUV$ flux indicating recent quenching. We find the number densities of the entire quiescent population (including post-starbursts) more than doubles from $3.5\pm2.2\times10^{-6}$ Mpc$^{-3}$ at $4<z<5$ to $1.4\pm0.4\times10^{-5}$ Mpc$^{-3}$ at $3<z<4$, confirming that the onset of massive galaxy quenching occurs as early as $3<z<5$.

It is well known that supermassive black holes (SMBHs) and their host galaxies co-evolve. A manifestation of this co-evolution is the correlation that has been found between the SMBH mass, M$_{BH}$, and the galaxy bulge or stellar mass, M$_*$. The cosmic evolution of this relation, though, is still a matter of debate. In this work, we examine the M$_{BH}-$M$_*$ relation, using 687 X-ray luminous (median $\rm log\,[L_{X,2-10keV}(ergs^{-1})]=44.3$), broad line AGN, at $\rm 0.2<z<4.0$ (median $\rm z\approx 1.4$) that lie in the XMM-{\it{XXL}} field. Their M$_{BH}$ and M$_*$ range from $\rm 7.5<log\,[M_{BH}\,(M_\odot)]<9.5$ and $\rm 10<log\,[M_*(M_\odot)]<12$, respectively. Most of the AGN live in star-forming galaxies and their Eddington ratios range from 0.01 to 1, with a median value of 0.06. Our results show that M$_{BH}$ and M$_*$ are correlated ($\rm r=0.47\pm0.21$, averaged over different redshift intervals). Our analysis also shows that the mean ratio of the M$_{BH}$ and M$_*$ does not evolve with redshift, at least up to $\rm z=2$ and has a value of $\rm log($M$_{BH}/$M$_*)=-2.44$. The majority of the AGN ($75\%$) are in a SMBH mass growth dominant phase. In these systems, the M$_{BH}-$M$_*$ correlation is weaker and their M$_*$ tends to be lower (for the same M$_{BH}$) compared to systems that are in a galaxy mass growth phase. Our findings suggest that the growth of black hole mass occurs first, while the early stellar mass assembly may not be so efficient.

Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass loss rates. In High Mass X-ray Binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform spectro-timing analysis of the HMXB OAO 1657-415 with non-simultaneous NuSTAR and NICER observations. We compute the hardness ratio from the energy-resolved light curves, and using an adaptive rebinning technique, we thus select appropriate time segments to search for rapid spectral variations on timescales of a few hundreds to thousands of seconds. Column density and intensity of Iron K$\alpha$ line were strongly correlated, and the recorded spectral variations were consistent with accretion from a clumpy wind. We also illustrate a novel framework to measure clump sizes, masses in HMXBs more accurately based on absorption measurements and orbital parameters of the source. We then discuss the limitations posed by current X-ray spacecrafts in such measurements and present prospects with future X-ray missions. We find that the source pulse profiles show a moderate dependence on energy. We identify a previously undetected dip in the pulse profile visible throughout the NuSTAR observation near spin phase 0.15 possibly caused by intrinsic changes in accretion geometry close to the neutron star. We do not find any evidence for the debated cyclotron line at $\sim$ 36\,keV in the time-averaged or the phase-resolved spectra with NuSTAR.

Using a single gravitational lens system observed at $\lesssim5$ milli-arcsecond resolution with very long baseline interferometry (VLBI), we place a lower bound on the mass of the fuzzy dark matter (FDM) particle, ruling out $m_\chi \leq 4.4\times10^{-21}~\mathrm{eV}$ with a 20:1 posterior odds ratio relative to a smooth lens model. We generalize our result to non-scalar and multiple-field models, such as vector FDM, with $m_{\chi,\mathrm{vec}} > 1.4 \times 10^{-21}~\mathrm{eV}$. Due to the extended source structure and high angular resolution of the observation, our analysis is directly sensitive to the presence of granule structures in the main dark matter halo of the lens, which is the most generic prediction of FDM theories. A model based on well-understood physics of ultra-light dark matter fields in a gravitational potential well makes our result robust to a wide range of assumed dark matter fractions and velocity dispersions in the lens galaxy. Our result is competitive with other lower bounds on $m_\chi$ from past analyses, which rely on intermediate modelling of structure formation and/or baryonic effects. Higher resolution observations taken at 10 to 100 GHz could improve our constraints by up to 2 orders of magnitude in the future.

We investigate the formation of brightest cluster galaxies (BCGs) in the TNG300 cosmological simulation of the IllustrisTNG project. Our cluster sample consists of 700 haloes with $M_{200} \geq 5 \times 10^{13} \, \mathrm{M}_{\odot}$ at $z=0$, along with their progenitors at earlier epochs. This includes 280 systems with $M_{200} \geq 10^{14} \, \mathrm{M}_{\odot}$ at $z=0$, as well as three haloes with $M_{200} \geq 10^{15} \, \mathrm{M}_{\odot}$. We find that the stellar masses and star formation rates of our simulated BCGs are in good agreement with observations at $z \lesssim 0.4$, and that they have experienced, on average, $\sim$2 ($\sim$3) major mergers since $z=1$ ($z=2$). Separating the BCG from the intracluster light (ICL) by means of a fixed 30 kpc aperture, we find that the fraction of stellar mass contributed by ex situ (i.e. accreted) stars at $z=0$ is approximately 70, 80, and 90 per cent for the BCG, BCG+ICL, and ICL, respectively. Tracking our simulated BCGs back in time using the merger trees, we find that they became dominated by ex situ stars at $z \sim $1-2, and that half of the stars that are part of the BCG at $z=0$ formed early ($z \sim 3$) in other galaxies, but `assembled' onto the BCG until later times ($z \approx 0.8$ for the whole sample, $z \approx 0.5$ for BCGs in $M_{200} \geq 5 \times 10^{14} \, \mathrm{M}_{\odot}$ haloes). Finally, we show that the stellar mass profiles of BCGs are often dominated by ex situ stars at all radii, with stars from major mergers being found closer to the centre, while stars that were tidally stripped from other galaxies dominate the outer regions.

Reflected sunlight observations from the Ultraviolet Spectrograph (UVS) on the Juno spacecraft were used to study the distribution of acetylene (C$_2$H$_2$) at Jupiter's south pole. We find that the shape of the C$_2$H$_2$ absorption feature varies significantly across the polar region, and this can be used to infer spatial variability in the C$_2$H$_2$ abundance. There is a localized region of enhanced C$_2$H$_2$ absorption which coincides with the location of Jupiter's southern polar aurora; the C$_2$H$_2$ abundance poleward of the auroral oval is a factor of 3 higher than adjacent quiescent, non-auroral longitudes. This builds on previous infrared studies which found enhanced C$_2$H$_2$ abundances within the northern auroral oval. This suggests that Jupiter's upper-atmosphere chemistry is being strongly influenced by the influx of charged auroral particles and demonstrates the necessity of developing ion-neutral photochemical models of Jupiter's polar regions.

Planetary magnetic fields can affect the predicted mass loss rate for close-in planets that experience large amounts of UV irradiation. In this work, we present a method to detect the magnetic fields of close-in exoplanets undergoing atmospheric escape using transit spectroscopy at the 10830 Angstrom line of helium. Motivated by previous work on hydrodynamic and magneto-hydrodynamic photoevaporation, we suggest that planets with magnetic fields that are too weak to control the outflow's topology lead to blue-shifted transits due to day-to-night-side flows. In contrast, strong magnetic fields prevent this day-to-night flow, as the gas is forced to follow the magnetic field's roughly dipolar topology. We post-process existing 2D photoevaporation simulations to test this concept, computing synthetic transit profiles in helium. As expected, we find that hydrodynamically dominated outflows lead to blue-shifted transits on the order of the sound speed of the gas. Strong surface magnetic fields lead to unshifted or slightly red-shifted transit profiles. High-resolution observations can distinguish between these profiles; however, eccentricity uncertainties generally mean that we cannot conclusively say velocity shifts are due to the outflow for individual planets. The majority of helium observations are blue-shifted, which could be a tentative indication that close-in planets generally have surface dipole magnetic field strengths $\lesssim 0.1$ gauss. More 3D hydrodynamic and magneto-hydrodynamic are needed to confirm this conclusion robustly.

We undertook a search for new dwarf galaxies in the Leo-I group using the data from the DECaLS digital sky survey. Five new presumed members of this group have been found in a wide vicinity of ${\rm M}\,105 ({\rm NGC}\,3379$). Currently, the group has a population of $83$ galaxies, $33$ of which have measured radial velocities. More than half of the group members belong to early types with no signs of ongoing star formation. About a quarter of the galaxies are outside the group's virial radius, $R_v = 385$~kpc. The presence of multiple systems with a size of about 15~kpc is evident in the group, but there are no noticeable global flat or filamentary substructures. The luminosity function of the group looks to be deficient in galaxies with absolute magnitudes in the interval $M_B = [-18, -15]$ mag. The ${\rm M}\,105$ group is characterized by a radial velocity dispersion of $136$~km~s$^{-1}$, orbital mass estimate $(5.76\pm1.32)\times 10^{12}~M_{\odot}$, and the total mass-to-K-band-luminosity ratio $(17.8\pm4.1) M_\odot/L_\odot$. The neighboring group of galaxies around ${\rm M}\,66 ({\rm NGC}\,3627$) has a similar virial radius, $390$~kpc, velocity dispersion, $135$~km~s$^{-1}$, and total mass-to-luminosity ratio, $(15.6\pm3.9) M_\odot/L_\odot$. Both groups in the Leo constellation are approaching the Local Group with a velocity of about 100~km~s$^{-1}$. In the background of the ${\rm M}\,105$ group, we noted a group of 6 galaxies with an unusually low virial mass-to-luminosity ratio, $M_T/L_K = (4.1\pm2.2) M_\odot/L_\odot$.

Vela X-1 is the archetypical eclipsing high-mass X-ray binary, composed of a neutron star (NS) accreting the B-star wind. It was observed by nearly all X-ray observatories, often multiple times, featuring a rich spectrum of variable emission lines. Yet, the precise origin of these lines in the binary system remains uncertain. We perform a systematic, orbital phase-dependent analysis of the reflected Fe K$\alpha$ fluorescence line at 6.4 keV using over 100 NICER observations. We resolve the line variability into 500s time bins and find that it is predominantly due to variation in the ionizing flux, with a moderate underlying phase dependence over the 9-day orbital period. Our analysis reveals a significant reflection component that cannot originate from the companion B-star alone. We also find that an appreciable portion of the B-star surface is obscured opposite the eclipse, and this obscuration is not symmetric around the mid-point (phase=0.5). We argue that an accretion stream, from the B-star to the NS and distorted by the orbital motion, is responsible for both the additional fluorescence emission component and for obscuring the B-star.

We present search results of 22 high latitude (b > 25 deg.) sightlines for OH 18-cm emission using the 305-m radio telescope at the Arecibo Observatory. These sightlines appear in neutral hydrogen emission at intermediate velocities (V_lsr values ranging from -90 to -20 km/s) and are predicted to have a sufficient molecular composition so as to be detectable in molecular emission. Such objects, known as Intermediate-Velocity Molecular Clouds (IVMCs), have historically been detected through 12CO emission. Recent studies indicate that IVMCs may be widespread in the Galaxy and have important implications for models of the interstellar medium and star formation. However, we report non-detections of OH emission toward the 22 sightlines and provide stringent upper limits on the OH column density. Using available HI and Av data in combination with existing state-of-the-art PDR models, we estimate H2 column densities and find that they are more than an order of magnitude lower than the predicted values. We also find that the hydrogen volume density of these clouds is less than roughly 25 per cubic centimeter. In addition, we discuss the known IVMCs with previous 12CO detections in the context of the PDR models. Our analysis of these clouds indicates that the structure of molecular material in IVMCs is morphologically clumpy. These results motivate the need for future sensitive, on-the-fly searches (rather than targeted searches) for CO emission from IVMCs with of order roughly 1' resolution. High angular resolution (1') HI and Av data will also be helpful to better constrain the structure and composition of IVMCs.

We present an overview of the capabilities and key algorithms employed in the so-called eMPT software suite developed for planning scientifically optimized, multi-object spectroscopic (MOS) observations with the Micro-Shutter Array (MSA) of the Near-Infrared Spectrograph (NIRSpec) instrument on board the James Webb Space Telescope (JWST), the first multi-object spectrograph to operate in space. NIRSpec MOS mode is enabled by a programmable MSA, a regular grid of ~250,000 individual apertures that projects to a static, semi-regular pattern of available slits on the sky and makes the planning and optimization of an MSA observation a rather complex task. As such, the eMPT package is offered to the NIRSpec user community as a supplement to the MSA Planning Tool (MPT) included in the STScI Astronomer's Proposal Tool (APT) to assist in the planning of NIRSpec MOS proposals requiring advanced functionality to meet ambitious science goals. The eMPT produces output that can readily be imported and incorporated into the user's observing program within the APT to generate a customized MPT MOS observation. Furthermore, its novel algorithms and modular approach make it highly flexible and customizable, providing users the option to finely control the workflow and even insert their own software modules to tune their MSA slit masks to the particular scientific objectives at hand.

Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Methods (abridged): We used a 500 Mpc/h large constrained simulation of the LU with initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalogue and follow galaxy formation physics directly in the hydro-dynamical simulations to base the comparison on stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 Mpc/h large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that haloes and galaxies in our constrained simulation trace the local DM density field very differently. Thereby, this simulation reproduces the observed 50% under-density of galaxy clusters and groups within the sphere of ~100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ~1.5$\sigma$ feature. At the same time, the simulation reproduces the observed over-density of massive galaxy clusters within the same sphere, which on its own also corresponds to a ~1.5$\sigma$ feature. Interestingly, we find that only 44 out of 15635 random realizations (i.e. 0.28%) are matching both anomalies, making the LU to be a ~3$\sigma$ environment. We finally compared a mock galaxy catalogue with the observed distribution of galaxies in the LU, finding also a match to the observed factor of two over-density at ~16 Mpc as well as the observed 15% under-density at ~40 Mpc distance.

Active comets have been detected in several exoplanetary systems, although so far only indirectly, when the dust or gas in the extended coma has transited in front of the stellar disk. The large optical surface and relatively high temperature of an active cometary coma also makes it suitable to study with direct imaging, but the angular separation is generally too small to be reachable with present-day facilities. However, future imaging facilities with the ability to detect terrestrial planets in the habitable zones of nearby systems will also be sensitive to exocomets in such systems. Here we examine several aspects of exocomet imaging, particularly in the context of the Large Interferometer for Exoplanets (LIFE), which is a proposed space mission for infrared imaging and spectroscopy through nulling interferometry. We study what capabilities LIFE would have for acquiring imaging and spectroscopy of exocomets, based on simulations of the LIFE performance as well as statistical properties of exocomets that have recently been deduced from transit surveys. We find that for systems with extreme cometary activities such as beta Pictoris, sufficiently bright comets may be so abundant that they overcrowd the LIFE inner field of view. More nearby and moderately active systems such as epsilon Eridani or Fomalhaut may turn out to be optimal targets. If the exocomets have strong silicate emission features, such as in comet Hale-Bopp, it may become possible to study the mineralogy of individual exocometary bodies. We also discuss the possibility of exocomets as false positives for planets, with recent deep imaging of alpha Centauri as one hypothetical example. Such contaminants could be common, primarily among young debris disk stars, but should be rare among the main sequence population. We discuss strategies to mitigate the risk of any such false positives.

Context: A detailed study of the close interacting binary V4142\,Sgr based on photometric and spectroscopic analysis is presented.This system belongs to the enigmatic class of Algol-like variables showing a long photometric cycle of unknown nature. Aims: Performing photometric data-mining and spectroscopic observations covering the orbital cycle, we obtain the orbital parameters and the stellar properties of the binary system, along with the physical properties of the accretion disk located around the hot star. Insights on the evolutive path of the system are obtained. Methods: The light curve was modeled through an inverse modeling method using a theoretical light curve of the binary system, considering the light curve contribution of both stars and the accretion disk of the hot star to obtain the fundamental parameters. To constrain the main stellar parameters the mass ratio was fixed, as well as the donor temperature using the obtained values from our spectroscopic analysis including deblending methods to isolate the spectral lines of the stellar components. The system parameters were compared with a grid of binary star evolutive models in order to get insights on the evolutionary history of the system. Results: The orbital period and the long cycle were re-calculated and found to be of $30.633 \pm 0.002 ~\mathrm{days}$ and $1201 \pm 14 ~\mathrm{days}$. The spectral analysis reveals H$\alpha$ double emission with a persistent $V \leq R$ asymmetry which is considered evidence of a possible wind emergin from the hotspot region...

Telescopes measuring cosmic microwave background (CMB) polarization on large angular scales require exquisite control of systematic errors to ensure the fidelity of the cosmological results. In particular, far-sidelobe contamination from wide angle scattering is a potentially prominent source of systematic error for large aperture microwave telescopes. Here we describe and demonstrate a ray-tracing-based modeling technique to predict far sidelobes for a Three Mirror Anistigmat (TMA) telescope designed to observe the CMB from the South Pole. Those sidelobes are produced by light scattered in the receiver optics subsequently interacting with the walls of the surrounding telescope enclosure. After comparing simulated sidelobe maps and angular power spectra for different enclosure wall treatments, we propose a highly scattering surface that would provide more than an order of magnitude reduction in the degree-scale far-sidelobe contrast compared to a typical reflective surface. We conclude by discussing the fabrication of a prototype scattering wall panel and presenting measurements of its angular scattering profile.

Context: The surface properties of rotating stars can vary from pole to equator, resulting in anisotropic stellar winds which are not included in the currently available evolutionary models. Aims: We develop a formalism to describe the mass and angular momentum loss of rotating stars which takes into account both the varying surface properties and distortion due to rotation. Methods: Adopting the mass-loss recipe for non-rotating stars, we assigned to each point on the surface of a rotating star an equivalent non-rotating star, for which the surface mass flux is given by the recipe. The global mass-loss and angular momentum loss rates are then given by integrating over the deformed stellar surface as appropriate. Evolutionary models were computed and our prescription is compared to the currently used simple mass-loss enhancement recipes for rotating stars. Results: We find that mass-loss rates are largely insensitive to rotation for models not affected by the bi-stability jump. For those affected by the bi-stability jump, the increase in mass-loss rates with respect to time is smoothed. As our prescription considers the variation of physical conditions over the stellar surface, the region affected by the bi-stability jump is able to grow gradually instead of the whole star suddenly being affected. Conclusion: We have provided an easy to implement and flexible, yet physically meaningful prescription for calculating mass and angular momentum loss rates of rotating stars in a one-dimensional stellar evolution code which compares favourably to more physically comprehensive models. The implementation of our scheme in the stellar evolution code MESA is available online: https://zenodo.org/record/7437006

The discovery of cosmic acceleration motivated extensive studies of dynamical dark energy and modified gravity models. Of particular interest are the scalar-tensor theories, with a scalar field dark energy non-minimally coupled to matter. Cosmological constraints on these models often employ the quasi-static approximation (QSA), in which the dynamics of the scalar field perturbations is proportional to the perturbation in the matter density. Using the QSA simplifies the physical interpretation of the phenomenology of scalar-tensor theories, and results in substantial savings of computing time when deriving parameter constraints. Focusing on the symmetron model, which is a well-motivated scalar-tensor theory with a screening mechanism, we compare the exact solution of the linearly perturbed field equations to those obtained under the QSA and identify the range of the model parameters for which the QSA is valid. We find that the evolution of background scalar field is most important, namely, whether it is dominated by the Hubble friction or the scalar field potential. This helps us derive a criterion for the symmetron model, but same argument can be applied to other scalar-tensor theories of generalized Brans-Dicke type. We consider two scenarios, one where the scalar field is only coupled to dark matter and where it couples to all of the matter.

The chemical compositions of Earth's core and mantle provide insight into the processes that led to their formation. N-body simulations, on the other hand, generally do not contain chemical information, and seek to only reproduce the masses and orbits of the terrestrial planets. These simulations can be grouped into four potentially viable scenarios of Solar System formation (Classical, Annulus, Grand Tack, and Early Instability) for which we compile a total of 433 N-body simulations. We relate the outputs of these simulations to the chemistry of Earth's core and mantle using a melt-scaling law combined with a multi-stage model of core formation. We find the compositions of Earth analogs to be largely governed by the fraction of equilibrating embryo cores and the initial embryo masses in N-body simulations. Simulation type may be important when considering magma ocean lifetimes, where Grand Tack simulations have the largest amounts of material accreted after the last giant impact. However, we cannot rule out any accretion scenarios or initial embryo masses due to the sensitivity of Earth's mantle composition to different parameters and the stochastic nature of N-body simulations. Comparing the last embryo impacts experienced by Earth analogs to specific Moon-forming scenarios, we find the characteristics of the Moon-forming impact are dependent on the initial conditions in N-body simulations where larger initial embryo masses promote larger and slower Moon-forming impactors. Mars-sized initial embryos are most consistent with the canonical hit-and-run scenario onto a solid mantle. Our results suggest that constraining the fraction of equilibrating impactor core and the initial embryo masses in N-body simulations could be significant for understanding both Earth's accretion history and characteristics of the Moon-forming impact.

Young $\delta$ Scuti stars have proven to be valuable asteroseismic targets but obtaining robust uncertainties on their inferred properties is challenging. We aim to quantify the random uncertainties in grid-based modelling of $\delta$ Sct stars. We apply Bayesian inference using nested sampling and a neural network emulator of stellar models, testing our method on both simulated and real stars. Based on results from simulated stars we demonstrate that our method can recover plausible posterior probability density estimates while accounting for both the random uncertainty from the observations and neural network emulation. We find that the posterior distributions of the fundamental parameters can be significantly non-Gaussian, multi-modal, and have strong covariance. We conclude that our method reliably estimates the random uncertainty in the modelling of $\delta$ Sct stars and paves the way for the investigation and quantification of the systematic uncertainty.

[abridged] ALMA observations of dust in protoplanetary disks are revealing the existence of sub-structures such as rings, gaps and cavities. Such morphology are expected to be the outcome of dynamical interaction between the disk and planets. However, other mechanisms are able to produce similar dust sub-structures. A solution is to look at the perturbation induced by the planet to the gas surface density and/or to the kinematics. In the case of the disk around AS 209, a prominent gap has been reported in the surface density of CO at $r \sim 100\,$au. Recently, Bae et al. (2022) detected a localized velocity perturbation in the $^{12}$CO $J=2-1$ emission along with a clump in $^{13}$CO $J=2-1$ at nearly 200 au, interpreted as a gaseous circumplanetary disk. We report a new analysis of ALMA archival observations of $^{12}$CO and $^{13}$CO J=2-1. A clear kinematics perturbation (kink) is detected in multiple channels and over a wide azimuth range in both dataset. We compared the observed perturbation with a semi-analytic model of velocity perturbations due to planet-disk interaction. The observed kink is not consistent with a planet at 200\,au as this would require a low gas disk scale height ($< 0.05$) in contradiction with previous estimate ($h/r \sim 0.118$ at $r = 100$ au). When we fix the disk scale height to 0.118 (at $r = 100$ au) we find instead that a planet of 3-5 M$_{\rm Jup}$ at 100 au induces a kinematics perturbation similar to the observed one. Thus, we conclude that a giant protoplanet orbiting at $r \sim 100\,$au is responsible of the large scale kink as well as of the perturbed dust and gas surface density previously detected. The position angle of the planet is constrained to be between 60$^{\circ}$-100$^{\circ}$. Future observations with high contrast imaging technique in the near- and mid- infrared are needed to confirm the presence and position of such a planet.

Superfluid vortex avalanches are one plausible cause of pulsar glitch activity. If they occur according to a state-dependent Poisson process, the measured long-term glitch rate is determined by the spin-down rate of the stellar crust, $\dot{\Omega}_{\rm c}$, and two phenomenological parameters quantifying the vortex-nucleus pinning force: a crust-superfluid angular velocity lag threshold, $X_{\rm cr}$, and a reference unpinning rate, $\lambda_0$. A Bayesian analysis of 541 glitches in 177 pulsars, with $N_{\rm g} \geq 1$ events per pulsar, yields $X_{\rm cr} = 0.15^{+0.09}_{-0.04} \, {\rm rad \, s^{-1}}$, $\lambda_{\rm ref} = 7.6^{+3.7}_{-2.6} \times 10^{-8} \, {\rm s^{-1}}$, and $a = -0.27^{+0.04}_{-0.03}$ assuming the phenomenological rate law $\lambda_0 = \lambda_{\rm ref} [\tau/(1 \, {\rm yr})]^a$, where $\tau$ denotes the characteristic spin-down age. The results are broadly similar, whether one includes or excludes quasiperiodic glitch activity, giant glitches, or pulsars with $N_{\rm g}=0$, up to uncertainties about the completeness of the sample and the total observation time per pulsar. The $X_{\rm cr}$ and $\lambda_0$ estimates are consistent with first-principles calculations based on nuclear theory, e.g. in the semiclassical local density approximation.

Using the MEGASIM, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and next-generation ground-based astronomical surveys. Our high-fidelity Earth Trojan Asteroid (ETA) distribution maps show never-before-seen high-resolution spatial features that evolve over timescales up to 1 Gyr. The simulation was synchronized to two surveys, 1) the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) and 2) the Zwicky Transient Facility (ZTF) public and partnership survey. Upper limits were calculated for the ETA population for ZTF (with assumed null detection) and LSST's cadence simulations. For LSST, a null detection is within three sigma of the simulation results out to H=19 for the twilight survey and H=20 for the baseline survey. Due to the Yarkovsky Effect, no ETAs are stable on billion year timescales are likely to be detected. ETAs large enough to remain stable on billion year timescales are very rare relative to the rest of the ETA population. A null detection in LSST will restrict that population to tens of objects larger than 100 meters. The null detection by ZTF to date has already done so.

We used data from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) to study the incidence of AGN in continuum-selected galaxies at $z\sim3$. From optical and infrared imaging in the 24 deg$^{2}$ Spitzer HETDEX Exploratory Large Area (SHELA) survey, we constructed a sample of photometric-redshift selected $z\sim3$ galaxies. We extracted HETDEX spectra at the position of 716 of these sources and used machine learning methods to identify those which exhibited AGN-like features. The dimensionality of the spectra was reduced using an autoencoder, and the latent space was visualized through t-distributed stochastic neighbor embedding (t-SNE). Gaussian mixture models were employed to cluster the encoded data and a labeled dataset was used to label each cluster as either AGN, stars, high-redshift galaxies, or low-redshift galaxies. Our photometric redshift (photo-z) sample was labeled with an estimated $92\%$ overall accuracy, an AGN accuracy of $83\%$, and an AGN contamination of $5\%$. The number of identified AGN was used to measure an AGN fraction for different magnitude bins. The UV absolute magnitude where the AGN fraction reaches $50\%$ is $M_{UV} = -23.8$. When combined with results in the literature, our measurements of AGN fraction imply that the bright end of the galaxy luminosity function exhibits a power-law rather than exponential decline, with a relatively shallow faint-end slope for the $z\sim3$ AGN luminosity function.

Solar eruptions often show the rotation of filaments, which is a manifestation of the rotation of erupting magnetic flux rope (MFR). Such a rotation of MFR can be induced by either the torque exerted by a background shear-field component (which is an external cause) or the relaxation of the magnetic twist of the MFR (an internal cause). For a given chirality of the erupting field, both the external and internal drivers cause the same rotation direction. Therefore, it remains elusive from direct observations which mechanism yields the dominant contribution to the rotation. In this paper, we exploit a full MHD simulation of solar eruption by tether-cutting magnetic reconnection to study the mechanism of MFR rotation. In the simulation, the MFR's height-rotation profile suggests that the force by the external shear-field component is a dominant contributor to the rotation. Furthermore, the torque analysis confirms that it is also the only factor in driving the counterclockwise rotation. On the contrary, the Lorentz torque inside the MFR makes a negative effect on this counterclockwise rotation.

Recently, the B.O.A.T. ("brightest of all time") gamma-ray burst, dubbed GRB 221009A, was detected by various instruments. Unprecedentedly, the GRB presented very-high-energy (VHE, energy above 0.1 TeV) gamma-ray emission with energy extending above 10 TeV, as reported by the Large High Altitude Air Shower Observatory (LHAASO). We here demonstrate that the VHE and especially >10 TeV emission may originate from the internal hadronic dissipation of the GRB, without the need of invoking any exotic processes as suggested by some previous studies. We also discuss the constraints on the properties of the GRB ejecta from multiwavelength and multi-messenger observations, which favors a magnetically dominated GRB ejecta. The suggested Poynting-flux-dominated GRB ejecta in this work supports the Blandford & Znajek (BZ) mechanism as the possible central engine model of GRB.

We present observations of the four 2 Pi 3/2 J = 3/2 ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitation temperature (Tex) and optical depth, while our ATCA observations give optical depth only. We perform Gaussian decomposition using the Automated Molecular Excitation Bayesian line-fitting Algorithm 'AMOEBA' (Petzler, Dawson, and Wardle 2021) fitting all four transitions simultaneously with shared centroid velocity and width. We identify 109 features across 38 sightlines (including 58 detections along 27 sightlines with excitation temperature measurements). While the main lines at 1665 and 1667 MHz tend to have similar excitation temperatures (median Tex(main) difference = 0.6 K, 84% show Tex(main) difference < 2 K), large differences in the 1612 and 1720 MHz satellite line excitation temperatures show that the gas is generally not in LTE. For a selection of sightlines we compare our OH features to associated (on-sky and in velocity) HI cold gas components (CNM) identified by Nguyen et al. (2019) and find no strong correlations. We speculate that this may indicate an effective decoupling of the molecular gas from the CNM once it accumulates.

Multi-wavelength images from the Hubble Space Telescope covering the wavelength range 0.27-1.6 $\mu$m show that the central area of the nearby dwarf galaxy NGC4449 contains several tens of compact sources that are emitting in the hydrogen recombination line Pa$\beta$ (1.2818 $\mu$m) but are only marginally detected in H$\alpha$ (0.6563 $\mu$m) and undetected at wavelengths $\lambda\le$0.55 $\mu$m. An analysis of the spectral energy distributions (SEDs) of these sources indicates that they are likely relatively young star clusters heavily attenuated by dust. The selection function used to identify the sources prevents meaningful statistical analyses of their age, mass, and dust extinction distributions. However, these cluster candidates have ages $\sim$5-6 Myr and A$_V>$6 mag, according to their SED fits, and are extremely compact, with typical deconvolved radii of 1 pc. The dusty clusters are located at the periphery of dark clouds within the galaxy and appear to be partially embedded. Density and pressure considerations indicate that the HII regions surrounding these clusters may be stalled, and that pre-supernova feedback has not been able to clear the clusters of their natal cocoons. These findings are in potential tension with existing models that regulate star formation with pre-supernova feedback, since pre-supernova feedback acts on short timescales, $\lesssim$4 Myr, for a standard Stellar Initial Mass function. The existence of a population of dusty star clusters with ages $>$4 Myr, if confirmed by future observations, paints a more complex picture for the role of stellar feedback in controlling star formation.

We have investigated the origin of a sub-class of carbon-polluted white dwarfs (DQ) originally identified as the ``hot DQ" white dwarfs. These objects are relatively hot (10 000 < T_eff < 25 000 K), have markedly higher carbon abundance (C-enriched), are more massive (M > 0.8 M_Sun) than ordinary DQs (M ~ 0.6 M_Sun), and display high space velocities. Hence, despite their young appearance their kinematic properties are those of an old white dwarf population. The way out of this dilemma is to assume that they formed via the merging of two white dwarfs. In this paper we examine the observed characteristics of this population of ``C-enriched" DQ white dwarfs and confirm that nearly half of the 63 known objects have kinematic properties consistent with those of the Galactic thick disc or halo. We have also conducted population synthesis studies and found that the merging hypothesis is indeed compatible with observations. Studies of this sub-class of white dwarfs have important implications for our understanding of Type Ia Supernovae (SNeIa), commonly used to determine the expansion history of the universe, since the same formation channel applies to both kinds of objects. Hence probing the properties of these white dwarfs that failed to explode may yield important constraints to the modelling of the mechanisms leading to a thermonuclear runaway.

In light of the recent measurement of the nonzero Cosmic Microwave Background (CMB) polarization rotation angle from the Planck 2018 data, we explore the possibility that such a cosmic birefringence effect is induced by coupling a fermionic current with photons via a Chern-Simons-like term. We begin our discussion by rederiving the general formulae of the cosmic birefringence angle with correcting a mistake in the previous study. We then identify the fermions in the current as the left-handed electron neutrinos and asymmetric dark matter (ADM) particles, since the rotation angle is sourced by the number density difference between particles and antiparticles. For the electron neutrino case, with the value of the degeneracy parameter $\xi_{\nu_e}$ recently measured by the EMPRESS survey, we find a large parameter space which can explain the CMB photon polarization rotations. On the other hand, for the ADM solution, we consider two benchmark cases with $M_\chi = 5$ GeV and 5 keV. The former is the natural value of the ADM mass if the observed ADM and baryon asymmetry in the Universe are produced by the same mechanism, while the latter provides a warm DM candidate. In addition, we explore the experimental constraints from the CMB power spectra and the DM direct detections.

Using Zwicky Transient Facility (ZTF) data, I noticed that MGAB-V240 = PS1-3PI J185529.82+323017.8 showed two different states: regularly outbursting state with a cycle length of 12 d and standstills. I found that the regularly outbursting state was in fact a sequence of superoutburst and intervening normal outbursts comprising a 12 d supercycle. During one of the superoutbursts, superhumps with a period of 0.015824(9) d (=22.79 min) were detected in the ZTF time-resolved data. This period and behavior have confirmed that MGAB-V240 is an AM CVn-type object with the shortest known supercycle and the second known AM CVn star showing genuine standstills. The standstills in this system were interrupted by short drops and the system often brightened after these drops. This phenomenon can be explained by the accumulation of the transferred matter in the outer part of the disk during the drops. This phenomenon favors a constant mass-transfer from the secondary combined with the difficulty in maintaining the hot state in a helium disk rather than a temporary decrease of the mass-transfer rate as the cause of these drops. MGAB-V240 should be close to the border of the thermal instability of a helium disk, and the observed superhump period agrees very well with the activity sequence expected by the disk instability theory and the evolutionary sequence of AM CVn stars.

The Reionization Cluster Survey (RELICS) imaged 41 galaxy clusters with the Hubble Space Telescope (HST), in order to detect lensed and high-redshift galaxies. Each cluster was imaged to about 26.5 AB mag in three optical and four near-infrared bands, taken in two distinct visits separated by varying time intervals. We make use of the multiple near-infrared epochs to search for transient sources in the cluster fields, with the primary motivation of building statistics for bright caustic crossing events in gravitational arcs. Over the whole sample, we do not find any significant ($\gtrsim5 \sigma$) caustic crossing events, in line with expectations from semi-analytic calculations but in contrast to what may be naively expected from previous detections of some bright events, or from deeper transient surveys that do find high rates of such events. Nevertheless, we find six prominent supernova (SN) candidates over the 41 fields: three of them were previously reported and three are new ones reported here for the first time. Out of the six candidates, four are likely core-collapse (CC) SNe -- three in cluster galaxies, and among which only one was known before, and one slightly behind the cluster at $z\sim0.6-0.7$. The other two are likely Ia -- both of them previously known, one probably in a cluster galaxy, and one behind it at $z\simeq2$. Our study supplies empirical bounds for the rate of caustic crossing events in galaxy cluster fields to typical HST magnitudes, and lays the groundwork for a future SN rate study.

Mass-to-flux ratios measured via the Zeeman effect suggest the existence of a transition from a magnetically sub-critical state in HI clouds to a super-critical state in molecular clouds. However, due to projection, chemical, and excitation effects, Zeeman measurements are subject to a number of biases, and may not reflect the true relations between gravitational and magnetic energies. In this paper, we carry out simulations of the formation of magnetised molecular clouds, zooming in from an entire galaxy to sub-pc scales, which we post-process to produce synthetic HI and OH Zeeman measurements. The mass-to-flux ratios we recover from the simulated observations show a transition in magnetic criticality that closely matches observations, but we find that the gravitational-magnetic energy ratios on corresponding scales are mostly super-critical, even in the HI regime. We conclude that HI clouds in the process of assembling to form molecular clouds are already super-critical even before H_2 forms, and that the apparent transition from sub- to super-criticality between HI and H_2 is primarily an illusion created by chemical and excitation biases affecting the Zeeman measurements.

Line-intensity mapping (LIM) is a promising technique to constrain the global distribution of galaxy properties. To combine LIM experiments probing different tracers with traditional galaxy surveys and fully exploit the scientific potential of these observations, it is necessary to have a physically motivated modeling framework. As part of developing such a framework, in this work we introduce and model the conditional galaxy property distribution (CGPD), i.e. the distribution of galaxy properties conditioned on the host halo mass and redshift. We consider five galaxy properties, including the galaxy stellar mass, molecular gas mass, galaxy radius, gas phase metallicity and star formation rate (SFR), which are important for predicting the emission lines of interest. The CGPD represents the full distribution of galaxies in the five dimensional property space; many important galaxy distribution functions and scaling relations, such as the stellar mass function and SFR main sequence, can be derived from integrating and projecting it. We utilize two different kinds of cosmological galaxy simulations, a semi-analytic model and the IllustrisTNG hydrodynamic simulation, to characterize the CGPD and explore how well it can be represented using a Gaussian mixture model (GMM). We find that with just a few ($\sim 3$) Gaussian components, a GMM can describe the CGPD of the simulated galaxies to high accuracy for both simulations. The CGPD can be mapped to LIM or other observables by constructing the appropriate relationship between galaxy properties and the relevant observable tracers.

We study how the nonlocal effects of the mean electromotive force affects the mean-field solar type dynamo model. Following suggestion of Rheinhardt and Brandenburg (2012) we approximate the integro-differential equation for the mean electromotive force by the reaction-diffusion type equation. This generalization alleviates the scale separation approximation. Solution of the eigenvalue problem reveals a few curious properties of the dynamo model with the nonlocal mean electromotive force. Beside a decrease of the critical dynamo instability threshold, reported in earlier studies, there is an increase the dynamo periods of the unstable modes. Simultaneously, the nonlocal model shows substantially lower growth rate of the unstable dynamo modes in vicinity above the critical threshold than the model which employ the scale separation approximation. Also, for the nonlocal model we find a number of the different oscillating and steady dynamo modes can be excited in the close vicinity of threshold of the first unstable dynamo mode. We verify these findings using the nonlinear dynamo model. The model shows the Parker's dynamo wave solutions with the wave propagating from the mid latitude at the bottom of the convection zone toward the solar equator at the surface. In the weakly nonlinear regime the interference of the dynamo mode of different spatial localization results into the Grand activity cycles of period about 300 years.

Very metal-poor (VMP, [Fe/H]<-2.0) stars offer a wealth of information on the nature and evolution of elemental production in the early galaxy and universe. The upcoming China Space Station Telescope (CSST) will provide us with a large amount of spectroscopic data that may contain plenty of VMP stars, and thus it is crucial to determine the stellar atmospheric parameters ($T_{eff}$, $\log g$, and [Fe/H]) for low-resolution spectra similar to the CSST spectra (R~200). In this paper, a two-dimensional Convolutional Neural Network (CNN) model with three convolutional layers and two fully connected layers is constructed. The principal aim of this work is to measure the ability of this model to estimate stellar parameters on low-resolution (R~200) spectra and to identify VMP stars so that we can better search for VMP stars in the spectra observed by CSST.We mainly use 10,008 observed spectra of VMP stars from LAMOST DR3, and 16,638 spectra of common stars ([Fe/H]>-2.0) from LAMOST DR8 for the experiment and make comparisons. All spectra are reduced to R~200 to match the resolution of the CSST and are preprocessed and collapsed into two-dimensional spectra for input to the CNN model. The results show that the MAE values are 99.40 K for $T_{eff}$, 0.22 dex for $\log g$, 0.14 dex for [Fe/H], and 0.26 dex for [C/Fe], respectively. Besides, the CNN model efficiently identifies VMP stars with a precision of 94.77%. The validation and practicality of this model are also tested on the MARCS synthetic spectra. This paper powerfully demonstrates the effectiveness of the proposed CNN model in estimating stellar parameters for low-resolution spectra (R~200) and recognizing VMP stars that are of interest for stellar population and galactic evolution work.

All the scientific instruments of the Italian National Telescope "Galileo" (TNG). as well as the tracking systems and the Shack-Hartmann wavefront systems use CCDs as detectors. The CCD procurement is crucial to optimize the instruments performance, and, of course, equal importance assumes the design and realization of controllers able to drive various kinds of CCDs with a very low readout noise. Detectors characterization is of fundamental importance when they have to be used in scientific instrumentation. Various CCDs have been assembled in cryostats and tested at our laboratory in order to select the suitable detector for the optical instruments of the TNG. The relevant phases of the group activity are here described, as well as the commissioning at the telescope and the work that is in progress.

The present paper describes the construction, the installation and the operation of the Optical Imager Galileo (OIG), a scientific instrument dedicated to the 'imaging' in the visible. OIG was the first instrument installed on the focal plane of the Telescopio Nazionale Galileo (TNG) and it has been extensively used for the functional verification of several parts of the telescope (as an example the optical quality, the rejection of spurious light, the active optics and the tracking), in the same way also several parts of the TNG informatics system (instrument commanding, telemetry and data archiving) have been verified making extensive use of OIG. This paper provides also a frame of work for a further development of the imaging dedicated instrumentation inside TNG. OIG, coupled with the first near-IR camera (ARNICA), has been the 'workhorse instrument' during the first period of telescope experimental and scientific scheduling.

We combine observations from different vantage points to perform a detailed study of a long duration eruptive C7 class flare that occurred on 17 April 2021 and was partially occulted from Earth view. The dynamics and thermal properties of the flare-related plasma flows, the flaring arcade, and the energy releases and particle acceleration are studied together with the kinematic evolution of the associated CME in order to place this long duration event in context of previous eruptive flare studies. The flare showed hard X-ray (HXR) bursts over the duration of an hour in two phases lasting from 16:04 UT to 17:05 UT. During the first phase, a strong increase in emission from hot plasma and impulsive acceleration of the CME was observed. The CME acceleration profile shows a three-part evolution of slow rise, acceleration, and propagation in line with the first STIX HXR burst phase, which is triggered by a rising hot (14 MK) plasmoid. During the CME acceleration phase, we find signatures of ongoing magnetic reconnection behind the erupting structure, in agreement with the standard eruptive flare scenario. The subsequent HXR bursts that occur about 30 minutes after the primary CME acceleration show a spectral hardening (from $\delta \approx $ 7 to $\delta \approx $ 4) but do not correspond to further CME acceleration and chromospheric evaporation. Therefore, the CME-flare feedback relationship may only be of significance within the first 25 minutes of the event under study, as thereafter the flare and the CME eruption evolve independently of each other.

The magnetic inclination angle $\chi$, namely the angle between the spin and magnetic axes of a neutron star (NS), plays a vital role in its observational characteristics. However, there are few systematic investigations on its long-term evolution, especially for accreting NSs in binary systems. Applying the model of \citet{2021MNRAS.505.1775B} and the binary evolution code \mesa{}, we simultaneously simulate the evolution of the accretion rate, spin period, magnetic field, and magnetic inclination angle of accreting NSs in intermediate/low X-ray binaries (I/LMXBs). We show that the evolution of $\chi$ depends not only on the initial parameters of the binary systems, but also on the mass transfer history and the efficiency of pulsar loss. Based on the calculated results we present the characteristic distribution of $\chi$ for various types of systems including ultracompact X-ray binaries, binary millisecond pulsars, and ultraluminous X-ray sources, and discuss their possible observational implications.

The halo mass-temperature relation for a sample of 216 galaxy clusters, groups, and individual galaxies observed by $Chandra$ X-ray Observatory is presented. Using accurate spectral measurements of their hot atmospheres, we derive the $M-T$ relation for systems with temperatures ranging between 0.4-15.0 keV. We measure the total mass of clusters, groups, and galaxies at radius $R_{2500}$, finding that the $M_{2500} \propto T^{\alpha}$ relation follows a power-law with $\alpha$ = 1.65$\pm$0.06. Our relation agrees with recent lensing studies of the $M-T$ relation at $R_{200}$ and is consistent with self-similar theoretical prediction and recent simulations. This agreement indicates that the $M-T$ relation is weakly affected by non-gravitational heating processes. Using lensing masses within $R_{200}$ we find $M_{200}-T$ follows a power-law with slope 1.61$\pm$0.19, consistent with the $M_{2500}-T$ relation. No evidence for a break or slope change is found in either relation. Potential biases associated with sample selection, evolution, and the assumption of hydrostatic equilibrium that may affect the scaling are examined. No significant impacts attributable to these biases are found. Non-cool-core clusters and early spirals produce higher scatter in the $M-T$ relation than cool-core clusters and elliptical galaxies.

MAXI J1820$+$070 is a transient black hole binary (BHB) discovered on 2018 March 11. The unprecedented rich statistics brought by the NICER X-ray telescope allows detailed timing analysis up to $\sim$1~kHz uncompromised by the photon shot noise. To estimate the time lags, the Fourier analysis was applied, which led to two different conclusions for the system configuration; one supporting the lamp-post configuration with a stable accretion disk extending close to the innermost stable circular orbit and the other supporting the truncated accretion disk contracting with time. Using the same data set, we present the results based on the cross-correlation function (CCF). The CCF is calculated between two different X-ray bands and one side is subtracted from the other side, which we call the differential CCF (dCCF). The soft and hard lags respectively of $\sim$0.03 and 3~s are clearly identified without being diluted by the spectral mixture, demonstrating the effectiveness of the dCCF analysis. The evolution of these lags is tracked, along with spectral changes for the first 120~days since the discovery. Both the dCCF and spectral fitting results are interpreted that the soft lag is a reverberation lag between the Comptonized emission and the soft excess emission and that the hard lag is between the disk black body emission and the Comptonized emission. The evolution of these lags is in line with the picture of the truncated disk contracting with time.

Understanding electron acceleration associated with magnetic energy release at sub-second scales presents a major challenges in solar physics. Solar radio spikes observed as sub-second, narrow bandwidth bursts with $\Delta{f}/f\sim10^{-3}-10^{-2}$ are indicative of sub-second evolution of the electron distribution. We present a statistical analysis of frequency, and time-resolved imaging of individual spikes and Type IIIb striae associated with a coronal mass ejection (CME). LOFAR imaging reveals that co-temporal ($<2$ s) spike and striae intensity contours almost completely overlap. On average, both burst types have similar source size with fast expansion at millisecond scales. The radio source centroid velocities are often superluminal, and independent of frequency over 30-45 MHz. The CME perturbs the field geometry, leading to increased spike emission likely due to frequent magnetic reconnection. As the field restores towards the prior configuration, the observed sky-plane emission locations drift to increased heights over tens of minutes. Combined with previous observations above 1 GHz, average decay time and source size estimates follow $\sim1/f$ dependency over three decades in frequency, similar to radio-wave scattering predictions. Both time and spatial characteristics of the bursts between 30-70 MHz are consistent with radio-wave scattering with strong anisotropy of the density fluctuation spectrum. Consequently, the site of radio-wave emission does not correspond to the observed burst locations and implies acceleration and emission near the CME flank. The bandwidths suggest intrinsic emission source sizes $<1$ arcsec at 30 MHz, and magnetic field strengths a factor of two larger than average in events that produce decameter spikes.

We collect a large sample with a reliable redshift detected by the Fermi satellite after 10 years of data (4FGL-DR2), including blazars, $\gamma$-ray Narrow-line Seyfert 1 galaxies ($\gamma$NLS1s), and radio galaxies. The spectral energy distributions (SEDs) of these Fermi sources are fitted by using a second-degree polynomial, and some important parameters including spectral curvature, synchrotron peak frequency, and peak luminosity are obtained. Based on those parameters, we discuss the Fermi blazar sequence and the particle acceleration mechanism. Our main results are as follows:(i) By studying the relationship between the synchrotron peak frequency and the synchrotron peak frequency luminosity, jet kinetic power, and $\gamma$-ray luminosity for jetted AGNs, we find an ``L'' shape in the Fermi blazar sequence. (ii) There is a significant anti-correlation between Compton dominance, black hole spin, and the synchrotron peak frequency for jetted AGNs, respectively. These results support that the $\gamma$NLS1s and radio galaxies belong to the Fermi blazar sequence. (iii) On the basis of previous work, statistical or stochastic acceleration mechanisms can be used to explain the relationship between synchrotron peak frequency and synchrotron curvature. For different subclasses, the correlation slopes are different, which implies that the Fermi sources of different subclasses have different acceleration mechanisms. (iv) The FSRQs and $\gamma$NLS1s have a higher median spin of a black hole than BL Lacs and radio galaxies.

The open-source Python package Gammapy, developed for the high-level analysis of gamma-ray data, requires gamma-like event lists combined with the corresponding instrument response functions. For a morphological analysis, these data have to include a background acceptance model. Here we report an approach to generate such a model for the MAGIC telescope data, accounting for the azimuth and zenith dependencies of the MAGIC background acceptance. We validate this method using observations of the Crab Nebula with different offsets from the pointing position.

Galaxy mergers are hugely important in our current dark matter cosmology. These powerful events cause the disruption of the merging galaxies, pushing the gas, stars and dust of the galaxies resulting in changes to morphologies. This disruption can also cause more extreme events inside the galaxies: periods of extreme star formation rates and the rapid increase in active galactic nuclei activity. Hence, to better understand what goes on in these rare events, we need to be able to identify statistically large samples. In the last few years, the growth of artificial intelligence techniques has seen application to identifying galaxy mergers. These techniques have shown to be highly accurate and their application has grown beyond academic studies of ``can we?'' to deeper scientific use. However, these classifications are not without their problems. In this proceedings, we will explore how galaxy merger classification can be improved by adding pre-extracted galaxy morphologies alongside the traditional imaging data. This demonstrates that current neural networks are not extracting all the information from the images they are given. It will also explore how the resulting samples of rare objects could be highly contaminated. This has a knock on impact on the upcoming large scale surveys like Euclid and Rubin-LSST.

We present the re-detection of a compact source in the face-on protoplanetary disc surrounding HD 169142, using VLT/SPHERE data in YJH bands. The source is found at a separation of 0.''319 ($\sim$37 au) from the star. Three lines of evidence argue in favour of the signal tracing a protoplanet: (i) it is found in the annular gap separating the two bright rings of the disc, as predicted by theory; (ii) it is moving at the expected Keplerian velocity for an object at $\sim$37 au in the 2015, 2017 and 2019 datasets; (iii) we also detect a spiral-shaped signal whose morphology is consistent with the expected outer spiral wake triggered by a planet in the gap, based on dedicated hydrodynamical simulations of the system. The YJH colours we extracted for the object are consistent with tracing scattered starlight, suggesting that the protoplanet is enshrouded in a significant amount of dust, as expected for a circumplanetary disc or envelope surrounding a gap-clearing Jovian-mass protoplanet.

We present 170 optical spectra of 35 low-redshift stripped-envelope core-collapse supernovae observed by the Carnegie Supernova Project-I between 2004 and 2009. The data extend from as early as -19 days (d) prior to the epoch of B-band maximum to +322 d, with the vast majority obtained during the so-called photospheric phase covering the weeks around peak luminosity. In addition to histogram plots characterizing the red-shift distribution, number of spectra per object, and the phase distribution of the sample, spectroscopic classification is also provided following standard criteria. The CSP-I spectra are electronically available and a detailed analysis of the data set is presented in a companion paper being the fifth and final paper of the series

An analysis leveraging 170 optical spectra of 35 stripped-envelope (SE) core-collapse supernovae observed by the Carnegie Supernova Project-I and published in a companion paper is presented. Mean template spectra are constructed for the SNe IIb, Ib and Ic sub-types and parent ions associated with designated spectral features are identified with the aid of the spectral synthesis code SYNAPPS. Our modeled mean spectra suggest the ~6150~\AA\ feature in SNe~IIb may have an underlying contribution due to silicon, while the same feature in some SNe Ib may have an underlying contribution due to hydrogen. Standard spectral line diagnostics consisting of pseudo-equivalent widths (pEW) and blue-shifted Doppler velocity are measured for each of the spectral features. Correlation matrices and rolling mean values of both spectral diagnostics are constructed. A Principle Component Analysis (PCA) is applied to various wavelength ranges of the entire data set and suggests clear separation among the different SE SN sub-types, which follows from trends previously identified in the literature. In addition, our finds reveal the presence of two SNe IIb sub-types, a handful of SNe Ib displaying signatures of weak, high-velocity hydrogen, and a single SN~Ic with evidence of weak helium features. Our PCA results can be leveraged to obtain robust sub-typing of SE SN based on a single spectrum taken during the so-called photospheric phase, separating SNe IIb from SNe Ib with ~80 percent completion.

We aim to study the alignment and rotational disruption of dust grains at the centre of our Galaxy using polarized thermal dust emission observed by SOFIA/HAWC+ and JCMT/SCUPOL at 53, 216, and 850 $\mu$m. We analyzed the relationship between the observed polarization degree with total intensity, dust temperature, column density, and polarization angle dispersion function. Polarization degree from this region follows the predictions of the RAdiative Torque (RAT) alignment theory, except at high temperatures and long wavelengths where we found evidence for the rotational disruption of grains as predicted by the RAdiative Torque Disruption mechanism. The alignment and disruption sizes for the grains were found to be around 0.1 $\mu$m and 1 $\mu$m respectively. The maximum polarization degree observed was around $p\sim13$% at 216 $\mu$m and comes from a region of high temperature, low column density, and ordered magnetic field. Magnetically Enhanced RAT alignment (MRAT) was found to play an important role in the grain alignment due to the presence of a strong magnetic field. MRAT can lead to perfect alignment of the grains if they are super-paramagnetic in nature ($N_{\rm cl}\geqslant20$). We estimated the mass fraction of aligned grains using a parametric model for the fraction of the grains at high-$J$ attractors and found it to correlate weakly with the observed polarization degree. We observe a change in the polarization ratio, from $p_{216\mu m}/p_{850\mu m}<1$ to $p_{216\mu m}/p_{850\mu m}>1$ at $T_{\rm d}\gtrsim35$ K, which suggests a change in the grain model from a composite to a separate population of carbon and silicate grains as implied by previous numerical modeling.

In the last decades an incredible amount of evidence for the existence of dark matter (DM) has been accumulating. At the same time, many efforts have been undertaken to try to identify what dark matter is made of. Indirect searches look at places in the Universe where dark matter is known to be abundant and seek for possible annihilation or decay signatures. Indirect searches with the Fermi Gamma ray Space Telescope and Imaging Atmospheric Cherenkov Telescopes (IACTs) are playing a crucial role in constraining the nature of the DM particle through the study of their annihilation into gamma rays from different astrophysical structures. In this talk I will review the status of the search with IACTs and I will describe the sensitivity projections for dark matter searches on the various targets taking into account the latest instrument response functions expected for the Cherenkov Telescope Array (CTA) together with estimations for the systematic uncertainties from diffuse astrophysical and cosmic-ray

We present the hot molecular and warm ionised gas kinematics for 33 nearby ($0.001\lesssim z\lesssim0.056$) X-ray selected active galaxies using the H$_2 2.1218 \mu$m and Br$\gamma$ emission lines observed in the K-band with the Gemini Near-Infrared Field Spectrograph (NIFS). The observations cover the inner 0.04$-$2 kpc of each AGN at spatial resolutions of 4$-$250 pc with a velocity resolution of $\sigma_{\rm inst}\approx$20 ${\rm km s^{-1}}$. We find that 31 objects (94 per cent) present a kinematically disturbed region (KDR) seen in ionised gas, while such regions are observed in hot molecular gas for 25 galaxies (76 per cent). We interpret the KDR as being due to outflows with masses of 10$^2-$10$^7$ M$_\odot$ and 10$^0-$10$^4$ M$_\odot$ for the ionised and hot molecular gas, respectively. The ranges of mass-outflow rates ($\dot{M}_{\rm out}$) and kinetic power ($\dot{E}_{\rm K}$) of the outflows are 10$^{-3}-$10$^{1}$ M$_\odot$yr$^{-1}$ and $\sim$10$^{37}$$-$10$^{43}$ erg s$^{-1}$ for the ionised gas outflows, and 10$^{-5}$$-$10$^{-2}$ M$_\odot$ yr$^{-1}$ and 10$^{35}$$-$10$^{39}$ erg s$^{-1}$ for the hot molecular gas outflows. The median coupling efficiency in our sample is $\dot{E}_{K}/L_{\rm bol}\approx1.8\times10^{-3}$ and the estimated momentum fluxes of the outflows suggest they are produced by radiation-pressure in low-density environment, with possible contribution from shocks.

Selk crater is an $\sim$ 80 km diameter impact crater on the Saturnian icy satellite, Titan. Melt pools associated with impact craters like Selk provide environments where liquid water and organics can mix and produce biomolecules like amino acids. It is partly for this reason that the Selk region has been selected as the area that NASA's Dragonfly mission will explore and address one of its primary goals: to search for biological signatures on Titan. Here we simulate Selk-sized impact craters on Titan to better understand the formation of Selk and its melt pool. We consider several structures for the icy target material by changing the thickness of the methane clathrate layer, which has a substantial effect on the target thermal structure and crater formation. Our numerical results show that a 4 km-diameter-impactor produces a Selk-sized crater when 5-15 km thick methane clathrate layers are considered. We confirm the production of melt pools in these cases and find that the melt volumes are similar regardless of methane clathrate layer thickness. The distribution of the melted material, however, is sensitive to the thickness of the methane clathrate layer. The melt pool appears as a torus-like shape with a few km depth in the case of 10-15 km thick methane clathrate layer, and as a shallower layer in the case of a 5 km thick clathrate layer. Melt pools of this thickness may take tens of thousands of years to freeze, allowing more time for complex organics to form.

This letter describes a new phenomenon of recurring plasma density enhancements having delta_n/n less than or equal to 10 percent that occur at a repetition rate of about 5 Hz. They were observed sporadically for about five hours between 14 and 15 solar radii on Parker Solar Probe orbit 12 and they are seen in the same radial range on both the inbound and the outbound orbit 11 and other orbits. Their apparently steady-state existence suggests that their pressure gradient is balanced by the radial electric field, which was not measured. However, the electric field components measured in the plane perpendicular to the Sun-satellite line had sharp features like those hypothesized for the radial component, which offers support for the pressure balance requirement. These structures do not correlate with any magnetic field so they are not associated with any electromagnetic wave and they are probably not the result of a three-wave or similar process.

We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with finite-temperature microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios as low as 10, i.e., close to the detectability threshold of post-merger signals.

The heterogeneity of the small body population complicates the prediction of the small body properties before the spacecraft's arrival. In the context of autonomous small body exploration, it is crucial to develop algorithms that estimate the small body characteristics before orbit insertion and close proximity operations. This paper develops a vision-based estimation of the small-body rotational state (i.e., the center of rotation and rotation axis direction) during the approach phase. In this mission phase, the spacecraft observes the celestial body rotating and tracks features in images. As feature tracks are the projection of landmarks' circular movement, the possible rotation axes are computed. Then, the rotation axis solution is chosen among the possible candidates by exploiting feature motion and a heuristic approach. Finally, the center of rotation is estimated from the center of brightness. The algorithm is tested on more than 800 test cases with two different asteroids (i.e., Bennu and Itokawa), three different lighting conditions, and more than 100 different rotation axis orientations. Results show that the rotation axis can be determined with limited error in most cases implying that the proposed algorithm is a valuable method for autonomous small body characterization.

The aim of this work is to investigate the applications of the neutral atom imaging to the environments of the Earth, Mars and Mercury. This innovative technique permits the study of energetic plasma by means of analysing the result of the interaction of this plasma with a neutral thermal population or with a surface. The main advantage, when compared to the direct ion detection, is that it is possible to have an instantaneous survey of the whole magnetosphere of a planet. An example could help. Before the first ENA data, most of the knowledge about the Earth magnetospheric plasma came from in situ measurements of ions, electrons and electromagnetic fields. Those measurements, of course, could not represent any real instantaneous situation, but only an averaged picture of it, since the temporal and spatial variation cannot be easily be distinguished. Some short time scale phenomena, such as substorms, have been found difficult to comprehend without a global and continuous imaging. Even if some information about the plasma may be extracted from other sources, such as UV imaging [like aurorae, e.g. Horwitz, 1987], some populations (for example, the ring current) remained invisible. Furthermore, neutral atom imaging gives information not only about the energetic plasma, but also about the thermal neutral population (in the case of charge-exchange) or about the surface composition (in the case of sputtering). Conversely, it is necessary to set up some dedicated unfolding techniques to recover the 3D plasma distributions from the 2D ENA images.

Wide-field time-domain photometric sky surveys are now finding hundreds of eclipsing white dwarf plus M dwarf binaries, a population encompassing a wealth of information and potential insight into white dwarf and close binary astrophysics. Precise follow-up observations are essential in order to fully constrain these systems and capitalise on the power of this sample. We present the first results from our program of high-speed, multi-band photometric follow-up. We develop a method to measure temperatures, (model-dependent) masses, and radii for both components from the eclipse photometry alone and characterize 34 white dwarf binaries, finding general agreement with independent estimates using an alternative approach while achieving around a factor of two increase in parameter precision. In addition to these parameter estimates, we discover a number of interesting systems -- finding four with sub-stellar secondaries, doubling the number of eclipsing examples, and at least six where we find the white dwarf to be strongly magnetic, making these the first eclipsing examples of such systems and key to investigating the mechanism of magnetic field generation in white dwarfs. We also discover the first two pulsating white dwarfs in detached and eclipsing post-common-envelope binaries -- one with a low-mass, likely helium core, and one with a relatively high mass, towards the upper end of the known sample of ZZ Cetis. Our results demonstrate the power of eclipse photometry, not only as a method of characterising the population, but as a way of discovering important systems that would have otherwise been missed by spectroscopic follow-up.

Since the gravitational wave event GW170817 and gamma-ray burst GW170817A there have been numerous studies constraining the burst properties through analysis of the afterglow light curves. Most agree that the burst was viewed off-axis with a ratio of the observer angle to the jet angle ($\theta_{obs}/\theta_j$) between 4 - 6. We use a parameterized model and broadband synchrotron data up to $\sim 800$ days post-merger to constrain parameters of the burst. To reproduce the hydrodynamics of a gamma-ray burst outflow we use a two-parameter "boosted fireball" model. The structure of a boosted fireball is determined by the specific internal energy, $\eta_0$, and the bulk Lorentz factor, $\gamma_B(\sim 1/\theta_j)$ with shapes varying smoothly from a quasi-spherical outflow for low values of $\gamma_B$ to a highly collimated jet for high values. We run simulations with $\gamma_B$ in the range $1-20$ and $\eta_0$ in the range $2-15$. To calculate light curves we use a synchrotron radiation model characterized by $F_{peak}$, $\nu_m$, and $\nu_c$ and calculate millions of spectra at different times and $\theta_{obs}$ values using the \texttt{boxfit} radiation code. We can tabulate the spectral parameter values from our spectra and rapidly generate arbitrary light curves for comparison to data in MCMC analysis. We find that our model prefers a gamma-ray burst with jet energy $E_j\sim10^{50}$ ergs and with an observer angle of $\theta_{obs}=0.65^{+0.13}_{-0.14}$ radians and ratio to jet opening angle of ($\theta_{obs}/\theta_j$) = 5.4$^{+0.53}_{-0.38}$.

With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SNe Ia color evolution (needed for distance estimation) is uncertainty in telescope transmission, both within and between surveys. We introduce here the Collimated Beam Projector (CBP), which is meant to measure a telescope transmission with collimated light. The collimated beam more closely mimics a stellar wavefront as compared to flat-field based instruments, allowing for more precise handling of systematic errors such as those from ghosting and filter angle-of-incidence dependence. As a proof of concept, we present CBP measurements of the StarDICE prototype telescope, achieving a standard (1 sigma) uncertainty of 3 % on average over the full wavelength range measured with a single beam illumination.

We present the most extensive catalog of exposures of volatiles on the 67P/Churyumov-Gerasimenko nucleus generated from observations acquired with the OSIRIS cameras on board the Rosetta mission. We identified more than 600 volatile exposures on the comet. Bright spots are found isolated on the nucleus or grouped in clusters, usually at the bottom of cliffs, and most of them are small, typically a few square meters or smaller. Several of them are clearly correlated with the cometary activity. We note a number of peculiar exposures of volatiles with negative spectral slope values in the high-resolution post-perihelion images, which we interpret as the presence of large ice grains ($>$ 1000 $\mu$m) or local frosts condensation. We observe a clear difference both in the spectral slope and in the area distributions of the bright spots pre- and post-perihelion, with these last having lower average spectral slope values and a smaller size, with a median surface of 0.7 m$^2$, even if the size difference is mainly due to the higher resolution achieved post-perihelion. The minimum duration of the bright spots shows three clusters: an area-independent cluster dominated by short-lifetime frosts; an area-independent cluster with lifetime of 0.5--2 days, probably associated with the seasonal fallout of dehydrated chunks; and an area-dependent cluster with lifetime longer than 2 days consistent with water-driven erosion of the nucleus. Even if numerous bright spots are detected, the total surface of exposed water ice is less than 0.1% of the total 67P nucleus surface, confirming that the 67P surface is dominated by refractory dark terrains, while exposed ice occupies only a tiny fraction. Moreover, the abundance of volatile exposures is six times less in the small lobe than in the big lobe, adding additional evidence to the hypothesis that comet 67P is composed of two distinct bodies.

We are in the era of the Big Data. In Astronomy and Astrophysics, the massive amounts of data generated are, as of today, in the Peta-scale if not already in the Exa-scale. In the near future, we will see the data collected size and complexity grow further constantly, setting new challenges for data processing, reduction and analysis. This will pose new needs in terms of software and hardware solutions but also in terms of new models for resource management, access and sharing. In Astronomy and Astrophysics, in the environment of the International Virtual Observatory Alliance (IVOA), a big work has already been done with regards to data, gaining complete data FAIRness. In this paper, a model is proposed, based on the IVOA architecture, for software and hardware solutions for data analysis. The goal of this model is to build a cloud to access Astronomy and Astrophysics resources following the FAIR principles.

ESCAPE (European Science Cluster of Astronomy and Particle physics ESFRI research infrastructures) is a project to set up a cluster of ESFRI (European Strategy Forum on Research Infrastructures) facilities for astronomy, astroparticle and particle physics to face the challenges emerging through the modern multi-disciplinary data driven science. One of the main goal of ESCAPE is the building of ESAP (ESFRI Science Analysis Platform), a science platform for the analysis of open access data available through the EOSC (European Open Science Cloud) environment. ESAP will allow EOSC researchers to identify and stage existing data collections for analysis, share data, share and run scientific workflows. For many of the concerned ESFRIs and RIs, the data scales involved require significant computational resources (storage and compute) to support processing and analysis. The EOSC-ESFRI science platform therefore must implement appropriate interfaces to an underlying HPC (High Performance Computing) or HTC (High Throughput Computing) infrastructure to take advantage of it. This poster describes the analysis done to identify the main requirements for the implementation of the interfaces enabling the ESAP data access and computation resources integration in HPC and HTC computation infrastructures in terms of authentication and authorization policies, data management, workflow deployment and run.

We examine the potential for using the LIGO-Virgo-KAGRA network of gravitational-wave detectors to provide constraints on the physical properties of core-collapse supernovae through the observation of their gravitational radiation. We use waveforms generated by 14 of the latest 3D hydrodynamic core-collapse supernova simulations, which are added to noise samples based on the predicted sensitivities of the GW detectors during the O5 observing run. Then we use the BayesWave algorithm to model-independently reconstruct the gravitational-wave waveforms, which are used as input for various machine learning algorithms. Our results demonstrate how these algorithms perform in terms of i) indicating the presence of specific features of the progenitor or the explosion, ii) predicting the explosion mechanism, and iii) estimating the mass and angular velocity of the progenitor, as a function of the signal-to-noise ratio of the observed supernova signal. The conclusions of our study highlight the potential for GW observations to complement electromagnetic detections of supernovae by providing unique information about the exact explosion mechanism and the dynamics of the collapse.

The morphology of a galaxy has been shown to encode the evolutionary history and correlates strongly with physical properties such as stellar mass, star formation rates and past merger events. While the majority of galaxies in the local universe can be classified on the Hubble sequence, little is known about the different types of galaxies we observe at high redshift. The irregular morphology of these galaxies makes visual classifications difficult, and with the future of astronomy consisting of many "Big Data" surveys we need an efficient, and unbiased classification system in place. In this work we explore the use of unsupervised machine learning techniques to preform feature extraction from galaxy images to separate high redshift galaxies into different morphological types based on the machine learning clusters. We expand on previous work by addressing observation biases such as the orientation, apparent size of the galaxies and noise before extracting features, thus reducing the number of clusters and forcing the network to learn meaningful features. We then compare the extracted clusters' physical properties, to investigate the separation between the groups.

Bright basal reflectors in radargram from the Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) of the Martian south polar layered deposits (SPLD) have been interpreted to be evidence of subglacial lakes. However, this interpretation is difficult to reconcile with the low Martian geothermal heat flow and the frigid surface temperature at the south pole. We conduct a comprehensive thermophysical evolution modeling of the SPLD and show that subglacial lakes may only form under exceptional circumstances. Subglacial lakes may form if the SPLD contains more than 60 % dust volumetrically or extremely porous ice (>30 %), which is unlikely. A thick (>100 m) layer of dirty ice (>90% dust) at the base of the SPLD may also enable basal melting, resembling a sludge instead of a lake. Other scenarios enabling subglacial lakes in the SPLD are equally unlikely, such as recent magmatic intrusions at shallow depths.

Based on our dedicated Swift monitoring program, MOMO, OJ 287 is one of the best-monitored blazars in the X-ray--UV--optical regime. Here, we report results from our accompanying, dense, multi-frequency (1.4--44 GHz) radio monitoring of OJ 287 between 2015 and 2022 covering a broad range of activity states. Fermi gamma-ray observations are added. We characterize the radio flux and spectral variability in detail, including DCF and other variability analyses, and discuss its connection with the multiwavelength emission. Deep fades of radio and optical--UV fluxes are found to occur every 1--2 years. Further, it is shown that a precursor flare of thermal bremsstrahlung predicted by one of the binary supermassive black hole (SMBH) models of OJ 287 was absent. We then focus on the nature of the extraordinary, nonthermal 2016/2017 outburst that we initially discovered with Swift. We interpret it as the latest of the famous optical double-peaked outbursts of OJ 287, favoring binary scenarios that do not require a highly precessing secondary SMBH.

We report transmission spectroscopy of the bloated hot Jupiter WASP-74b in the wavelength range from 4000 to 6200 \r{A}. We observe two transit events with the Very Large Telescope FOcal Reducer and Spectrograph (VLT FORS2) and present a new method to measure the exoplanet transit depth as a function of wavelength. The new method removes the need for a reference star in correcting the spectroscopic light curves for the impact of atmospheric extinction. It also provides improved precision, compared to other techniques, reaching an average transit depth uncertainty of 211 ppm for a solar-type star of V=9.8 mag and over wavelength bins of 80 \r{A}. The VLT transmission spectrum is analysed both individually and in combination with published data from Hubble Space Telescope (HST) and Spitzer. The spectrum is found to exhibit a mostly featureless slope and equilibrium chemistry retrievals with PLATON favour hazes in the upper atmosphere of the exoplanet. Free chemistry retrievals with AURA further support the presence of hazes. While additional constraints are possible depending on the choice of atmospheric model, they are not robust and may be influenced by residual systematics in the data sets. Our results demonstrate the utility of new techniques in the analysis of optical, ground-based spectroscopic data and can be highly complementary to follow-up observations in the infrared with JWST.

Magnetic reconnection in the deep solar atmosphere can give rise to enhanced emission in the Balmer hydrogen lines, a phenomenon referred to as Ellerman bombs. To effectively trace magnetic reconnection below the canopy of chromospheric fibrils, we analyzed unique spectroscopic observations of Ellerman bombs in the H-alpha line. We analyzed a 10 min dataset of a young emerging active region observed with the prototype of the Microlensed Hyperspectral Imager (MiHI) at the Swedish 1-m Solar Telescope (SST). The MiHI instrument is an integral field spectrograph that is capable of achieving simultaneous ultra-high resolution in the spatial, temporal and spectral domains. With the combination of the SST adaptive optics system and image restoration techniques, MiHI can deliver diffraction limited observations if the atmospheric seeing conditions allow. The dataset samples the H-alpha line over 4.5 A with 10 mA/pix, with 0.065"/pix over a field of view of 8.6" x 7.7", and at a temporal cadence of 1.33s. This constitutes a hyperspectral data cube that measures 132 x 118 spatial pixels, 456 spectral pixels, and 455 time steps. There were multiple sites with Ellerman bomb activity associated with strong magnetic flux emergence. The Ellerman bomb activity is very dynamic, showing rapid variability and small-scale substructure. We found a number of plasmoid-like blobs with full-width-half-maximum sizes between 0.1" - 0.4" and moving with apparent velocities between 14 and 77 km/s. Some of these blobs have Ellerman bomb spectral profiles with a single peak at a Doppler offset between 47 and 57 km/s. Our observations support the idea that fast magnetic reconnection in Ellerman bombs is mediated by the formation of plasmoids. These MiHI observations demonstrate that a micro-lens based integral field spectrograph is capable of probing fundamental physical processes in the solar atmosphere.

We present a survey strategy to detect the neutral hydrogen (HI) power spectrum at $5<z<6$ using the SKA-Low radio telescope in presence of foregrounds and instrumental effects. We simulate observations of the inherently weak HI signal post-reionization with varying levels of noise and contamination with foreground amplitudes equivalent to residuals after sky model subtraction. We find that blind signal separation methods on imaged data are required in order to recover the HI signal at large cosmological scales. Comparing different methods of foreground cleaning, we find that Gaussian Process Regression (GPR) performs better than Principle Component Analysis (PCA), with the key difference being that GPR uses smooth kernels for the total data covariance. The integration time of one field needs to be larger than $\sim 250$ hours to provide large enough signal-to-noise ratio to accurately model the data covariance for foreground cleaning. Images within the primary beam field-of-view give measurements of the HI power spectrum at scales $k\sim 0.02\,{\rm Mpc^{-1}}-0.3\,{\rm Mpc^{-1} }$ with signal-to-noise ratio $\sim 2-5$ in $\Delta[{\rm log}( k/{\rm Mpc^{-1}})] = 0.25$ bins assuming an integration time of 600 hours. Systematic effects, which introduce small-scale fluctuations across frequency channels, need to be $\lesssim 1\%$ to enable unbiased measurements outside the foreground wedge. Our results provide an important validation towards using the SKA-Low array for measuring the HI power spectrum in the post-reionization Universe.

GRB 221009A is the brightest gamma-ray burst (GRB) ever detected and occurred at low Galactic latitude. Owing to this exceptional combination, its prompt X-ray emission could be detected for weeks in the form of expanding X-ray rings produced by scattering in Galactic dust clouds. We report on the analysis of 20 rings, generated by dust at distances ranging from 0.3 to 18.6 kpc, detected during two XMM-Newton observations performed about 2 and 5 days after the GRB. By fitting the spectra of the rings with different models for the dust composition and grain size distribution, we reconstructed the spectrum of the GRB prompt emission in the 0.7-4 keV energy range as an absorbed power law with photon index 1-1.4 and absorption in the host galaxy nHz=(4.1-5.3)E21 cm-2. Taking into account the systematic uncertainties on the column density of dust contained in the clouds producing the rings, the 0.5-5 keV fluence of GRB 221009A can be constrained between 1E-3 and 7E-3 erg cm-2.

In this contribution, we present the results of a study on the high abundance discrepancy factor (ADF $\sim$ 10) planetary nebula (PN) NGC 6153 with MUSE. We have constructed flux maps for dozens of emission lines, that allowed us to build spatially resolved maps of extinction, electron temperature ($T_{\rm e}$), electron density ($n_{\rm e}$), and ionic abundances. We have simultaneously constructed ADF maps for O$^+$ and O$^{2+}$ and found that they centrally peak in this PN, with a remarkable spatial coincidence with the low $T_{\rm e}$ found from recombination line diagnostics. This finding strongly supports the hypothesis that two distinct gas phases co-exist: one cold and metal-rich, and a second warm and with ``normal'' metal content. We show that to build $T_{\rm e}$([N II]) and ionic abundance maps of low-ionization species for these objects, recombination contribution to the auroral [N II] and [O II] lines must be properly evaluated and corrected.

Solar flares are among the most powerful and disruptive events in our solar system, however the physical mechanisms driving and transporting this energetic release are not fully understood. An important signature associated with flare energy release is highly variable emission on timescales of sub-seconds to minutes which often exhibit oscillatory behaviour, features collectively known as quasi-periodic pulsations (QPPs). To fully identify the driving mechanism of QPPs, exploit their potential as a diagnostic tool, and incorporate them into our understanding of solar and stellar flares, new observational capabilities and initiatives are required. There is a clear community need for flare-focused, rapid cadence, high resolution, multi-wavelength imaging of the Sun, with high enough sensitivity and dynamic range to observe small fluctuations in intensity in the presence of a large overall intensity. Furthermore, multidisciplinary funding and initiatives are required to narrow the gap between numerical models and observations. QPPs are direct signatures of the physics occurring in flare magnetic reconnection and energy release sites and hence are critical to include in a unified flare model. Despite significant modelling and theoretical work, no single mechanism or model can fully explain the presence of QPPs in flares. Moreover, it is also likely that QPPs fall into different categories that are produced by different mechanisms. At present we have insufficient information to observationally distinguish between mechanisms. The motivation to understand QPPs is strengthened by the geo-effectiveness of flares on the Earth's ionosphere, and by the fact that stellar flares exhibit similar QPP signatures. QPPs present a golden opportunity to better understand flare physics and exploit the solar-stellary analogy, benefiting both astrophysics, heliophysics, and the solar-terrestrial connection.

Climate modeling has shown that tidally influenced terrestrial exoplanets, particularly those orbiting M-dwarfs, have unique atmospheric dynamics and surface conditions that may enhance their likelihood to host viable habitats. However, sporadic libration and rotation induced by planetary interactions, such as that due to mean motion resonances (MMRs) in compact planetary systems may destabilize attendant exoplanets away from synchronized states (or 1:1 spin-orbit ratio). Here, we use a three-dimensional N-Rigid-Body integrator and an intermediately-complex general circulation model to simulate the evolving climates of TRAPPIST-1 e and f with different orbital and spin evolution pathways. Planet f perturbed by MMR effects with chaotic spin-variations are colder and dryer compared to their synchronized counterparts due to the zonal drift of the substellar point away from open ocean basins of their initial eyeball states. On the other hand, the differences between perturbed and synchronized planet e are minor due to higher instellation, warmer surfaces, and reduced climate hysteresis. This is the first study to incorporate the time-dependent outcomes of direct gravitational N-Rigid-Body simulations into 3D climate modeling of extrasolar planets and our results show that planets at the outer edge of the habitable zones in compact multiplanet systems are vulnerable to rapid global glaciations. In the absence of external mechanisms such as orbital forcing or tidal heating, these planets could be trapped in permanent snowball states.

The spatial distribution, Galactic model parameters and luminosity function of cataclysmic variables (CVs) are established using re-estimated trigonometric parallaxes of {\it Gaia} DR3. The data sample of 1,587 CVs in this study is claimed to be suitable for Galactic model parameter estimation as the distances are based on trigonometric parallaxes and the {\it Gaia} DR3 photometric completeness limits were taken into account when the sample was created. According to the analysis, the scale height of All CVs increases from 248$\pm$2 to 430$\pm$4 pc towards shorter periods near the lower limit of the period gap and suddenly drops to 300$\pm$2 pc for the shortest orbital period CVs. The exponential scale heights of All CVs and magnetic systems are found to be 375$\pm$2 and 281$\pm$3 pc, respectively, considerably larger than those suggested in previous observational studies. The local space density of All CVs and magnetic systems in the sample are $6.8^{+1.3}_{-1.1}\times$10$^{-6}$ and $2.1^{+0.5}_{-0.4}\times10^{-6}$ pc$^{-3}$, respectively. Our measurements strengthen the 1-2 order of magnitude discrepancy between CV space densities predicted by population synthesis models and observations. It is likely that this discrepancy is due to objects undetected by CV surveys, such as the systems with very low $\dot{M}$ and the ones in the period gap. The comparisons of the luminosity function of white dwarfs with the luminosity function of All CVs in this study show that 500 times the luminosity function of CVs fits very well to the luminosity function of white dwarfs. We conclude that the estimations and data sample in this study can be confidently used in further analysis of CVs.

Models where symmetries are predominantly broken (and masses are then generated) through radiative corrections typically produce strong first-order phase transitions with a period of supercooling, when the temperature dropped by several orders of magnitude. Here it is shown that a model-independent description of these phenomena and the consequent production of potentially observable gravitational waves is possible in terms of few parameters (which are computable once the model is specified) if enough supercooling occurred. It is explicitly found how large the supercooling should be in terms of those parameters, in order for the model-independent description to be valid. It is also explained how to systematically improve the accuracy of such description by computing higher-order corrections in an expansion in powers of a small quantity, which is a function of the above-mentioned parameters. Furthermore, the corresponding gravitational wave spectrum is compared with the existing experimental results from the latest observing run of LIGO and VIRGO and the expected sensitivities of future gravitational wave experiments to find regions of the parameter space that are either ruled out or can lead to a future detection.

Solar activities may disturb the geomagnetic field and impact the power grid via geomagnetically induced currents. We study active and reactive powers as well as the power factor of Iran's power grid transformers (230 kV and 400 kV) and their correlations with geomagnetic disturbances indices (SYM-H < -30 nT and sizable horizontal geomagnetic field fluctuation) from 19 March 2018 to 20 March 2020. Out of 128,627 cases with a transformer power factor of less than 0.7, we observe that 12,112 samples correlated with SYM-H. Our investigation shows that about 4 percent of two years, the SYM-H has values less than -30 nT. Analysis of high-performance transformers (a power factor greater than 0.7 at 95 percent of working time) shows at least a 55 percent correlation of power factor less than 0.7 and SYM-H less than -30 nT. We observe that the transformers' power factor of Rafsanjan-Kerman and Sarcheshmeh-Kerman (wye connection on 230 kV side) substations decreased to less than 0.7 and correlated with SYM-H. We show that the reactive power of the Sefidrood-Guilan and Shahid Beheshti-Guilan transformers (wye configurations) increased considerably on 9 January 2020 and positively correlated with SYM-H may produce large GICs at this part of the grid. We observe that the increase in reactive power at the Shahid Beheshti-Guilan substation correlated with the sizable changes in the horizontal field recorded by the Jaipur station. However, more details (temperature and current of transformers) need records to estimate the impact of geomagnetic induction current on the power grid.

We scrutinize the hypothesis that gauge singlet fermions -- sterile neutrinos -- interact with Standard Model particles through the transition magnetic moment portal. These interactions lead to the production of sterile neutrinos in supernovae followed by their decay into photons and active neutrinos which can be detected at $\gamma$-ray telescopes and neutrino detectors, respectively. We find that the non-observation of active neutrinos and photons from sterile-neutrino decay associated to SN1987A yields the strongest constraints to date on magnetic-moment-coupled sterile neutrinos if their masses are inside a $0.1-100$ MeV window. Assuming a near-future galactic supernova explosion, we estimate the sensitivity of several present and near-future experiments, including Fermi-LAT, e-ASTROGAM, DUNE, and Hyper-Kamiokande, to magnetic-moment-coupled sterile neutrinos. We also study the diffuse photon and neutrino fluxes produced in the decay of magnetic-moment coupled sterile neutrinos produced in all past supernova explosions and find that the absence of these decay daughters yields the strongest constraints to date for sterile neutrino masses inside a $1-100$ keV window.

Interaction of the standard model particles (electrons, quarks, gluons, photons) with dark matter field $\phi$ produces oscillating shift of atomic energy levels. To resolve these oscillations with frequency $m_{\phi} c^2/\hbar$, the measurements of the interaction strength have been done for very small masses of dark matter particles. However, for interaction proportional to $\phi^2$ one may measure energy shifts averaged over oscillations and probe dark matter particles with bigger mass. The problem is that only time dependence of the energy shifts produced by new interactions can be measured accurately. The constant part is hidden by the uncertainty of the theoretical values of the energies. However, amplitude of scalar or pseudoscalar (axion) field $\phi_0$ fluctuates on the time scale $\tau \sim 10^6 \hbar/m_{\phi} c^2$ and this causes fluctuations of the energy shift ($E \propto \phi_0^2$) of atomic, molecular and nuclear transition energies. Measurements of the shift fluctuation variance $\delta E^2 =\overline{(E - \overline{E})^2}=(\overline {E})^2$ will give us interaction strength constants for scalars and pseudoscalars. In the case of linear in $\phi$ interaction similar effects appear in the second order of perturbation theory which may be enhanced by small energy denominators.

Eccentric compact binary mergers are significant scientific targets for current and future gravitational wave observatories. To detect and analyze eccentric signals, there is an increasing effort to develop waveform models, numerical relativity simulations, and parameter estimation frameworks for eccentric binaries. Unfortunately, current models and simulations adopt different internal parameterisations of eccentricity in the absence of a unique natural definition of eccentricity in general relativity, which can result in incompatible eccentricity measurements. In this paper, we present a standard definition of eccentricity and mean anomaly based solely on waveform quantities. This definition is free of gauge ambiguities, has the correct Newtonian limit, and can be applied as a postprocessing step when comparing eccentricity measurements from different models. This standardization puts all models and simulations on the same footing and enables direct comparisons between eccentricity estimates from gravitational wave observations and astrophysical predictions. We demonstrate the applicability of our definition for waveforms of different origins, including post-Newtonian theory, effective one body, extreme mass ratio inspirals, and numerical relativity simulations. We focus on binaries without spin-precession in this work, but possible generalizations to spin-precessing binaries are discussed. We make our implementation publicly available through an easy-to-use Python package, gw_eccentricity.

Gravitational waves (GWs) from binary neutron stars (NSs) have opened unique opportunities to constrain the nuclear equation of state by measuring tidal effects associated with the excitation of characteristic modes of the NSs. This includes gravitomagnetic modes associated with the Coriolis effect, whose frequencies are proportional to the NS's spin frequency, and for which the spin orientation determines the subclass of modes that are predominantly excited. We advance the GW models for these effects that are needed for data analysis by first developing a description for the adiabatic signatures from gravitomagnetic modes in slowly rotating NSs. We show that they can be encapsulated in an effective Love number which differs before and after a mode resonance. Combining this with a known generic model for abrupt changes in the GWs at the mode resonance and a point-mass baseline leads to an efficient description which we use to perform case studies of the impacts of gravitomagnetic effects for measurements with Cosmic Explorer, an envisioned next-generation GW detector. We quantify the extent to which neglecting (including) the effect of gravitomagnetic modes induces biases (significantly reduces statistical errors) in the measured tidal deformability parameters, which depend on the equation of state. Our results substantiate the importance of dynamical gravitomagnetic tidal effects for measurements with third generation detectors.

We study the equation of state (EOS) dependence of polar $f$- and $p$- oscillation modes frequencies of cold and hot compact stars. Correlations between oscillation frequencies of cold purely nucleonic neutron stars (NSs), their EOS and global parameters and properties of nuclear matter (NM) are investigated by considering a set of covariant density functional (CDF) models generated within a Bayesian approach, where a minimal number of constraints on the saturation properties of NM, the behavior of pure neutron matter (PNM) as a function of density and the lower bound of maximal NS mass were imposed. Then the role of finite temperature effects is addressed by assuming various idealized profiles of temperature or entropy per baryon and charge fraction for models accounting, in addition to nucleons, also for hyperons, $\Delta$-resonances, anti-kaon condensates or a hadron to quark phase transition. We find that finite temperature effects reduce the oscillation frequencies of nucleonic stars while the opposite effect is obtained for stars with exotic particle degrees of freedom. When the $\Gamma$-law is employed to build finite temperature EOS errors in estimating oscillation modes frequencies are of the order of 10\% to 30\%, depending on the mass. Throughout this work the Cowling approximation is used.

We extract the Hubble law by the frequency shift considerations of test particles revolving the Kerr black hole in asymptotically de Sitter spacetime. To this end, we take into account massive geodesic particles circularly orbiting the Kerr-de Sitter black holes that emit redshifted photons towards a distant observer whose going away from the emitter--black hole system. By considering this configuration, we obtain an expression for redshift in terms of the spacetime parameters, such as mass, angular momentum, and the cosmological constant. Then, we find the frequency shift of photons versus the Hubble constant with the help of some physically motivated approximations. Finally, some exact formulas for the Schwarzschild black hole mass and the Hubble constant in terms of the observational redshift of massive bodies circularly orbiting this black hole are extracted. Our results suggest a new independent general relativistic approach to obtaining the late--time Hubble constant in terms of observable quantities.