Papers for Friday, Apr 28 2023

We present Gemini-N/'Alopeke speckle observations of KIC 9832227, a system originally predicted to become a red nova. The diffraction limited observations do not find an outer companion between 11 and 678 AU that could be responsible for the period changes of the system.

In order to test active galactic nucleus (AGN) unification and evolutionary models, we measured the AGN clustering properties as a function of AGN obscuration defined in terms of hydrogen column density, $N_{\rm H}$. In addition to measuring the clustering of unobscured ($N_{\rm H} < 10^{22}\,{\rm cm}^{-2}$) and moderately obscured ($10^{22} \leq N_{\rm H} < 10^{23.5}$) AGNs, we also targeted highly obscured sources ($N_{\rm H}\geq 10^{23.5}$) up to redshifts of $z=3$. We have compiled one of the largest samples of X-ray-selected AGNs from a total of eight deep XMM/Chandra surveys. We measured the clustering as a function of both AGN obscuration and redshift using the projected two-point correlation function, $w_{\rm p}(r_{\rm p})$. We modeled the large-scale clustering signal, measured the AGN bias, $b(z, N_{\rm H})$, and interpreted it in terms of the typical AGN host dark matter halo, $M_{\rm halo}(z, N_{\rm H}$). We find no significant dependence of AGN clustering on obscuration, suggesting similar typical masses of the hosting halos as a function of $N_{\rm H}$. This result matches expectations of AGN unification models, in which AGN obscuration depends mainly on the viewing angle of the obscuring torus. We measured, for the first time, the clustering of highly obscured AGNs and find that these objects reside in halos with typical mass $\log M_{\rm halo} = 12.98_{-0.22}^{+0.17} [h^{-1} M_\odot]$ ($12.28_{-0.19}^{+0.13}$) at low $z \sim 0.7$ (high $z \sim 1.8$) redshifts. We find that irrespective of obscuration, an increase in AGN bias with redshift is slower than the expectation for a constant halo mass and instead follows the growth rate of halos, known as the passive evolution track. This implies that for those AGNs the clustering is mainly driven by the mass growth rate of the hosting halos and galaxies across cosmic time.

Exascale super-computers now becoming available rely on hybrid energy-efficient architectures that involve an accelerator such as Graphics Processing Units (GPU). Leveraging the computational power of these machines often means a significant rewrite of the numerical tools each time a new architecture becomes available. To address these issues, we present Idefix, a new code for astrophysical flows that relies on the Kokkos meta-programming library to guarantee performance portability on a wide variety of architectures while keeping the code as simple as possible for the user. Idefix is based on a Godunov finite-volume method that solves the non-relativistic HD and MHD equations on various grid geometries. Idefix includes a wide choice of solvers and several additional modules (constrained transport, orbital advection, non-ideal MHD) allowing users to address complex astrophysical problems. Idefix has been successfully tested on Intel and AMD CPUs (up to 131 072 CPU cores on Irene-Rome at TGCC) as well as NVidia and AMD GPUs (up to 1024 GPUs on Adastra at CINES). Idefix achieves more than 1e8 cell/s in MHD on a single NVidia V100 GPU and 3e11 cell/s on 256 Adastra nodes (1024 GPUs) with 95% parallelization efficiency (compared to a single node). For the same problem, Idefix is up to 6 times more energy efficient on GPUs compared to Intel Cascade Lake CPUs. Idefix is now a mature exascale-ready open-source code that can be used on a large variety of astrophysical and fluid dynamics applications.

Ground-based optical astronomical observations supported by or in the vicinity of laser guide-star systems can be contaminated by Raman-scattered laser photons. Anticipating, alleviating, and correcting for the impact of this self-inflicted contamination requires a detailed knowledge of the pure-rotational and rotational-vibrational spectrum of the molecules in our atmosphere. We present a 15.3hr-deep combined spectrum of the 4LGSF's 589nm $\approx$ 509THz sodium laser beams of Paranal observatory, acquired with the ESPRESSO spectrograph at a resolution $\lambda/\Delta\lambda\cong140'000\approx0.12$ cm$^{-1}$ and an altitude of 23 km above mean sea level. We identify 865 Raman lines over the spectral range of [3770; 7900]{\AA}$\approx$[+9540; -4315] cm$^{-1}$, with relative intensities spanning ~5 orders of magnitudes. These lines are associated to the most abundant molecules of dry air, including their isotopes: 14N14N, 14N15N, 16O16O, 16O17O, 16O18O, and 12C16O16O. The signal-to-noise of these observations implies that professional observatories can treat the resulting catalogue of Raman lines as exhaustive (for the detected molecules, over the observed Raman shift range) for the purpose of predicting/correcting/exploiting Raman lines in astronomical data. Our observations also reveal that the four laser units of the 4LGSF do not all lase at the same central wavelength. [...] The [measured] offsets [...] are larger than the observed 4LGSF spectral stability of $\pm$3 MHz over hours. They remain well within the operational requirements for creating artificial laser guide-stars, but hinder the assessment of the radial velocity accuracy of ESPRESSO at the required level of 10 m/s. Altogether, our observations demonstrate how Raman lines can be exploited by professional observatories as highly-accurate, on-sky wavelength references.

We present follow-up spectroscopy of 21 cataclysmic variables (CVs) with evolved secondaries and ongoing or recently-terminated mass transfer. Evolutionary models predict that the secondaries should have anomalous surface abundances owing to nuclear burning in their cores during their main-sequence evolution and subsequent envelope stripping by their companion white dwarfs. To test these models, we measure sodium (Na) abundances of the donors from the Fraunhofer "D" doublet. Accounting for interstellar absorption, we find that {\it all} objects in our sample have enhanced Na abundances. We measure 0.3 $\lesssim$ [Na/H] $\lesssim$ 1.5 dex across the sample, with a median [Na/H] = 0.956 dex, i.e., about an order of magnitude enhancement over solar values. To interpret these values, we run MESA binary evolution models of CVs in which mass transfer begins just as the donor leaves the main sequence. These generically predict Na enhancement in donors with initial donor masses $\gtrsim 1\,M_{\odot}$, consistent with our observations. In the models, Na enrichment occurs in the donors' cores via the NeNa cycle near the end of their main-sequence evolution. Na-enhanced material is exposed when the binaries reach orbital periods of a few hours. Donors with higher initial masses are predicted to have higher Na abundances at fixed orbital period owing to their higher core temperatures during main-sequence evolution. The observed [Na/H] values are on average $\approx$0.3 dex higher than predicted by the models. Surface abundances of evolved CV donors provide a unique opportunity to study nuclear burning products in the cores of intermediate-mass stars.

In the past few years, we built a Hubble diagram of quasars up to redshift z$\sim$7, based on the nonlinear relation between quasars' x-ray and UV luminosities. Such a Hubble diagram shows a >4$\sigma$ deviation from the standard flat $\Lambda$CDM model at z>1.5. Given the important consequences of this result, it is fundamental to rule out any systematic effect in the selection of the sample and/or in the flux measurements, and to investigate possible redshift dependences of the relation, that would invalidate the use of quasars as standard candles. Here we review all the observational results supporting our method: the match of the Hubble diagram of quasars with that of supernovae in the common redshift range, the constant slope of the relation at all redshifts, the redshift non-evolution of the spectral properties of our sources both in the x-rays and in the UV. An independent test of our results requires the observation of other standard candles at high redshift. In particular, we expect that future observations of supernovas at z>2 will confirm the deviation from the concordance model found with the Hubble diagram of quasars.

Irregular satellites are the minor bodies found orbiting all four Solar System giant planets, with large semi-major axes, eccentricities, and inclinations. Previous studies have determined that the Solar System's irregular satellites are extremely collisionally evolved populations today, having lost $\sim$99 per cent of their initial mass over the course of hundreds of Myr. Such an evolution implies that the irregular satellites must have produced a population of dusty collisional debris in the past, which is potentially observable due to the resulting reprocessing of stellar light. In this paper we examine the signatures of the debris discs produced by extrasolar analogues of this process. Radiation pressure, quantified by the parameter $\beta$, is the driving force behind the liberation of dust grains from the planetary Hill sphere, and results in the formation of circumstellar dust rings, even in the absence of an underlying belt of asteroids in the system. Our simulated discs reproduce many of the same features seen in some classes of observed debris discs, such as thin ring morphology, a large blowout size, and azimuthal symmetry. We compare our simulated discs' radial profiles to those of the narrow dust rings observed around Fomalhaut and HR 4796A, and show that they can broadly reproduce the observed radial distribution of dust.

Recent observations by JWST have uncovered galaxies in the very early universe via the JADES and CEERS surveys. These galaxies have been measured to have very high stellar masses with substantial star formation rates. There are concerns that these observations are in tension with the $\Lambda$CDM model of the universe, as the stellar masses of the galaxies are relatively high for their respective redshifts. Recent studies have compared the JWST observations with large-scale cosmological simulations. While they were successful in reproducing the galaxies seen in JADES and CEERS, the mass and spatial resolution of these simulations were insufficient to fully capture the early assembly history of the simulated galaxies. In this study, we use results from the Renaissance simulations, which are a suite of high resolution simulations designed to model galaxy formation in the early universe. We find that the most massive galaxies in Renaissance have stellar masses and star formation rates that are entirely consistent with the observations from the JADES and CEERS surveys. The exquisite resolution afforded by Renaissance allows us to model the build-up of early galaxies from stellar masses as low as 10$^4$ M$_\odot$ up to a maximum stellar mass of a few times 10$^{7}$ M$_\odot$. Within this galaxy formation paradigm, we find excellent agreement with JADES and CEERS. We find no tension between the $\Lambda$CDM model and current JWST measurements. As JWST continues to explore the high redshift universe, high resolution simulations, such as Renaissance, will continue to be crucial in understanding the formation history of early embryonic galaxies.

