Papers for Friday, Jul 21 2023

BINA-3 has been the third workshop of this series involving scientists from India and Belgium aimed at fostering future joint research in the view of cutting-edge observatories and advances in theory. BINA-3 was held at the Graphic Era Hill University, 22-24 March 2023 at Bhimtal (near Nainital), Uttarakhand, India. A major event was the inauguration of the International Liquid-Mirror Telescope (ILMT), the first liquid mirror telescope devoted exclusively to astronomy. BINA-3 provided impressive highlights encompassing topics of both general astrophysics and solar physics. Research results and future projects have been featured through invited and contributed talks, and poster presentations.

We present the Lyman-$a$ (Lya) Luminosity Function (LF) at $2.05<z<3.75$, estimated from a sample of 67 Lya-emitter (LAE) candidates in the J-PAS Pathfinder surveys: miniJPAS and J-NEP. These two surveys cover a total effective area of $\sim 1.14$ deg$^2$ with 54 Narrow Band (NB) filters across the optical range, with typical limiting magnitudes of $\sim 23$. This set of NBs allows to probe Lya emission in a wide and continuous range of redshifts. We develop a method for detecting Lya emission for the estimation of the Lya LF using the whole J-PAS filter set. We test this method by applying it to the miniJPAS and J-NEP data. In order to compute the corrections needed to estimate the Lya LF and to test the performance of the candidates selection method, we build mock catalogs. These include representative populations of Lya Emitters at $1.9<z<4.5$ as well as their expected contaminants, namely low-$z$ galaxies and $z<2$ QSOs. We show that our method is able to provide the Lya LF at the intermediate-bright range of luminosity ($\rm 10^{43.5} erg\,s^{-1} \lesssim L_{Lya} \lesssim 10^{44.5} erg\,s^{-1}$). The photometric information provided by these surveys suggests that our samples are dominated by bright, Lya-emitting Active Galactic Nuclei. At $L_{{\rm Ly}a}<10^{44.5}$ erg\,s$^{-1}$, we fit our Lya LF to a power-law with slope $A=0.70\pm0.25$. We also fit a Schechter function to our data, obtaining: Log$(\Phi^* / \text{Mpc$^{-3}$})=-6.30^{+0.48}_{-0.70}$, Log$(L^*/ \rm erg\,s^{-1})=44.85^{+0.50}_{-0.32}$, $a=-1.65^{+0.29}_{-0.27}$. Overall, our results confirm the presence of an AGN component at the bright-end of the Lya LF. In particular, we find no significant contribution of star-forming LAEs to the Lya LF at Log$(L_{\rm Lya}$ / erg\,s$^{-1}$)>43.5. This work serves as a proof-of-concept for the results that can be obtained with the upcoming data releases of the J-PAS survey.

It has been reported that the application of convolutional neural-network techniques to infer the Dark Matter distribution in the local IGM has revealed how it follows the hierarchical distribution of galaxies in the locality, rather than exhibiting homogeneity. This result makes it natural to consider the possibility that, on scales at least as big as $20 \sim 30\,Mpc$, the distribution of all material comprising the local IGM is hierarchically distributed. Given this possibility, any model of galaxy formation must then involve a process in which all of the hierarchically distributed material $M_0$ within a sphere $R_0$ coalesces about a unique center so that hierarchical symmetry is broken on the scale $(M_0,R_0)$. In the particular case of the hierarchical distribution being quasi-fractal $D \approx 2$ in the local cosmos then, for circular velocity $V_0$ on $R_0$, the scaling relation $V_0^4 \sim M_0$ emerges automatically when the condition that such a galaxy formation process must be gravitationally stable in the Newtonian sense is applied. In other words, subject to the caveat that the analysis applies to a highly idealized model, the Baryonic Tully-Fisher Relationship (BTFR) is shown to arise as a consequence of Newtonian gravitation acting in a hierarchical Universe. We discuss the ramifications of this result, which are significant and non-trivial.

We investigate the ionised gas morphology, excitation properties, and kinematics in 19 nearby star-forming galaxies from the PHANGS-MUSE survey. We directly compare the kinetic energy of expanding superbubbles and the turbulent motions in the interstellar medium with the mechanical energy deposited by massive stars in the form of winds and supernovae, with the aim to answer whether the stellar feedback is responsible for the observed turbulent motions and to quantify the fraction of mechanical energy retained in the superbubbles. Based on the distribution of the flux and velocity dispersion in the H$\alpha$ line, we select 1484 regions of locally elevated velocity dispersion ($\sigma$(H$\alpha$)>45 km/s), including at least 171 expanding superbubbles. We analyse these regions and relate their properties to those of the young stellar associations and star clusters identified in PHANGS-HST data. We find a good correlation between the kinetic energy of the ionised gas and the total mechanical energy input from supernovae and stellar winds from the stellar associations, with a typical efficiency of 10-20%. The contribution of mechanical energy by the supernovae alone is not sufficient to explain the measured kinetic energy of the ionised gas, which implies that pre-supernova feedback in the form of radiation/thermal pressure and winds is necessary. We find that the gas kinetic energy decreases with metallicity for our sample covering Z=0.5-1.0 Zsun, reflecting the lower impact of stellar feedback. For the sample of superbubbles, we find that about 40% of the young stellar associations are preferentially located in their rims. We also find a slightly higher (by ~15%) fraction of the youngest (1-2.5 Myr) stellar associations in the rims of the superbubbles than in the centres, and the opposite for older associations, which implies possible propagation or triggering of star formation.

In galactic nuclei, stellar densities are so high that stars can physically collide with each other. In this work we focus on the collision of red giants and in particular on the formation of non-thermal processes through collisions and their properties. We analytically address these points by evaluating head-on collisions but also take into account scenarios with a deviation from the radial orbit, which we treat in a perturbative fashion. The collisions produce internal shocks with supersonic Mach numbers. Almost immediately, jet-like structures with important Lorentz factors form. The debris from the collision produces another shock wave which, when interacting with the interstellar medium of a galactic nucleus, leads to particle acceleration. We estimate the background flux in X- and gamma rays created by the background of these collisions by deriving the spectral index within a radius of 100 Mpc and find that they are high. Additionally, we make an estimate of the neutrino production and find about $10^{11}$ neutrinos per square meter per second for a collision at 100 Mpc from Earth. Also, we derive that there is a non-negligible chance to ignite fusion during the collision, due to the squeezing of the material. We investigate the possibility that the degenerate cores collide with each other, leading to a high afterglow luminosity, and find that it is non-negligible, although this should be addressed with dedicated numerical simulations. Colliding red giants in galactic nuclei trigger a plethora of high-energy phenomena, and have a particular gravitational wave emission associated, as shown by us, so that their detection will allow us to rule out alternatives.

In a previous work we analysed the white-light coronal brightness as a function of elongation and time from Wide-Field Imager (WISPR) observations on board the Parker Solar Probe (PSP) mission when PSP reached a minimum heliocentric distance of ~ 28 Rs. We found 4-5 transient outflows per day over a narrow wedge in the PSP orbital plane, which is close to the solar equatorial plane. However, the elongation versus time map (J-map) analysis supplied only lower limits on the number of released density structures due to the small spatial-scales of the transient outflows and line-of-sight integration effects. In this work we place constraints on the properties of slow solar wind transient mass release from the entire solar equatorial plane. We simulated the release and propagation of transient density structures in the solar equatorial plane for four scenarios: (1) periodic release in time and longitude with random speeds; (2) corotating release in longitude, periodic release in time with random speeds; (3) random release in longitude, periodic release in time and speed; and (4) random release in longitude, time, and speed. The simulations were used in the construction of synthetic J-maps, which are similar to the observed J-map. The four considered scenarios have similar ranges (35-45 for the minimum values and 96-127 for the maximum values) of released density structures per day from the solar equatorial plane and consequently from the streamer belt, given its proximity to the solar equatorial plane during the WISPR observation. Our results also predict that density structures with sizes in the range 2-8 Rs, covering 1-20 % of the perihelion could have been detectable by PSP in situ observations during that interval.

GRB 220627A is a rare burst with two distinct gamma-ray emission episodes separated by almost 1000 s that triggered the Fermi Gamma-ray Burst Monitor twice. High-energy GeV emission was detected by the Fermi Large Area Telescope coincident with the first emission episode but not the second. The discovery of the optical afterglow with MeerLICHT led to MUSE observations which secured the burst redshift to z=3.08, making this the most distant ultra-long gamma-ray burst (GRB) detected to date. The progenitors of some ultra-long GRBs have been suggested in the literature to be different to those of normal long GRBs. Our aim is to determine whether the afterglow and host properties of GRB 220627A agree with this interpretation. We performed empirical and theoretical modelling of the afterglow data within the external forward shock framework, and determined the metallicity of the GRB environment through modelling the absorption lines in the MUSE spectrum. Our optical data show evidence for a jet break in the light curve at ~1.2 days, while our theoretical modelling shows a preference for a homogeneous circumburst medium. Our forward shock parameters are typical for the wider GRB population, and we find that the environment of the burst is characterised by a sub-solar metallicity. Our observations and modelling of GRB 220627A do not suggest that a different progenitor compared to the progenitor of normal long GRBs is required. We find that more observations of ultra-long GRBs are needed to determine if they form a separate population with distinct prompt and afterglow features, and possibly distinct progenitors.

