Papers for Wednesday, Mar 30 2022

The SPT0311-58 system at z=6.900 is an extremely massive structure within the reionization epoch, and offers a chance to understand the formation of galaxies in an extreme peak in the primordial density field. We present 70mas Atacama Large Millimeter/submillimeter Array observations of the dust continuum and CII 158um emission in the central pair of galaxies and reach physical resolution ~100-350pc, among the most detailed views of any reionization-era system to date. The observations resolve the source into at least a dozen kiloparsec-size clumps. The global kinematics and high turbulent velocity dispersion within the galaxies present a striking contrast to recent claims of dynamically cold thin-disk kinematics in some dusty galaxies just 800Myr later at z~4. We speculate that both gravitational interactions and fragmentation from massive parent disks have likely played a role in the overall dynamics and formation of clumps in the system. Each clump individually is comparable in mass to other 6<z<8 galaxies identified in rest-UV/optical deep field surveys, but with star formation rates elevated by ~3-5x. Internally, the clumps themselves bear close resemblance to greatly scaled-up versions of virialized cloud-scale structures identified in low-redshift galaxies. Our observations are qualitatively similar to the chaotic and clumpy assembly within massive halos seen in simulations of high-redshift galaxies.

A complete census of dusty star-forming galaxies (DSFGs) at early epochs is necessary to constrain the obscured contribution to the cosmic star formation rate density (CSFRD), however DSFGs beyond $z \sim 4$ are both rare and hard to identify from photometric data alone due to degeneracies in submillimeter photometry with redshift. Here, we present a pilot study obtaining follow-up Atacama Large Millimeter Array (ALMA) $2\,$mm observations of a complete sample of 39 $850\,\rm\mu m$-bright dusty galaxies in the SSA22 field. Empirical modeling suggests $2\,$mm imaging of existing samples of DSFGs selected at $850\,\rm\mu m - 1\,$mm can quickly and easily isolate the "needle in a haystack" DSFGs that sit at $z>4$ or beyond. Combining archival submillimeter imaging with our measured ALMA $2\,$mm photometry ($1\sigma \sim 0.08\,$mJy$\,$beam$^{-1}$ rms), we characterize the galaxies' IR SEDs and use them to constrain redshifts. With available redshift constraints fit via the combination of six submillimeter bands, we identify 6/39 high-$z$ candidates each with $>50\%$ likelihood to sit at $z > 4$, and find a positive correlation between redshift and $2\,$mm flux density. Specifically, our models suggest the addition of $2\,$mm to a moderately constrained IR SED will improve the accuracy of a millimeter-derived redshift from $\Delta z/(1+z) = 0.3$ to $\Delta z/(1+z) = 0.2$. Our IR SED characterizations provide evidence for relatively high emissivity spectral indices ($\langle \beta \rangle = 2.4\pm0.3$) in the sample. We measure that especially bright ($S_{850\rm\mu m}>5.55\,$mJy) DSFGs contribute $\sim10$% to the cosmic-averaged CSFRD from $2<z<5$, confirming findings from previous work with similar samples.

Determining the sound speed $c_s$ in compact stars is an important open question with numerous implications on the behaviour of matter at large densities and hence on gravitational-wave emission from neutron stars. To this scope, we construct more than $10^7$ equations of state (EOSs) with continuous sound speed and build more than $10^8$ nonrotating stellar models consistent not only with nuclear theory and perturbative QCD, but also with astronomical observations. In this way, we find that EOSs with sub-conformal sound speeds, i.e. with $c^2_s < 1/3$ within the stars, are possible in principle but very unlikely in practice, being only $0.03\%$ of our sample. Hence, it is natural to expect that $c^2_s > 1/3$ somewhere in the stellar interior. Using our large sample, we obtain estimates at $95\%$ credibility of neutron-star radii for representative stars with $1.4$ and $2.0$ solar masses, $R_{1.4}=12.42^{+0.52}_{-0.99}\,{\rm km}$, $R_{2.0}=12.12^{+1.11}_{-1.23}\,{\rm km}$, and for the binary tidal deformability of the GW170817 event, $\tilde\Lambda_{1.186}=485^{+225}_{-211}$. Interestingly, our lower-bounds on the radii are in very good agreement with the prediction derived from very different arguments, namely, the threshold mass. Finally, we provide simple analytic expressions to determine the minimum and maximum values of $\tilde\Lambda$ as a function of the chirp mass.

We assemble a large comprehensive sample of 2534 z~2, 3, 4, 5, 6, 7, 8, and 9 galaxies lensed by the six clusters from the Hubble Frontier Fields (HFF) program. Making use of the availability of multiple independent magnification models for each of the HFF clusters and alternatively treating one of the models as the "truth," we show that the median magnification factors from the v4 parametric models are typically reliable to values of 30 to 50, and in one case to 100. Using the median magnification factor from the latest v4 models, we estimate the UV luminosities of the 2534 lensed z~2-9 galaxies, finding sources as faint as -12.4 mag at z~3 and -12.9 mag at z~7. We explicitly demonstrate the power of the surface density-magnification relations Sigma(z) vs. mu in the HFF clusters to constrain both distant galaxy properties and cluster lensing properties. Based on the Sigma(z) vs. mu relations, we show that the median magnification estimates from existing public models must be reliable predictors of the true magnification mu to mu<15 (95% confidence). We also use the observed Sigma(z) vs. mu relations to derive constraints on the evolution of the luminosity function faint-end slope from z~7 to z~2, showing that faint-end slope results can be consistent with blank-field studies if, and only if, the selection efficiency shows no strong dependence on the magnification factor mu. This can only be the case if very low luminosity galaxies are very small, being unresolved in deep lensing probes.

We study gas inflows onto supermassive black holes using hydrodynamics simulations of isolated galaxies and idealized galaxy mergers with an explicit, multiphase interstellar medium (ISM). Our simulations use the recently developed ISM and stellar evolution model called Stars and MUltiphase Gas in GaLaxiEs (SMUGGLE). We implement a novel super-Lagrangian refinement scheme that increases the gas mass resolution in the immediate neighborhood of the black holes (BHs) to accurately resolve gas accretion. We do not include black hole feedback in our simulations. We find that the complex and turbulent nature of the SMUGGLE ISM leads to highly variable BH accretion. BH growth in SMUGGLE converges at gas mass resolutions $\lesssim3\times10^3{\rm M_\odot}$. We show that the low resolution simulations combined with the super-Lagrangian refinement scheme are able to produce central gas dynamics and BH accretion rates very similar to that of the uniform high resolution simulations. We further explore BH fueling by simulating galaxy mergers. The interaction between the galaxies causes an inflow of gas towards the galactic centres and results in elevated and bursty star formation. The peak gas densities near the BHs increase by orders of magnitude resulting in enhanced accretion. Our results support the idea that galaxy mergers can trigger AGN activity, although the instantaneous accretion rate depends strongly on the local ISM. We also show that the level of merger-induced enhancement of BH fueling predicted by the SMUGGLE model is much smaller compared to the predictions by simulations using an effective equation of state model of the ISM.