Radio used as a star formation rate (SFR) tracer presents enormous advantages by being unaffected by dust and radio sources being pinpointed at the sub-arc-second level. The interpretation of the low frequency 1.4 GHz luminosity is hampered by the difficulty in modeling the cosmic ray paths in the interstellar medium, and their interactions with the magnetic field. In this work, we compare the SFR derived from radio observations, and the ones derived from spectral energy distribution (SED) modeling. We aim at better understand the behavior of the SFR radio tracer, with a specific emphasis on the link with star-formation histories. We used the SED modeling code Code Investigating GALaxy Emission, CIGALE, with a non-parametric star formation history model (SFH) and fit the data over the wavelength range from the ultraviolet (UV) up to the mid-infrared (mid-IR). We interpret the difference between radio and SED-based SFR tracers in the light of recent gradients in the derived SFH. To validate the robustness of the results, we checked for any remaining active galaxy nuclei (AGN) contribution and tested the impact of our SFH modeling approach. Approximately 27% our galaxies present a radio SFR (SFR$_{\rm radio}$) at least ten times larger than the instantaneous SFR from SED-fitting (SFR$_{\rm SED}$). This trend affects primarily the galaxies that show a declining SFH activity over the last 300 Myr. Both SFR indicators converge toward a consistent value, when the SFHs are averaged over a period larger than 150 Myr to derive SFR$_{\rm SED}$. Although the radio at low frequency 1.4 GHz is a good tracer of the star formation activity of galaxies with constant or increasing SFH, our results indicate that this is not the case for galaxies that are quenching. Our analysis suggests that the star formation time sensitivity of the radio low frequency could be longer than 150 Myr.

We characterise the planet-occulted line distortions (POLDs) in absorption spectra of transiting planets, that arise from proxies used for the occulted stellar lines and investigate the impact of stellar rotation, centre-to-limb variations, and broadband limb-darkening. We used the EVaporating Exoplanets (EVE) code to generate realistic stellar spectra during the transit of exoplanets, accounting for the 3D geometry of the system's architecture and atmospheric transit, as well as for spectral variations over the stellar disc. The absorption spectra were calculated using approaches drawn from the literature and compared to the expected signal. The POLDs from stellar rotation are dominant for moderate to fast rotating stars, reaching amplitudes comparable to atmospheric signals, but they can be mitigated by shifting the stellar line proxies to the radial velocity of the planet-occulted region. Centre-to-limb variations become dominant for slow rotators and are more easily mitigated at the stellar limb. We re-interpret the ESPRESSO data of two iconic systems and confirm that the sodium signature from HD209458b mainly arises from POLDs. However, we unveil a possible contribution from the planetary atmosphere that warrants further observations. For MASCARA-1b, we did not find evidence for atmospheric sodium absorption and we can fully explain the observed signature by a POLD for super-solar stellar sodium abundance. We studied POLDs dependency on star and planet properties, and on the proxy used for planet-occulted lines. Distinguishing planetary absorption signatures from POLDs is challenging without access to accurate estimates of the local stellar spectrum and system orbital parameters. We propose a way to mitigate POLDs and improve atmospheric characterisation, by using simultaneous forward modelling of both the star and the planet to simulate the global observed signatures.

We present the wavelength-resolved reverberation mapping (RM) of combined Mg II and UV Fe II broad-line emissions for two intermediate redshifts (z$\sim$1), luminous quasars - HE 0413-4031 and HE 0435-4312, monitored by the SALT and 1-m class telescopes between 2012-2022. Through this technique, we aim to disentangle the Mg II and Fe II emission regions and to build a radius-luminosity relation for UV Fe II emission, which has so far remained unconstrained. Several methodologies have been applied to constrain the time delays for total Mg II and Fe II emissions. In addition, this technique is performed to quantify the inflow or outflow of broad-line region gas around the supermassive black hole and to disentangle the emission/emitting regions from lines produced in proximity to each other. The mean total FeII time delay is nearly equal to the mean total MgII time delay for HE 0435-4312 suggesting a co-spatiality of their emissions. However, in HE 0413-4031, the mean FeII time delay is found to be longer than the mean MgII time delay, suggesting that FeII is produced at longer distances from the black hole. The UV FeII R-L relation is updated with these two quasars and compared with the optical FeII relation, which suggests that the optical FeII region is located further than the UV FeII by a factor of 1.7-1.9 i.e. $R_{\rm FeII-opt}\sim(1.7-1.9)R_{\rm FeII-UV}$. We detected a weak pattern in the time delay vs. wavelength relation, suggesting that the MgII broad-line originates a bit closer to the SMBH than the UV FeII, however, the difference is not very significant. Comparison of MgII, UV, and optical FeII R-L relations suggests that the difference may be larger for lower-luminosity sources, possibly with the MgII emission originating further from the SMBH. In the future, more RM data will be acquired to put better constraints on these trends, in particular the UV FeII R-L relation.

In dwarf nov{\ae} and low-mass X-ray binaries, the tidal potential excites spiral waves in the accretion disc. Spiral wave driven accretion may be important in quiescent discs, where the angular momentum transport mechanism has yet to be identified. Previous studies were limited to unrealistically high temperatures for numerical studies or to specific regimes for analytical studies. We perform the first numerical simulation of spiral wave driven accretion in the cold temperature regime appropriate to quiescent discs, which have Mach numbers > 100. We use the new GPU-accelerated finite volume code Idefix to produce global hydrodynamics 2D simulations of the accretion discs of dwarf nov{\ae} systems with a fine-enough spatial resolution to capture the short scale-height of cold, quiescent discs with Mach numbers ranging from 80 to 370. Running the simulations on timescales of tens of binary orbits shows transient angular momentum transport that decays as the disc relaxes from its initial conditions. We find the angular momentum parameter {\alpha} drops to values << 0.01 , too weak to drive accretion in quiescence.

We present a detailed analysis of the point spread function (PSF) of JWST NIRCam imaging in eight filters: F070W, F115W, F150W, F200W, F277W, F356W, F444W, and F480M, using publicly available data. Spatial variations in the PSF FWHM generally decrease with wavelength: the maximum and RMS fractional variations are $\sim20\%$ and $5\%$ in F070W, reduced to $\sim3\%$ and $0.6\%$ in F444W. We compare three commonly-used methods (SWarp, photutils, and PSFEx) to construct model PSFs and conclude that PSFEx delivers the best performance. Using simulated images of broad-line AGNs, we evaluate the impact of PSF mismatches on the recoverability of host galaxy properties. Host fluxes are generally overestimated when adopting mismatched PSF models, with larger overestimation for more AGN-dominated systems. Broader PSFs tend to produce less concentrated hosts while narrower PSFs tend to produce more concentrated and compact hosts. Systematic uncertainties in host measurements from PSF and model mismatches are generally larger than the formal fitting uncertainties for high signal-to-noise ratio data. Image decomposition can also lead to an artificial offset between the AGN and host centroids, which is common (e.g., $>1\sigma$ [$3\sigma$] detection in $\sim 80%$ [$\sim 20-30\%$] of systems), and scales with the mean host surface brightness. Near the surface brightness limit, this artificial offset can reach as large as $\sim80\%$, $26\%$, and $7\%$ of $R_e$ in systems with $R_e=$0.12", 0.48", and 1.92", respectively. We demonstrate our PSF construction and image decomposition methods with an example broad-line quasar at $z=1.646$ in the CEERS field.

Mobile communication towers represent a relatively new but growing contributor to the total radio leakage associated with planet Earth. We investigate the overall power contribution of mobile communication towers to the Earth\'s radio leakage budget, as seen from a selection of different nearby stellar systems. We created a model of this leakage using publicly available data of mobile tower locations. The model grids the planet's surface into small, computationally manageable regions, assuming a simple integrated transmission pattern for the mobile antennas. In this model, these mobile tower regions rise and set as the Earth rotates. In this way, a dynamic power spectrum of the Earth was determined, summed over all cellular frequency bands. We calculated this dynamic power spectrum from three different viewing points, HD 95735, Barnard star, and Alpha Centauri A. Our preliminary results demonstrate that the peak power leaking into space from mobile towers is $\sim 4$GW. This is associated with LTE mobile tower technology emanating from the East Coast of China as viewed from HD 95735. We demonstrate that the mobile tower leakage is periodic, direction dependent, and could not currently be detected by a nearby civilization located within 10 light years of the Earth, using instrumentation with a sensitivity similar to the Green Bank Telescope. We plan to extend our model to include more powerful 5G mobile systems, radar installations, ground based uplinks (including the Deep Space Network), and various types of satellite services, including low Earth orbit constellations such as Starlink and OneWeb.

Ground-based characterization of spacecraft targets prior to mission operations is critical to properly plan and execute measurements. Understanding surface properties, like mineralogical composition and phase curves (expected brightness at different viewing geometries) informs data acquisition during the flybys. Binary near-Earth asteroids (NEA) (35107) 1991 VH and (175706) 1996 FG3 were selected as potential targets of the National Aeronautics and Space Administration's (NASA) dual spacecraft Janus mission. We observed 1991 VH using the 3-m NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawaii, on July 26, 2008. 1996 FG3 was observed with the IRTF for seven nights during the spring of 2022. Compositional analysis of 1991 VH revealed that this NEA is classified as an Sq-type in the Bus-DeMeo taxonomy classification, with a composition consistent with LL ordinary chondrites. Using thermal modeling, we computed the thermally corrected spectra for 1996 FG3 and the corresponding best fit albedo of about 2-3% for the best spectra averaged for each night. Our spectral analysis indicates that this NEA is a Ch-type. The best possible meteorite analogs for 1996 FG3, based on curve matching, are two carbonaceous chondrites, Y-86789 and Murchison. No rotational variation was detected in the spectra of 1996 FG3, which means there may not be any heterogeneities on the surface of the primary. However, a clear phase reddening effect was observed in our data, confirming findings from previous ground-based studies.