The nature of two recently discovered radio emitters with unusually long periods of 18min (GLEAM-X J1627-52) and 21min (GPM J1839-10) is highly debated. Their bright radio emission resembles that of radio magnetars, but their long periodicities and lack of detection at other wavelengths challenge the neutron-star interpretation. In contrast, long rotational periods are common in white dwarfs but, although predicted, dipolar radio emission from isolated magnetic white dwarfs has never been unambiguously observed. In this work, we investigate these long-period objects as potential isolated neutron-star or white-dwarf dipolar radio emitters and find that both scenarios pose significant challenges to our understanding of radio emission via pair production in dipolar magnetospheres. We also perform population-synthesis simulations based on dipolar spin-down in both pictures, assuming different initial-period distributions, masses, radii, beaming fractions, and magnetic-field prescriptions, to assess their impact on the ultra-long pulsar population. In the neutron-star scenario, we cannot reproduce the large number of expected ultra-long period pulsars under any physically motivated (or even extreme) assumptions. Thus, if GLEAM-X J1627-52 and GPM J1839-10 are confirmed as neutron-star pulsars (even if they are magnetars), this would necessarily call for a significant revision of our understanding of birth parameters at the population level. On the other hand, in the white-dwarf scenario, no mechanism can explain the production of such a bright coherent radio emission in isolated magnetic white dwarf systems (binaries with low mass companions are still viable), although we can easily accommodate a large population of long-period radio emitters.

The exact evolution of elements in the universe, from primordial hydrogen and helium to heavier elements like gold and platinum produced via the r-process, is still under scrutiny. The supernova deaths of the very first stars to form in the universe led to the enrichment of their local environments with new metals, and can leave behind neutron stars as remnants. These remnants can end up in binary systems with other neutron stars, and eventually merge, allowing for the r-process to occur. In this work, we study the scenario where a single neutron star merger (NSM) enriches a halo early in its evolution to understand the impact on the second generation of stars and their metal abundances. We perform a suite of high resolution cosmological zoom-in simulations using Enzo where we have implemented a new NSM model varying the explosion energy and the delay time. In general, a NSM leads to a significant r-process enhancement in the second generation of stars. A high explosion energy leads to almost all enhanced r-process stars being highly enhanced, while a lower explosion energy leads to a higher mass fraction of stars being r-process enhanced, but not as many being highly enhanced. When a NSM has a short delay time, there is a higher mass fraction of stars being r-process enhanced, but a smaller mass fraction being highly enhanced compared to longer delay times. This work represents a stepping stone towards understanding how NSMs impact their environments and metal abundances of descendant generations of stars.

Convolution Neural Networks trained for the task of lens finding with similar architecture and training data as is commonly found in the literature are biased classifiers. An understanding of the selection function of lens finding neural networks will be key to fully realising the potential of the large samples of strong gravitational lens systems that will be found in upcoming wide-field surveys. We use three training datasets, representative of those used to train galaxy-galaxy and galaxy-quasar lens finding neural networks. The networks preferentially select systems with larger Einstein radii and larger sources with more concentrated source-light distributions. Increasing the detection significance threshold to 12$\sigma$ from 8$\sigma$ results in 50 per cent of the selected strong lens systems having Einstein radii $\theta_\mathrm{E}$ $\ge$ 1.04 arcsec from $\theta_\mathrm{E}$ $\ge$ 0.879 arcsec, source radii $R_S$ $\ge$ 0.194 arcsec from $R_S$ $\ge$ 0.178 arcsec and source S\'ersic indices $n_{\mathrm{Sc}}^{\mathrm{S}}$ $\ge$ 2.62 from $n_{\mathrm{Sc}}^{\mathrm{S}}$ $\ge$ 2.55. The model trained to find lensed quasars shows a stronger preference for higher lens ellipticities than those trained to find lensed galaxies. The selection function is independent of the slope of the power-law of the mass profiles, hence measurements of this quantity will be unaffected. The lens finder selection function reinforces that of the lensing cross-section, and thus we expect our findings to be a general result for all galaxy-galaxy and galaxy-quasar lens finding neural networks.

Identifying the anisotropies in a cosmologically sourced stochastic gravitational wave background (SGWB) would be of significance in shedding light on the nature of primordial inhomogeneities. For example, if SGWB carries isocurvature fluctuations, it would provide evidence for a multi-field inflationary origin of these inhomogeneities. However, this is challenging in practice due to finite detector sensitivity and also the presence of the astrophysical foregrounds that can compete with the cosmological signal. In this work, we explore the prospects for measuring cosmological SGWB anisotropies in the presence of an astrophysical counterpart and detector noise. To illustrate the main idea, we perform a Fisher analysis using a well-motivated cosmological SGWB template corresponding to a first order phase transition, and an astrophysical SGWB template corresponding to extra-galactic binary mergers, and compute the uncertainty with which various parameters characterizing the isotropic and anisotropic components can be extracted. We also discuss some subtleties and caveats involving shot noise in the astrophysical foreground. Overall, we show that upcoming experiments, e.g., LISA, Taiji, Einstein Telescope, Cosmic Explorer, and BBO, can all be effective in discovering plausible anisotropic cosmological SGWBs.

We report the discovery of Mothra, an extremely magnified monster star, likely a binary system of two supergiant stars, in one of the strongly lensed galaxies behind the galaxy cluster MACS0416. The star is in a galaxy with spectroscopic redshift $z=2.091$ in a portion of the galaxy that is parsecs away from the cluster caustic. The binary star is observed only on the side of the critical curve with negative parity but has been detectable for at least eight years, implying the presence of a small lensing perturber. Microlenses alone cannot explain the earlier observations of this object made with the Hubble Space Telescope. A larger perturber with a mass of at least $10^4$\,\Msun\ offers a more satisfactory explanation. Based on the lack of perturbation on other nearby sources in the same arc, the maximum mass of the perturber is $M< 2.5\times10^6$\,\Msun, making it the smallest substructure constrained by lensing above redshift 0.3. The existence of this millilens is fully consistent with the expectations from the standard cold dark matter model. On the other hand, the existence of such small substructure in a cluster environment has implications for other dark matter models. In particular, warm dark matter models with particle masses below 8.7\,keV are excluded by our observations. Similarly, axion dark matter models are consistent with the observations only if the axion mass is in the range $0.5\times10^{-22}\, {\rm eV} < m_a < 5\times10^{-22}\, {\rm eV}$.

In core-collapse supernovae, the neutrino density is so large that neutrino flavor instabilities, leading to flavor conversion, can be triggered by the forward scattering of neutrinos among each other, if a crossing between the angular distributions of electron neutrino and antineutrinos exists (fast instability) or in the presence of perturbations induced by the neutrino vacuum frequency (slow instability). Recently, it has been advanced the conjecture that neutrino collisions with the medium could be another mean to kickstart flavor change (collisional instability). We rely on a spherically symmetric core-collapse supernova model with mass $18.6\ M_\odot$, compute the neutrino angular distributions solving the kinetic equations and investigate the occurrence of flavor instabilities at different post-bounce times, ranging from the accretion phase to the early cooling phase. We find that fast and slow flavor instabilities largely dominate over the collisional ones in the decoupling region for all post-bounce times. While more work is needed to assess the relevance of collisional instabilities in neutrino-dense environments, our findings suggest that neutrino collisions with matter affect the flavor evolution in the decoupling region, but are not responsible for triggering flavor conversion.

Aims: We aim to compute the binary fraction of "classical" dwarf spheroidal galaxies (dSphs) that are satellites of the Milky Way (MW). This value can offer insights into the binary fraction in environments that are less dense and more metal-poor than our own galaxy. Additionally, knowledge of the binary fraction in dwarf galaxies is important with respect to avoiding overestimations of their dark matter content, inferred from stellar kinematics. Methods: We refined an existing method from the literature, placing an emphasis on providing robust uncertainties on the value of the binary fraction. We applied this modified method to a VLT/FLAMES dataset for Sculptor, specifically acquired for the purpose of velocity monitoring of individual stars, as well as to literature datasets for other six MW "classical" dSphs. In all cases, the targeted stars were mainly red giant branch stars (RGBs), with expected masses of around 0.8 M$_{\odot}$. The VLT/FLAMES dataset offers the most precise binary fractions compared to literature datasets, due to its time baseline of 12 years, along with at least nine repeated observations for each star. Results: We found that the binary fraction of Sculptor is 0.55$^{+0.17}_{-0.19}$. We find that it is important to take into account the Roche lobe overflow for constraining the period distribution of binary stars. In contrast to what has recently been proposed in the literature, our analysis indicates that there is no evidence to support varying the properties of the binary stellar population or their deviations from those established for the solar neighborhood, based on the sample of MW dSphs analyzed here.