Supersonic motions are common in molecular clouds. (Sub)sonic turbulence is usually detected toward dense cores and filaments. However, it remains unknown whether (sub)sonic motions at larger scales ($\gtrsim$1~pc) can be present in different environments or not. Located at a distance of about 110 pc, Ophiuchus North 1 (Oph N1) is one of the nearest molecular clouds that allows in-depth investigation of its turbulence properties by large-scale mapping observations of single-dish telescopes. We carried out the $^{12}$CO ($J=1-0$) and C$^{18}$O ($J=1-0$) imaging observations toward Oph N1 with the Purple Mountain Observatory 13.7 m telescope. The observations have an angular resolution of $\sim$55\arcsec (i.e., 0.03~pc). Most of the whole C$^{18}$O emitting regions have Mach numbers of $\lesssim$1, demonstrating the large-scale (sub)sonic turbulence across Oph N1. Based on the polarization measurements, we estimate the magnetic field strength of the plane-of-sky component to be $\gtrsim$9~$\mu$G. We infer that Oph N1 is globally sub-Alfv{\'e}nic, and is supported against gravity mainly by the magnetic field. The steep velocity structure function can be caused by the expansion of the Sh~2-27 H{\scriptsize II} region or the dissipative range of incompressible turbulence. Our observations reveal a surprising case of clouds characterised by widespread subsonic turbulence and steep size-linewidth relationship. This cloud is magnetized where ion-neutral friction should play an important role.

The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $1.05\,M_\odot$ and a $0.7\,M_\odot$ carbon-oxygen white dwarf with $0.03\,M_\odot$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $v>15000\,\mathrm{km\,s^{-1}}$ are indistinguishable. Light curves and spectra are very similar until $\sim 40$d after explosion and the ejecta are much more spherical than for violent merger models. The inner ejecta differ significantly which slows down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by about 20 per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model.

The chemical makeup of a star provides the fossil information of the environment where it formed. Under this premise, it should be possible to use chemical abundances to tag stars that formed within the same stellar association. This idea - known as chemical tagging - has not produced the expected results, especially within the thin disk where open stellar clusters have chemical patterns that are difficult to disentangle. The ultimate goal of this study is to probe the feasibility of chemical tagging within the thin disk population using high-quality data from a controlled sample of stars. We also aim at improving the existing techniques of chemical tagging and giving guidance on different strategies of clustering analysis in the elemental abundance space. Here we develop the first blind search of open clusters' members through clustering analysis in the elemental abundance space using the OPTICS algorithm applied to data from the Gaia-ESO survey. First, we evaluate different strategies of analysis, determining which ones are more performing. Second, we apply these methods to a data set including both field stars and open clusters attempting a blind recover of as many open clusters as possible. We show how specific strategies of data analysis can improve the final results. Specifically, we demonstrate that open clusters can be more efficaciously recovered with the Manhattan metric and on a space whose dimensions are carefully selected. Using these (and other) prescriptions we are able to recover open clusters hidden in our data set and find new members of these stellar associations. Our results indicate that there are chances of recovering open clusters' members via clustering analysis in the elemental abundance space. Presumably, the performances of chemical tagging will further increase with higher quality data and more sophisticated clustering algorithms.

We overview angle-based models to study planar symmetric central configurations of four bodies. We present models to determine the masses given the shape of the configuration and the shape-type of the configuration given the masses. We also describe a diagram-based method for the counting of the shape-types in the case of the concave configurations. As an application, we determine planar symmetric central configurations containing bodies of Earth's and Moon's masses.

We have analyzed the kinematics of OB2 stars with proper motions and parallaxes selected by Xu et al. from the Gaia EDR3 catalog. The relative parallax errors for all the stars in this sample do not exceed 10\%. Based on a sample of 9750 stars, the group velocity components $(U,V,W)_\odot=(7.21,7.46,8.52)\pm(0.13,0.20,0.10)$~km/s were obtained and the parameters of the angular velocity of rotation of the Galaxy: $\Omega_0 =29.712\pm0.062$~km/s/kpc, $\Omega^{'}_0=-4.014\pm0.018$~km/s/kpc$^{2}$ and $\Omega^{''}_0=0.674\pm0.009$~km/s/kpc$^{3}$. The circular velocity of rotation of the solar neighborhood around the center of the Galaxy is $V_0=240.7\pm3.0$~km/s for the assumed distance of the Sun to the galactic center $R_0=8.1\pm0.1$~kpc. It is shown that the influence of the systematic correction to the trigonometric parallaxes of the Gaia EDR3 catalog with the value $\Delta\pi=-0.040$~mas does not exceed the $\sim1\sigma$ level of the errors of the sought-for kinematic parameters of the model. Based on the proper motions of OB stars, the following variances of the residual velocities were found: $(\sigma_1,\sigma_2,\sigma_3)=(11.79,9.66,7.21)\pm(0.06,0.05,0.04)$~km/s. It is shown that the first axis of this ellipsoid slightly deviates from the direction to the center of the Galaxy $L_1=12.4\pm0.1^\circ$, and the third axis is oriented almost exactly to the north pole of the Galaxy, $B_3=87.7\pm0.1^\circ.$

We compute the linear order, general relativistic corrections to angular redshift fluctuations (ARF), a new cosmological observable built upon density-weighted two-dimensional (2D) maps of galaxy redshifts. We start with an existing approach for galaxy/source counts developed in the Newtonian gauge, and generalize it to ARF, modifying for this purpose a standard Boltzmann code. Our calculations allow us identifying the velocity terms as the leading corrections on large scales, emphasizing the sensitivity of ARF to peculiar, cosmological velocity fields. Just like for standard 2D clustering, the impact of gravitational lensing on ARF is dominant on small angular scales and for wide redshift shells, while the signatures associated to gravitational potentials are extremely small and hardly detectable. The ARF also present interesting correlation properties to anisotropies of the Cosmic Microwave Background (CMB): they are highly correlated to CMB lensing potential fluctuations, while also exhibiting a significant (S/N$\sim 4$-$5$) {\em anti-}correlation with the Integrated Sachs-Wolfe effect (ISW). This negative ARF$\times$ISW signal is quite complementary to the standard 2D clustering$\times$ISW correlation, since the former appears mostly at higher redshift ($z\sim 2$) than the latter ($z\lesssim 1)$, and the combination of the two observables significantly increases the $\chi^2$ statistics testing the null (no ISW) hypothesis. We conclude that ARF constitute a novel, alternative, and potentially powerful tool to constrain the nature of Dark Energy component that gives rise to the ISW.

Determining the frequency and duration of changing--look (CL) active galactic nuclei (AGNs) phenomena, where the optical broad emission lines appear or disappear, is crucial to understand the evolution of the accretion flow around supermassive black holes. We present a strategy to select new CL candidates starting from a spectroscopic type 2 AGNs sample and searching for current type 1 photometric variability. We use the publicly available Zwicky Transient Facility (ZTF) alert stream and the Automatic Learning for the Rapid Classification of Events (ALeRCE) light curve classifier to produce a list of CL candidates with a highly automated algorithm, resulting in 60 candidates. Visual inspection reduced the sample to 30. We performed new spectroscopic observations of six candidates of our clean sample, without further refinement, finding the appearance of clear broad Balmer lines in four of them and tentative evidence of type changes in the remaining two, which suggests a promising success rate of $\geq66$ per cent for this CL selection method.

We examine the physical conditions, environments, and statistical properties of intergalactic Ovi, Ovii and Oviii absorbers in the Simba cosmological hydrodynamic simulation suite. The goal is to understand the nature of these high ionisation absorbers, and test Simba's surprising prediction that $\sim 70\%$ of cosmic baryons at $z=0$ are in the Warm-Hot Intergalactic Medium (WHIM) driven by jet feedback from active galactic nuclei (AGN). By comparing a full-physics Simba run versus one with jets turned off, we find that jet feedback causes widespread heating that impacts the absorption morphology particularly of the higher ions. However, the distribution of the physical properties of detectable absorbers are not dramatically affected. Higher ionisation absorbers probe hotter gas as expected, but in Simba all ions arise at similar overdensities (typically $\delta\sim20-30$), similar environments (predominantly filaments), and similar nearest-halo distances (typically $\sim2-3r_{200c}$). Simba matches the observed Ovi column density distribution function (CDDF) fairly well, but under-predicts the CDDF preliminarily derived from two detected intergalactic Ovii absorbers. Predicted CDDFs are very similar at $z=1$ with or without jets, but show differences by $z=0$ particularly at the high-column end. Despite some discrepancies, Simba reproduces available observations as well as or better than other comparable simulations, suggesting that Simba's widespread jet heating cannot be ruled out by these data. These results offer hope that future X-ray and ultraviolet facilities could provide significant constraints on galactic feedback models from high-ionisation IGM metal absorbers.

With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as well as the reanalysis of the 4.5 $\micron$ Qatar-1b phase curve due to the similar equilibrium temperature. In total, five 4.5 $\micron$ phase curves and three 3.6 $\micron$ phase curves are analyzed here with a uniform approach. Using these new results, in combination with literature values for the entire population of published Spitzer phase curves of hot Jupiters, we present evidence for a linear trend of increasing hot spot offset with increasing orbital period, as well as observational evidence for two classes of planets in apparent redistribution vs. equilibrium temperature parameter space, and tentative evidence for a dependence of hot spot offset on planetary surface gravity in our $\sim$ 1300 K sample. We do not find trends in apparent heat redistribution with orbital period or gravity. Non-uniformity in literature Spitzer data analysis techniques precludes a definitive determination of the sources or lack of trends.

In the present paper, we report the photometric and spectroscopic observations obtained by the 1.88 m telescope at the Kottamia astronomical observatory of the pulsating star BL Cam. Fourier analysis of the light curves reveals that the fundamental mode has two harmonics. The O-C method is used to establish the period changes. So far, the analysis has been very successful in mapping the pulsation amplitude of the star across the instability strip. By using the formalism of Eddington and Plakidis (1929), we found significant results and strong indications of the evolutionary period change. A total of 55 new maximum light timings are reported. New values of (1/P) dP/dt are estimated using the O-C diagram based on all newly obtained times of maximum light combined with those taken from the literature, assuming the periods are decreasing and changing smoothly. To compute the effective temperature and surface gravity of the star, we performed model atmosphere analysis on its spectra. The physical parameters of the star are calculated and compared with the evolutionary models.

We present the first measurement of pulse scattering close to the eclipse region of PSR B1957+20, which is in a compact binary system with a low-mass star. We measured pulse scattering time-scales up to 0.2 ms close to the eclipse and showed that it scales with the dispersion measure (DM) excess roughly as $\tau\propto\Delta{\rm DM}^{2}$. Our observations provide the first evidence of strong scattering due to multi-path propagation effects in the eclipsing material. We show that Kolmogorov turbulence in the eclipsing material with an inner scale of $\sim100$ m and an outer scale of the size of the eclipse region can naturally explain the observation. Our results show that the eclipsing material in such systems can be highly turbulent and suggest that scattering is one of the main eclipsing mechanisms at around 1.4 GHz.

We report 18 dust continuum detections ($\geq 3.3\sigma$) at $\sim88{\rm \mu m}$ and $158{\rm \mu m}$ out of 49 ultraviolet(UV)-bright galaxies ($M_{\rm UV} < -21.3$ mag) at $z>6.5$, observed by the Cycle-7 ALMA Large Program, REBELS and its pilot programs. This has more than tripled the number of dust continuum detections known at $z>6.5$. Out of these 18 detections, 12 are reported for the first time as part of REBELS. In addition, 15 of the dust continuum detected galaxies also show a [CII]$_{\rm 158{\rm \mu m}}$ emission line, providing us with accurate redshifts. We anticipate more line emission detections from six targets (including three continuum detected targets) where observations are still ongoing. The dust continuum detected sources in our sample tend to have a redder UV spectral slope than the ones without a dust continuum detection. We estimate that all of the sources have an infrared (IR) luminosity ($L_{\rm IR}$) in a range of $3-8 \times 10^{11} L_\odot$, except for one with $L_{\rm IR} = 1.5^{+0.8}_{-0.5} \times 10^{12}\,L_{\odot}$. Their fraction of obscured star formation is significant at $\gtrsim 50\%$. Some of the dust continuum detected galaxies show spatial offsets ($\sim 0.5-1.5''$) between the rest-UV and far-IR emission peaks. These separations appear to have an increasing trend against an indicator that suggests spatially decoupled phases of obscured and unobscured star formation. REBELS offers the best available statistical constraints on obscured star formation in UV-bright, massive galaxies at $z > 6.5$.

Multiple pieces of evidence suggest that neutron stars receive large kicks when formed from the remnant of a collapsing star. However, the evidence for whether black holes (BH) receive natal kicks is less clear, reliant on weak constraints from the analysis of BH X-ray binaries and massive runaway and walkaway stars. Here we show for the first time that recent microlensing detections offer a new method for measuring the kicks BHs receive at birth. When a BH is identified through both photometric and astrometric microlensing and when the lensed star has a known distance and proper motion, the mass, distance and proper motion of the BH can be determined. We study the runaway velocities for components of eccentric binaries disrupted during a supernova, finding the peculiar velocity correlates strongly with the kick a BH received a birth, typically within 20\%, even when the natal kick is smaller than the orbital velocity. Therefore, by measuring the peculiar velocity of a BH or other compact object that formed from a binary that disrupted during core collapse, we are in effect measuring the natal kick that object received. We focus on MOA-2011-BLG-191/OGLE-2011-BLG-0462, an isolated, single BH detected by microlensing and consider a range of possible formation scenarios, including its formation from the disruption of a binary during a supernova event. We determine that MOA-2011-BLG-191/OGLE-2011-BLG-0462 has a Milky Way orbit consistent with a thick disk population, but if it were formed within the kinematic thin disk it received a natal kick $\lesssim$100 km s$^{-1}$.

A particle orbiting a misaligned eccentric orbit binary undergoes nodal precession either around the binary angular momentum vector (a circulating orbit) or around a stationary inclination (a librating orbit). In the absence of general relativity, the stationary inclination is inclined by 90 degrees to the binary angular momentum vector (aligned with the binary eccentricity vector) and does not depend on the particle semi-major axis. General relativity causes apsidal precession of the binary orbit. Close to the binary, the behaviour of the particle is not significantly affected, a librating particle precesses with the binary. However, we find that the stationary inclination and the minimum inclination required for libration both increase with the particle semi-major axis. There is a critical radius beyond which there are no librating orbits, only circulating orbits, and therefore there is a maximum orbital radius for a stationary polar orbiting body. The critical radius is within planet forming regions around binaries with semi-major axis <= 1 au. This has implications for the search for misaligned circumbinary planets and the radial extent of polar circumbinary disks.

Stellar active regions, including spots and faculae, can create radial velocity (RV) signals that interfere with the detection and mass measurements of low mass exoplanets. In doing so, these active regions affect each spectral line differently, but the origin of these differences is not fully understood. Here we explore how spectral line variability correlated with S-index (Ca H & K emission) is related to the atomic properties of each spectral line. Next we develop a simple analytic stellar atmosphere model that can account for the largest sources of line variability with S-index. Then we apply this model to HARPS spectra of {\alpha} Cen B to explain Fe I line depth changes in terms of a disk-averaged temperature difference between active and quiet regions on the visible hemisphere of the star. This work helps establish a physical basis for understanding how stellar activity manifests differently in each spectral line, and may help future work mitigating the impact of stellar activity on exoplanet RV surveys.

The nonthermal particle acceleration during magnetic reconnection remains a fundamental topic in several astrophysical phenomena, such as solar flares, pulsar wind, magnetars, etc, for more than half a century, and one of the unresolved questions is its efficiency. Recently, nonthermal particle acceleration mechanisms during reconnection have been extensively studied by particle-in-cell simulations, yet it is an intriguing enigma as to how the magnetic field energy is divided into thermally heated plasmas and nonthermal particles. Here we study both non-relativistic and relativistic magnetic reconnections using large-scale particle-in-cell simulation for a pair plasma, and indicate that the production of the nonthermal particle becomes efficient with increasing the plasma temperature. In the relativistic hot plasma case, we determine that the heated plasmas by reconnection can be approximated by a kappa distribution function with the kappa index of approximately 3 or less (equivalent to 2 or less for the power-law index), and the nonthermal energy density of reconnection is approximately over 95% of the total internal energy in the downstream exhaust.

Supernovae emit large fluxes of neutrinos which can be detected by detectors on Earth. Future tonne-scale detectors will be sensitive to several neutrino interaction channels, with thousands of events expected if a supernova emerges in the galaxy neighborhood. There is a limited number of tools to study the interaction rates of supernova neutrinos, although a plethora of available supernova models exists. EstrellaNueva is an open-source software to calculate expected rates of supernova neutrinos in detectors using target materials with typical compositions, and additional compositions can be easily added. This software considers the flavor transformation of neutrinos in the supernova through the adiabatic Mikheyev--Smirnov--Wolfenstein effect, and their interaction in detectors through several channels. Most of the interaction cross sections have been analytically implemented, such as neutrino-electron and neutrino-proton elastic scattering, inverse beta decay, and coherent elastic neutrino-nucleus scattering. This software provides a link between supernova simulations and the expected events in detectors by calculating fluences and event rates to ease any comparison between theory and observation. It provides a simple and standalone tool to explore many physics scenarios offering an option to add analytical cross sections and define any target material.

We present a high-cadence short term photometric and spectroscopic monitoring campaign of a type Ibn SN 2019wep, which is one of the rare SN Ibn after SNe 2010al and 2019uo to display signatures of flash ionization (\ion{He}{2}, \ion{C}{3}, \ion{N}{3}). We compare the decline rates and rise time of SN 2019wep with other SNe Ibn and fast transients. The post-peak decline in all bands (0.1 mag d$^{-1}$) are consistent with SNe Ibn but less than the fast transients. On the other hand, the $\Delta$m$_{15}$ values are slightly lower than the average values for SNe Ibn but consistent with the fast transients. The rise time is typically shorter than SNe Ibn but longer than fast transients. SN 2019wep lies at the fainter end of SNe Ibn but possesses an average luminosity amongst the fast transients sample. The peculiar color evolution places it between SNe Ib and the most extreme SNe Ibn. The bolometric light curve modelling shows resemblance with SN 2019uo with ejecta masses consistent with SNe Ib. SN 2019wep belongs to the "P cygni" sub-class of SNe Ibn and shows faster evolution in line velocities as compared to the "emission" sub-class. The post-maximum spectra show close resemblance with ASASSN-15ed hinting it to be of SN Ib nature. The low \ion{He}{1} CSM velocities and residual H$\alpha$ further justifies it and gives evidence of an intermittent progenitor between WR and LBV star.

Recently, a tight correlation between the dynamical radial acceleration and the baryonic radial acceleration in galaxies - the radial acceleration relation - has been discovered. This has been claimed as an indirect support of the modified gravity theories. However, whether the radial acceleration relation could also be found in galaxy clusters is controversial. In this article, we derive and present an analytic radial acceleration relation for the central region of galaxy clusters. We examine the data of some large galaxy clusters and we find that the resulting radial acceleration relation has a very large scatter. Moreover, although the radial acceleration relation for galaxy clusters shows some agreement with the one discovered in galaxies for a certain range of baryonic radial acceleration, their functional forms are somewhat different from each other. This suggests that the radial acceleration relation may not be a universal relation in general.

We estimate the amount of negative feedback energy injected into the ISM of the host galaxy of 3C273, a prototypical radio loud quasar. We obtained 93, 233 and 343 GHz continuum images with the Atacama Large Millimeter/Sub-millimeter Array (ALMA). After self calibration and point source subtraction, we reach an image dynamic range of $\sim 85000$ at 93\ GHz, $\sim 39000$ at 233\ GHz and $\sim 2500$ at 343\ GHz. These are currently the highest image dynamic range obtained using ALMA. We detect spatially extended millimeter emission associated with the host galaxy, cospatial with the Extended Emission Line Region (EELR) observed in the optical. The millimeter spectral energy distribution and comparison with centimeter data show that the extended emission cannot be explained by dust thermal emission, synchrotron or thermal bremsstrahlung arising from massive star formation. We interpret the extended millimeter emission as thermal bremsstrahlung from gas directly ionized by the central source. The extended flux indicates that at least $\sim 7\%$ of the bolometric flux of the nuclear source was used to ionize atomic hydrogen in the host galaxy. The ionized gas is estimated to be as massive as $10^{10}$ to $10^{11}\ \mathrm{M_\odot}$, but the molecular gas fraction with respect to the stellar mass is consistent with other ellipticals, suggesting that direct ionization ISM by the QSO may not be sufficient to suppress star formation, or we are witnessing a short timescale before negative feedback becomes observable. The discovery of a radio counterpart to EELRs provides a new pathway to studying the QSO-host ISM interaction.

Motivated by previous findings that the magnitude gap between certain satellite galaxy and the central galaxy can be used to improve the estimation of halo mass, we carry out a systematic study of the information content of different member galaxies in the modelling of the host halo mass using a machine learning approach. We employ data from the hydrodynamical simulation IllustrisTNG and train a Random Forest (RF) algorithm to predict a halo mass from the stellar masses of its member galaxies. Exhaustive feature selection is adopted to disentangle the importances of different galaxy members. We confirm that an additional satellite does improve the halo mass estimation compared to that estimated by the central alone. However, the magnitude of this improvement does not differ significantly using different satellite galaxies. When three galaxies are used in the halo mass prediction, the best combination is always that of the central galaxy with the most massive satellite and the smallest satellite. Furthermore, among the top 7 galaxies, the combination of a central galaxy and two or three satellite galaxies gives a near-optimal estimation of halo mass, and further addition of galaxies does not raise the precision of the prediction. We demonstrate that these dependences can be understood from the shape variation of the conditional satellite distribution, with different member galaxies accounting for distinct halo-dependent features in different parts of the cumulative stellar mass function.

Efforts to discover and characterize habitable zone planets have primarily focused on Sun-like stars and M dwarfs. K stars, however, provide an appealing compromise between these two alternatives that has been relatively unexplored. Understanding the ultraviolet (UV) environment around such stars is critical to our understanding of their planets, as the UV can drastically alter the photochemistry of a planet's atmosphere. Here we present near-UV and far-UV \textit{Hubble Space Telescope}'s Cosmic Origins Spectrograph observations of 39 K stars at three distinct ages: 40 Myr, 650 Myr, and $\approx$5 Gyr. We find that the K star (0.6 -- 0.8 M$_{\odot}$) UV flux remains constant beyond 650 Myr before falling off by an order of magnitude by field age. This is distinct from early M stars (0.3 -- 0.6 M$_{\odot}$), which begin to decline after only a few hundred Myr. However, the rotation-UV activity relation for K stars is nearly identical to that of early M stars. These results may be a consequence of the spin-down stalling effect recently reported for K dwarfs, in which the spin-down of K stars halts for over a Gyr when their rotation periods reach $\approx$10 d, rather than the continuous spin down that G stars experience. These results imply that exoplanets orbiting K dwarfs may experience a stronger UV environment than thought, weakening the case for K stars as hosts of potential "super-habitable" planets.

We investigate the importances of various dynamical features in predicting the dynamical state (DS) of galaxy clusters, based on the Random Forest (RF) machine learning approach. We use a large sample of galaxy clusters from the Three Hundred Project of hydrodynamical zoomed-in simulations, and construct dynamical features from the raw data as well as from the corresponding mock maps in the optical, X-ray, and Sunyaev-Zel'dovich (SZ) channels. Instead of relying on the impurity based feature importance of the RF algorithm, we directly use the out-of-bag (OOB) scores to evaluate the importances of individual features and different feature combinations. Among all the features studied, we find the virial ratio, $\eta$, to be the most important single feature. The features calculated directly from the simulations and in 3-dimensions carry more information on the DS than those constructed from the mock maps. Compared with the features based on X-ray or SZ maps, features related to the centroid positions are more important. Despite the large number of investigated features, a combination of up to three features of different types can already saturate the score of the prediction. Lastly, we show that the most sensitive feature $\eta$ is strongly correlated with the well-known half-mass bias in dynamical modelling. Without a selection in DS, cluster halos have an asymmetric distribution in $\eta$, corresponding to an overall positive half-mass bias. Our work provides a quantitative reference for selecting the best features to discriminate the DS of galaxy clusters in both simulations and observations.

Chromo-natural inflation (CNI) is an inflationary model where an axion coupled with SU$(2)$ gauge fields acts as the inflaton. In CNI, the gauge fields have nonzero vacuum expectation values (VEVs), which results in the enhancement of gravitational waves (GWs). The original CNI is ruled out by the Planck observations due to the overproduction of GWs. In this work, we consider an inflationary model where the gauge fields acquire nonzero VEVs after the CMB modes exit the horizon. Moreover, we add to the model another field that dominates the universe and drives inflation after the axion starts to oscillate and the gauge field VEVs vanish. By performing numerical simulations, we find a parameter space where the enhanced GWs do not violate the CMB constraints and can be detected by the future GWs observations such as LISA, BBO, and ET.

We input solar wind parameters responsible for the main phases of 15 great geomagnetic storms (GGSs: $\Delta$SYM-H$\le-$200 nT) into the empirical formulae created by \cite{Burton1975}(hereafter Burton equation), and by \cite{OBrien2000}(hereafter OM equation) to evaluate whether \textbf{two equations} can correctly estimate the intensities of GGSs. The results show that the intensities of most GGSs estimated by OM equation are much smaller than the observed intensities. The RMS error between the intensities estimated by OM equation and the observed intensities is \textbf{203} nT, implying that the estimated storm intensity deviates significantly from the observed one. The RMS error between the intensities estimated by Burton equation and the observed intensities is 130.8 nT. The relative error caused by Burton equation for the storms with intensities $\Delta$SYM-H$<$-400 nT is larger than 27\%, implying that the absolute error will be large for the storms with $\Delta$SYM-H$<$-400 nT. The results indicate that the two equations cannot work effectively in the estimation of GGSs. On the contrary, the intensity of a GGS estimated by the empirical formula created by \cite{WangCB2003} can always be very close to the observed one if we select the right weight for solar wind dynamic pressure, proving that solar wind dynamic pressure is an important factor for GGS intensity, but it is overlooked in the ring current injection terms of Burton equation or OM equation. This is the reason why the two equations cannot work effectively in the estimation of GGSs.

We revisit the role of longitudinal waves in driving the solar wind. We study how the the $p$-mode-like vertical oscillation on the photosphere affects the properties of solar winds under the framework of Alfv\'en-wave-driven winds. We perform a series of one-dimensional magnetohydrodynamical numerical simulations from the photosphere to beyond several tens of solar radii. We find that the mass-loss rate drastically increases with the longitudinal wave amplitude at the photosphere up to $\sim 4$ times, in contrast to the classical understanding that the acoustic wave hardly affects the energetics of the solar wind. The addition of the longitudinal fluctuation induces the longitudinal-to-transverse wave mode conversion in the chromosphere, which results in the enhanced Alfv\'enic Poynting flux in the corona. Consequently, the coronal heating is promoted to give higher coronal density by the chromospheric evaporation, leading to the increased mass-loss rate. This study clearly shows the importance of the longnitudinal oscillation in the photosphere and the mode conversion in the chromosphere in determining the basic properties of the wind from solar-like stars.

We present the baseline conceptual design of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope. CUBES will provide unprecedented sensitivity for spectroscopy on a 8 - 10 m class telescope in the ground ultraviolet (UV), spanning a bandwidth of > 100 nm that starts at 300 nm, the shortest wavelength accessible from the ground. The design has been optimized for end-to-end efficiency and provides a spectral resolving power of R > 20000, that will unlock a broad range of new topics across solar system, Galactic and extraglactic astronomy. The design also features a second, lower-resolution (R \sim 7000) mode and has the option of a fiberlink to the UVES instrument for simultaneous observations at longer wavelengths. Here we present the optical, mechanical and software design of the various subsystems of the instrument after the Phase A study of the project. We discuss the expected performances for the layout choices and highlight some of the performance trade-offs considered to best meet the instrument top-level requirements. We also introduce the model-based system engineering approach used to organize and manage the project activities and interfaces, in the context that it is increasingly necessary to integrate such tools in the development of complex astronomical projects.

Understanding structures and evolutions of the magnetic fields and plasma in multiple layers on the Sun is very important. A force-free magnetic field which is an accurate approximation of the solar corona due to the low plasma $\beta$ has been widely studied and used to model the coronal magnetic structure. While the force-freeness assumption is well satisfied in the solar corona, the lower atmosphere is not force-free given the high plasma $\beta$. Therefore, a magnetohydrostatic (MHS) equilibrium which takes into account plasma forces, such as pressure gradient and gravitational force, is considered to be more appropriate to describe the lower atmosphere. This paper reviews both analytical and numerical extrapolation methods based on the MHS assumption for calculating the magnetic fields and plasma in the solar atmosphere from measured magnetograms.

Pulsations and binarity are both common features of massive stars. The study of pulsating massive stars in eclipsing binary systems hold great potential for constraining stellar structure and evolution theory. However, prior to the all-sky Transiting Exoplanet Survey Satellite (TESS) mission, few such systems had been discovered or studied in detail. We have inspected the TESS light curves of a large number of eclipsing binaries known to contain high-mass stars, and compiled a list of 18 objects which show intrinsic variability. The light curves were modelled both to determine the physical properties of the systems, and to remove the effects of binarity in order to leave residual light curves suitable for asteroseismic analysis. Precise mass and radius measurements were obtained for delta Cir, CC Cas, SZ Cam, V436 Per and V539 Ara. We searched the residual light curves for pulsation signatures and, within our sample of 18 objects, we find six definite and eight possible cases of beta Cephei pulsation, seven cases of stochastic low-frequency (SLF) variability, and eight instances of possible slowly pulsating B (SPB) star pulsation. The large number of pulsating eclipsing systems we have identified makes asteroseismology of high-mass stars in eclipsing binaries a feasible avenue to constrain the interior physics of a large sample of massive stars for the first time.

A comprehensive study of Nitric oxide (NO) chemistry in the Martian upper atmosphere is restricted due to the lack of requisite measurements. NO is an abundant form of odd nitrogen species in the Martian lower atmosphere and its density depends on several photochemical processes. We have developed a photochemical model to study the NO density in the dayside of Martian upper atmosphere by accounting for various production and loss mechanisms. By utilizing the Neutral Gas and Ion Mass Spectrometer (NGIMS) on-board Mars Atmosphere and Volatile Evolution (MAVEN) mission measured neutral and ion densities during deep dip 8 and 9 campaigns, we modelled NO number density in the Martian sunlit upper atmosphere for the altitudes between 120 and 200 km. The modelled NO densities are employed to calculate NO (1,0) gamma band emission intensity profiles in the dayside upper atmosphere of Mars. The calculated NO density and its gamma band intensity profiles are found to be consistent with Imaging Ultraviolet Spectrograph (IUVS) onboard MAVEN observations and also with other modelling studies. We found that the local CO2 and N2 density variations can lead to a change in NO density and consequently its dayglow intensity by a factor of 2 to 5. Since NO is a trace constituent and also its dayglow emissions are strongly obscured by CO Cameron band emissions, we suggest that the derivation of NO number density based on our approach can constrain its abundance in the dayside upper atmosphere of Mars. More observations of (1-0) gamma band emission along with modelling will help to study the global distribution of NO in the Martian atmosphere.

Extragalactic water megamaser emissions at 22 GHz have been playing vital roles in astrophysics. The limited detection rate of these masers has been motivating researchers to find clues that can help characterize them. The physical environments masers formed in is still ambiguous, accordingly, statistical studies have been thoroughly used to resolve these favorable environments. In this work, we go through the most essential parameter of Active Galactic Nuclei (AGN), namely, the mass of the central supermassive black hole (MBH) of the maser host galaxy. We study the correlation between maser luminosity (LH2O) and MBH. The regression line of the relation is also calculated. Additionally, sub-samples of megamasers (MMs), kilomasers (KMs), and disc masers are studied. Our results show a very significant LH2O-MBH correlation for our 68 galaxy sample. Unlike the results of previous works, dividing the sample into MMs and KMs gives no privilege to MM galaxies. In opposite to expectation, KMs have weak and low significant LH2O - MBH correlation, while MMs show no correlation. The positive correlation in KMs can be explained by the role of AGN therein, while the diversity of MMs types, with some of which are not strongly related to AGN, may explain the correlation missing. The 28 disc maser sample, where tight correlation is expected, surprisingly shows a very weak and low significant LH2O - MBH correlation. Future VLBI studies will eventually lead to a certain classification of a good number of maser galaxies, which is essential to clearly establish the LH2O-MBH relation.

We present the simulation tools developed to aid the design phase of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope (VLT), exploring aspects of the system design and evaluating the performance for different design configurations. CUBES aims to be the 'ultimate' ultraviolet (UV) instrument at the European Southern Observatory (ESO) in terms of throughput, with the goal to cover the bluest part of the spectrum accessible from the ground (300 nm to 400 nm) with the highest possible efficiency. Here we introduce the End-to-End (E2E) and the Exposure Time Calculator (ETC) tools. The E2E simulator has been developed with different versions to meet the needs of different users, including a version that can be accessed for use by the broader scientific community using a Jupyter notebook. The E2E tool was used by the system team to help define the Phase A baseline design of the instrument, as well as in scientific evaluation of a possible low-resolution mode. The ETC is a web-based tool through which the science community are able to test a range of science cases for CUBES, demonstrating its potential to push the limiting magnitude for the detection of specific UV-features, such as abundance estimates of beryllium in main sequence stars.

One of the most intriguing properties of the GD-1 stellar stream is the existence of three gaps. If these gaps were formed by close encounters with dark matter subhalos, the GD-1 stream opens an exciting window through which we can see the size, mass, and velocity distributions of the dark matter subhalos in the Milky Way. However, in order to use the GD-1 stream as a probe of the dark matter substructure, we need to disprove that these gaps are not due to the perturbations from baryonic components of the Milky Way. Here we ran a large number of test particle simulations to investigate the probability that each of the known globular clusters (GCs) can form a GD-1-like gap, by using the kinematical data of the GD-1 stream and GCs from Gaia EDR3 and by fully taking into account the observational uncertainty. We found that the probability that all of the three gaps were formed by GCs is as low as $1.2\times10^{-5}$ and the expected number of gaps formed by GCs is only $0.057$ in our fiducial model. Our result highly disfavors a scenario in which GCs form the gaps. Given that other baryonic perturbers (e.g., giant molecular clouds) are even less likely to form a gap in the retrograde-moving GD-1 stream, we conclude that at least one of the gaps in the GD-1 stream was formed by dark matter subhalos if the gaps were formed by flyby perturbations.

To evaluate the role of radio activity in galaxy evolution, we designed a large archival CO survey of radio galaxies (RGs) to determine their molecular gas masses at different epochs. We used a sample of 120 RGs representative of the NVSS 1.4GHz survey, when flux limited at 0.4Jy. Of those, 66 galaxies belonged to the ALMA Radio-source Catalogue (ARC) of calibrators and had spectral window tunings around CO (1-0), (2-1), (3-2), or (4-3). We reduced their ALMA data, determined their H2 mass contents, and combined the results with similar results for the remaining 54 galaxies from the literature. We found that, while at all epochs the majority of RGs have undetectable reservoirs, there is a rapid increase in the H2 mass content of the CO-detected RGs with z. At 1<z<2.5, 1/4 of the RGs have at least as much molecular gas as simulations would indicate for a typical halo mass of that epoch. These galaxies plausibly have 'normal' or even starbursty hosts. Taking into account the completeness correction of the sample, we created the corresponding H2 mass functions at 0.005<z<0.3 and 1<z<2.5. The local mass function reveals that the number density of low-z RGs with detectable molecular gas reservoirs is only a little lower (a factor of ~4) than that of type 1 and 2 AGN in simulations. At 1<z<2.5, there is a significant decrease in the number density of high-z RGs due to the rarity of bright radio galaxies. An estimate for the missing faint RGs would, nonetheless, bring populations close again. Finally, we find that the volume density of molecular gas locked up in the brightest 1/5000-1/7000 RGs is similar in the examined z bins. This result likely indicates that the inflow rate on one hand and the star-formation depletion rate plus the jet-driven expulsion rate on the other hand counteract each other in the most luminous RGs of each epoch.

We propose to use an elongated rectangular waveguide near its cutoff frequency for axionic dark matter searches. The detector's large surface area allows for significant signal power, while its narrow transverse dimension and tapered-waveguide coupling suppress parasitic modes. The proposed system can fit inside a solenoid magnet and is sensitive to the QCD-axion in the axion mass $40-400\,\mu$eV. We describe the theoretical principles of the new design, present simulation results, and discuss the implementation.

We train neural models to represent both the optimal policy (i.e. the optimal thrust direction) and the value function (i.e. the time of flight) for a time optimal, constant acceleration low-thrust rendezvous. In both cases we develop and make use of the data augmentation technique we call backward generation of optimal examples. We are thus able to produce and work with large dataset and to fully exploit the benefit of employing a deep learning framework. We achieve, in all cases, accuracies resulting in successful rendezvous (simulated following the learned policy) and time of flight predictions (using the learned value function). We find that residuals as small as a few m/s, thus well within the possibility of a spacecraft navigation $\Delta V$ budget, are achievable for the velocity at rendezvous. We also find that, on average, the absolute error to predict the optimal time of flight to rendezvous from any orbit in the asteroid belt to an Earth-like orbit is small (less than 4\%) and thus also of interest for practical uses, for example, during preliminary mission design phases.

Solar extreme ultraviolet (EUV) waves are spectacular propagating disturbances with EUV enhancements in annular shapes in the solar corona. These EUV waves carry critical information about the coronal magnetised plasma that can shed light on the elusive physical parameters (e.g. the magnetic field strength) by global solar coronal magneto-seismology. EUV waves are closely associated with a wide range of solar atmospheric eruptions, from violent flares and coronal mass ejections (CMEs) to less energetic plasma jets or mini-filament eruptions. However, the physical nature and driving mechanism of EUV waves is still controversial. Here, we report the unique discovery of twin EUV waves (TEWs) that were formed in a single eruption with observations from two different perspectives. In all earlier studies, a single eruption was associated at most with a single EUV wave. The newly found TEWs urge to re-visit our theoretical understanding about the underlying formation mechanism(s) of coronal EUV waves. Two distinct scenarios of TEWs were found. In the first scenario, the two waves were separately associated with a filament eruption and a precursor jet, while in another scenario the two waves were successively associated with a filament eruption. Hence, we label these distinguished scenarios as "fraternal TEWs" and "identical TEWs", respectively. Further, we also suggest that impulsive lateral expansions of two distinct groups of coronal loops are critical to the formation of TEWs in a single eruption.

Gamma-Ray Bursts (GRBs), observed up to $z=9.4$, can be employed as standardized candles, extending the distance ladder beyond Supernovae Type Ia (SNe Ia, $z=2.26$). We standardize GRBs using the 3D fundamental plane relation among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity. Combining SNe Ia and GRBs, we constrain the matter content $\Omega_{\text{M}}= 0.299 \pm 0.009$ assuming a flat $\Lambda$CDM cosmology with and without correcting for selection biases and redshift evolution. Using a similar 3D correlation in optical, we find that our optical GRB sample is as efficacious in the determination of $\Omega_{\text{M}}$ as the X-ray sample. We trimmed our GRB samples to achieve tighter planes to be used to simulate additional GRBs. We determined how many GRBs are needed to be used as standalone probes to achieve a comparable precision on $\Omega_{\text{M}}$ to the one obtained by SNe Ia only. {\bf We achieve the error measurements derived using SNe Ia in (Conley et al. 2011) with 142 GRBs in optical considering the errorbars on the variables halved and in Betoule et al. (2014) with 284 optical GRBs, by considering that the error bars on the variables halved. Using a doubled sample (obtained by future machine learning approaches that allow both a lightcurve reconstruction and the estimates of GRBs for which the GRB is unknown) compared to our current one, with errorbars halved we will reach the same precision of (Scolnic et al. 2018) with 390 GRBs which will be uses as efficient standalone probes as SNe Ia by 2054.

The search for pulsars in a sample of pulsar candidates found based on a multi-year survey conducted with low (6 channels; sampling 0.1s) time-frequency resolution on declinations -9^o < \delta < +42^o was carried out. A Large Phased Array (LPA) transit telescope operating at 111 MHz in the 2.5 MHz band was used. Search, analysis and evidence of pulsar detection were carried out using a visualization program of summed up power spectra obtained from the survey data with high 32 channels; sampling 12.5ms) time-frequency resolution. 11 new pulsars with periods P0 = 0.41-3.75 s and dispersion measure DM=15-154 pc/cm^3 have been discovered. In total, in the survey with a low time-frequency resolution for the period 2016-2021 in a blind search, 208 pulsars were found, of them 42 new and 166 known pulsars. It is shown that in the search on the data with high time-frequency resolution accumulated over a time interval of 7 years, pulsars with a flux density of 0.1 - 0.2 mJy at the frequency of 111 MHz can be detected. When searching for pulsars with regular (periodic) emission at declinations +21^o < \delta < +42^o, all pulsars located outside the galactic plane, having P0 \ge 0.5 s, DM \le 100 pc/cm^3 and the flux density S \ge 0.5 mJy, can be detected

The tangential YORP effect is a radiation pressure torque produced by asymmetric thermal emission by structures on the asteroid surface. As such structures, previous works considered boulders of different shapes lying on the surface of the asteroid. We study the tangential YORP produced by the rough interface of the asteroid's regolith. We create an approximate analytic theory of heat conduction on a slightly non-flat sinusoidal surface. We analyze the published data on the small-scale shape of the asteroid (162173) Ryugu and estimate its tangential YORP due to the surface roughness. We derive an analytic formula that expresses the TYORP of a sinusoidal surface in terms of its geometric and thermal properties. TYORP is maximal at the thermal parameter of the order of unity and for the shape irregularities of the order of the thermal wavelength. Application of this equation to Ryugu predicts TYORP, which is 5-70 times greater than its normal YORP effect. The contribution of the small-scale regolith roughness to the YORP effect of the asteroid can be comparable to the normal YORP and the tangential YORP produced by boulders. The same theory can describe roughness on the asteroid boulders, thus adding a new term to the previously considered TYORP created by boulders.

The proposed CUBES spectrograph for ESO's Very Large Telescope will be an exceptionally powerful instrument for the study of comets. The gas coma of a comet contains a large number of emission features in the near-UV range covered by CUBES (305-400 nm), which are diagnostic of the composition of the ices in its nucleus and the chemistry in the coma. Production rates and relative ratios between different species reveal how much ice is present and inform models of the conditions in the early solar system. In particular, CUBES will lead to advances in detection of water from very faint comets, revealing how much ice may be hidden in the main asteroid belt, and in measuring isotopic and molecular composition ratios in a much wider range of comets than currently possible, provide constraints on their formation temperatures. CUBES will also be sensitive to emissions from gaseous metals (e.g., FeI and NiI), which have recently been identified in comets and offer an entirely new area of investigation to understand these enigmatic objects.

The science case on studies of accretion and outflows in low-mass ($<$1.5 $M_{\odot}$) young stellar objects (YSOs) with the new CUBES instrument is presented. We show the need for a high-sensitivity, near-ultraviolet (NUV) spectrograph like CUBES, with a resolving power at least four times that of X-Shooter and combined with UVES via a fibrelink for simultaneous observations. Simulations with the CUBES exposure time calculator and the end-to-end software show that a significant gain in signal-to -noise can be achieved compared to current instruments, for both the spectral continuum and emission lines, including for relatively embedded YSOs. Our simulations also show that the low-resolution mode of CUBES will be able to observe much fainter YSOs (V $\sim$22 mag) in the NUV than we can today, allowing us extend studies to YSOs with background-limited magnitudes. The performance of CUBES in terms of sensitivity in the NUV will provide important new insights into the evolution of circumstellar disks, by studying the accretion, jets/winds and photo-evaporation processes, down to the low-mass brown dwarf regime. CUBES will also open-up new science as it will be able to observe targets that are several magnitudes fainter than those reachable with current instruments, facilitating studies of YSOs at distances of $\sim$ kpc scale. This means a step-change in the field of low-mass star formation, as it will be possible to expand the science case from relatively local star-forming regions to a large swathe of distances within the Milky Way.

Absorption lines from molecular hydrogen ($\rm H_2$) in the spectra of background sources are a powerful probe of the physical conditions in intervening cold neutral medium. At high redshift, $z>2$, $\rm H_2$ lines are conveniently shifted in the optical domain, allowing the use of ground-based telescopes to perform high-resolution spectroscopy, which is essential for a proper analysis of the cold gas. We describe recent observational progress, based on the development of efficient pre-selection techniques in low-resolution spectroscopic surveys such as the Sloan Digital Sky Survey (SDSS). The next generation of spectrographs with high blue-throughput, such as CUBES, will certainly significantly boost the efficiency and outcome of follow-up observations. In this paper, we discuss high priority science cases for CUBES, building on recent $\rm H_2$ observations at high-z: probing the physical conditions in the cold phase of regular galaxies and outflowing gas from active galactic nucleus.

We study the signature of primordial non-Gaussianity imprinted on the power spectrum of the 21-cm line differential brightness temperature during dark ages. Employing the perturbative treatment of gravitational clustering, we quantitatively estimate the effects of the non-Gaussian and one-loop corrections on the 21-cm power spectrum. The potential impact of the use of the 21-cm power spectrum for the constraint on local-type primordial non-Gaussianity is investigated based on the Fisher matrix analysis. Our results show that the 21-cm power spectrum for an array with a baseline of several tens of kilometers can constrain the primordial non-Gaussianity to a level severer than that from cosmic microwave background measurements and its constraining power is stronger than that of the 21-cm bispectrum, while in the ultimate situation the 21-cm bispectrum eventually becomes more powerful.

Globular clusters host multiple stellar populations that display star-to-star variation of light elements that are affected by hot hydrogen burning (e.g., He, C, N, O). Several scenarios have been suggested to explain these variations. Most involve multiple star formation episodes, where later generations are born from material contaminated by the nucleosynthetic products of the previous stellar generation(s). One difficulty in the modelling of such scenarios is knowing the extent to which processed and pristine material are mixed. In this context, beryllium abundances measured in turn-off stars of different generations can provide new information. Beryllium originates from cosmic-ray spallation and can only be destroyed inside stars. Beryllium abundances can thus directly measure the degree of pollution of the material that formed stars in globular clusters. Turn-off stars in globular clusters are however faint and such studies are beyond the capabilities of current instrumentation. In this work, we show the progress that the CUBES spectrograph will bring to this area. Our simulations indicate that CUBES will enable the detection of variations of about 0.6 dex in the Be abundances between stars from different generations, in several nearby globular clusters with turn-off magnitude down to $V$ = 18 mag.

We extend our previous work on the enhancement of the curvature spectrum during inflation to the two-field case. We identify the slow-roll parameter $\eta$ as the quantity that can trigger the rapid growth of perturbations. Its two components, $\eta_\parallel$ along the background trajectory and $\eta_\perp$ perpendicular to it, remain small during most of the evolution, apart from short intervals during which they take large, positive or negative, values. The typical reason for the appearance of strong features in $\eta_\parallel$ is sharp steps or inflection points in the inflaton potential, while $\eta_\perp$ grows large during sharp turns in field space. We focus on the additive effect of several features leading to the resonant growth of the curvature spectrum. Three or four features in the evolution of $\eta$ are sufficient in order to induce an enhancement of the power spectrum by six or seven orders of magnitude, which can lead to the significant production of primordial black holes and stochastic gravitational waves. A big part of our study focuses on understanding the evolution of the perturbations and the resulting spectra through analytic means. The presence of multiple features in the background evolution points to a more complex inflationary paradigm, which is also more natural in the multi-field case. The critical examination of this possibility is within the reach of experiment.

We study the impact of the COVID-19 pandemic on astronomy using public records of astronomical publications. We show that COVID-19 has had both positive and negative impacts on research in astronomy. We find that the overall output of the field, measured by the yearly paper count, has increased. This is mainly driven by boosted individual productivity seen across most countries, possibly the result of cultural and technological changes in the scientific community during COVID. However, a decreasing number of incoming new researchers is seen in most of the countries we studied, indicating larger barriers for new researchers to enter the field or for junior researchers to complete their first project during COVID. Unfortunately, the overall improvement in productivity seen in the field is not equally shared by female astronomers. By fraction, fewer papers are written by women and fewer women are among incoming new researchers in most countries. Even though female astronomers also became more productive during COVID, the level of improvement is smaller than for men. Pre-COVID, female astronomers in the Netherlands, Australia, Switzerland were equally as or even more productive than their male colleagues. During COVID, no single country's female astronomers were able to be more productive than their male colleagues on average.

Diffusive shock acceleration at collisionless shocks remains the most likely process for accelerating particles in a variety of astrophysical sources. While the standard prediction for strong shocks is that the spectrum of accelerated particles is universal, $f(p)\propto p^{-4}$, numerous phenomena affect this simple conclusion. In general, the non-linear dynamical reaction of accelerated particles leads to a concave spectrum, steeper than $p^{-4}$ at momenta below a few tens of GeV/c and harder than the standard prediction at high energies. However, the non-linear effects become important in the presence of magnetic field amplification, which in turn leads to higher values of the maximum momentum $p_{max}$. It was recently discovered that the self-generated perturbations that enhance particle scattering, when advected downstream, move in the same direction as the background plasma, so that the effective compression factor at the shock decreases and the spectrum becomes steeper. We investigate the implications of the excitation of the non-resonant streaming instability on these spectral deformations, the dependence of the spectral steepening on the shock velocity and the role played by the injection momentum.

T Aurigae is an eclipsing old nova which exploded in 1891. At a Gaia EDR3 distance of 815-871 pc, it is a relatively nearby old nova. Through ultraviolet spectral modeling and using the new precise Gaia distance, we find that the HST/STIS spectrum of T Aurigae is consistent with an accretion disk with a mass transfer rate $\dot{M}$ of the order of $10^{-8}M_{\odot}$/yr, for a white dwarf mass of $M_{\rm wd} \approx 0.7 \pm 0.2 M_{\odot}$, an inclination of $i \sim 60^{\circ}$, and a Gaia distance of of $840_{-25}^{+31}$~pc. The sharp absorption lines of metals cannot form in the disk and are likely forming in material above the disk (e.g. due stream disk overflow), in circumbinary material, and/or in material associated with the ejected shell from the 1891 nova explosion. The saturated hydrogen Ly$\alpha$ absorption feature is attributed to a large interstellar medium hydrogen column density of the order of $10^{21}$cm$^{-2}$ towards T Aur, as corroborated by the value of its reddening $E(B-V)=0.42 \pm 0.08$.

CUBES is the Cassegrain U-Band Efficient Spectrograph, a high-efficiency instrument operating in the UV spectral range between 300nm and 400nm with a resolution not less than 20000. CUBES is to be installed at a Cassegrain focus of the Very Large Telescope of the European Southern Observatory. The paper briefly reviews various types of devices used as dispersing elements in astronomical spectrographs to achieve high resolution, before identifying binary transmission gratings produced by microlithography as the best candidate technology for the CUBES instrument. We describe the lithographic fabrication technology in general, two different design considerations to achieve the required high-resolution transmission grating, its prototyping by a direct-write lithographic fabrication technology, and the characterization of the achieved optical performance. An outlook to the realization of the grating for the final instrument, taking the most recent developments of lithographic writing capabilities into consideration is given.

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Ly{\alpha}, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

Context: As a ubiquitous phenomenon, large-amplitude longitudinal filament oscillations usually decay in 1--4 periods. Recently, we observed a decayless case of such oscillations in the corona. Aims: We try to understand the physical process that maintains the decayless oscillation of the filament. Methods: Multi-wavelength imaging observations and magnetograms are collected to study the dynamics of the filament oscillation and its associated phenomena. To explain the decayless oscillations, we also perform one-dimensional hydrodynamic numerical simulations using the MPI-AMRVAC code. Results: In observations, the filament oscillates decaylessly with a period of $36.4 \pm 0.3$ min for almost 4 hours before eruption. During oscillations, four quasi-periodic jets emanate from a magnetic cancellation site near the filament. The time interval between neighboring jets is $\sim 68.9 \pm 1.0$ min. Numerical simulations constrained by the observations reproduced the decayless longitudinal oscillations. However, it is surprising to find that the period of the decayless oscillations is not consistent with the pendulum model. Conclusions: We propose that the decayless longitudinal oscillations of the filament are maintained by quasi-periodic jets, which is verified by the hydrodynamic simulations. More importantly, it is found that, when driven by quasi-periodic jets, the period of the filament longitudinal oscillations depends also on the driving period of the jets, not simply the pendulum period. With a parameter survey in simulations, we derived a formula, by which one can derive the pendulum oscillation period using the observed period of decayless filament oscillations and the driving periods of jets.

As the Universe evolves, it develops a web of filamentary structure of matter. This cosmic web is filled with gas, with the most diffuse gas lying in the intergalactic regions. At low redshift, the gas is predominantly warm-hot, and one of its best tracers is X-ray absorption in sightlines to background quasars. In this Chapter, we present the theoretical background for the formation of the warm-hot intergalactic medium (WHIM) and present the physical properties of the WHIM from cosmological hydro-dynamical simulations. We discuss the feasibility of detecting the WHIM with X-ray absorption lines, with high-resolution and high signal-to-noise spectra. We present detailed discussion of observing techniques, including the WHIM ionization balance, observable lines, the curve of growth, and the diagnostics using the X-ray lines. We present the current efforts of detecting the WHIM with gratings on-board Chandra and XMM-Newton observatories. We discuss the criticality of WHIM detections reported in literature, where robust detections are likely from the circumgalactic medium of intervening galaxies, or intra-group medium, rather than truly diffuse gas in the intergalactic medium. Secure detections of the most diffuse gas in the low redshift large scale structure may have to await next generation of X-ray telescopes. We end our Chapter with the discussion of future missions carrying dispersive and non-dispersive spectrometers. We present figure-of-merit parameters for line detectibility as well as for the number of WHIM systems that can be detected with future missions. These will define our ability to account for the missing low-redshift baryons and to understand the evolution of the Universe over half of its life.

The most recent detections of LIGO/Virgo and NICER have placed strong constraints on neutron stars' properties. In this work, we study neutron stars modeling them as hybrid stars, compact objects with a quark matter core surrounded by layers of hadronic matter. In addition, we consider the presence of strong magnetic fields, motivated by magnetars detection, like the recently discovered Swift J1818.0-1607. We incorporate the effects of the anomalous magnetic moment of the constituent particles into the equation of state of dense matter and analyze the implications of different hadron-quark phase transitions on the dynamic stability of these compact objects. This study shows that the constraints on the mass, radius, and tidal deformability, imposed by observations of massive pulsars and gravitational waves can be satisfied within the framework of our model.

We derive a simple analytical model for the propagation of HSSs from the Sun to Earth and thereby show how the area of coronal holes and the size of their boundary regions affect the HSS velocity, temperature, and density near Earth. We presume that velocity, temperature, and density profiles form across the HSS cross section close to the Sun and that these spatial profiles translate into corresponding temporal profiles in a given radial direction due to the solar rotation. These temporal distributions drive the stream interface to the preceding slow solar wind plasma and disperse with distance from the Sun. The HSS properties at 1 AU are then given by all HSS plasma parcels launched from the Sun that did not run into the stream interface at Earth distance. We show that the velocity plateau region of HSSs as seen at 1 AU, if apparent, originates from the center region of the HSS close to the Sun, whereas the velocity tail at 1 AU originates from the trailing boundary region. The peak velocity of HSSs at Earth further depends on the longitudinal width of the HSS close to the Sun. The temperature and density of HSS plasma parcels at Earth depend on their radial expansion from the Sun to Earth. The radial expansion is determined by the velocity gradient across the HSS boundary region close to the Sun and gives the velocity-temperature and density-temperature relationships at Earth their specific shape. When considering a large number of HSSs, the presumed correlation between the HSS velocities and temperatures close to the Sun degrades only slightly up to 1 AU, but the correlation between the velocities and densities is strongly disrupted up to 1 AU due to the radial expansion. Finally, we show how the number of particles of the piled-up slow solar wind in the stream interaction region depends on the velocities and densities of the HSS and preceding slow solar wind plasma.

We investigate the feasibility of robust abundances for selected neutron-capture elements (Ge, Bi, Hf, U) from near-UV spectroscopy with the CUBES instrument now in development for the Very Large Telescope. We use the CUBES end-to-end simulator to synthesise observations of the Ge I 3039 {\AA} and Hf II 3400 and 3719 {\AA} lines in a very metal-poor star, using the well-studied star CS 31082-001 as a template. From simulated 4 hr exposures, we recover estimated abundances to $\pm$0.1 dex for Ge for U $\sim$ 14.25 mag., and for Hf for U = 18 mag. These performances neatly highlight the powerful gain of CUBES for near-UV observations of targets that are two-to-three magnitudes fainter than the existing observations of CS 31082-001 (U = 12.5 mag.). We also investigate the weak Bi I 3025 {\AA} and U II 3860 {\AA} lines (for U $\sim$ 14.25 and 16mag., respectively), finding that simulated 4hr exposures should provide upper limits to these observationally challenging lines.

The first measurements of the 21 cm brightness temperature power spectrum from the epoch of reionization will be very likely achieved in the near future by radio interferometric array experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the precursors of the Square Kilometre Array (SKA). Standard MCMC analyses use an explicit likelihood approximation to infer the reionization and astrophysical parameters from the 21 cm power spectrum. In this paper, we present a new Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). Realistic effects including thermal noise and foreground avoidance are also applied to the mock observations from the HERA and SKA. We demonstrate that this method recovers accurate posterior distributions for the reionization parameters, and outperforms the standard MCMC analysis in terms of the location and size of credible parameter regions. With the minutes-level processing time once the network is trained, this technique is a promising approach for the scientific interpretation of future 21 cm power spectrum observation data. Our code 21cmDELFI-PS is publicly available at this link.

The Tibet AS$\gamma$ and LHAASO collaborations recently reported the observation of a $\gamma$-ray diffuse emission with energy up to the PeV from the Galactic plane.} {We discuss the relevance of non-uniform cosmic-ray transport scenarios and the implications of these results for cosmic-ray physics.} {We use the {\tt DRAGON} and {\tt HERMES} codes to build high-resolution maps and spectral distributions of that emission for several representative models under the condition that they reproduce a wide set of local cosmic-ray data up to 100 PeV.} {We show that the energy spectra measured by Tibet AS$\gamma$, LHAASO, ARGO-YBJ and Fermi-LAT in several regions of interest in the sky can all be consistently described in terms of the emission arising by the Galactic cosmic-ray "sea". We also show that all our models are compatible with IceTop $\gamma$-ray upper limits.} {Our results favor transport models characterized by spatial-dependent diffusion although some degeneracy remains between the choice of the transport scenario and that of the cosmic-ray spectral shape above 10 TeV. We discuss the role of forthcoming measurements in resolving that ambiguity.

Large observational surveys, like those that will be conducted at the Vera C. Rubin Observatory, are expected to discover up to one million new asteroids in the first year of operation. This will more than double the database of known asteroids in a very short time. New methods and techniques will be needed to handle the large influx of data. Here, we tested some of these new methods by studying the population of asteroids on stable orbits inside the ${\nu}_6$ secular resonance. This resonance is one of the strongest mechanisms for destabilizing the orbits of main-belt bodies and producing Near-Earth Asteroids (NEAs). Yet, stable orbital configurations where the asteroid pericenter is either aligned or anti-aligned with that of Saturn exist inside the resonance. The population of stable ${\nu}_6$ resonators, first discovered in the early 2010s, is now the largest population of asteroids in stable orbits inside a secular resonance. Here we obtained the largest sample of asteroids' proper elements for numbered and multi-opposition objects most likely to be affected by the resonance. Artificial Neural Networks (ANN) were then used to identify the images of resonant angles of asteroids in stable orbits, more than doubling their number. Clustering methods and the use of machine learning algorithms allowed the identification of the known asteroid families crossed by the ${\nu}_6$ resonance, the Tina, Euphrosyne, and Svea clusters, and of two entirely new groups: the Tiffanykapler and the 138605 QW177 families.

We present the first results of our survey of asteroid polarization-phase curves in the near-infrared J and H bands using the WIRC+Pol instrument on the Palomar 200-inch telescope. We confirm through observations of standard stars that WIRC+Pol can reach the 0.1% precision needed for asteroid phase curve characterization, and show that C-complex asteroids could act as an alternate calibration source as they show less wavelength variation than stellar polarized standards. Initial polarization-phase curve results for S-complex asteroids show a shift in behavior as a function of wavelength from visible to near-infrared bands, extending previously observed trends. Full near-infrared polarization-phase curve characterization of individual asteroids will provide a unique constraint on surface composition of these objects by probing the wavelength dependence of albedo and index of refraction of the surface material.

Primordial black holes (PBHs) formed in the early Universe are sources of neutrinos emitted via Hawking radiation. Such astrophysical neutrinos could be detected at Earth and constraints on the abundance of comet-mass PBHs could be derived from the null observation of this neutrino flux. Here, we consider non-rotating PBHs and improve constraints using Super-Kamiokande neutrino data, as well as we perform forecasts for next-generation neutrino (Hyper-Kamiokande, JUNO, DUNE) and dark matter (DARWIN, ARGO) detectors, which we compare. For PBHs less massive than $\sim \textrm{few} \times 10^{14}$ g, PBHs would have already evaporated by now, whereas more massive PBHs would still be present and would constitute a fraction of the dark matter of the Universe. We consider monochromatic and extended (log-normal) mass distributions, and a PBH mass range spanning from $10^{12}$ g to $\sim 10^{16}$ g. Finally, we also compare our results with previous ones in the literature.

This work analyzes the phase correlation of the three lunar cycles and the Saros/Exeligmos cycle, after the study of the chapter About Exeligmos in Introduction to the Phenomena by Geminus. Geminus, refers that each Exeligmos cycle began on very specific and rare dates, when the Moon positioned at the points of the three lunar cycles beginning: New moon at Apogee and at the Node. The extremely large duration of the Annular Solar eclipse occurred on December 22 178BC (Saros series 58), marks the start of the Prominent Saros Cycle Apokatastasis. The next day, 23 December 178BC, the Winter Solstice started. During these two neighboring dates, the celebration of the religious festival of Isia started in Egypt and the Hellenistic Greece. After the analysis of the Mechanism's Parapegma events specific position, 22/23 December 178BC is an ideal, functional and representative initial date, in order to calibrate the initial position of the Mechanism's pointers.

We investigate whether true physical observables associated with the measurements of large scale structure in the universe are frame-independent. In particular, we study if cosmological observables such as the galaxy number counts and weak lensing observables are invariant under the disformal transformations. In a previous work, it was shown that this frame-invariance holds true for the case of conformal transformations. In this work, we find that although the cosmological number counts remain invariant under the disformal transformations, convergence and cosmic shear associated with the weak lensing potential are generally not invariant. Since the lightcone structure does not remain causal under disformal transformations, photon geodesics do not remain null anymore, and as a result, weak lensing observables are indeed affected. We also briefly comment on the disformal invariance of other cosmological observables.

Several pulsar-timing array (PTA) collaborations are finding tantalising hints for a stochastic gravitational wave background signal in the nano-Hertz regime. So far, though, no convincing evidence for the expected Hellings-Downs quadrupolar correlations has been found. While this issue might get fixed at the light of more accurate, forthcoming data, it is important to keep an eye open on different possibilities, and explore scenarios able to produce different types of PTA angular correlations. We point out that a stationary non-Gaussian component to the gravitational wave background can modulate the 2-point PTA overlap reduction function, adding contributions that can help in fitting the angular distribution of PTA data. We discuss possible sources for such non-Gaussian signal in terms of cosmological processes occurring after inflation ends, and we investigate further tests for this idea.

Dark matter and inflation are two key elements to understand the origin of cosmic structures in modern cosmology, and yet their exact physical models remain largely uncertain. The Weyl scaling invariant theory of gravity may provide a feasible scheme to solve these two puzzles jointly, which contains a massive gauge boson playing the role of dark matter candidate, and allows the quadratic scalar curvature term, namely $R^2$, to realize a viable inflationary mechanism in agreement with current observations. We ponder on the production of dark matters in the Weyl $R^2$ model, including the contribution from the non-perturbative production due to the quantum fluctuations from inflationary vacuum and perturbative ones from scattering. We demonstrate that there are generally three parameter ranges for viable dark matter production: (1) If the reheating temperature is larger than $10^4~\mathrm{GeV}$, the Weyl gauge boson as dark matter can be produced abundantly with mass larger than the inflation scale $\sim 10^{13}~\mathrm{GeV}$. (2) Small mass region with $2\times10^{-11}~\mathrm{GeV}$ for a higher reheating temperature. (3) Annihilation channel becomes important in the case of higher reheating temperature, which enables the Weyl gauge boson with mass up to $4\times10^{16}~\mathrm{GeV}$ to be produced through freeze-in.

Stochastic gravitational wave background from the early Universe has a cut-off frequency close to 100 MHz, due to the horizon of the inflationary phase. To detect gravitational waves at such frequencies, resonant electromagnetic cavities are very suitable. In this work, we study the frequency sensitivity of such detectors, and show how we could use them to probe this cut-off frequency and also the energy density per frequency of this stochastic background. This paper paves the way for further experimental studies to probe the most ancient relic of the Universe.