Extragalactic fast X-ray transients (FXTs) are short flashes of X-ray photons of unknown origin that last a few minutes to hours. We extend the search for extragalactic FXTs from Quirola et al. 2022 (Paper I; based on sources in the Chandra Source Catalog 2.0, CSC2) to further Chandra archival data between 2014-2022. We extract X-ray data using a method similar to that employed by CSC2 and apply identical search criteria as in Paper I. We report the detection of eight FXT candidates, with peak 0.3-10 keV fluxes between 1$\times$10$^{-13}$ to 1$\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ and $T_{90}$ values from 0.3 to 12.1 ks. This sample of FXTs has likely redshifts between 0.7 to 1.8. Three FXT candidates exhibit light curves with a plateau (${\approx}$1-3 ks duration) followed by a power-law decay and X-ray spectral softening, similar to what was observed for a few previously reported FXTs in Paper I. In light of the new, expanded source lists (eight FXTs with known redshifts from Paper I and this work), we update the event sky rates derived in Paper I, finding 36.9$_{-8.3}^{+9.7}$ deg$^{-2}$ yr$^{-1}$ for the extragalactic samples for a limiting flux of ${\gtrsim}$1${\times}$10$^{-13}$ erg cm$^{-2}$ s$^{-1}$, calculate the first FXT X-ray luminosity function, and compare the volumetric density rate between FXTs and other transient classes. Our latest Chandra-detected extragalactic FXT candidates boost the total Chandra sample by $\sim$50 %, and appear to have a similar diversity of possible progenitors.

Software is a critical aspect of large-scale science, providing essential capabilities for making scientific discoveries. Large-scale scientific projects are vast in scope, with lifespans measured in decades and costs exceeding hundreds of millions of dollars. Successfully designing software that can exist for that span of time, at that scale, is challenging for even the most capable software companies. Yet scientific endeavors face challenges with funding, staffing, and operate in complex, poorly understood software settings. In this paper we discuss the practice of early-phase software architecture in the Square Kilometre Array Observatory's Science Data Processor. The Science Data Processor is a critical software component in this next-generation radio astronomy instrument. We customized an existing set of processes for software architecture analysis and design to this project's unique circumstances. We report on the series of comprehensive software architecture plans that were the result. The plans were used to obtain construction approval in a critical design review with outside stakeholders. We conclude with implications for other long-lived software architectures in the scientific domain, including potential risks and mitigations.

Within the $\Lambda$CDM cosmology, dark matter haloes are comprised of both a smooth component and a population of smaller, gravitationally bound subhaloes. These components are often treated as a single halo when halo properties, such as density profiles, are extracted from simulations. Recent work has shown that density profiles change substantially when subhalo mass is excluded. In this paper, we expand on this result by analysing the change in three specific host halo properties -- concentration ($c_{\rm{NFW}}$), spin ($\lambda_{\rm Bullock}$), and shape ($c/a$), -- when calculated only from the smooth component of the halo. This analysis is performed on both Milky Way-mass haloes and cluster-mass haloes in high-resolution, zoom-in, $N$-body simulations. We find that when subhaloes are excluded the median value of (1) $c_{\rm{NFW}}$ is enhanced by $\approx 38 \pm 12\%$ and $\approx 88 \pm 7.7\%$ for Milky Way mass ($10^{12.1}\,\text{M}_\odot$) and cluster mass ($10^{14.8}\,\text{M}_\odot$) haloes respectively, (2) $\lambda_{\rm Bullock}$ is reduced for Milky Way mass by $\approx 16 \pm 6.8\%$ and cluster mass haloes by $\approx 32 \pm 8.9\%$. Additionally, with the removal of subhaloes, cluster mass haloes tend to become more spherical as the ratio of minor-to-major axis, $c/a$, increases by $\approx 12 \pm 4\%$, whereas Milky Way mass haloes remain approximately the same shape with $c/a$ changed by $\approx 1.2 \pm 5.6\%$. The fractional change of each of these properties depends primarily on the amount of mass that is removed from the halo system and, to a lesser extent, mass accretion history. Our findings demonstrate that the properties of the smooth components of dark matter haloes are biased relative to the total mass of the halo including subhaloes.

Ground-based, all-sky astronomical surveys are imposed with an inevitable day-night cadence that can introduce aliases in period-finding methods. We examined four different methods -- three from the literature and a new one that we developed -- that remove aliases to improve the accuracy of period-finding algorithms. We investigate the effectiveness of these methods in decreasing the fraction of aliased period solutions by applying them to the Zwicky Transient Facility (ZTF) and the LSST Solar System Products Data Base (SSPDB) asteroid datasets. We find that the VanderPlas method had the worst accuracy for each survey. The mask and our newly proposed window method yields the highest accuracy when averaged across both datasets. However, the Monte Carlo method had the highest accuracy for the ZTF dataset, while for SSPDB, it had lower accuracy than the baseline where none of these methods are applied. Where possible, detailed de-aliasing studies should be carried out for every survey with a unique cadence.

In this paper, we revisit the hydrogen recombination history from a novel perspective: the evolution of chemical potentials. We derive expressions for the chemical potentials, which depend on the thermal bath temperature and the ionization degree of the universe. Our main finding reveals a constraint between the chemical potentials of hydrogen and proton at $z\approx 1200$ when the free electron fraction is $X_e\approx 1/3$. Furthermore, we present important data on the chemical potentials during recombination, highlighting the differences between the predictions of the Peebles' and CosmoRec code solutions. Finally, we discuss a particular case related to the chemical potential of hydrogen.

We present an analytical model that describes the response of companion stars after being impacted by a supernova in a close binary system. This model can be used to constrain the pre-supernova binary properties using photometry of the companion star several years after the explosion in a relatively simple manner. The derived binary parameters are useful in constraining the evolutionary scenario for the progenitors and the physics of binary interactions. We apply our model to the observed photometry of some known stripped-envelope supernova companions (SN1993J, SN2001ig, SN2006jc, SN2011dh, SN2013ge). Combined with other observational constraints such as the pre-supernova progenitor photometry, we find that SN1993J and SN2011dh likely had relatively massive companions on wide orbits, while SN2006jc may have had a relatively low-mass companion on a tight orbit. This trend suggests that type IIb supernova progenitors evolved from stable mass transfer channels and type Ibc progenitors may have formed from common-envelope channels. The constraints on orbital separation helps us probe the highly uncertain common-envelope physics for massive stars, especially with multiple epochs of companion observations.

We present Atacama Large Millimeter/submillimeter Array Band 3 data toward five massive young stellar objects (MYSOs), and investigate relationships between unsaturated carbon-chain species and saturated complex organic molecules (COMs). An HC$_{5}$N ($J=35-34$) line has been detected from three MYSOs, where nitrogen(N)-bearing COMs (CH$_{2}$CHCN and CH$_{3}$CH$_{2}$CN) have been detected. The HC$_{5}$N spatial distributions show compact features and match with a methanol (CH$_{3}$OH) line with an upper-state energy around 300 K, which should trace hot cores. The hot regions are more extended around the MYSOs where N-bearing COMs and HC$_{5}$N have been detected compared to two MYSOs without these molecular lines, while there are no clear differences in the bolometric luminosity and temperature. We run chemical simulations of hot-core models with a warm-up stage, and compare with the observational results. The observed abundances of HC$_{5}$N and COMs show good agreements with the model at the hot-core stage with temperatures above 160 K. These results indicate that carbon-chain chemistry around the MYSOs cannot be reproduced by warm carbon-chain chemistry, and a new type of carbon-chain chemistry occurs in hot regions around MYSOs.

We model the effect of grain size distribution in a galaxy on the evolution of CO and H$_2$ abundances. The formation and dissociation of CO and H$_2$ in typical dense clouds are modelled in a manner consistent with the grain size distribution. The evolution of grain size distribution is calculated based on our previous model, which treats the galaxy as a one-zone object but includes various dust processing mechanisms in the interstellar medium (ISM). We find that typical dense clouds become fully molecular (H$_2$) when the dust surface area increases by shattering while an increase of dust abundance by dust growth in the ISM is necessary for a significant rise of the CO abundance. Accordingly, the metallicity dependence of the CO-to-H$_2$ conversion factor, $X_\mathrm{CO}$, is predominantly driven by dust growth. We also examine the effect of grain size distribution in the galaxy by changing the dense gas fraction, which controls the balance between coagulation and shattering, clarifying that the difference in the grain size distribution significantly affects $X_\mathrm{CO}$ even if the dust-to-gas ratio is the same. The star formation time-scale, which controls the speed of metal enrichment also affects the metallicity at which the CO abundance rapidly increases (or $X_\mathrm{CO}$ drops). We also propose dust-based formulae for $X_\mathrm{CO}$, which need further tests for establishing their usefulness.