Strong gravitational lensing is a powerful tool to provide constraints on galaxy mass distributions and cosmological parameters, such as the Hubble constant, $H_0$. Nevertheless, inference of such parameters from images of lensing systems is not trivial as parameter degeneracies can limit the precision in the measured lens mass and cosmological results. External information on the mass of the lens, in the form of kinematic measurements, is needed to ensure a precise and unbiased inference. Traditionally, such kinematic information has been included in the inference after the image modeling, using spherical Jeans approximations to match the measured velocity dispersion integrated within an aperture. However, as spatially resolved kinematic measurements become available via IFU data, more sophisticated dynamical modeling is necessary. Such kinematic modeling is expensive, and constitutes a computational bottleneck which we aim to overcome with our Stellar Kinematics Neural Network (SKiNN). SKiNN emulates axisymmetric modeling using a neural network, quickly synthesizing from a given mass model a kinematic map which can be compared to the observations to evaluate a likelihood. With a joint lensing plus kinematic framework, this likelihood constrains the mass model at the same time as the imaging data. We show that SKiNN's emulation of a kinematic map is accurate to considerably better precision than can be measured (better than $1\%$ in almost all cases). Using SKiNN speeds up the likelihood evaluation by a factor of $\sim 200$. This speedup makes dynamical modeling economical, and enables lens modelers to make effective use of modern data quality in the JWST era.

To accurately characterize the planets a star may be hosting, stellar parameters must first be well-determined. $\tau$ Ceti is a nearby solar analog and often a target for exoplanet searches. Uncertainties in the observed rotational velocities have made constraining $\tau$ Ceti's inclination difficult. For planet candidates from radial velocity (RV) observations, this leads to substantial uncertainties in the planetary masses, as only the minimum mass ($m \sin i$) can be constrained with RV. In this paper, we used new long-baseline optical interferometric data from the CHARA Array with the MIRC-X beam combiner and extreme precision spectroscopic data from the Lowell Discovery Telescope with EXPRES to improve constraints on the stellar parameters of $\tau$ Ceti. Additional archival data were obtained from a Tennessee State University Automatic Photometric Telescope and the Mount Wilson Observatory HK project. These new and archival data sets led to improved stellar parameter determinations, including a limb-darkened angular diameter of $2.019 \pm 0.012$ mas and rotation period of $46 \pm 4$ days. By combining parameters from our data sets, we obtained an estimate for the stellar inclination of $7\pm7^\circ$. This nearly-pole-on orientation has implications for the previously-reported exoplanets. An analysis of the system dynamics suggests that the planetary architecture described by Feng et al. (2017) may not retain long-term stability for low orbital inclinations. Additionally, the inclination of $\tau$ Ceti reveals a misalignment between the inclinations of the stellar rotation axis and the previously-measured debris disk rotation axis ($i_\mathrm{disk} = 35 \pm 10^\circ$).

We present RDSim, a fast and comprehensive framework for the simulation of the radio emission and detection of downgoing air showers. It can handle any downgoing shower that can be simulated with ZHAireS including those induced by CC and NC neutrino interactions and $\tau$ decays. RDSim is based on a superposition toymodel that disentangles the Askaryan and geomagnetic components of the shower emission. By using full ZHAireS simulations as input, it is able to estimate the full radio footprint on the ground. A single input simulation at a given energy and arrival direction can be scaled in energy and rotated in azimuth by taking into account all relevant effects. This makes it possible to simulate a huge number of geometries and energies using just a few ZHAireS input simulations. The framework takes into account the main characteristics of the detector, such as trigger setups, thresholds and antenna patterns. To accommodate arrays that use particle detectors for triggering, such as the Auger RD extension, it also features a second toymodel to estimate the muon density at ground level, which is used to perform simple particle trigger simulations. It's speed makes it possible to investigate in detail events with a very low trigger probability, as well as many geometrical effects due to the array layout. In case more detailed studies of the radio detection are needed, RDSim can also be used to sweep the phase-space for the efficient creation of dedicated full simulation sets. This is particularly important in the case of neutrino events, that have extra variables that greatly impact shower characteristics, such as interaction or $\tau$ decay depth as well as the type of interaction and it's fluctuations.

The standard siren method using gravitational-wave observations has great potential to resolve the tension in measurements of the Hubble constant from different experiments. To realize this goal, we must thoroughly understand the sources of potential systematic bias. Among the known sources of systematic uncertainties, selection effects originating from electromagnetic counterpart observations of gravitational-wave sources may dominate the measurements and no method to mitigate this effect is currently established. In this Letter, we develop a new formalism to mitigate the counterpart selection effect. With realistic examples, we show that our formalism can reduce the systematic uncertainty of standard siren Hubble constant measurement to less than 0.6%. We conclude with how to apply our formalism to different electromagnetic emissions and observing scenarios.

We present the first spatially resolved maps of gas-phase metallicity for dust-obscured star-forming galaxies (DSFGs) at $z$ $\sim$ 4, from the JWST TEMPLATES Early Release Science program, derived from NIRSpec integral field unit spectroscopy of the H$\alpha$ and [NII] emission lines. Empirically derived literature optical line calibrations are used to determine that the sources are highly metal rich, with both appearing to display regions of supersolar metallicity, particularly in SPT2147-50. While we cannot rule out shocks or AGN in these regions, we suggest that the two systems have already undergone significant enrichment as a result of their extremely high star-formation rates. Utilising ALMA rest-frame 380$\mu$m continuum and [CI]($^3$P$_2$-$^3$P$_1$) line maps we compare metallicity and gas-to-dust ratio variations in the two galaxies, finding the two to be anticorrelated on highly resolved spatial scales, consistent with various literature studies of $z$ $\sim$ 0 galaxies. The data are indicative of the enormous potential of JWST to probe the enrichment of the interstellar medium on $\sim$kpc scales in extremely dust-obscured systems at $z$ $\sim$ 4 and beyond.