Recent observations made by the JWST have revealed a number of massive galaxies at high redshift ($z$). The presence of these galaxies appears at odds with the current $\Lambda$CDM cosmology. Here we investigate the possibility of alleviating the tension by incorporating uncertainties from three sources in counting massive galaxies at high $z$: cosmic variance, error in stellar mass estimate, and contribution by backsplash. We find that each of the sources can significantly increase the cumulative stellar mass density $\rho_*(>M_*)$ at the high-mass end, and the combination of them can boost the density by more than one order of magnitude. Assuming a star formation efficiency of $\epsilon_* \sim 0.5$, cosmic variance alone can reduce the tension to $2\sigma$ level, except the most massive galaxy at $z=8$. Including in addition a lognormal dispersion with a width of 0.3 dex in the stellar mass can bring the observed stellar mass density at $z \sim 7 - 10$ to the $2\sigma$ range of the cosmic variance. The tension is completely eliminated when gas stripped from backsplash halos is also taken into account. Our results highlight the importance of fully modeling uncertainties when interpreting observational data of rare objects. We use the constrained simulation, ELUCID, to investigate the descendants of high $z$ massive galaxies. We find that a significant portion of these galaxies end up in massive halos with mass $M_{\rm halo} > 10^{13} h^{-1}M_\odot $ at $z=0$. A large fraction of central galaxies in $M_{\rm halo} \geqslant 10^{14.5} h^{-1}M_\odot$ halos today are predicted to contain significant amounts of ancient stars formed in massive galaxies at $z\sim 8$. This prediction can be tested by studying the structure and stellar population of central galaxies in present-day massive clusters.

The proton synchrotron radiation is considered as the origin of high-energy emission of blazars at times. However, extreme physical parameters are often required. In this work, we propose an analytical method to study the parameter space when applying the proton synchrotron radiation to fit the keV, GeV, and very-high-energy emission of blazar jets. We find that proton synchrotron radiation can fit the high-energy hump when it peaks beyond tens GeV without violating basic observations and theories. For the high-energy hump peaked around GeV band, extreme parameters, such as a super-Eddington jet power and a very strong magnetic field, are required. For the high-energy hump peaked around keV band, if an acceptable parameter space can be found depends on the object's keV luminosity.

We analysed the nuclear region of all 56 early-type galaxies from the DIVING$^\mathrm{3D}$ Project, which is a statistically complete sample of objects that contains all 170 galaxies of the Southern Hemisphere with B < 12.0 mag and galactic latitude |b| < 15$^{\circ}$. Observations were performed with the Integral Field Unit of the Gemini Multi-Object Spectrograph. Emission lines were detected in the nucleus of 86$\pm$5% of the objects. Diagnostic diagrams were used to classify 52$\pm$7% of the objects as LINERs or Seyferts, while the other 34$\pm$6% galaxies without H$\beta$ or [O III] lines in their spectra were classified as weak emission line objects. Transition Objects are not seen in the sample, possibly because the seeing-limited data cubes of the objects allow one to isolate the nuclei of the galaxies from their circumnuclear regions, avoiding contamination from H II regions. A broad line region is seen in 29$\pm$6% of the galaxies. Of the 48 galaxies with emission-line nuclei, 41 have signs of AGNs. Some objects have also indications of shocks in their nuclei. Lenticular galaxies are more likely to have emission lines than ellipticals. Also, more luminous objects have higher [N II]/H$\alpha$ ratios, which may be associated with the mass-metalicity relation of galaxies. A direct comparison of our results with the Palomar Survey indicates that the detection rates of emission lines and also of type 1 AGNs are higher in the DIVING$^\mathrm{3D}$ objects. This is a consequence of using a more modern instrument with a better spatial resolution than the Palomar Survey observations.

Fast radio bursts (FRBs) are mysterious astronomical phenomena, and it is still uncertain whether they consist of multiple types. In this study we use two nonlinear dimensionality reduction algorithms - Uniform Manifold Approximation and Projection (UMAP) and t-distributed stochastic neighbour embedding (t-SNE) - to differentiate repeaters from apparently non-repeaters in FRBs. Based on the first Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB catalogue, these two methods are applied to standardized parameter data and image data from a sample of 594 sub-bursts and 535 FRBs, respectively. Both methods are able to differentiate repeaters from apparently non-repeaters. The UMAP algorithm using image data produces more accurate results and is a more model-independent method. Our result shows that in general repeater clusters tend to be narrowband, which implies a difference in burst morphology between repeaters and apparently non-repeaters. We also compared our UMAP predictions with the CHIME/FRB discovery of 6 new repeaters, the performance was generally good except for one outlier. Finally, we highlight the need for a larger and more complete sample of FRBs.

In principle, the local cosmic void can be simply modeled by the spherically symmetric Lemaitre-Tolman-Bondi (LTB) metric. In practice, the real local cosmic void is probably not spherically symmetric. In this paper, to reconstruct the realistic profile of the local cosmic void, we divide it into several segments. Each segment with certain solid angle is modeled by its own LTB metric. Meanwhile, we divide the 1048 type Ia supernovae (SNIa) of Pantheon into corresponding subsets according to their distribution in the galactic coordinate system. Obviously, each SNIa subset can only be used to reconstruct the profile of one segment. Finally, we can patch together an irregular profile for the local cosmic void with the whole Pantheon sample. But our constraints are too weak to challenge the cosmic homogeneity and the cosmic isotropy.

Stellar chemical abundances are crucial and fundamental in astrophysics. However, they could suffer from substantial systematic errors according to several investigations but still lack calibrations in bulk. By using Gaia wide binaries, we find the temperature-dependent bias between the two binary components for [Fe/H] and [alpha/Fe] measurements from the LAMOST low-resolution spectra and Gaia RVS spectra. At Teff=4000 K, the LAMOST [Fe/H] is significantly underestimated by approximately 0.4 dex, compared with its typical uncertainty of 0.1 dex. Its [alpha/Fe] is overestimated by about 0.2 dex. For Gaia, the underestimation of [M/H] and overestimation of [alpha/Fe] becomes pronounced near 7000 K with smaller magnitudes. We perform an internal calibration by minimizing the differences between binary components and provide the correction curves. After corrections, the standard deviations of the residuals compared to the PASTEL catalog decrease from about 0.045/0.1 to 0.02/0.043 for LAMOST and Gaia, respectively. The chemical homogeneity of the open cluster M 44 is also improved by a factor of two. We stress that the underestimation of [Fe/H] could lead to an overestimation of binary fractions when selecting binary stars by the excess of luminosity. The method of this work could be applied to other data-sets in the future. Our results will benefit statistic studies that use LAMOST and Gaia samples with a wide temperature range.

We report the discovery of a 100 kpc HI tail in the merging galaxy pair NGC 4490/85 detected by the Five-Hundred-meter Aperture Spherical radio Telescope (FAST). The tidal tails extended in both the south and north directions, and they are much longer than that reported previously based on the VLA interferometric maps. The NGC 4490/85 is surrounded by a large gas envelope, and a starburst low metallicity dwarf galaxy MAPS 1231+42 is found to be connected with the gas envelope, indicating that galaxy interaction trigged the intense star formation in it. Based on the fact that the metallicity in MAPS 1231+42 is one order of magnitude lower than that in the two disks of NGC 4490 and NGC 4485, we speculate that the gas near this galaxy should be primordial and could be due to gas inflow from the circum-galactic medium (CGM). We also found a collimated gas component pointing at a nearby dwarf galaxy KK 149, suggesting that this galaxy might also be interacting with the NGC 4490 pair. We discuss the possible origin of the long tidal tails and the extended gas envelope in this merging system based on the new data from FAST.

The joint detection of GW signals by a network of instruments will increase the detecting ability of faint and far GW signals with higher signal-to-noise ratios (SNRs), which could improve the ability of detecting the lensed GWs as well, especially for the 3rd generation detectors, e.g. Einstein Telescope (ET) and Cosmic Explorer (CE). However, identifying Strongly Lensed Gravitational Waves (SLGWs) is still challenging. We focus on the identification ability of 3G detectors in this article. We predict and analyze the SNR distribution of SLGW signals and prove only 50.6\% of SLGW pairs detected by ET alone can be identified by Lens Bayes factor (LBF), which is a popular method at present to identify SLGWs. For SLGW pairs detected by CE\&ET network, owing to the superior spatial resolution, this number rises to 87.3\%. Moreover, we get an approximate analytical relation between SNR and LBF. We give clear SNR limits to identify SLGWs and estimate the expected yearly detection rates of galaxy-scale lensed GWs that can get identified with 3G detector network.

Time-delay interferometry (TDI) is a crucial technology for space-based gravitational wave detectors. Previous studies have identified the optimal TDI configuration for the first-generation. In this research, we used an Algebraic approach theory to describe the TDI space and employed a method to maximize the signal-to-noise ratio (SNR) to derive the optimal TDI combination for the second-generation. When this combination is used in the sensitivity curve, we observed enhancements of up to 1.91 times in the low-frequency domain and 2 to 3.5 times in the high-frequency domain compared to the Michelson combination. Furthermore, changes in the detector index significantly affect the optimization effect. We also present detection scenarios for several low-frequency gravitational wave sources. Compared to the first-generation TDI optimization, the SNR value for verification double white dwarfs (DWD) and the detection rate for DWD increase by 16.5%.

The vector Apodizing Phase Plate (vAPP) is a pupil plane coronagraph that suppresses starlight by forming a dark hole in its point spread function (PSF). The unconventional and non-axisymmetrical PSF arising from the phase modification applied by this coronagraph presents a special challenge to post-processing techniques. We aim to implement a recently developed post-processing algorithm, temporal reference analysis of planets (TRAP) on vAPP coronagraphic data. The property of TRAP that uses non-local training pixels, combined with the unconventional PSF of vAPP, allows for more flexibility than previous spatial algorithms in selecting reference pixels to model systematic noise. Datasets from two types of vAPPs are analysed: a double grating-vAPP (dgvAPP360) that produces a single symmetric PSF and a grating-vAPP (gvAPP180) that produces two D-shaped PSFs. We explore how to choose reference pixels to build temporal systematic noise models in TRAP for them. We then compare the performance of TRAP with previously implemented algorithms that produced the best signal-to-noise ratio (S/N) in companion detections in these datasets. We find that the systematic noise between the two D-shaped PSFs is not as temporally associated as expected. Conversely, there is still a significant number of systematic noise sources that are shared by the dark hole and the bright side in the same PSF. We should choose reference pixels from the same PSF when reducing the dgvAPP360 dataset or the gvAPP180 dataset with TRAP. In these datasets, TRAP achieves results consistent with previous best detections, with an improved S/N for the gvAPP180 dataset.