Galaxy clusters enable unique opportunities to study cosmology, dark matter, galaxy evolution, and strongly-lensed transients. We here present a new cluster-finding algorithm, CluMPR (Clusters from Masses and Photometric Redshifts), that exploits photometric redshifts (photo-z's) as well as photometric stellar mass measurements. CluMPR uses a 2-dimensional binary search tree to search for overdensities of massive galaxies with similar redshifts on the sky and then probabilistically assigns cluster membership by accounting for photo-z uncertainties. We leverage the deep DESI Legacy Surveys grzW1W2 imaging over one-third of the sky to create a catalogue of ~ 300,000 galaxy cluster candidates out to z = 1, including tabulations of member galaxies and estimates of each cluster's total stellar mass. Compared to other methods, CluMPR is particularly effective at identifying clusters at the high end of the redshift range considered (z = 0.75-1), with minimal contamination from low-mass groups. These characteristics make it ideal for identifying strongly lensed high-redshift supernovae and quasars that are powerful probes of cosmology, dark matter, and stellar astrophysics. As an example application of this cluster catalogue, we present a catalogue of candidate wide-angle strongly-lensed quasars in Appendix C. The five best candidates identified from this sample include two known lensed quasar systems and a possible changing-look lensed QSO with SDSS spectroscopy. All code and catalogues produced in this work are publicly available (see Data Availability).

We exploit a sample of 80,000 SDSS central galaxies to investigate the effect of AGN feedback on their evolution. We trace the demographics of optically-selected AGN (Seyferts) as a function of their internal properties and environment. We find that the preeminence of AGN as the dominant ionising mechanism increases with stellar mass, overtaking star-formation for galaxies with $M_\text{stellar} \geq 10^{11}M_\odot$. The AGN fraction changes systematically with the galaxies' star-formation activity. Within the blue cloud, this fraction increases as star-formation activity declines, reaching a maximum near the green valley ($\sim 17 \pm 4\%$), followed by a decrease as the galaxies transition into the red sequence. This systematic trend provides evidence that AGN feedback plays a key role in regulating and suppressing star formation. In general, Seyfert central galaxies achieve an early-type morphology while they still host residual star formation. This suggests that, in all environments, the morphology of Seyfert galaxies evolves from late- to early-type before their star formation is fully quenched. Stellar mass plays an important role in this morphological transformation: while low mass systems tend to emerge from the green valley with an elliptical morphology (T-Type $\sim -2.5 \pm 0.7$), their high-mass counterparts maintain a spiral morphology deeper into the red sequence. In high-stellar-mass centrals, the fraction of Seyferts increases from early- to late-type galaxies, indicating that AGN feedback may be linked with the morphology and its transformation. Our analysis further suggests that AGN are fuelled by their own host halo gas reservoir, but when in group centrals can also increase their gas reservoir via interactions with satellite galaxies.

We present analysis of OI 63 micron and CO $J$ = 5-4 and 8-7 multi-position data in the W3A region and use it to develop a model for the extended low-density foreground gas that produces absorption features in the OI and $J$ = 5-4 CO lines. We employ the extinction to the exciting stars of the background HII region to constrain the total column density of the foreground gas. We have used the Meudon PDR code to model the physical conditions and chemistry in the region employing a two-component model with high density layer near the HII region responsible for the fine structure line emission, and an extended low density foreground layer. The best-fitting total proton density, constrained largely by the CO lines, is $n$(H) = 250 cm$^{-3}$ in the foreground gas, and 5$\times$10$^5$ cm$^{-3}$ in the material near the HII region. The absorption is distributed over the region mapped in W3A, and is not restricted to the foreground of either the embedded exciting stars of the HII region or the protostar W3 IRS5. The low-density material associated with regions of massive star formation, based on an earlier study by Goldsmith et al. (2021), is quite common, and we now see that it is extended over a significant portion of W3A. It thus should be included in modeling of fine structure line emission, including interpreting low-velocity resolution observations made with incoherent spectrometer systems, in order to use these lines as accurate tracers of massive star formation.

Intense mass loss through cool, low-velocity winds is a defining characteristic of low-to-intermediate mass stars during the asymptotic giant branch (AGB) evolutionary stage. Such winds return up ~80% of the initial stellar mass to the interstellar medium and play a major role in enriching it with dust and heavy elements. A challenge to understanding the physics underlying AGB mass loss is its dependence on an interplay between complex and highly dynamic processes, including pulsations, convective flows, shocks, magnetic fields, and opacity changes resulting from dust and molecule formation. I highlight some examples of recent advances in our understanding of late-stage stellar mass loss that are emerging from radio and (sub)millimeter observations, with a particular focus on those that resolve the surfaces and extended atmospheres of evolved stars in space, time, and frequency.

In 2018, the Greenland Telescope (GLT) started scientific observation in Greenland. Since then, we have completed several significant improvements and added new capabilities to the telescope system. This paper presents a full review of the GLT system, a summary of our observation activities since 2018, the lessons learned from the operations in the Arctic regions, and the prospect of the telescope.

We study solar wind magnetic turbulence with Parker Solar Probe during its first perihelion (at 0.17 au), from MHD to kinetic plasma scales. Using Morlet wavelet decomposition, we detect intermittent coherent structures, appearing as localized in time energetic events and covering a large range of scales. This implies the presence of embedded coherent structures from MHD down to sub-ion scales. For example, we observe a current sheet at MHD scales ($\sim 200$ s) with magnetic fluctuations inherent for a magnetic vortex at ion scales ($\sim 1$ s) and at sub-ion scales ($\sim 0.1$ s). The amplitude anisotropy of magnetic fluctuations is analyzed within nearly (i) $200$ structures at MHD scales, (ii) $\sim 2 \cdot 10^3$ events at ion scales and (iii) $\sim 10^4$ events at sub-ion scales. We compare it with crossings of model structures, such as Alfv\'en vortices, current sheets and magnetic holes, along the trajectories with various impact parameters. From this comparison, we conclude that at MHD and ion scales the majority of the structures are incompressible and represent mainly dipole Alfv\'en vortices (>80 %), monopole Alfv\'en vortices (<10 %) and current sheets (<10 %). On subion scales coherent structures represent monopole vortices (7 %), dipole vortices (49 %) current sheets (5 %) and magnetic holes (0.4 %). Around 40 % of structures at sub-ion scales do not fit any of the considered models. These events might represent compressible vortices.

This thesis explores the missing baryon problem in a computational context. An overview of the problem is given, along with a discussion regarding the relevance of the Circumgalactic Medium (CMG) and cosmological Zoom-in simulations. The mechanisms underlying the N-body code ChaNGa (H. Menon, et al., Computational Astrophysics and Cosmology 2, 1 (2015), arXiv:1409.1929), as well as the data visualization and analysis tools yt (M. J. Turk, et al., 192, 9 (2011), arXiv:1011.3514) and trident (Hummels, et al., 847, 59 (2017), arXiv:1612.03935) are presented at a conceptual level. Finally, a series of synthetic quasar absorption spectra produced by using trident on a ChaNGa dataset from (S. Roca-Fabrega, et al., 917, 64 (2021), arXiv:2106.09738) at redshift of $z\sim4$ are shown. The low relative flux exhibited by these spectra render absorption features indistinguishable from background noise, and possible explanations for this phenomena such as high redshift are discussed. Though the resulting spectra exhibit serious obstacles for both qualitative and quantitative interpretation, they provide a "proof-of-concept" for future work, demonstrating trident's compatibility with ChaNGa's data format. Future prospects for using trident to analyze the CGM as simulated by ChaNGa are discussed, as well as possible extensions of this project.

A sidereal rotation counting approach is demonstrated by resolving an ambiguity in the synodic rotation period of Koronis family member (3032) Evans, whose rotation lightcurves' features did not easily distinguish between doubly- and quadruply-periodic. It confirms that Evans's spin rate does not exceed the rubble-pile spin barrier and thus presents no inconsistency with being a ~14-km reaccumulated object. The full spin vector solution for Evans is comparable to those for the known prograde low-obliquity comparably-fast rotators in the Koronis family, consistent with having been spun up by YORP thermal radiation torques.

In 2018, Jewitt identified the "The Trojan Color Conundrum", namely that Neptune's Trojan asteroids (NTs) had no ultra-red members, unlike the the nearby Kuiper Belt. Since then, numerous ultra-red NTs have been discovered, seemingly resolving this conundrum (Lin et al. 2019; Bolin et al.12 2023). However, it is still unclear whether or not the Kuiper Belt has a color distribution consistent with the NT population, as would be expected if it were the source population. In this work, we present a new photometric survey of 15 out of 31 NTs. We utilized the Sloan g'r'i'z' filters on the IMACS f/4 instrument which is mounted on the 6.5m Baade telescope. In this survey, we identify four NTs as being ultra-red using a Principal Component Analysis (PCA). This result brings the ratio of red to ultra-red NTs to 7.75:1, more consistent with the corresponding Trans-Neptunian Object (TNO) ratio of 4-11:1. We also identify three targets as being blue (nearly Solar) in color. Such objects may be C-type surfaces, but we see more of these blue NTs than has been observed in the Kuiper Belt (Seccull et al. 2018). Finally, we show that there are hints of a color-absolute magnitude (H) correlation, with larger H (smaller sized, lower albedo) tending to be more red, but more data is needed to confirm this result. The origin of such a correlation remains an open question which will be addressed by future observations of the surface composition of these targets and their rotational properties.

Some interaction-powered supernovae have long rise times of more than 100 days. We show that such long rise times are naturally expected if circumstellar matters (CSM) have a flat density structure (s <~ 1.5, where rho_CSM ~ r^{-s}). In such cases, bolometric luminosities from the CSM interaction keep increasing as long as the CSM interacts with the outer layers of the SN ejecta. Thus, the rise time is determined by the dynamical timescale in which the reverse shock propagates the outer layers of the SN ejecta, not by the timescales in which photons diffuse in the CSM as often considered. Interaction-powered supernovae with very long rise times can be an important probe of extensive non-steady mass loss in massive stars.

To demonstrate the magnetic energy dissipation via relativistic shocks, we carry out spherically symmetrical one-dimensional special relativistic magneto-hydrodynamic simulations of highly magnetised outflows with an adaptive mesh refinement method. We first investigate the detail of the dynamical energy dissipation via interaction between a single ejecta and an external medium. The energy dissipation timescales, which affect the early behaviour of the afterglow emission in gamma-ray bursts, are estimated for a wide range of magnetisation. In addition, we demonstrate the internal shock dissipation in multiple interactions between magnetically dominated relativistic ejecta and kinetically dominated non-relativistic winds. Our numerical results show that almost 10% of the magnetic energy in the ejecta can be converted into the thermal energy of the relativistic and low-magnetised outflows via shocks in the rarefaction waves or the winds. Such hot and less magnetised outflows are relevant for observed non-thermal emissions in blazars or gamma-ray bursts.

X-ray photons can penetrate deep into the intergalactic medium (IGM), leading to preheating of the IGM prior to cosmic reionization. X-ray preheating wipes out some of the small-scale structures that would otherwise be present prior to the passage of an ionization front. Accurate modeling of the small-scale structure is vital to the post-reionization IGM since the small-scale structure is ultimately the dominant source of long-lasting relics from hydrogen reionization. However, the precise impact of X-ray preheating in the fossils from hydrogen reionization is highly uncertain. In this work, we explore and establish for the first time, the long-lasting impact of X-ray preheating in the memory of reionization present in the IGM via hydrodynamic simulations with high-mass resolution. We find that the addition of X-ray preheating astrophysics leads to an overall lesser impact of the memory of reionization on the Ly$\alpha$ forest -- depending on specific X-ray prescription -- at low redshifts ($z \sim 2$) with respect to a model with no X-ray preheating. However, at high redshifts ($z \sim 4$), our results indicate a strengthening of the memory of reionization in the Ly$\alpha$ forest because the IGM becomes more transparent compared to the scenario with no preheating. Thus, the absence of X-ray preheating in Ly$\alpha$ modeling can lead to a biased inference of cosmological parameters. Nevertheless, optimistically, the inclusion of X-ray preheating emerges as a promising novel avenue to probe the astrophysics of cosmic dawn.

Observations of elliptical galaxies suggest that black holes (BHs) might serve as dark energy candidates, coupled to the expansion of the Universe. According to this hypothesis, the mass of a BH could increase as the Universe expands. BH low-mass X-ray binaries (LMXBs) in the Galactic disk were born several gigayears ago, making the coupling effect potentially significant. In this work, we calculate the evolution of BH binaries with a binary population synthesis method to examine the possible influence of cosmologically-coupled growth of BHs, if it really exists. The measured masses of the compact objects in LMXBs show a gap around $\sim 2.5-5~{\rm M_\odot}$, separating the most massive neutron stars from the least massive BHs. Our calculated results indicate that, considering the mass growth seem to (partially) account for the mass gap and the formation of compact BH LMXBs, alleviating the challenges in modeling the formation and evolution of BH LMXBs with traditional theory. However, critical observational evidence like the detection of intermediate-mass black hole binaries is required to test this hypothesis.

Previous theoretical works using the pre-merger orbital evolution of coalescing neutron stars to constrain properties of dense nuclear matter assume a gravitational wave phase uncertainty of a few radians, or about a half cycle. However, recent studies of the signal from GW170817 and next generation detector sensitivities indicate actual phase uncertainties at least twenty times better. Using these refined estimates, we show that future observations of nearby sources like GW170817 may be able to reveal neutron star properties beyond just radius and tidal deformability, such as the matter composition and/or presence of a superfluid inside neutron stars, via tidal excitation of g-mode oscillations. Data from GW170817 already limits the amount of orbital energy that is transferred to the neutron star to <2x10^47 erg and the g-mode tidal coupling to Qmode<10^-3 at 50 Hz (5x10^48 erg and 4x10^-3 at 200 Hz), and future observations and detectors will greatly improve upon these constraints. In addition, analysis using general parameterization models that have been applied to the so-called p-g instability show that the instability already appears to be restricted to regimes where the mechanism is likely to be inconsequential; in particular, we show that the number of unstable modes is <<100 at <~100 Hz, and next generation detectors will essentially rule out this mechanism (assuming that the instability remains undetected). Finally, we illustrate that measurements of tidal excitation of r-mode oscillations in nearby rapidly rotating neutron stars are within reach of current detectors and note that even non-detections will limit the inferred inspiralling neutron star spin rate to <20 Hz, which will be useful when determining other parameters such as neutron star mass and tidal deformability.

The Extragalactic Background Light (EBL) that spans the UV-IR band originates from direct and dust-reprocessed starlight integrated over the history of the Universe. EBL measurements are very challenging due to foreground emission like the zodiacal light and interplanetary dust emission. Indeed, some optical/NIR direct measurements overpredict EBL models based on galaxy counts. On the other hand, there is some debate on possible additional components of the Optical-NIR photon density: e.g., population-III stars, axion-photon decay, direct collapse of black holes, intra-halo light etc. Owing to the absorption of Very High Energy (VHE) $\gamma$ rays by interaction with EBL photons, we study the prospects of accommodating an additional population of EBL sources in the Optical-NIR band on top of the standard galaxy-count--based component. To this aim we use 105 VHE spectra of 37 blazars with known redshifts, $0.03<z<0.94$. We correct the observed spectra for absorption by our model EBL. By requiring the intrinsic spectra to be non-concave and with a VHE spectral index $>$1.5, we estimate, at different wavelengths, upper limits to the additional low energy photon fields which would contribute to the absorption of $\gamma$-rays. Considering these limits we suggest that there is room for photons from Pop III stars and axion-like particle (ALP) annihilation. However, these additional hypothetical photon fields are bound to fall significantly below direct published EBL measurements by several instruments, and therefore our limits are either in tension or even inconsistent with such measurements.

We study two possible cosmological consequences of a first-order phase transition in the temperature range of $1$ GeV to $10^3$ TeV: the generation of a stochastic gravitational wave background (SGWB) within the sensitivity of the Laser Interferometer Space Antenna (LISA) and, simultaneously, primordial magnetic fields that would evolve through the Universe's history and could be compatible with the lower bound from $\gamma$-ray telescopes on intergalactic magnetic fields (IGMF) at present time. We find that, if even a small fraction of the kinetic energy in sound waves is converted into MHD turbulence, a first order phase transition occurring at temperature between $1$ and $10^6$ GeV can give rise to an observable SGWB signal in LISA and, at the same time, an IGMF compatible with the lower bound from the $\gamma$-ray telescope MAGIC, for all proposed evolutionary paths of the magnetic fields throughout the radiation dominated era. For two values of the fraction of energy density converted into turbulence, $\varepsilon_{\rm turb}=0.1$ and $1$, we provide the range of first-order phase transition parameters (strength $\alpha$, duration $\beta^{-1}$, bubbles wall speed $v_w$, and temperature $T_*$), together with the corresponding range of magnetic field strength $B$ and correlation length $\lambda$, that would lead to the SGWB and IGMF observable with LISA and MAGIC. The resulting magnetic field strength at recombination can also correspond to the one that has been proposed to induce baryon clumping, previously suggested as a possible way to ease the Hubble tension. In the limiting case $\varepsilon_{\rm turb} \ll 1$, the SGWB is only sourced by sound waves, however, an IGMF is still generated. We find that values as small as $\varepsilon_{\rm turb} \sim O(10^{-13})$ (helical) and $O (10^{-9})$ (non-helical) can provide IGMF compatible with MAGIC's lower bound.

Forming planets around young, fast-rotating solar-like stars are exposed to an intense X-ray/extreme ultraviolet radiation field and strongly magnetized stellar winds, as a consequence of the high magnetic activity of these stars. Under these conditions, Earth-like exoplanets may experience a rapid loss of their primordial hydrogen atmospheres, resulting in atmosphere-less rocky obstacles for the stellar winds. The interaction of stellar winds with those planets leads to the formation of potentially observable structures due to the formation of large-scale magnetic field and density disturbances in the vicinity of these planets, such as bow shocks, induced magnetospheres and comet-like tails. In this work, we study the interaction between the stellar winds of active, fast-rotating solar-like stars in the superfast-magnetosonic regime with Earth-like, unmagnetized, tenuous atmosphere, planetary obstacles through numerical 3D simulations using the PLUTO magnetohydrodynamical code. The properties of AB Doradus, a nearby young star with a small rotation period (0.51 days) and a strong flaring activity, have been used to parameterize this early wind state. Bow shock and induced magnetosphere formation are characterized through the alfv\'enic Mach number MA of the wind, for different stellar wind configurations. Large bow shocks, up to an extension of ~7.0 planetary radii are found for low-MA winds. The general increase of density, temperature and magnetic field in these large-scale structures formed around planets may result in potentially detectable spectral signatures.

Within the framework of investigating the link between central super massive black holes in the core of galaxies and the galaxies themselves, we detected a variable X-ray source in the center of CGCG 077-102 NED02, member of the CGCG 077-102 galaxy pair within the Abell 2063 galaxy cluster. Our goal was then to combine X-ray and optical data to demonstrate that this object harbors an active super massive black hole in its core, and to relate this to the dynamical status of the galaxy pair within the Abell 2063 cluster. We used Chandra and XMM-Newton archival data to derive the X-ray spectral shape and variability. We also obtained optical spectroscopy to detect the expected emission lines that are typically found in Active Galactic Nuclei. And we finally used public ZTF imaging data to investigate the optical variability. There is no evidence of multiple X-ray sources or extended component within CGCG 077-102 NED02. Single X-ray spectral models fit well the source. Non-random significant X-ray flux inter-observation X-ray variabilities were detected, between ~4days for short term variations and up to ~700days for long term variations. Optical spectroscopy points toward a passive galaxy for CGCG 077-102 NED01 and a Seyfert for CGCG 077-102 NED02. We did not detect short-term variability in the optical ZTF light curves. However, we found a significant long-term stochastic variability in the g-band that can be well described by the damped random walk model. Finally, the CGCG 077-102 galaxy pair is deeply embedded within the Abell 2063 potential, and has underwent the cluster influence for a long time. Our observations point toward a moderatly massive black hole in the center of CGCG 077-102 NED02, of ~10^6 Msol. CGCG 077-102 NED02 is not heavily obscured, perhaps due to surrounding intra cluster medium ram pressure stripping.

The observations and comprehensive study of intermediate initial mass stars at the late stages of evolution, and after the asymptotic giant branch (AGB) in particular, are of crucial importance to identify the common properties for the stars of given group and to reveal binaries among them. This work aims to investigate photometric and spectral peculiarities of a poorly studied post-AGB candidate and infrared source IRAS 07253-2001. We present the new multicolour $UBVR_{C}I_{C}YJHK$ photometry obtained with the telescopes of the Caucasian mountain observatory and analyse it together with the data acquired by the All Sky Automated Survey for SuperNovae. We report on the detection of multiperiod brightness variability caused by pulsations. A beating of close periods, the main one of 73 days and additional ones of 68 and 70 days, leads to amplitude variations. We have also detected a long-term sine trend in brightness with a period of nearly 1800 days. We suppose it to be orbital and IRAS 07253-2001 to be binary. Based on new low-resolution spectroscopic data obtained with the 2.5-m telescope of the Caucasian mountain observatory in 2020 and 2023 in the $\lambda$3500-7500 wavelength range we have identified spectral lines and compiled a spectral atlas. We have found the [N II], [Ni II] and [S II] forbidden emission lines in the spectrum and discuss their origin. The H$\alpha$ line has a variable double-peaked emission component. We have derived preliminary estimates of the star's parameters and detected a variation of radial velocity with a peak-to-peak amplitude of about 30 km s$^{-1}$.

Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosphere. This Letter focuses on CX Tau, a ${\sim}0.4\,M_\odot$ star with an extended gas emission and a compact and apparently structureless dust disk, with an average millimeter flux when compared to Class II sources in Taurus. We present Karl G. Jansky Very Large Array (VLA) observations in 4 bands (between 9.0 mm and 6.0 cm) and combine them with archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), the Submillimeter Array (SMA) and the Plateau de Bure Interferometer (PdBI). Such a multiwavelength approach allows to separate the dust continuum from other emissions. After isolating the dust thermal emission, we derived an upper limit of the dust disk extent at 1.3 cm which is consistent with theoretical predictions of a radial drift-dominated disk. Centimeter data show a peculiar behavior: deep observations at 6.0 cm did not detect the source, while at 1.3 cm the flux density is anomalously higher than adjacent bands. Intraband spectral indices suggest a dominant contribution from free-free emission, whereas gyro-synchrotron emission is excluded. To explain these observations, we propose strong variability of the free-free emission with timescales shorter than a month. Another possible interpretation is the presence of anomalous microwave emission from spinning dust grains.

Binary systems constitute a valuable astrophysics tool for testing our understanding of stellar structure and evolution. Systems containing a oscillating component are interesting as asteroseismology offers independent parameters for the oscillating component that aid the analysis. About 150 of such systems are known in the literature. To enlarge the sample of these benchmark objects, we crossmatch the Two-Body-Orbit Catalogue (TBO) of Gaia DR3, with catalogs of confirmed solar-like oscillators on the main-sequence and red-giant phase from NASA Kepler and TESS. We obtain 954 new binary system candidates hosting solar-like oscillators, of which 45 and 909 stars are on the main sequence and red-giant, resp., including 2 new red giants in eclipsing systems. 918 oscillators in potentially long-periodic systems are reported. We increase the sample size of known solar-like oscillators in binary systems by an order of magnitude. We present the seismic properties of the full sample and conclude that the grand majority of the orbital elements in the TBO is physically reasonable. 82% of all TBO binary candidates with multiple times with APOGEE are confirmed from radial-velocity measurement. However, we suggest that due to instrumental noise of the TESS satellite the seismically inferred masses and radii of stars with $\nu_\textrm{max}$$\lesssim$30$\mu$Hz could be significantly overestimated. For 146 giants the seismically inferred evolutionary state has been determined and shows clear differences in their distribution in the orbital parameters, which are accounted the accumulative effect of the equilibrium tide acting in these evolved binary systems. For other 146 systems hosting oscillating stars values for the orbital inclination were found in the TBO. From testing the TBO on the SB9 catalogue, we obtain a completeness factor of 1/3.