We report the discovery of the isolated neutron star (INS) candidates eRASSU J065715.3+260428 and eRASSU J131716.9-402647 from the Spectrum Roentgen Gamma (SRG) eROSITA All-Sky Survey. Selected for their soft X-ray emission and absence of catalogued counterparts, both objects were recently targeted with the Large Binocular Telescope and the Southern African Large Telescope. The absence of counterparts down to deep optical limits (25 mag, 5$\sigma$) and, as a result, large X-ray-to-optical flux ratios in both cases strongly suggest an INS nature. The X-ray spectra of both sources are well described by a simple absorbed blackbody, whereas other thermal and non-thermal models (e.g. a hot-plasma emission spectrum or power law) are disfavoured by the spectral analysis. Within the current observational limits, and as expected for cooling INSs, no significant variation ($>2\sigma$) has been identified over the first two-year time span of the survey. Upcoming dedicated follow-up observations will help us to confirm the candidates' nature.

The statistics of the magnetic field line separation provide insight into how a bundle of field lines spreads out and the dispersion of non-thermal particles in a turbulent environment, which underlies various astrophysical phenomena. Its diffusive character depends on the distance along the field line, the initial separation, and the characteristics of the magnetic turbulence. This work considers the separation of two magnetic field lines in general transverse turbulence in terms of the magnetic power spectrum in three-dimensional wavenumber space. We apply non-perturbative methods using Corrsin's hypothesis and assume random ballistic decorrelation to calculate the ensemble average field line separation for general transverse magnetic turbulence. For 2D+slab power spectra, our analytic formulae and computer simulations give similar results, especially at low slab fraction. Our analytical expression also demonstrates several features of field line separation that are verified by computer simulations.

X-ray appearance of normal galaxies is mainly determined by X-ray binaries powered by accretion onto a neutron star or a stellar mass black hole. Their populations scale with the star-formation rate and stellar mass of the host galaxy and their X-ray luminosity distributions show a significant split between star-forming and passive galaxies, both facts being consequences of the dichotomy between high- and low-mass X-ray binaries. Metallicity, IMF and stellar age dependencies, and dynamical formation channels add complexity to this picture. The numbers of high-mass X-ray binaries observed in star-forming galaxies indicate quite high probability for a massive star to become an accretion powered X-ray source once upon its lifetime. This explains the unexpectedly high contribution of X-ray binaries to the Cosmic X-ray Background, of the order of $\sim 10\%$, mostly via X-ray emission of faint star-forming galaxies located at moderate redshifts which may account for the unresolved part of the CXB. Cosmological evolution of the $L_X-{\rm SFR}$ relation can make high-mass X-ray binaries a potentially significant factor in (pre)heating of intergalactic medium in the early Universe.

In this paper, we analyzed 12 years of Fermi LAT gamma-ray data towards three nearby giant molecular clouds (GMCs), i.e., R~CrA, Chamaeleon, and Lupus. We calibrated the gas column density of these regions by using the Planck dust opacity map as well as the Gaia extinction map. With both the gamma-ray observations and gas column density maps, we derived the cosmic ray densities in the three GMCs. We found the derived CR spectra have almost the same shape but significantly different normalizations, which may reflect that the distributions of CRs in the vicinity of solar systems are inhomogeneous.

M dwarfs are the most abundant stars in the Solar Neighborhood and they are prime targets for searching for rocky planets in habitable zones. Consequently, a detailed characterization of these stars is in demand. The spectral sub-type is one of the parameters that is used for the characterization and it is traditionally derived from the observed spectra. However, obtaining the spectra of M dwarfs is expensive in terms of observation time and resources due to their intrinsic faintness. We study the performance of four machine-learning (ML) models: K-Nearest Neighbor (KNN), Random Forest (RF), Probabilistic Random Forest (PRF), and Multilayer Perceptron (MLP), in identifying the spectral sub-types of M dwarfs at a grand scale by deploying broadband photometry in the optical and near-infrared. We trained the ML models by using the spectroscopically identified M dwarfs from the Sloan Digital Sky Survey Data Release (SDSS) 7, together with their photometric colors that were derived from the SDSS, Two-Micron All-Sky Survey, and Wide-field Infrared Survey Explorer. We found that the RF, PRF, and MLP give a comparable prediction accuracy, 74%, while the KNN provides slightly lower accuracy, 71%. We also found that these models can predict the spectral sub-type of M dwarfs with ~99% accuracy within +/-1 sub-type. The five most useful features for the prediction are r-z, r-i, r-J, r-H, and g-z, and hence lacking data in all SDSS bands substantially reduces the prediction accuracy. However, we can achieve an accuracy of over 70% when the r and i magnitudes are available. Since the stars in this study are nearby (d~1300 pc for 95% of the stars), the dust extinction can reduce the prediction accuracy by only 3%. Finally, we used our optimized RF models to predict the spectral sub-types of M dwarfs from the Catalog of Cool Dwarf Targets for TESS, and we provide the optimized RF models for public use.

The near-Earth asteroid, Kamo`oalewa (469219), is one of a small number of known quasi-satellites of Earth. Numerical simulations show that it transitions between quasi-satellite and horseshoe orbital states on centennial timescales, maintaining this dynamics over megayears. Its reflectance spectrum suggest a similarity to lunar silicates. Considering its Earth-like orbit and its physical resemblance to lunar surface materials, we explore the hypothesis that it might have originated as a debris-fragment from a meteoroidal impact with the lunar surface. We carry out numerical simulations of the dynamical evolution of particles launched from different locations on the lunar surface with a range of ejection velocities. As these ejecta escape the Earth-Moon environment and evolve into heliocentric orbits, we find that a small fraction of launch conditions yield outcomes that are compatible with Kamo`oalewa's dynamical behavior. The most favored conditions are launch velocities slightly above the escape velocity from the trailing lunar hemisphere.