4U 1543-47 underwent its brightest outburst in 2021 after two decades of inactivity. During its decay phase, AstroSat conducted nine observations of the source spanning from July $1^{st}$ to September $26^{th}$, 2021. The first three observations were performed with an offset of 40 arcmin with AstroSat/LAXPC, while the remaining six were on-axis observations. In this report, we present a comprehensive spectral analysis of the source as it was in the High/Soft state during the entire observation period. The source exhibited a disk-dominated spectra with a weak high-energy tail (power-law index $\geq2.5$) and a high inner disk temperature ($\sim$0.84 keV). Modelling the disk continuum with non-relativistic and relativistic models, we find inner radius to be significantly truncated at $>$$10~R_g$. Alternatively, to model the spectral evolution with the assumption that the inner disk is at the ISCO, it is necessary to introduce variation in the spectral hardening in the range $\sim$1.5-1.9.

Gyrochronology enables the derivation of ages of late-type main sequence stars based on their rotation periods and a mass proxy, such as color. It has been explored in open clusters, but a connection to field stars has yet to be successfully established. We explore the rotation rates of wide binaries, representing enlightening intermediaries between clusters and field stars, and their overlap with those of open cluster stars. We investigated a recently created catalog of wide binaries, matched the cataloged binaries to observations by the Kepler mission (and its K2 extension), validated or re-derived their rotation periods, identified 283 systems where both stars are on the main sequence and have vetted rotation periods, and compared the systems with open cluster data. We find that the vast majority of these wide binaries (236) line up directly along the curvilinear ribs defined by open clusters in color-period diagrams or along the equivalent interstitial gaps between successive open clusters. The parallelism in shape is remarkable. Twelve additional systems are clearly rotationally older. The deviant systems, a minority, are mostly demonstrably hierarchical. Furthermore, the position of the evolved component in the color-magnitude diagram for the additional wide binary systems that contain one is consistent with the main sequence component's rotational age. We conclude that wide binaries, despite their diversity, follow the same spindown relationship as observed in open clusters, and we find that rotation-based age estimates yield the same ages for both components in a wide binary. This suggests that cluster and field stars spin down in the same way and that gyrochronology can be applied to field stars to determine their ages, provided that they are sufficiently distant from any companions to be considered effectively single.