M31 has experienced a recent tumultuous merger history as evidenced from the many substructures that are still present in its inner halo, particularly the G1-Clump, NE- and W- shelves, and the Giant Stream (GS). We present planetary nebulae (PNe) line-of-sight velocity (LOSV) measurements covering the entire spatial extent of these four substructures. We further use predictions for the satellite and host stellar particle phase space distributions for a major merger (mass ratio = 1:4) simulation to help interpret the data. The measured PN LOSVs for the two shelves and GS are consistent with those from red giant branch stars. Their projected radius vs. LOSV phase space, links the formation of these substructures in a single unique event, consistent with a major merger. We find the G1-clump to be dynamically cold compared to the M31 disc ($\rm\sigma_{LOS, PN}=27$ km s$^{-1}$), consistent with pre-merger disc material. Such a structure can not form in a minor merger (mass ratio $\sim$1:20), and is therefore a smoking gun for the recent major merger event in M31. The simulation also predicts the formation of a predominantly in-situ halo from splashed-out pre-merger disc material, in qualitative agreement with observations of a metal-rich inner halo in M31. Juxtaposed with previous results for its discs, we conclude that M31 has had a recent (2.5 - 4 Gyr ago) `wet' major merger with the satellite falling along the GS, heating the pre-merger disc to form the M31 thicker disc, rebuilding the M31 thin disc, and creating the aforementioned inner-halo substructures.

In November 2021, Solar Orbiter started its nominal mission phase. The remote-sensing instruments on board the spacecraft acquired scientific data during three observing windows surrounding the perihelion of the first orbit of this phase. The aim of the analysis is the detection of magnetohydrodynamic (MHD) wave modes in an active region by exploiting the capabilities of spectropolarimetric measurements. The High Resolution Telescope (HRT) of the Polarimetric and Helioseismic Imager (SO/PHI) on board the Solar Orbiter acquired a high-cadence data set of an active region. This is studied in the paper. B-$\omega$ and phase-difference analyses are applied on line-of-sight velocity and circular polarization maps and other averaged quantities. We find that several MHD modes at different frequencies are excited in all analysed structures. The leading sunspot shows a linear dependence of the phase lag on the angle between the magnetic field and the line of sight of the observer in its penumbra. The magnetic pore exhibits global resonances at several frequencies, which are also excited by different wave modes. The SO/PHI measurements clearly confirm the presence of magnetic and velocity oscillations that are compatible with one or more MHD wave modes in pores and a sunspot. Improvements in modelling are still necessary to interpret the relation between the fluctuations of different diagnostics.

The newest generation of radio telescopes are able to survey large areas with high sensitivity and cadence, producing data volumes that require new methods to better understand the transient sky. Here we describe the results from the first citizen science project dedicated to commensal radio transients, using data from the MeerKAT telescope with weekly cadence. Bursts from Space: MeerKAT was launched late in 2021 and received ~89000 classifications from over 1000 volunteers in 3 months. Our volunteers discovered 142 new variable sources which, along with the known transients in our fields, allowed us to estimate that at least 2.1 per cent of radio sources are varying at 1.28 GHz at the sampled cadence and sensitivity, in line with previous work. We provide the full catalogue of these sources, the largest of candidate radio variables to date. Transient sources found with archival counterparts include a pulsar (B1845-01) and an OH maser star (OH 30.1-0.7), in addition to the recovery of known stellar flares and X-ray binary jets in our observations. Data from the MeerLICHT optical telescope, along with estimates of long time-scale variability induced by scintillation, imply that the majority of the new variables are active galactic nuclei. This tells us that citizen scientists can discover phenomena varying on time-scales from weeks to several years. The success both in terms of volunteer engagement and scientific merit warrants the continued development of the project, whilst we use the classifications from volunteers to develop machine learning techniques for finding transients.

The hypothesis that tidal forces on the Sun are related to the modulations of the solar-activity cycle has gained increasing attention. The works proposing physical mechanisms of planetary action via tidal forcing have in common that quasi-alignments between Venus, Earth, and Jupiter (V-E-J configurations) would provide a basic periodicity of $\approx 11.0$ years able to synchronize the operation of solar dynamo with these planetary configurations. Nevertheless, the evidence behind this particular tidal forcing is still controversial. In this context we develop, for the first time, the complete Sun's tide-generating potential (STGP) in terms of a harmonic series, where the effects of different planets on the STGP are clearly separated and identified. We use a modification of the spectral analysis method devised by Kudryavtsev (J. Geodesy. 77, 829, 2004; Astron. Astrophys. 471, 1069, 2007b) that permits to expand any function of planetary coordinates to a harmonic series over long time intervals. We build a catalog of 713 harmonic terms able to represent the STGP with a high degree of precision. We look for tidal forcings related to V-E-J configurations and specifically the existence of periodicities around $11.0$ years. Although the obtained tidal periods range from $\approx$ 1000 years to 1 week, we do not find any $\approx$ 11.0 years period. The V-E-J configurations do not produce any significant tidal term at this or other periods. The Venus tidal interaction is absent in the 11-year spectral band, which is dominated by Jupiter's orbital motion. The planet that contributes the most to the STGP in three planets configurations, along with Venus and Earth, is Saturn. An $\approx 11.0$ years tidal period with a direct physical relevance on the 11-year-like solar-activity cycle is highly improbable.

The Epoch of Reionization (EoR) neutral Hydrogen (HI) 21-cm signal evolves significantly along the line-of-sight (LoS) due to the light-cone (LC) effect. It is important to accurately incorporate this in simulations in order to correctly interpret the signal. 21-cm LC simulations are typically produced by stitching together slices from a finite number $(N_{\rm RS})$ of ''reionization snapshot'', each corresponding to a different stage of reionization. In this paper, we have quantified the errors in the 21-cm LC simulation due to the finite value of $N_{\rm RS}$. We show that this can introduce large discontinuities $(> 200 \%)$ at the stitching boundaries when $N_{\rm RS}$ is small $(= 2,4)$ and the mean neutral fraction jumps by $\delta \bar{x}_{\rm HI} = 0.2,0.1$ respectively at the stitching boundaries. This drops to $17 \%$ for $N_{\rm RS} = 13$ where $\delta \bar{x}_{\rm HI}=0.02$. We present and also validate a method for mitigating this error by increasing $N_{\rm RS}$ without a proportional increase in the computational costs which are mainly incurred in generating the dark matter and halo density fields. Our method generates these fields only at a few redshifts, and interpolates them to generate reionization snapshots at closely spaced redshifts. We use this to generate 21-cm LC simulations with $N_{\rm RS} = 26,51,101$ and $201$, and show that the errors go down as $N_{\rm RS}^{-1}$.

Understanding the temperature structure of protoplanetary disks is crucial for answering the fundamental question of when and where in the disks rocky planets like our own form. However, the thermal structure of the inner few au of the disks is poorly understood not only because of lack of observational constraints but also because of the uncertainty of accretion heating processes. Here, we propose thermal tomography of the inner regions of protoplanetary disks with the ngVLA and ALMA. The proposed approach is based on the assumption that the inner disk regions are optically thick at submillimeter wavelengths but are marginally optically thin at longer millimeter wavelengths. By combining high-resolution millimeter continuum images from the ngVLA with submillimeter images at comparable resolutions from ALMA, we will be able to reconstruct the radial and vertical structure of the inner few au disk regions. We demonstrate that the thermal tomography we propose can be used to constrain the efficiency of midplane accretion heating, a process that controls the timing of snow-line migration to the rocky planet-forming region, in the few au regions of protoplanetary disks at a distance of 140 pc.

Context: The conventional approach to direct imaging has been the use of a single aperture coronagraph with wavefront correction via extreme adaptive optics. Such systems are limited to observing beyond an inner working (IWA) of a few {\lambda}/D. Nulling interferometry with two or more apertures will enable detections of companions at separations at and beyond the formal diffraction limit. Aims: This paper evaluates the astrophysical potential of a kernel-nuller as the prime high-contrast imaging mode of the Very Large Telescope Interferometer (VLTI). Methods: By taking into account baseline projection effects which are induced by Earth rotation, we introduce some diversity in the response of the nuller as a function of time. This response is depicted by transmission maps. We also determine whether we can extract the astrometric parameters of a companion from the kernel outputs, which are the primary intended observable quantities of the kernel-nuller. This then leads us to comment on the characteristics of a possible observing program for the discovery of exoplanets. Results: We present transmission maps for both the raw nuller outputs and their subsequent kernel outputs. To further examine the properties of the kernel-nuller, we introduce maps of the absolute value of the kernel output. We also identify 38 targets for the direct detection of exoplanets with a kernel-nuller at the focus of the VLTI. Conclusions: With continued upgrades of the VLTI infrastructure that will reduce fringe tracking residuals, a kernel-nuller would enable the detection of young giant exoplanets at separations < 10 AU, where radial velocity and transit methods are more sensitive.

Stellar winds of massive stars are known to be driven by line absorption of UV photons, a mechanism which is prone to instabilities, causing the wind to be clumpy. The clumpy structure hampers wind mass-loss estimates, limiting our understanding of massive star evolution. The wind structure also impacts accretion in high-mass X-ray binary (HMXB) systems. We analyse the wavelength-dependent variability of X-ray absorption in the wind to study its structure. Such an approach is possible in HMXBs, where the compact object serves as an X-ray backlight. We probe different parts of the wind by analysing data taken at superior and inferior conjunction. We apply excess variance spectroscopy to study the wavelength-dependent soft X-ray variability of the HMXB Cygnus X-1 in the low/hard spectral state. Excess variance spectroscopy quantifies the variability of an object above the statistical noise as a function of wavelength, which allows us to study the variability of individual spectral lines. As one of the first studies, we apply this technique to high-resolution gratings spectra provided by Chandra, accounting for various systematic effects. The frequency dependence is investigated by changing the time binning. The strong orbital phase dependence we observe in the excess variance is consistent with column density variations predicted by a simple model for a clumpy wind. We identify spikes of increased variability with spectral features found by previous spectroscopic analyses of the same data set, most notably from silicon in over-dense clumps in the wind. In the silicon line region, the variability power is redistributed towards lower frequencies, hinting at increased line variability in large clumps. In prospect of the microcalorimetry missions that are scheduled to launch within the next decade, excess variance spectra present a promising approach to constrain the wind structure.

Relativistic jets launched by Active Galactic Nuclei are among the most powerful particle accelerators in the Universe. The emission over the entire electromagnetic spectrum of these relativistic jets can be extremely variable with scales of variability from less than few minutes up to several years. These variability patterns, which can be very complex, contain information about the acceleration processes of the particles and the area(s) of emission. Thanks to its sensitivity, five-to twenty-times better than the current generation of Imaging Atmospheric Cherenkov Telescopes depending on energy, the Cherenkov Telescope Array will be able to follow the emission from these objects with a very accurate time sampling and over a wide spectral coverage from 20 GeV to > 20 TeV and thus reveal the nature of the acceleration processes at work in these objects. We will show the first results of our lightcurve simulations and long-term behavior of AGN as will be observed by CTA, based on state-of-art particle acceleration models.

Quasars emission is highly variable, and this variability gives us clues to understand the accretion process onto supermassive black holes. We can expect variability properties to correlate with the main physical properties of the accreting black hole, i.e., its mass and accretion rate. It has been established that the relative amplitude of variability anti-correlates with the accretion rate.The dependence of the variance on black hole mass has remained elusive, and contradicting results, including positive, negative, or no correlation, have been reported. In this work, we show that the key to these contradictions lies in the timescales of variability studied (e.g., the length of the light curves available). By isolating the variance on different timescales as well as mass and accretion rate bins we show that there is indeed a negative correlation between black hole mass and variance and that this anti-correlation is stronger for shorter timescale fluctuations. The behavior can be explained in terms of a universal variability power spectrum for all quasars, resembling a broken power law where the variance is constant at low temporal frequencies and then drops continuously for frequencies higher than a characteristic frequency $f_b$, where $f_b$ correlates with the black hole mass. Furthermore, to explain all the variance results presented here, not only the normalization of this power spectrum must anti-correlate with the accretion rate, but also the shape of the power spectra at short timescales must depend on this parameter as well.

CONTEXT: The unparalleled characteristics of Gaia photometry make it an excellent choice to study stellar variability. AIMS: To measure the phot. dispersion in G+G_BP+G_RP of the 145 677 450 Gaia DR3 5-parameter sources with G <= 17 mag and G_BP-G_RP with -1.0 to 8.0 mag. To use that unbiased sample to analyze stellar variability in the Milky Way, LMC, and SMC. METHODS: We convert from magnitude uncertainties to the observed phot. dispersions, calculate the instrumental component as a function of apparent magnitude and color, and use it to transform the observed dispersions into the astrophysical ones. We give variability indices in the three bands for the whole sample. We use the subsample of Rimoldini et al. that includes light curves and variability types to calibrate our results and establish their limitations. RESULTS: We use information from the MW, LMC, and SMC CAMDs to discuss variability across the HRD. Most WDs and sdBs are variable and follow a distribution in s_G peaking around 12 mmag but variability decreases for the former with age. The MS region in the Gaia CAMD has an s_G distribution peaks at low values (~1-2 mmag) and has a large tail dominated by EBs, RR Lyr stars, and YSOs. RC stars are characterized by little variability, with their s_G distribution peaking at 1 mmag or less. The stars in the PMS region are highly variable, with a power law distribution in s_G with slope 2.75 and a cutoff for values lower than 7 mmag. The luminous red stars region of the Gaia CAMD has the highest variability, with its extreme dominated by AGB stars and with a power law in s_G with a slope of ~2.2 that extends from there to a cutoff of 7 mmag. We show that our method can be used to search for LMC Cepheids. We analyze four stellar clusters with O stars and detect a strong difference in s_G between stars that are already in the MS and those that are still in the PMS. [ABRIDGED]

Next generation surveys will provide us with an unprecedented number of detections of supernovae Type Ia and gravitational wave merger events. Cross-correlations of such objects offer novel and powerful insights into the large-scale distribution of matter in the universe. Both of these sources carry information on their luminosity distance, but remain uninformative about their redshifts; hence their clustering analyses and cross-correlations need to be carried out in luminosity distance space, as opposed to redshift space. In this paper, we calculate the full expression for the number count fluctuation in terms of a perturbation to the observed luminosity distance. We find the expression to differ significantly from the one commonly used in redshift space. Furthermore, we present a comparison of the number count angular power spectra between luminosity distance and redshift spaces. We see a wide divergence between the two at large scales, and we note that lensing is the main contribution to such differences. On such scales and at higher redshifts the difference between the angular power spectra in luminosity distance and redshift spaces can be roughly 50$\%$. We also investigate cross-correlating different redshift bins using different tracers, i.e. one in luminosity distance space and one in redshift, simulating the cross-correlation angular power spectrum between background gravitational waves/supernovae and foreground galaxies. Finally, we show that in a cosmic variance limited survey, the relativistic corrections to the density-only term ought to be included.

Recently multi-field inflation models that can produce large scalar fluctuations on small scales have drawn a lot of attention, primarily because they could lead to primordial black hole production and generation of large second-order gravitational waves. In this work, we focus on models where the scalar fields responsible for inflation live on a hyperbolic field space. In this case, geometrical destabilisation and non-geodesic motion are responsible for the peak in the scalar power spectrum. We present new results for scalar non-Gaussianity and discuss its dependence on the model's parameters. On scales around the peak, we typically find that the non-Gaussianity is large and close to local in form. We validate our results by employing two different numerical techniques, utilising the transport approach, based on full cosmological perturbation theory, and the $\delta N$ formalism, based on the separate universe approximation. We discuss implications of our results for the perturbativity of the underlying theory, focusing in particular on versions of these models with potentially relevant phenomenology at interferometer scales.

Most active galactic nuclei (AGNs) at the center of the nearby galaxies have super-massive black holes accreting at sub-Eddington rates through hot accretion flows or radiatively inefficient accretion flows, which efficiently produce jets. The association of radio and X-ray flares with the knot ejection from M81* inspires us to model its multi-wavelength spectral energy distribution during these flares to constrain the physical parameters of the jet. Moreover, we construct a long-term light curve in X-rays to identify the flares in the available data and constrain the jet parameters during those periods. The jet activity may vary on short and long time scales, which may produce flares in different frequency bands. The spectral energy distributions from radio to X-ray during the quiescent as well as flaring states are found to be satisfactorily explained by synchrotron emission of relativistic electrons from a single zone. The variation in the values of the jet parameters during the different states is shown and compared with high synchrotron peaked blazars.

We present an investigation of the open cluster Trumpler 2 using Gaia DR3 photometric, astrometric and spectroscopic data. 92 stars were identified as likely members of the cluster, with membership probabilities greater than 0.5. The mean proper-motion components of the cluster are derived as ($\mu_{\alpha}\cos \delta$, $\mu_{\delta}$)=($1.494 \pm 0.004$, $-5.386 \pm 0.005$) mas yr$^{-1}$. By comparing the Gaia based colour-magnitude diagram with the PARSEC isochrones scaled to $z=0.0088$, age, distance modulus and reddening are simultaneously estimated as $t=110 \pm 10$ Myr, $\mu=10.027 \pm0.149$ mag and $E(G_{\rm BP}-G_{\rm RP})=0.452\pm 0.019$ mag, respectively. The total mass of the cluster is estimated as 162 $M/M_{\odot}$ based on the stars with membership probabilities $P > 0$. The Mass function slope is derived to be $\Gamma = 1.33 \pm 0.13$ for Trumpler 2. This value is in a good agreement with that of of Salpeter. Galactic orbit analyses show that the Trumpler 2 orbits in a boxy pattern outside the solar circle and belongs to the young thin-disc component of the Galaxy.

Current space-based missions, such as the Transiting Exoplanet Survey Satellite (TESS), provide a large database of light curves that must be analysed efficiently and systematically. In recent years, deep learning (DL) methods, particularly convolutional neural networks (CNN), have been used to classify transit signals of candidate exoplanets automatically. However, CNNs have some drawbacks; for example, they require many layers to capture dependencies on sequential data, such as light curves, making the network so large that it eventually becomes impractical. The self-attention mechanism is a DL technique that attempts to mimic the action of selectively focusing on some relevant things while ignoring others. Models, such as the Transformer architecture, were recently proposed for sequential data with successful results. Based on these successful models, we present a new architecture for the automatic classification of transit signals. Our proposed architecture is designed to capture the most significant features of a transit signal and stellar parameters through the self-attention mechanism. In addition to model prediction, we take advantage of attention map inspection, obtaining a more interpretable DL approach. Thus, we can identify the relevance of each element to differentiate a transit signal from false positives, simplifying the manual examination of candidates. We show that our architecture achieves competitive results concerning the CNNs applied for recognizing exoplanetary transit signals in data from the TESS telescope. Based on these results, we demonstrate that applying this state-of-the-art DL model to light curves can be a powerful technique for transit signal detection while offering a level of interpretability.

Mean plane measurements of the Kuiper Belt from observational data are of interest for their potential to test dynamical models of the solar system. Recent measurements have yielded inconsistent results. Here we report a measurement of the Kuiper Belt's mean plane with a sample size more than twice as large as in previous measurements. The sample of interest is the non-resonant Kuiper belt objects, which we identify by using machine learning on the observed Kuiper Belt population whose orbits are well-determined. We estimate the measurement error with a Monte Carlo procedure. We find that the overall mean plane of the non-resonant Kuiper Belt (semimajor axis range 35-150 au) and also that of the classical Kuiper Belt (semimajor axis range 42-48 au) are both close to (within about 0.7 degrees) but distinguishable from the invariable plane of the solar system to greater than 99.7% confidence. When binning the sample into smaller semimajor axis bins, we find the measured mean plane mostly consistent with both the invariable plane and the theoretically expected Laplace surface forced by the known planets. Statistically significant discrepancies are found only in the semimajor axis ranges 40.3-42 au and 45-50 au; these ranges are in proximity to a secular resonance and Neptune's 2:1 mean motion resonance where the theory for the Laplace surface is likely to be inaccurate. These results do not support a previously reported anomalous warp at semimajor axes above 50 au.

The extraordinary 2021 September-October outburst of Centaur 29P/Schwassmann-Wachmann 1 afforded an opportunity to test the composition of primitive Kuiper disk material at high sensitivity. We conducted nearly simultaneous multi-wavelength spectroscopic observations of 29P/Schwassmann-Wachmann 1 using iSHELL at the NASA Infrared Telescope Facility and nFLASH at the Atacama Pathfinder EXperiment (APEX) on 2021 October 6, with follow-up APEX/nFLASH observations on 2021 October 7 and 2022 April 3. This coordinated campaign between near-infrared and radio wavelengths enabled us to sample molecular emission from a wealth of coma molecules and to perform measurements that cannot be accomplished with either wavelength alone. We securely detected CO emission on all dates with both facilities, including velocity-resolved spectra of the CO (J=2-1) transition with APEX/nFLASH and multiple CO (v=1-0) rovibrational transitions with IRTF/iSHELL. We report rotational temperatures, coma kinematics, and production rates for CO and stringent (3-sigma) upper limits on abundance ratios relative to CO for CH4, C2H6, CH3OH, H2CO, CS, and OCS. Our upper limits for CS/CO and OCS/CO represent their first values in the literature for this Centaur. Upper limits for CH4, C2H6, CH3OH, and H2CO are the most stringent reported to date, and are most similar to values found in ultra CO-rich Oort cloud comet C/2016 R2 (PanSTARRS), which may have implications for how ices are preserved in cometary nuclei. We demonstrate the superb synergy of coordinated radio and near-infrared measurements, and advocate for future small body studies that jointly leverage the capabilities of each wavelength.

The afterglows of exceptionally bright gamma-ray bursts (GRBs) can reveal the angular structure of their ultra-relativistic jets after they emerge from the confining medium, e.g. the progenitor's stellar envelope in long-soft GRBs. These jets appear to have a narrow core (of half-opening angle $\theta_c$), beyond which their kinetic energy drops as a power-law with the angle $\theta$ from the jet's symmetry axis, $E_{k,\rm iso}(\theta)\propto[1+(\theta/\theta_c)^2]^{-a/2}$. The power-law index $a$ reflects the amount of mixing between the shocked jet and confining medium, which depends on the jet's inital magnetization. Weakly magnetized jets undergo significant mixing, leading to shallow ($a\lesssim2$) angular profiles. Here we use the exquisite multi-waveband afterglow observations of GRB 221009A to constrain the jet angular structure using a dynamical model that accounts for both the forward and reverse shocks, for a power-law external density radial profile, $n_{\rm{}ext}\propto{}R^{-k}$. Both the forward shock emission, that dominates the optical and X-ray flux, and the reverse shock emission, that produces the radio afterglow, require a jet with a narrow core ($\theta_c\approx0.021$) and a shallow angular structure ($a\approx0.8$) expanding into a stellar wind ($k\approx2$). In addition, the fraction of shock-heated electrons forming a relativistic power-law energy distribution is limited to $\xi_e\approx10^{-2}$ in both shocks.

The latest generation of radio surveys are now producing sky survey images containing many millions of radio sources. In this context it is highly desirable to understand the performance of radio image source finder (SF) software and to identify an approach that optimises source detection capabilities. We have created Hydra to be an extensible multi-SF and cataloguing tool that can be used to compare and evaluate different SFs. Hydra, which currently includes the SFs Aegean, Caesar, ProFound, PyBDSF, and Selavy, provides for the addition of new SFs through containerisation and configuration files. The SF input RMS noise and island parameters are optimised to a 90\% ''percentage real detections'' threshold (calculated from the difference between detections in the real and inverted images), to enable comparison between SFs. Hydra provides completeness and reliability diagnostics through observed-deep ($\mathcal{D}$) and generated-shallow ($\mathcal{S}$) images, as well as other statistics. In addition, it has a visual inspection tool for comparing residual images through various selection filters, such as S/N bins in completeness or reliability. The tool allows the user to easily compare and evaluate different SFs in order to choose their desired SF, or a combination thereof. This paper is part one of a two part series. In this paper we introduce the Hydra software suite and validate its $\mathcal{D/S}$ metrics using simulated data. The companion paper demonstrates the utility of Hydra by comparing the performance of SFs using both simulated and real images.

We present a comparison between the performance of a selection of source finders using a new software tool called Hydra. The companion paper, Paper~I, introduced the Hydra tool and demonstrated its performance using simulated data. Here we apply Hydra to assess the performance of different source finders by analysing real observational data taken from the Evolutionary Map of the Universe (EMU) Pilot Survey. EMU is a wide-field radio continuum survey whose primary goal is to make a deep ($20\mu$Jy/beam RMS noise), intermediate angular resolution ($15^{\prime\prime}$), 1\,GHz survey of the entire sky south of $+30^{\circ}$ declination, and expecting to detect and catalogue up to 40 million sources. With the main EMU survey expected to begin in 2022 it is highly desirable to understand the performance of radio image source finder software and to identify an approach that optimises source detection capabilities. Hydra has been developed to refine this process, as well as to deliver a range of metrics and source finding data products from multiple source finders. We present the performance of the five source finders tested here in terms of their completeness and reliability statistics, their flux density and source size measurements, and an exploration of case studies to highlight finder-specific limitations.

We present a kinematic analysis of the Small Magellanic Cloud using 3700 spectra extracted from the European Southern Observatory archive. We used data from Gaia and near-infrared photometry to select stellar populations and discard Galactic foreground stars. The sample includes main-sequence, red giant branch and red clump stars, observed with the Fibre Large Array Multi Wavelength Spectrograph. The spectra have a resolving power lambda/Delta(lambda) from 6500 to 38000. We derive radial velocities by employing a full spectrum fitting method using a penalised pixel fitting routine. We obtain a mean radial velocity for the galaxy of 159+/-2 km/s, with a velocity dispersion of 33+/-2 km/s. Our velocities agree with literature estimates for similar (young or old) stellar populations. The radial velocity of stars in the Wing and bar-like structure differ as a consequence of the dynamical interaction with the Large Magellanic Cloud. The higher radial velocity of young main-sequence stars in the bar compared to that of supergiants can be attributed to star formation around 40 Myr ago from gas already influenced by tidal stripping. Similarly, young main-sequence stars in the northern part of the bar, resulting from a prominent episode 25 Myr ago, have a higher radial velocity than stars in the southern part. Radial velocity differences between the northern and southern bar over densities are also traced by giant stars. They are corroborated by studies of the cold gas and proper motion indicating stretching/tidal stripping of the galaxy.

Voids possess a very complex internal structure and dynamics. Using $N$-body simulations we study the hierarchical nature of sub-structures present in the cosmic web (CW). We use the SpineWeb method which provides a complete characterization of the CW into its primary constituents: voids, walls, filaments, and nodes. We aim to characterize the inner compositions of voids by detecting their internal filamentary structure and explore the impact of this on the properties of void galaxies. Using a semi-analytical galaxy evolution model we explore the impact of the CW on several galaxies' properties. We find the fraction of haloes living in various CW components to be a function of their mass, with the majority of the haloes of mass below $10^{12}M_{\odot}/h$, residing in voids and haloes of higher masses distributed mostly in walls. Similarly, in the Stellar-to-Halo mass relationship, we observe an environmental dependence for haloes of masses below $10^{12}M_{\odot}/h$, showing an increased stellar mass fraction for the densest environments. The spin is lower for galaxies in the densest environments for the mass range of $10^{10}-10^{12}M_{\odot}/h$. Finally, we found a strong trend of higher metallicity fractions for filaments and node galaxies, with respect to the full sample, in the range of $M_*<10^{10}M_{\odot}/h$. Our results show that cosmic voids possess an intricate internal network of substructures. This in turn makes them a complex environment for galaxy formation, impacting in an unique way the properties and evolution of the chosen few galaxies that form inside them.

Axions mix with neutral pions after the QCD phase transition through their common coupling to the radiation bath via a Chern-Simons term, as a consequence of the $U(1)$ anomaly. The non-equilibrium effective action that describes this mixing phenomenon is obtained to second order in the coupling of neutral pions and axions to photons. We show that a misaligned axion condensate induces a neutral pion condensate after the QCD phase transition. The dynamics of the pion condensate displays long and short time scales and decays on the longer time scale exhibiting a phenomenon akin to the ``purification'' in a Kaon beam. On the intermediate time scales the macroscopic pion condensate is proportional to a condensate of the abelian Chern-Simons term induced by the axion. We argue that the coupling to the common bath also induces kinetic mixing. We obtain the axion and pion populations, and these exhibit thermalization with the bath. The mutual coupling to the bath induces long-lived axion- neutral pion coherence independent of initial conditions. The framework of the effective action and many of the consequences are more broadly general and applicable to scalar or pseudoscalar particles mixing in a medium.

The electroweak sphaleron process breaks the baryon number conservation within the realms of the Standard Model of particle physics (SM). Recently, it is pointed out that its decoupling may provide the out-of-equilibrium condition required for baryogenesis. In this paper, we study such a scenario taking into account the baryon-number wash-out effect of the sphaleron itself to improve the estimate. We clarify the amount of CP violation required for this scenario to explain the observed asymmetry.

The advent of gravitational wave astronomy has seen a huge influx of new predictions for potential discoveries of beyond the Standard Model fields. The coupling of all fundamental fields to gravity, together with its dominance on large scales, makes gravitational physics a rich laboratory to study fundamental physics. This holds especially true for the search for the elusive dark photon, a promising dark matter candidate. The dark photon is predicted to generate instabilities in a rotating black hole spacetime, birthing a macroscopic Bose-Einstein condensate. These condensates can especially form around super massive black holes, modifying the dynamical inspiralling process. This then opens another window to leverage future space-borne gravitational wave antennas to join the hunt for the elusive dark matter particle. This study builds a preliminary model for the gravitational waveform emitted by such a dressed extreme mass-ratio inspiral. Comparing these waveforms to the vacuum scenario allows projections to the potential constrainability on the dark photon mass by space-borne gravitational wave antennas. The superradiant instability of a massive vector field on a Kerr background is calculated and the modification to the dynamics of an inspiralling solar mass-scale compact object is determined with approximations on the backreaction effect of the cloud on the compact object. The end result is the projection that the LISA mission should be able to constrain the dark photon mass using extreme mass ratio inspirals in the range $[1.8 \times 10^{-17}, 4.47 \times 10^{-16}]$ eV.

This article provides a detailed investigation into the motion of the surrounding particles around a polymer black hole in loop quantum gravity (LQG). Using effective potential, the critical bound orbits and innermost stable circular orbits (ISCO) are analyzed. The study finds that the radii and angular momentum of the critical bound orbits decrease with an increase in the parameter $A_\lambda$ which labels the LQG effects, while the energy and angular momentum of the ISCO also decreases with an increase in $A_\lambda$. Based on these findings, we then explore the periodic orbits of the polymer black hole in LQG using rational numbers composed of three integers. Our results show that the rational numbers increase with the energy of particles and decrease with the increase of angular momentum based on a classification scheme. Moreover, compared to a Schwarzschild black hole, the periodic orbits in a polymer black hole in LQG consistently have lower energy, providing a potential method for distinguishing a polymer black hole in LQG from a Schwarzschild black hole. Finally, we also examine the gravitational wave radiations of the periodic orbits of a test object which orbits a supermassive polymer black hole in LQG, which generates intricate GW waveforms that can aid in exhibiting the gravitational structure of the system.

This paper discusses the leading-order correction induced by cosmological perturbations on the average expansion rate of an expanding spacetime, containing one or many perfect fluids. The calculation is carried out up to the second order in the perturbations, and is kept as general as possible. In particular, no approximation such as a long-wavelength or a short-wavelength limit is invoked, and all three types of perturbations (scalar, vector, and tensor) are considered. First, the average value of the expansion rate is computed over a three-dimensional space-like surface where the total density of the fluids is constant. Then, a formula is derived relating that average value to the one over any other surface, on which a different scalar property of the fluids is constant. Moreover, the general formulas giving the correction to the average expansion rate are applied, in particular, to the case of a spacetime containing a single fluid with a constant equation of state. The sign and the effective equation of state of the corresponding back-reaction effect in the first Friedmann equation are examined.

From a classical analysis, we show that gravitational waves in a cosmological medium with equation of state $\omega=-1/3$ can follow a London-like equation, implying that some gravitational wave solutions present a decay for certain wavelengths. This scenario, corresponding to a cosmic string cosmology, induces an attenuation temporal scale on the gravitational wave propagation. We discuss on how these solutions impose a limit on the wavelength of the waves that can propagate, which depends on the type of spatial curvature and the energy density content of this type of cosmology.