A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at $\sim$7 MeV cm$^{-2}$ (time period 1000~s) and $\sim$35 MeV cm$^{-2}$ (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total gravitational wave energy. These are the first limits on UHE photons from GW sources.

We present a comprehensive analysis of the $\Lambda_{\rm s}$CDM model, which explores the recent conjecture suggesting a rapid transition of the Universe from anti-de Sitter vacua to de Sitter vacua (viz., the cosmological constant switches sign from negative to positive) at redshift ${z_\dagger\sim 2}$, inspired by the graduated dark energy (gDE) model. Our analysis shows that, predicting $z_\dagger\approx1.7$, $\Lambda_{\rm s}$CDM simultaneously addresses the major cosmological tensions of the standard $\Lambda$CDM model, viz., the Hubble constant $H_0$, the Type Ia Supernovae absolute magnitude $M_{\rm B}$, and the growth parameter $S_8$ tensions, along with other less significant tensions such as the BAO Lyman-$\alpha$ discrepancy.

We present numerical results from a parameter study of the standing accretion shock instability (SASI), investigating the impact of general relativity (GR) on the dynamics. Using GR hydrodynamics and gravity, and non-relativistic (NR) hydrodynamics and gravity, in an idealized model setting, we vary the initial radius of the shock and, by varying its mass and radius in concert, the proto-neutron star (PNS) compactness. We investigate two regimes expected in a post-bounce core-collapse supernova (CCSN): one meant to resemble a relatively low-compactness configuration and one meant to resemble a relatively high-compactness configuration. We find that GR leads to a longer SASI oscillation period, with ratios between the GR and NR cases as large as 1.29 for the high-compactness suite. We also find that GR leads to a slower SASI growth rate, with ratios between the GR and NR cases as low as 0.47 for the high-compactness suite. We discuss implications of our results for CCSN simulations.

The Sunyaev-Zeldovich thermal (tSZ) and kinetic (kSZ) effects can be used to constrain the thermodynamic properties of pressure and density, respectively, of galaxies and their surrounding regions. As SZ observations continue to improve, it is important to understand any modeling systematics when inferring properties from the data. Thus, a pipeline to forward model observed SZ profiles was developed called Mop-c-GT. Previous studies have used this repository to create modeled SZ profiles by selecting halos from the IllustrisTNG simulation and found significant differences between the simulated profiles and those recently measured by the Atacama Cosmology Telescope. There are many uncertainties involved in modeling observed samples and in the forward modeling process, so in this study, we explore methods implemented in Mop-c-GT and in the selection of the simulated halos to test the effects on the resulting modeled profiles. After testing several methods within the forward modeling process and varying the halo selection from the simulation, we find minimal differences between the simulated tSZ profiles of the original calculation and the updated methods, indicating that the observations still pose a challenge for the numerical methods used to describe the astrophysics of these systems.

The smoothed particle hydrodynamics technique strongly relies on the proper choice of interpolating functions. In this work, we revisit and extend the main properties of a family of interpolators called $Sinc~kernels$ and compare them with those of the widely used family of Wendland kernels. We show that a linear combination of low and high-order Sinc kernels generates good quality interpolators, which are resistant to the pairing instability while keeping good sampling properties in a wide range of neighbor interpolating points, $60\le n_b\le 400$. We show that a particular case of this linear mix of Sincs produces a well-balanced and robust kernel, that improves previous results in the Gresho-Chan vortex experiment even when the number of neighbors is not large, while yielding a good convergence rate. Although such a mixing technique is ideally suited for Sinc kernels owing to their excellent flexibility, it can be easily applied to other interpolating families such as the B-splines and Wendland kernels.

Transmission spectroscopy supports the presence of uncharacterised, light-scattering and -absorbing aerosols in the atmospheres of many exoplanets. The complexity of factors influencing the formation, 3-D transport, radiative impact, and removal of aerosols makes it challenging to match theoretical models to the existing data. Our study simplifies these factors to focus on the interaction between planetary general circulation and haze distribution at the planetary limb. We use an intermediate complexity general circulation model, ExoPlaSim, to simulate idealised organic haze particles as radiatively active tracers in the atmospheres of tidally locked terrestrial planets for 32 rotation rates. We find three distinct 3-D spatial haze distributions, corresponding to three circulation regimes, each with a different haze profile at the limb. All regimes display significant terminator asymmetry. In our parameter space, super-Earth-sized planets with rotation periods greater than 13 days have the lowest haze optical depths at the terminator, supporting the choice of slower rotators as observing targets.

Asteroseismology provides a unique opportunity to probe the interiors of evolved stars and constrain their internal rotation. The correct reproduction of the core rotation evolution is key to understanding the internal processes involved in low-mass stars. We explore the efficiency required to reproduce the behaviour of the transport of angular momentum (AM) in view of asteroseismic constraints. We computed a series of models and investigated an updated AM transport by including a time-dependent extra viscosity related to the AMRI. We compared our predictions to the asteroseismic measurements of the core and surface rotation of a sample of SGB and RGB stars. We confirm that a time-dependent additional viscosity is required to reproduce the general behaviour of the core rotation rate along evolution. We show that it results in stronger Li and Be depletions for low-mass stars over evolution. We confirm that predicted Li abundances at the RGB bump by classical models, commonly used as references, cannot reproduce the Li depletion along the MS and evolved phases of stellar evolution. We show that the observed amount of Li of stars less massive than 1Msun leads to a discrepancy between model predictions and observations at the RGB bump. We show that a semi-parametric model can reproduce the rotational behaviour along the first phases of evolution well, with the exception of the sharp transition observed during the SGB phase. This suggests that two distinct transport processes may be involved. The processes required to transport chemicals during the MS, and AM until the RGB phase impact the Li depletion all along the evolutionary duration. A good prediction of the Li abundance at young phases places strong constraints on the predicted one at more evolved phases. It also impacts the threshold that defines Li-rich giant stars, showing that classical models tend to overestimate its value.

The prediction of spin magnitudes in binary black hole and neutron star mergers is crucial for understanding the astrophysical processes and gravitational wave (GW) signals emitted during these cataclysmic events. In this paper, we present a novel Neuro-Symbolic Architecture (NSA) that combines the power of neural networks and symbolic regression to accurately predict spin magnitudes of black hole and neutron star mergers. Our approach utilizes GW waveform data obtained from numerical relativity simulations in the SXS Waveform catalog. By combining these two approaches, we leverage the strengths of both paradigms, enabling a comprehensive and accurate prediction of spin magnitudes. Our experiments demonstrate that the proposed architecture achieves an impressive root-mean-squared-error (RMSE) of 0.05 and mean-squared-error (MSE) of 0.03 for the NSA model and an RMSE of 0.12 for the symbolic regression model alone. We train this model to handle higher-order multipole waveforms, with a specific focus on eccentric candidates, which are known to exhibit unique characteristics. Our results provide a robust and interpretable framework for predicting spin magnitudes in mergers. This has implications for understanding the astrophysical properties of black holes and deciphering the physics underlying the GW signals.

The orbits of some warm Jupiters are highly inclined (20$^\circ$-50$^\circ$) to those of their exterior companions. Comparable misalignments are inferred between the outer and inner portions of some transition discs. These large inclinations may originate from planet-planet and planet-disc secular resonances that sweep across interplanetary space as parent discs disperse. The maximum factor by which a seed mutual inclination can be amplified is of order the square root of the angular momentum ratio of the resonant pair. We identify those giant planet systems (e.g. Kepler-448 and Kepler-693) which may have crossed a secular resonance, and estimate the required planet masses and semimajor axes in transition discs needed to warp their innermost portions (e.g. in CQ Tau). Passage through an inclination secular resonance could also explain the hypothesized large mutual inclinations in apsidally-orthogonal warm Jupiter systems (e.g. HD 147018).

The non-thermal nature of the stochastic gravitational-wave background of cosmological origin (CGWB) poses a challenge in defining the initial conditions for the graviton overdensity. Specifically, the adiabatic initial condition, which holds for Cosmic Microwave Background (CMB) photons, is not guaranteed a priori for the primordial GWs. In this letter, we compute the initial conditions for the cosmological background generated by quantum fluctuations of the metric during inflation. Our analysis reveals that adiabatic initial conditions are no longer valid. The violation of adiabaticity arises from the presence of independent tensor perturbations during inflation, which behave as two extra fields that affect the standard single-clock argument. Since the energy density of the CGWB is subdominant compared to ordinary matter, gravitational radiation plays a negligible role in Einstein's equations. Therefore, the only way to compute the initial conditions is to perturb the energy-momentum tensor defined in terms of the gravitational strain. A direct consequence of our finding is that the initial conditions from inflation enhance the total CGWB angular power-spectrum by an order of magnitude compared to the standard adiabatic case.

The discovery of ultra-high-energy (UHE) neutrinos has the potential to offer unique insight into fundamental questions. To capitalize on the upcoming opportunity provided by new UHE neutrino telescopes, we provide state-of-the-art forecasts of the discovery of a diffuse flux of UHE neutrinos over the next 10-20 years, focusing on neutrino radio-detection in the planned IceCube-Gen2 detector. We use state-of-the-art flux predictions and detector modeling. We find that, even under conservative analysis choices, most benchmark UHE neutrino flux models from the literature may be discovered within 10 years of detector exposure, with many discoverable sooner, and may be distinguished from each other. Our results demonstrate the transformative potential of next-generation UHE neutrino telescopes.

GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta.

The Effective Field Theory of Large-Scale Structure (EFTofLSS) attempts to amend some of the shortcomings of the traditional perturbative methods used in cosmology. It models the evolution of long-wavelength perturbations above a cutoff scale without the need for a detailed description of the short-wavelength ones. Short-scale physics is encoded in the coefficients of a series of operators composed of the long-wavelength fields, and ordered in a systematic expansion. As applied in the literature, the EFTofLSS corrects a summary statistic (such as the power spectrum) calculated from standard perturbation theory by matching it to $N$-body simulations or observations. This `bottom-up' construction is remarkably successful in extending the range of validity of perturbation theory. In this work, we compare this framework to a `top-down' approach, which estimates the EFT coefficients from the stress tensor of an $N$-body simulation, and propagates the corrections to the summary statistic. We consider simple initial conditions, viz. two sinusoidal, plane-parallel density perturbations with substantially different frequencies and amplitudes. We find that the leading EFT correction to the power spectrum in the top-down model is in excellent agreement with that inferred from the bottom-up approach which, by construction, provides an exact match to the numerical data. This result is robust to changes in the wavelength separation between the two linear perturbations. However, in our setup, the leading EFT coefficient does not always grow linearly with the cosmic expansion factor as assumed in the literature based on perturbative considerations. Instead, it decreases after orbit crossing takes place.

In standard inflationary cosmology, scalar and tensor perturbations grew as the Universe expanded and froze when their wavelengths exceeded the Hubble horizon, producing a tell-tale signature in the fluctuation spectrum and amplitude of the cosmic microwave background (CMB). But there are now very good reasons to examine whether structure formation could also have begun via the seeding of quantum fluctuations in a non-inflationary field. In this Letter, we study and compare the scalar and tensor modes produced in these two scenarios, and demonstrate that upcoming observations to measure the B-mode polarization of the CMB may be able to differentiate between them. Whereas both scalar and tensor modes should be observable if the field was inflationary, only scalar modes would be present in the CMB if it were not. Should gravity be purely classical, however, the tensor modes would have avoided canonical quantization in all cases, resulting in unmeasurably weak gravitational waves. A non-detection of B-mode polarization would thus not completely rule out inflation.

We investigate the late-time cosmological dynamics in a simple case of explicit spacetime-symmetry breaking. By expanding in a small symmetry-breaking coefficient we are able to write the Friedmann equations as $\Lambda$CDM + dynamical dark energy, which we show contains logarithmic dependence of the scale factor. We find that the dark energy equation of state displays divergences and phantom behaviour for certain values of the symmetry-breaking coefficient, where the Null Energy Condition is also broken. We also discuss the adiabatic sound speed of dark energy and compare the model to current constraints using the Chevallier-Polarski-Linder parametrisation.

General consistency requirements of Quantum Gravity demand the existence of one axion per Yang-Mills group. In this work, we consider theories with dark Yang-Mills sectors and investigate general phenomenological implications of these necessary axions. We carry out computations for two simple models, namely a pure Yang-Mills sector and $N$ exact Standard Model copies. For the former, the misalignment mechanism results in a minimal dark confinement scale $\Lambda_{\rm conf} \gtrsim 1 \, {\rm eV}$ if the dark sector axion is supposed to make up the dark matter. For the latter, the misalignment mechanism without fine-tuning of the initial misalignment angle places an upper bound on $N$ below the species bound. When the PQ symmetries are broken during inflation, the collective isocurvature fluctuations do not necessarily tighten the bound on the inflationary Hubble scale arising from a single axion. We also point out that axion stars collectively made from axions of different dark sectors with a suppressed mass spectrum are not possible. Lastly, for the two models at hand, intersector interaction through axion kinetic mixing leads to the existence of two distinct axion states. For a single dark YM sector, the upper bound $\Lambda_{\rm conf} \lesssim 10^{12} \, {\rm GeV}$ emerges from the stability requirement of the dark sector axion. For $N$ exact SM copies, the mass and photon coupling of the second state is completely determined after a potential measurement of the analogous parameters of the first axion.

We demonstrate that enhanced early galaxy formation can generically arise in axion-like particle (ALP) dark matter (DM) theories with the kinetic misalignment mechanism, potentially addressing the excess recently observed by the James Webb Space Telescope (JWST), while being consistent with existing constraints. We identify viable parameter regions with the ALP mass in the range of $10^{-22}~{\rm eV}<m_a<10^{-19}~\rm eV$. In addition, we show that the ALP parameter regions of interest can lead to intriguing complementary signatures in the small-scale structure of DM halos and existing experimental searches for ALPs.

Linear axion monodromy models modulated with higher powers of fields naturally realize the quantum-mechanical flux discharge mechanism for relaxing the cosmological constant toward zero. Working with multiple copies of superposed linear and quadratic flux monodromies, each copy spanned by a pair of fluxes, we show that when the axion is very massive and so effectively decoupled, the membrane discharges relax the cosmological constant toward an attractor $0 < \Lambda/\mpl^4 \ll 1$. If we restrict the flux variations and the intermediate flux values to never venture beyond a finite flux range, the terminal value of the cosmological constant will be tiny but finite. We show how it can reproduce the observed scale of dark energy, and explain how to incorporate matter sector phase transitions.

We propose a novel explanation for the 18 TeV gamma ray from GRB 221009A observed by LHAASO. High-energy neutrinos are converted into axion-like particles (ALPs) via their interaction with the cosmic neutrino background. Subsequently, ALPs are converted into high-energy photons in the magnetic field of our galaxy. We compute the fluxes of neutrinos, ALPs, and photons reaching Earth. IceCube's constraints on the neutrino flux from GRB 221009A translate into a severe upper bound on the photon flux. We find a range of parameters where all existing bounds are satisfied and the 18 TeV LHAASO photon can be explained. In the future, the specific correlation between the photon and neutrino flux reaching Earth from powerful neutrino sources with energies larger than 10 TeV such as GRBs or AGNs, can be used as a tool to differentiate our explanation from the alternatives suggested in the literature. We discuss how the interactions of our scenario can be embedded within electroweak gauge-invariant models, avoiding various cosmological and terrestrial bounds. We comment on the possibility of explaining the 251 TeV photon observed by the Carpet-2 detector, taking into account the bounds from the observation of high-energy neutrinos from TXS 0506+056.

We match scattering amplitudes in point particle effective field theory (EFT) and general relativity to extract low frequency dynamical tidal responses of rotating (Kerr) black holes to all orders in spin. In the conservative sector, we study local worldline couplings that correspond to the time-derivative expansion of the black hole tidal response function. These are dynamical (frequency-dependent) generalizations of the static Love numbers. We identify and extract couplings of three types of subleading local worldline operators: the curvature time derivative terms, the spin - curvature time derivative couplings, and quadrupole - octupole mixing operators that arise due to the violation of spherical symmetry. The first two subleading couplings are non-zero and exhibit a classical renormalization group running; we explicitly present their scheme-independent beta functions. The conservative mixing terms, however, vanish as a consequence of vanishing static Love numbers. In the non-conservative sector, we match the dissipation numbers at next-to-leading and next-to-next-to leading orders in frequency. In passing, we identify terms in the general relativity absorption probabilities that originate from tails and short-scale logarithmic corrections to the lowest order dissipation contributions.

Gravitational-wave observations by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo have provided us a new tool to explore the universe on all scales from nuclear physics to the cosmos and have the massive potential to further impact fundamental physics, astrophysics, and cosmology for decades to come. In this paper we have studied the science capabilities of a network of LIGO detectors when they reach their best possible sensitivity, called A#, and a new generation of observatories that are factor of 10 to 100 times more sensitive (depending on the frequency), in particular a pair of L-shaped Cosmic Explorer observatories (one 40 km and one 20 km arm length) in the US and the triangular Einstein Telescope with 10 km arms in Europe. We use a set of science metrics derived from the top priorities of several funding agencies to characterize the science capabilities of different networks. The presence of one or two A# observatories in a network containing two or one next generation observatories, respectively, will provide good localization capabilities for facilitating multimessenger astronomy and precision measurement of the Hubble parameter. A network of two Cosmic Explorer observatories and the Einstein Telescope is critical for accomplishing all the identified science metrics including the nuclear equation of state, cosmological parameters, growth of black holes through cosmic history, and make new discoveries such as the presence of dark matter within or around neutron stars and black holes, continuous gravitational waves from rotating neutron stars, transient signals from supernovae, and the production of stellar-mass black holes in the early universe. For most metrics the triple network of next generation terrestrial observatories are a factor 100 better than what can be accomplished by a network of three A# observatories.

In the violent post-merger of binary neutron-star mergers strong oscillations are present that impact the emitted gravitational-wave (GW) signal. The frequencies, temperatures and densities involved in these oscillations allow for violations of the chemical equilibrium promoted by weak-interactions, thus leading to a nonzero bulk viscosity that can impact dynamics and GW signals. We present the first simulations of binary neutron-star mergers employing the self-consistent and second-order formulation of the equations of relativistic hydrodynamics for dissipative fluids proposed by M\"uller, Israel and Stewart. With the spirit of obtaining a first assessment of the impact of bulk viscosity on the structure and radiative efficiency of the merger remnant we adopt a simplified approach for the viscosity, which we assume to be constant within the stars, but which we vary in strength for different binaries, thus exploring the possible behaviours and obtaining strict upper limits. In this way, we find that large bulk viscosities are very effective at damping the collision-and-bounce oscillations that characterize the dynamics of the stellar cores right after the merger. As a result, the $m=2$ deformations and the gravitational-radiation efficiency of the remnant are considerably reduced, with qualitative and quantitative changes in the post-merger spectrum that can be large in the case of the most extreme configurations. Overall, our crude but self-consistent results indicate that bulk viscosity reduces the energy radiated in GWs by $\lesssim 1\%$ in the (realistic) scenario of small viscosity, and by $\lesssim 15\%$ in the (unrealistic) scenario of large viscosity.

Since the beginning of the 20th century, a continuous evolution and perfection of what we today call the standard cosmological model has been produced, although some authors like to distinguish separate periods within this evolution. A possible historical division of the development of cosmology into six periods is: (1) the initial period (1917-1927); (2) the period of development (1927-1945); (3) the period of consolidation (1945-1965); (4) the period of acceptance (1965-1980); (5) the period of enlargement (1980-1998); and (6) the period of high-precision experimental cosmology (1998-now). The last period started with a epistemological optimism that has declined with time, and the expression "crisis in cosmology" is now stubbornly reverberating in the media. The initial expectation of removing the pending minor problems arising from the increased accuracy of measurements has backfired: the higher the precision with which the standard model tries to fit the data, the greater the number of tensions that arise, the problems proliferating rather than diminishing.

The Bardeen black hole holds historical significance as the first model of a regular black hole. Recently, there have been proposed interpretations of the Bardeen spacetime as quantum corrections to the Schwarzschild solution. Our study focuses on investigating the quasinormal modes and Hawking radiation of the Bardeen black hole. We have observed that previous studies on the quasinormal modes for the Bardeen black hole suffer from inaccuracies that cannot be neglected. Therefore, we propose accurate calculations of the quasinormal modes for scalar, electromagnetic, and neutrino fields in the Bardeen spacetime. Additionally, we have computed the grey-body factors and analyzed the emission rates of Hawking radiation. Even when the quantum correction is small and the fundamental mode only slightly differs from its Schwarzschild value, the first several overtones deviate at an increasingly stronger rate. This deviation leads to the appearance of overtones with very small real oscillation frequencies. This outburst of overtones is closely linked to the fact that the quantum-corrected black hole differs from its classical limit primarily near the event horizon. Moreover, the intensity of the Hawking radiation is significantly suppressed (up to three orders of magnitude) by the quantum correction.

This paper introduces the Circular Restricted n-Body Problem (CRNBP), an extension of the bicircular restricted four-body problem (BCR4BP) designed to describe the dynamics of an n-body system. In the CRNBP, each massive body in the system is constrained to follow a Keplerian motion, similar to the BCR4BP's artificial constraint. The CRNBP is an efficient alternative for trajectory design in multiple-body systems, particularly for outer planetary systems, as it requires integrating only six first-order ordinary differential equations compared to the 6N equations in an ephemerides model. By reproducing complex dynamical behaviors observed in ephemerides n-body problems, we demonstrate the structural stability of the CRNBP. Additionally, we propose a straightforward approach to relate the CRNBP with ephemerides, enabling the exploration of trajectory design possibilities before committing to a dedicated ephemerides analysis. This allows for the identification of general dynamical behaviors and provides valuable insights into the dynamics of multiple body systems. Finally, illustrative examples highlight the richness of trajectories and potential advantages of using the CRNBP for designing complex trajectories in outer planetary systems. The CRNBP proves to be a valuable tool for preliminary trajectory design, facilitating the identification of low-energy trajectories and providing a foundation for further exploration in future dedicated studies.

The purpose of this work is to investigate the formation and evaporation of the primordial black holes in the inflationary scenarios. Thermodynamic parameters such as mass, temperature and entropy are expressed in terms of NANOGrav frequency. By numerical calculations we show that the constraint on the mass range $10^{-5}kg-10^{50}kg$ is well confirmed. We discuss the relation between the redshift and the probability for gravitational wave source populations. A new parameter associated with the frequency and Hubble rate is presented, by which for the spectral index $n_{s}\approx 0.996$ and the Hubble constant $H_{0}\approx 67.27km.s^{-1}.Mpc^{-1}$, the effective Hubble constant is calculated to be $H_{eff,0}\approx 73.24km.s^{-1}.Mpc^{-1} $ which is compatible with the observational data. We make a comparison between the Hubble tension and the primordial perturbations and the expression of the mass loss rate, chemical potential and central charge needed to describe the Hawking evaporation will be established.

The COSINUS direct dark matter experiment situated at Laboratori Nazionali del Gran Sasso in Italy is set to investigate the nature of the annually modulating signal detected by the DAMA/LIBRA experiment. COSINUS has already demonstrated that sodium iodide crystals can be operated at mK temperature as cryogenic scintillating calorimeters using transition edge sensors, despite the complication of handling a hygroscopic and low melting point material. With results from a new COSINUS prototype, we show that particle discrimination on an event-by-event basis in NaI is feasible using the dual-channel readout of both phonons and scintillation light. The detector was mounted in the novel remoTES design and operated in an above-ground facility for 9.06 g$\cdot$d of exposure. With a 3.7 g NaI crystal, e$^-$/$\gamma$ events could be clearly distinguished from nuclear recoils down to the nuclear recoil energy threshold of 15 keV.