Papers for Tuesday, Jun 20 2023

Deep learning techniques are required for the analysis of synoptic (multi-band and multi-epoch) light curves in massive data of quasars, as expected from the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this follow-up study, we introduced an upgraded version of a conditional neural process (CNP) embedded in a multistep approach for analysis of large data of quasars in the LSST Active Galactic Nuclei Scientific Collaboration data challenge database. We present a case study of a stratified set of the u-band light curves for 283 quasars with very low variability $\sim 0.03$. In this sample, CNP average mean square error is found to be $\sim 5\% $($\sim 0.5$ mag). Interestingly, beside similar level of variability there are indications that individual light curves show flare like features. According to preliminary structure function analysis, these occurrences may be associated to microlensing events with larger time scales $5-10$ years.

The green valley represents the population of galaxies that are transitioning from the actively star-forming blue cloud to the passively evolving red sequence. Studying the properties of the green valley galaxies is crucial for our understanding of the exact mechanisms and processes that drive this transition. The green valley does not have a universally accepted definition. The boundaries of the green valley are often determined by empirical lines that are subjective and vary across studies. We present an unambiguous definition of the green valley in the colour-stellar mass plane using the entropic thresholding. We first divide the galaxy population into the blue cloud and the red sequence based on a colour threshold that minimizes the intra-class variance and maximizes the inter-class variance. Our method splits the region between the mean colours of the blue cloud and the red sequence into three parts by maximizing the total entropy of that region. We repeat our analysis in a number of independent stellar mass bins to define the boundaries of the green valley in the colour-mass diagram. Our method provides a robust and natural definition of the green valley.

Recent results from spectroscopic and astrometric surveys of nearby stars suggest that the stellar disk of our Milky Way (MW) was formed quite early, within the first few billion years of its evolution. Chemo-kinematic signatures of disk formation in cosmological zoom-in simulations appear to be in tension with these data, suggesting that MW-like disk formation is delayed in simulations. We investigate the formation of galactic disks using a representative sample of MW-like galaxies from the cosmological-volume simulation TNG50. We find that on average MW-mass disks indeed form later than the local data suggest. However, their formation time and metallicity exhibit a substantial scatter, such that $\sim$10% of MW-mass galaxies form disks early, similar to the MW. Thus, although the MW is unusual, it is consistent with the overall population of MW-mass disk galaxies. The direct MW analogs assemble most of their mass early, $\gtrsim 10$ Gyr ago, and are not affected by destructive mergers after that. In addition, these galaxies form their disks during the early enrichment stage when the ISM metallicity increases rapidly, with only $\sim$25% of early-forming disks being as metal-poor as the MW was at the onset of disk formation, [Fe/H] $\approx -1.0$. In contrast, most MW-mass galaxies either form disks from already enriched material or experience late destructive mergers that reset the signatures of galactic disk formation to later times and higher metallicities. Finally, we also show that the earlier disk formation leads to more dominant rotationally-supported stellar disks at redshift zero.

So far, stellar population studies have mainly focused on the evolution of single and binary stars. Recent observations show that triple and higher order multiple star systems are common, especially among massive stars. Introducing three-body dynamical effects can influence the evolution of an individual stellar system and hence affect the predicted rates of astrophysical sources that are a product of stellar evolution. We aim to constrain the main evolutionary pathways of massive hierarchical triple star systems and to quantify the effect of the third star on the evolution of the system. We model the massive triple star population by performing simulations of triple star evolution with the TRES code, which combines stellar evolution with secular evolution of triple systems, and explore how robust the predictions of these simulations are under variations of uncertain initial conditions. We find that interactions are common among massive triple stars. The majority of systems (65-77\%) experience a phase of mass transfer in the inner binary, often with an unevolved donor star. This differs significantly from isolated binary evolution, where mass transfer is less frequent (52.3\% instead of 67\% for our fiducial model) and the donors are typically post-main sequence stars. Initial constraints for dynamical stability as well as eccentricity oscillations driven by the third body facilitate the occurrence of interactions, for instance mass transfer. The requirement of dynamical stability at formation places quite stringent constraints on allowed orbital properties, reducing uncertainties in triple evolution that resort from these initial conditions. Ignoring three-body dynamics during evolution of non-interacting triples leads to triple compact-object systems with stronger eccentricity oscillations and thereby likely over-predicts the merger rate of compact objects in such systems.

The high-redshift galaxy UV luminosity function (UVLF) has become essential for understanding the formation and evolution of the first galaxies. Yet, UVLFs only measure galaxy abundances, giving rise to a degeneracy between the mean galaxy luminosity and its stochasticity. Here, we show that upcoming clustering measurements with the James Webb Space Telescope (JWST), as well as with Roman, will be able to break this degeneracy, even at redshifts $z \gtrsim 10$. First, we demonstrate that current Subaru Hyper Suprime-Cam (HSC) measurements of the galaxy bias at $z\sim 4-6$ point to a relatively tight halo-galaxy connection, with low stochasticity. Then, we show that the larger UVLFs observed by JWST at $z\gtrsim 10$ can be explained with either a boosted average UV emission or an enhanced stochasticity. These two models, however, predict different galaxy biases, which are potentially distinguishable in JWST and Roman surveys. Galaxy-clustering measurements, therefore, will provide crucial insights into the connection between the first galaxies and their dark-matter halos, and identify the root cause of the enhanced abundance of $z \gtrsim 10$ galaxies revealed with JWST during its first year of operations.

We explore models of satellite quenching in massive clusters at $z\gtrsim1$ using an MCMC framework, focusing on two primary parameters: $R_{\rm quench}$ (the host-centric radius at which quenching begins) and $\tau_{\rm quench}$ (the timescale upon which a satellite quenches after crossing $R_{\rm quench}$). Our MCMC analysis shows two local maxima in the 1D posterior probability distribution of $R_{\rm quench}$ at approximately $0.25$ and $1.0~R_{\rm{200}}$. Analyzing four distinct solutions in the $\tau_{\rm quench}$-$R_{\rm quench}$ parameter space, nearly all of which yield quiescent fractions consistent with observational data from the GOGREEN survey, we investigate whether these solutions represent distinct quenching pathways and find that they can be separated between "starvation" and "core quenching" scenarios. The starvation pathway is characterized by quenching timescales that are roughly consistent with the total cold gas (H$_{2}$+H${\scriptsize \rm I}$) depletion timescale at intermediate $z$, while core quenching is characterized by satellites with relatively high line-of-sight velocities that quench on short timescales ($\sim 0.25$ Gyr) after reaching the inner region of the cluster ($\lt 0.30~R_{\rm{200}}$). Lastly, we break the degeneracy between these solutions by comparing the observed properties of transition galaxies from the GOGREEN survey. We conclude that only the "starvation" pathway is consistent with the projected phase-space distribution and relative abundance of transition galaxies at $z \sim 1$. However, we acknowledge that ram pressure might contribute as a secondary quenching mechanism.

Large-scale surveys allow us to construct cosmological probes such as galaxy clustering, weak gravitational lensing, the luminosity distance, and cosmic microwave background anisotropies. The gauge-invariant descriptions of these cosmological probes reveal the presence of numerous relativistic effects in the cosmological probes, and they are sensitive (or even divergent) to the long wave-length fluctuations in the initial conditions. In the standard $\Lambda$CDM model, this infrared sensitivity is absent due to subtle cancellations among the relativistic contributions, once the Einstein equation is used. Here we derive the most general conditions for the absence of infrared sensitivity in the cosmological probes without committing to general relativity. We discuss the implications of our results for gravity theories beyond general relativity.

The Baryon Acoustic Oscillations (BAO) feature in the clustering of galaxies or quasars provides a ``standard ruler" for distance measurements in cosmology. In this work, we report a $2\sim3\sigma$ signal of the BAO dip feature in the galaxy density-ellipticity (GI) cross-correlation functions using the spectroscopic sample of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS, combined with the deep DESI Legacy Imaging Surveys for precise galaxy shape measurements. We measure the GI correlation functions and model them using the linear alignment model. We constrain the distance $D_V/r_{\mathrm{d}}$ to redshift $0.57$ to a precision of $3\sim5\%$, depending on the details of modeling. The GI measurement reduces the uncertainty of distance measurement by $\sim10\%$ on top of that derived from the galaxy-galaxy (GG) correlation. More importantly, for future large and deep galaxy surveys, the independent GI measurements can help sort out the systematics in the BAO studies.

Gravitational waves from a first-order cosmological phase transition, at temperatures at the MeV-scale, would arguably be the most exciting explanation of the common red spectrum reported by the NANOGrav collaboration, not the least because this would be direct evidence of physics beyond the standard model. Here we perform a detailed analysis of whether such an interpretation is consistent with constraints on the released energy deriving from big bang nucleosynthesis and the cosmic microwave background. We find that a phase transition in a completely secluded dark sector is strongly disfavoured with respect to the more conventional astrophysical explanation of the putative gravitational wave signal in terms of supermassive black hole binaries. On the other hand, a phase transition in a dark sector that subsequently decays, before the time of neutrino decoupling, remains an intriguing possibility to explain the data. From the model-building perspective, such an option is easily satisfied for couplings with the visible sector that are small enough to evade current collider and astrophysical constraints. The first indication that could eventually corroborate such an interpretation, once the observed common red spectrum is confirmed as a nHz gravitational wave background, could be the spectral tilt of the signal. In fact, the current data already show a very slight preference for a spectrum that is softer than what is expected from the leading astrophysical explanation.

Despite the increasing number of Gravitational Wave (GW) detections, the astrophysical origin of Binary Black Hole (BBH) mergers remains elusive. A promising formation channel for BBHs is inside accretion discs around supermassive black holes, that power Active Galactic Nuclei (AGN). In this paper, we test for the first time the spatial correlation between observed GW events and AGN. To this end, we assemble all sky catalogues with 1,412 (242) AGN with a bolometric luminosity greater than $10^{45.5} {\rm erg\ s}^{-1}$ ($10^{46}\,{\rm erg\,s}^{-1}$) with spectroscopic redshift of $z\leq0.3$ from the Milliquas catalogue, version 7.7b. These AGN are cross-matched with localisation volumes of BBH mergers observed in the same redshift range by the LIGO and Virgo interferometers during their third observing run. We find that the fraction of the detected mergers originated in AGN brighter than $10^{45.5}\,{\rm erg\,s}^{-1}$ ($10^{46}\,{\rm erg\,s}^{-1}$) cannot be higher than $0.74$ ($0.33$) at a 95 per cent credibility level. Our upper limits imply a limited BBH merger production efficiency of the brightest AGN, while most or all GW events may still come from lower luminosity ones. Alternatively, the AGN formation path for merging stellar-mass BBHs may be actually overall subdominant in the local Universe. To our knowledge, ours are the first observational constraints on the fractional contribution of the AGN channel to the observed BBH mergers.

Betelgeuse has fascinated people since they first looked at the sky. Here we present a contemporary summary of the observations and theory that lead to our understanding of Betelgeuse as a massive red supergiant doomed to collapse and explosion. At only ~200 parsecs from Earth, Betelgeuse can be spatially resolved yet uncertainties in its distance remain a critical impediment to deeper understanding. The surface of Betelgeuse is rent with a complex structure as deep convective eddies arise to the surface affecting most of its measured physical properties. Determination of the equatorial rotation velocity is critical since some current estimates indicate that Betelgeuse is rotating anomalously rapidly, a property that cannot be explained by single-star evolutionary models. Betelgeuse is also moving through space at relatively high velocity that indicates that it received a boost, likely via collective interaction with other stars in its birth cluster. A bow shock and other structure in the direction of the star's motion suggest that it has affected the organization of the circumstellar and interstellar medium. Betelgeuse varies in brightness on a variety of time scales with 200, 400 and 2000 days being prominent. Betelgeuse is probable to have been born in a binary system, and the high space velocity and apparent rotation have been related to binary star evolution. One possibility is that Betelgeuse underwent common envelope evolution culminating in a final merger with the core of a massive primary. Such merger models have been invoked to account for the anomalous rotation velocity. Betelgeuse underwent a Great Dimming in 2020 that received widespread attention. Explanations have focused on large cool spots on the surface and the expulsion of a cloud of dust that obscured the surface. We sketch the nature of the explosion to come and discuss perspectives for further research.

Isocurvature perturbations with a blue power spectrum are one of the natural targets for the future large scale structure observations which are probing shorter length scales with greater accuracy. We present a Fisher forecast for the Euclid and MegaMapper (MM) experiments in their ability to detect blue isocurvature perturbations. We construct the theoretical predictions in the EFTofLSS and bias expansion formalisms at quartic order in overdensities which allows us to compute the power spectrum at one loop order and bispectrum at tree level and further include theoretical error at the next to leading order for the covariance determination. We find that Euclid is expected to provide at least a factor of few improvement on the isocurvature spectral amplitude compared to the existing Planck constraints for large spectral indices while MM is expected to provide about 1 to 1.5 order of magnitude im provement for a broad range of spectral indices. We find features that are specific to the blue isocurvature scenario including the leading parametric degeneracy being with the Laplacian bias and a UV sensitive bare sound speed parameter.

Many studies have found that neutron star mergers leave a fraction of the stars' mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site of $r-$ process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars' interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star.

IT Cas is a detached eclipsing binary system containing two F3 V stars in an orbit of period 3.90 d and eccentricity 0.089. Light curves are available from three sectors of observations from the Transiting Exoplanet Survey Satellite (TESS), and extensive radial velocity measurements have been published by Lacy et al (1997). We model these data using the JKTEBOP code to determine the physical properties of the system. We find masses of 1.324 +/- 0.009 and 1.322 +/- 0.008 Msun, and radii of 1.555 +/- 0.004 and 1.551 +/- 0.005 Rsun. The two stars are identical to within the uncertainties, and the depths of the primary and secondary eclipses are also indistinguishable. Using the effective temperature of 6740 +/- 105 K from Lacy et al. (for both stars) gives a distance to the system of 505.5 +/- 8.3 pc, in good agreement with the value of 515.0 +/- 4.4 pc from the Gaia DR3 parallax. The properties of the stars are consistent with theoretical predictions for a solar chemical composition and an age of 2 Gyr. No pulsations are apparent in the TESS photometry.

We acquired ALMA ground state absorption profiles of HCO+ and other molecules toward 33 extragalactic continuum sources seen toward the Galactic anticenter, deriving N(H2) = N(HCO+)/3x10^{-9}. We observed J=1-0 CO emission with the IRAM 30m in directions where HCO+ was newly detected. HCO+ absorption was detected in 28 of 33 new directions and CO emission along 19 of those 28. The 5 sightlines lacking detectable HCO+ have 3 times lower mean EBV and N(DNM). Binned in EBV, N(H2) and N(DNM) are strongly correlated and vary by factors of 50-100 over the observed range EBV~0.05-1 mag, while N(HI) varies by factors of only 2-3. On average N(DNM) and N(H2) are well matched, and detecting HCO+ absorption adds little/no H2 in excess of the previously inferred DNM. There are 5 cases where 2N(H2) < N(DNM)/2 indicates saturation of the HI emission. For sightlines with \WCO > 1 K-\kms the CO-H2 conversion factor N(H2)/\WCO\ = 2-3x10^{20}\pcc/K-\kms is higher than derived from studies of resolved clouds in gamma-rays. Our work sampled primarily atomic gas with a mean H2 fraction ~1/3, but the DNM is almost entirely molecular. CO fulfills its role as an H2 tracer when its emission is strong, but large-scale CO surveys are not sensitive to H2 columns associated with typical values N(DNM) = 2-6x10^{20}\pcc. Lower \XCO\ values from $\gamma$-ray studies arise in part from different definitions and usage. Sightlines with \WCO\ \ge 1 K-\kms\ represent 2/3 of the H2 detected in HCO+ and detecting 90% of the H2 would require detecting CO at levels \WCO\~0.2-0.3 K-\kms For full abstract see the paper

Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the Solar System. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the $\nu_6$ and 3:1 resonances producing $\simeq 30$\% of NEOs with absolute magnitudes $H = 15$ and $\simeq 80$\% of NEOs with $H = 25$. Hence, the large and small NEOs have different orbital distributions. The inferred flux of $H<18$ bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate ($\simeq 2 \times 10^{-4}$ au Myr$^{-1}$ for diameter $D=1$ km and semimajor axis $a=2.5$~au). This implies obliquities $\theta \simeq 0^\circ$ for $a<2.5$~au and $\theta \simeq 180^\circ$ for $a>2.5$~au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model.

We study the ionised and highly ionised gas phases in the Seyfert~2 galaxy IC~5063 by means of VLT/MUSE integral field spectroscopy. Our analysis allowed us to detect a high-ionisation gas outflow traced by the coronal lines [\ion{Fe}{vii}]~$\lambda$6087 and [\ion{Fe}{x}]~$\lambda$6375. Both emissions are found to be extended. The former up to 1.2~kpc and 700~pc NW and SE from the nucleus, respectively. The latter reaches 700~pc NW of the nucleus. This is the first time that [\ion{Fe}{x}] emission is observed at such distances from the central engine in an active galactic nucleus. The [\ion{Fe}{vii}]~$\lambda$6087 emission peaks at the nucleus, with two secondary peaks at the position of the NW and SE radio-lobes. The gas kinematics is complex, with the coronal emission displaying split line profiles along the radio jet and line widths of several hundreds km~s$^{-1}$. Velocity shifts of up to 600~km~s$^{-1}$ in excess of the systemic velocity of the galaxy are found very close to the radio lobes and along the jet propagation. The extended coronal gas is characterised by temperatures reaching 20000~K and electron densities $>10^2$~cm$^{-3}$, with the larger values associated to the regions of larger turbulence, likely due to the passage of the radio jet. This hypothesis is supported by photoionisation models that combine the effects of the central engine and shocks. Our work highlights the strong relationship between extended coronal emission and the radio jet, with the former suitably tracing the latter, which in the case of IC~5063, propagates very close to the galaxy disc.

We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameterised by the combined mass density slope $\gamma$ (baryons $+$ dark matter), which is the logarithmic change of density with radius. We use a MAGPI sample of 28 galaxies from low-to-mid density environments and compare to density slopes derived from galaxies in the high density Frontier Fields clusters in the redshift range $0.29 <z < 0.55$, corresponding to a lookback time of $\sim 5$ Gyr. We find a median density slope of $\gamma = -2.22 \pm 0.05$ for the MAGPI sample, which is significantly steeper than the Frontier Fields median slope ($\gamma = -2.01 \pm 0.04$), implying the cluster galaxies are less centrally concentrated in their mass distribution than MAGPI galaxies. We also compare to the distribution of density slopes from galaxies in Atlas3D at $z \sim 0$, because the sample probes a similar environmental range as MAGPI. The Atlas3D median total slope is $\gamma = -2.25 \pm 0.02$, consistent with the MAGPI median. Our results indicate environment plays a role in the internal mass distribution of galaxies, with no evolution of the slope in the last 3-4 Gyr. These results are in agreement with the predictions of cosmological simulations.

We analyze the spectrum of gravitational waves generated by the induced spectrum of tensor fluctuation during warm natural inflation. In our previous work it has been demonstrated that an epoch of warm natural inflation can lead to cosmologically relevant dark matter production in the form of primordial black holes. Here we show that models which solve the dark-matter production also produce a contribution to the cosmic gravitational wave background that satisfies current constraints from pulsar timing and big bang nucleosynthesis. More importantly, this gravitational wave background may be observable in the next generation of space-based and ground-based gravitational wave interferometers.

So far no designated mission to either of the two ice giants, Uranus and Neptune, exists. Almost all of our gathered information on these planets comes from remote sensing. In recent years, NASA and ESA have started planning for future mission to Uranus and Neptune, with both agencies focusing their attention on orbiters and atmospheric probes. Whereas information provided by remote sensing is undoubtedly highly valuable, remote sensing of planetary atmospheres also presents some shortcomings, most of which can be overcome by mass spectrometers. In most studies presented to date a mass spectrometer experiment is thus a favored science instrument for in situ composition measurements on an atmospheric probe. Mass spectrometric measurements can provide unique scientific data, i.e., sensitive and quantitative measurements of the chemical composition of the atmosphere, including isotopic, elemental, and molecular abundances. In this review paper we present the technical aspects of mass spectrometry relevant to atmospheric probes. This includes the individual components that make up mass spectrometers and possible implementation choices for each of these components. We then give an overview of mass spectrometers that were sent to space with the intent of probing planetary atmospheres, and discuss three instruments, the heritage of which is especially relevant to Uranus and Neptune probes, in detail. The main part of this paper presents the current state-of-art in mass spectrometry intended for atmospheric probe. Finally, we present a possible descent probe implementation in detail, including measurement phases and associated expected accuracies for selected species.

Type IIn supernovae occur when stellar explosions are surrounded by dense hydrogen-rich circumstellar matter. The dense circumstellar matter is likely formed by extreme mass loss from their progenitors shortly before they explode. The nature of Type IIn supernova progenitors and the mass-loss mechanism forming the dense circumstellar matter are still unknown. In this work, we investigate if there are any correlations between Type IIn supernova properties and their local environments. We use Type IIn supernovae with well-observed light-curves and host-galaxy integral field spectroscopic data so that we can estimate both supernova and environmental properties. We find that Type IIn supernovae with a higher peak luminosity tend to occur in environments with lower metallicity and/or younger stellar populations. The circumstellar matter density around Type IIn supernovae is not significantly correlated with metallicity, so the mass-loss mechanism forming the dense circumstellar matter around Type IIn supernovae might be insensitive to metallicity.

Galaxy clusters, being the most massive objects in the Universe, exhibit the strongest alignment with the large-scale structure. However, mis-identification of members due to projection effects from the large scale structure can occur. We studied the impact of projection effects on the measurement of the intrinsic alignment of galaxy clusters, using galaxy cluster mock catalogs. Our findings showed that projection effects result in a decrease of the large scale intrinsic alignment signal of the cluster and produce a bump at $r_p\sim 1h^{-1}/Mpc$, most likely due to interlopers and missed member galaxies. This decrease in signal explains the observed similar alignment strength between bright central galaxies and clusters in the SDSS redMaPPer cluster catalog. The projection effect and cluster intrinsic alignment signal are coupled, with clusters having lower fractions of missing members or having higher fraction of interlopers exhibiting higher alignment signals in their projected shapes. We aim to use these findings to determine the impact of projection effects on galaxy cluster cosmology in future studies.

Regular black holes are singularity-free black hole spacetimes proposed to solve the problem of the presence of spacetime singularities that plagues the black holes of general relativity and most theories of gravity. In this work, we consider the regular black holes recently proposed by Mazza, Franzin \& Liberati and we extend previous studies to get a more stringent observational constraint on the regularization parameter $l$. We study simultaneous observations of \textit{NuSTAR} and \textit{Swift} of the Galactic black hole in GX~339--4 during its outburst in 2015. The quality of the \textit{NuSTAR} data is exceptionally good and the spectrum of the source presents both a strong thermal component and prominent relativistically blurred reflection features. This permits us to measure the regularization parameter $l$ from the simultaneous analysis of the thermal spectrum and the reflection features. From our analysis, we find the constraint $l/M < 0.39$ (90\% CL), which is stronger than previous constraints inferred with X-ray and gravitational wave data.

An algorithm for robust initial orbit determination (IOD) under perturbed orbital dynamics is presented. By leveraging map inversion techniques defined in the algebra of Taylor polynomials, this tool is capable of not only returning an highly accurate solution to the IOD problem, but also estimating a range of validity for the aforementioned solution in which the true orbit state should lie. Automatic domain splitting is then used on top of the IOD routines to ensure the local truncation error introduced by a polynomial representation of the state estimate remains below a predefined threshold to meet the specified requirements in accuracy. The algorithm is adapted to three types of ground based sensors, namely range radars, Doppler-only radars and optical telescopes by taking into account their different constraints in terms of available measurements and sensor noise. Its improved performance with respect to a Keplerian based IOD solution is finally demonstrated with large scale numerical simulations over a subset of tracked objects in low Earth orbit.

Constraining the vertical and radial structure of debris discs is crucial to understanding their formation, evolution and dynamics. To measure both the radial and vertical structure, a disc must be sufficiently inclined. However, if a disc is too close to edge-on, deprojecting its emission becomes non-trivial. In this paper we show how Frankenstein, a non-parametric tool to extract the radial brightness profile of circumstellar discs, can be used to deproject their emission at any inclination as long as they are optically thin and axisymmetric. Furthermore, we extend Frankenstein to account for the vertical thickness of an optically thin disc ($H(r)$) and show how it can be constrained by sampling its posterior probability distribution and assuming a functional form (e.g. constant $h=H/r$), while fitting the radial profile non-parametrically. We use this new method to determine the radial and vertical structure of 16 highly inclined debris discs observed by ALMA. We find a wide range of vertical aspect ratios, $h$, ranging from $0.020\pm0.002$ (AU Mic) to $0.20\pm0.03$ (HD 110058), which are consistent with parametric models. We find a tentative correlation between $h$ and the disc fractional width, as expected if wide discs were more stirred. Assuming discs are self-stirred, the thinnest discs would require the presence of at least 500 km-sized planetesimals. The thickest discs would likely require the presence of planets. We also recover previously inferred and new radial structures, including a potential gap in the radial distribution of HD 61005. Finally, our new extension of Frankenstein also allows constraining how $h$ varies as a function of radius, which we test on 49 Ceti, finding that $h$ is consistent with being constant.

We report X-ray observations of the newly discovered pulsating white dwarf eRASSU J191213.9-441044 with Spectrum Roentgen Gamma and eROSITA (SRG/eROSITA) and XMM-Newton. The new source was discovered during the first eROSITA all-sky survey at a flux level of fX (0.2 - 2.3 keV) = 3.3 e-13 erg cm-2 s-1 and found to be spatially coincident with a G = 17.1 stellar Gaia-source at a distance of 237 pc. The flux dropped to about fX = 1 e-13 erg cm-2 s-1 during the three following eROSITA all-sky surveys and remained at this lower level during dedicated XMM-Newton observations performed in September 2022. With XMM-Newton, pulsations with a period of 319 s were found at X-ray and ultraviolet wavelengths occurring simultaneously in time, thus confirming the nature of eRASSU J191213.9-441044 as the second white-dwarf pulsar. The X-ray and UV-pulses correspond to broad optical pulses. Narrow optical pulses that occurred occasionally during simultaneous XMM-Newton/ULTRACAM observations have no X-ray counterpart. The orbital variability of the X-ray signal with a roughly sinusoidal shape was observed with a pulsed fraction of ~28% and maximum emission at orbital phase ~0.25. The ultraviolet light curve peaks at around binary phase 0.45. The X-ray spectrum can be described with the sum of a power law spectrum and a thermal component with a mean X-ray luminosity of Lx(0.2-10 keV) = 1.4 e30 erg s-1. The spectral and variability properties could indicate some residual accretion, in contrast to the case of the prototypical object AR Sco.

NGC 5024 (M53) is the most distant globular cluster (GC) with known pulsars. In this study, we report the discovery of a new binary millisecond pulsar PSR J1312+1810E (M53E) and present the new timing solutions for M53B to M53E, based on 22 observations from the Five-hundred-meter Aperture Spherical radio Telescope (FAST).These discoveries and timing work benefit from FAST's high sensitivity. We find that M53C is the only isolated millisecond pulsar known in this distant globular cluster, with a spin period of 12.53 ms and spin period derivative of $5.26 \times 10^{-20} \, \rm s \; s^{-1}$. Our results reveal the orbital periods of 47.7, 5.8, and 2.4 days for M53B, D, and E, respectively. The companions, with a mass of 0.25, 0.27, and 0.18 ${\rm M}_\odot$, respectively, are likely to be white dwarf stars; if they are extended objects, they don't eclipse the pulsars. We find no X-ray counterparts for these millisecond pulsars in archival $Chandra$ images in the band of 0.3-8 keV. The characteristics of this pulsar population are similar to the population of millisecond pulsars in the Galactic disk, as expected from the low stellar density of M53.

The third realization of the International Celestial Reference Frame (ICRF3) was adopted in August 2018 and includes positions of extragalactic objects at three frequencies: 8.4 GHz, 24 GHz, and 32 GHz. In this paper, we present celestial reference frames estimated from Very Long Baseline Interferometry measurements at K-band (24 GHz) including data until June 2022. The data set starts in May 2002 and currently consists of more than 120 24h observing sessions performed over the past 20 years. Since the publication of ICRF3, the additional observations of the sources during the last four years allow maintenance of the celestial reference frame and more than 200 additional radio sources ensure an expansion of the frame. A study of the presented solutions is carried out helping us to understand systematic differences between the astrometric catalogs and moving us towards a better next ICRF solution. We compare K-band solutions (VIE-K-2022b and USNO-K-2022July05) computed by two analysts with two independent software packages (VieVS and Calc/Solve) and describe the differences in the solution strategy. We assess the systematic differences using vector spherical harmonics and describe the reasons for the most prominent ones.

Fundamental stellar parameters such as mass and radius are some of the most important building blocks in astronomy, both when it comes to understanding the star itself and when deriving the properties of any exoplanet(s) they may host. Asteroseismology of solar-like oscillations allows us to determine these parameters with high precision. We investigate the solar-like oscillations of the red-giant-branch star $\gamma$ Cep A, which harbours a giant planet on a wide orbit. We did this by utilising both ground-based radial velocities from the SONG network and space-borne photometry from the NASA TESS mission. From the radial velocities and photometric observations, we created a combined power spectrum, which we used in an asteroseismic analysis to extract individual frequencies. We clearly identify several radial and quadrupole modes as well as multiple mixed, dipole modes. We used these frequencies along with spectroscopic and astrometric constraints to model the star, and we find a mass of $1.27^{+0.05}_{-0.07}$ M$_\odot$, a radius of $4.74^{+0.07}_{-0.08}$ R$_\odot$, and an age of $5.7^{+0.8}_{-0.9}$ Gyr. We then used the mass of $\gamma$ Cep A and our SONG radial velocities to derive masses for $\gamma$ Cep B and $\gamma$ Cep Ab of $0.328^{+0.009}_{-0.012}$ M$_\odot$ and $6.6^{+2.3}_{-2.8}$ M$_{\rm Jup}$, respectively.

SN 2023emq is a fast-evolving transient initially classified as a rare Type Icn supernova (SN), interacting with a H- and He-free circumstellar medium (CSM) around maximum light. Subsequent spectroscopy revealed the unambiguous emergence of narrow He lines, confidently placing SN 2023emq in the more common Type Ibn class. Photometrically SN 2023emq has several uncommon properties regardless of its class, including its extreme initial decay (faster than > 90% of Ibn/Icn SNe) and sharp transition in the decline rate from 0.18 mag/d to 0.05 mag/d at +20 d. The bolometric light curve can be modelled as CSM interaction with 0.31M_Sun of ejecta and 0.13M_Sun of CSM, with 0.009M_Sun of nickel, as expected of fast interacting SN. Furthermore, broad-band polarimetry at +8.7 days (P = 0.55+/-0.30%) is consistent with high spherical symmetry. A discovery of a transitional Icn/Ibn SN would be unprecedented and would give valuable insights into the nature of mass loss suffered by the progenitor just before death, but we favour an interpretation that the emission lines in the classification spectrum are flash ionisation features commonly seen in young SNe the first days after the explosion. However, one of the features (5700 {\AA}) is significantly more prominent in SN 2023emq than in the few flash-ionised Type Ibn SNe and in that regard the SN is more similar to Icn SNe possibly implying continuum of properties between the two classes.

We study active galactic nucleus (AGN) feedback in nearby (z<0.35) galaxy clusters from the Planck Sunyaev-Zeldovich (SZ) sample using Chandra observations. This nearly unbiased mass-selected sample includes both relaxed and disturbed clusters and may reflect the entire AGN feedback cycle. We find that relaxed clusters better follow the one-to-one relation of cavity power versus cooling luminosity, while disturbed clusters display higher cavity power for a given cooling luminosity, likely reflecting a difference in cooling and feedback efficiency. Disturbed clusters are also found to contain asymmetric cavities when compared to relaxed clusters, hinting toward the influence of the intracluster medium (ICM) weather on the distribution and morphology of the cavities. Disturbed clusters do not have fewer cavities than relaxed clusters, suggesting that cavities are difficult to disrupt. Thus, multiple cavities are a natural outcome of recurrent AGN outbursts. As in previous studies, we confirm that clusters with short central cooling times, tcool, and low central entropy values, K0, contain warm ionized (10000 K) or cold molecular (<100 K) gas, consistent with ICM cooling and a precipitation/chaotic cold accretion (CCA) scenario. We analyzed archival MUSE observations that are available for 18 clusters. In 11/18 of the cases, the projected optical line emission filaments appear to be located beneath or around the cavity rims, indicating that AGN feedback plays an important role in forming the warm filaments by likely enhancing turbulence or uplift. In the remaining cases (7/18), the clusters either lack cavities or their association of filaments with cavities is vague, suggesting alternative turbulence-driven mechanisms (sloshing/mergers) or physical time delays are involved.

We use new observations from the Canada-France-Hawaii Telescope to study the white dwarf cooling sequence of the open cluster M37, a cluster that displays an extended main sequence turn-off and, according to a recent photometric analysis, also a spread of initial chemical composition. By taking advantage of a first epoch collected in 1999 with the same telescope, we have been able to calculate proper motions for sources as faint as g ~ 26 (about ~ 6 magnitudes fainter than the Gaia limit), allowing us to separate cluster members from field stars. This has enabled us to isolate a sample of the white dwarf population of M37, reaching the end of the cooling sequence (at g ~ 23.5). The here-derived atlas and calibrated catalogue of the sources in the field of view is publicly released as supplementary on-line material. Finally, we present an exhaustive comparison of the white dwarf luminosity function with theoretical models, which has allowed us to exclude the age-spread scenario as the main responsible for the extended turnoff seen in the cluster colour-magnitude-diagram.

We discuss the applicability of the Melnikov and Landau-Teller theories in obtaining semi-analytical estimates of the speed of chaotic diffusion in systems driven by the separatrix-like stochastic layers of a resonance belonging to the `second fundamental model' (SFM)\cite{henrard1983second}. Stemming from the analytic solution for the SFM in terms of Weierstrass elliptic functions, we introduce stochastic Melnikov and Landau-Teller models allowing to locally approximate chaotic diffusion as a sequence of uncorrelated `jumps' observed in the time series yielding the slow evolution of an ensemble of trajectories in the space of the adiabatic actions of the system. Such jumps occur in steps of one per homoclinic loop. We show how a semi-analytical determination of the probability distribution of the size of the jumps can be arrived at by the Melnikov and Landau-Teller approximate theories. Computing also the mean time required per homoclinic loop, we arrive at estimates of the chaotic diffusion coefficient in such systems. As a concrete example, we refer to the long-term diffusion of a small object (e.g. Earth navigation satellite or space debris) within the chaotic layers of the so-called $2g+h$ lunisolar resonance, which is of the SFM type. After a suitable normal form reduction of the Hamiltonian, we compute estimates of the speed of diffusion of these objects, which compare well with the results of numerical experiments.

We search for Milky Way-like galaxies among a sample of approximately 500 galaxies. The characteristics we considered of the candidate galaxies are the following: stellar mass M_star, optical radius R_25, rotation velocity V_rot, central oxygen abundance (O/H)_0, and abundance at the optical radius (O/H)_R25. If the values of R_25 and M_star of the galaxy were close to that of the Milky Way, then the galaxy was referred to as a structural Milky Way analogue (sMWA). The oxygen abundance at a given radius of a galaxy is defined by the evolution of that region, and we then assumed that the similarity of (O/H)_0 and (O/H)_R25 in two galaxies suggests a similarity in their evolution. If the values of (O/H)_0 and (O/H)_R25 in the galaxy were close to that of the Milky Way, then the galaxy was referred to as an evolutionary Milky Way analogue (eMWA). If the galaxy was simultaneously an eMWA and sMWA, then the galaxy was considered a Milky Way twin. We find that the position of the Milky Way on the (O/H)_0 - (O/H)_R25 diagram shows a large deviation from the general trend in the sense that the (O/H)_R25 in the Milky Way is appreciably lower than in other galaxies of similar (O/H)_0. This feature of the Milky Way evidences that its (chemical) evolution is not typical. We identify four galaxies (NGC~3521, NGC~4651, NGC~2903, and MaNGA galaxy M-8341-09101) that are simultaneously sMWA and eMWA and can therefore be considered as Milky Way twins. In previous studies, Milky Way-like galaxies were selected using structural and morphological characteristics, that is, sMWAs were selected. We find that the abundances at the centre and at the optical radius (evolutionary characteristics) provide a stricter criterion for selecting real Milky Way twins

Observations of the effect of microlensing in gravitationally lensed quasars can be used to study the structure of active galactic nuclei on distance scales down to the sizes of the supermassive black holes powering source activity. We search for the microlensing in the signal from a gravitationally lensed blazar B0218+357 in very-high-energy gamma-ray band. We combine observations of a bright flare of the source in 2014 with Fermi/LAT and MAGIC telescopes in 0.1-100 GeV energy range. Using the time-delayed leading and trailing signals from two gravitationally lensed images of the source, we measure magnification factor at the moment of the flare. We use the scaling of the maximal magnification factor with the source size to constrain the size of gamma-ray emission region in the entire 0.1-100 GeV energy range. The magnification factor in the very-high-energy band derived from a comparison of Fermi/LAT and MAGIC data is $\mu_{VHE} = 36^{+40}_{-26}$, substantially larger than that in the radio band. This suggests that one of the source images is strongly affected by microlensing at the moment of the flare. Assuming that the microlensing is produced by a stellar mass object in the lens galaxy, we constrain the size of the emission region in the $E>100$ GeV band to be $\mathrm{R_{VHE} = 2.2^{+27}_{-1.7} \times 10^{13}~cm}$. We note that the spectrum of the microlensed source was unusually hard at the moment of the flare and speculate that this hardening may be due to the energy dependent microlensing effect. This interpretation suggests that the source size decreases with energy in 0.1-100 GeV energy range studied.

We investigate the Stochastic Gravitational Wave Background (SGWB) produced by merging binary black holes (BBHs) and binary neutron stars (BNSs) in the frequency ranges of LIGO/Virgo/Kagra and LISA. We develop three analytical models, that are calibrated to the measured local merger rates, and complement them with three population synthesis models based on the COSMIC code. We discuss the uncertainties, focusing on the impact of the BBH mass distribution, the effect of the metallicity of the progenitor stars and the time delay distribution between star formation and compact binary merger. We also explore the effect of uncertainties in binary stellar evolution on the background. For BBHs, our analytical models predict $\Omega_{GW}$ in the range $[4.10^{-10}-1.10^{-9}]$ (25 Hz) and $[1.10^{-12}-4.10^{-12}]$ (3 mHz), and between $[2.10^{-10}-2.10^{-9}]$ (25 Hz) and $[7.10^{-13}- 7.10^{-12}]$ (3 mHz) for our population synthesis models. This background is unlikely to be detected during the LIGO/Virgo/Kagra O4 run, but could be detectable with LISA. We predict about 10 BBH and no BNS mergers that could be individually detectable by LISA for a period of observation of 4 years. Our study provides new insights into the population of compact binaries and the main sources of uncertainty in the astrophysical SGWB.

We present $\gamma$-ray optical-depth calculations from a recently published extragalactic background light (EBL) model built from multiwavelength galaxy data from the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (HST/CANDELS). CANDELS gathers one of the deepest and most complete observations of stellar and dust emissions in galaxies. This model resulted in a robust derivation of the evolving EBL spectral energy distribution up to $z\sim 6$, including the far-infrared peak. Therefore, the optical depths derived from this model will be useful for determining the attenuation of $\gamma$-ray photons coming from high-redshift sources, such as those detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, and for multi-TeV photons that will be detected from nearby sources by the future Cherenkov Telescope Array. From these newly calculated optical depths, we derive the cosmic $\gamma$-ray horizon and also measure the expansion rate and matter content of the Universe including an assessment of the impact of the EBL uncertainties. We find $H_{0}=61.9$ $^{+2.9}_{-2.4}$ km s$^{-1}$ Mpc$^{-1}$ when fixing $\Omega_{m}=0.32$, and $H_{0}=65.6$ $^{+5.6}_{-5.0}$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_{m}=0.19\pm 0.07$, when exploring these two parameters simultaneously.

Galaxy evolution is an important topic, and our physical understanding must be complete to establish a correct picture. This includes a thorough treatment of feedback. The effects of thermal-mechanical and radiative feedback have been widely considered, however cosmic rays (CRs) are also powerful energy carriers in galactic ecosystems. Resolving the capability of CRs to operate as a feedback agent is therefore essential to advance our understanding of the processes regulating galaxies. The effects of CRs are yet to be fully understood, and their complex multi-channel feedback mechanisms operating across the hierarchy of galaxy structures pose a significant technical challenge. This review examines the role of CRs in galaxies, from the scale of molecular clouds to the circum-galactic medium. An overview of their interaction processes, their implications for galaxy evolution, and their observable signatures is provided and their capability to modify the thermal and hydrodynamic configuration of galactic ecosystems is discussed. We present recent advancements in our understanding of CR processes and interpretation of their signatures, and highlight where technical challenges and unresolved questions persist. We discuss how these may be addressed with upcoming opportunities.

We present the results of current observations of the young compact cluster of massive stars Westerlund 2 with the Mikhail Pavlinsky ART-XC telescope aboard the Spectrum-Roentgen-Gamma (SRG) observatory which we analysed together with the archival Chandra data. In general, Westerlund 2 was detected over the whole electromagnetic spectrum including high-energy gamma rays, which revealed a cosmic ray acceleration in this object to the energies up to tens of TeV. The detection of Westerlund 2 with ART-XC allowed us to perform a joint spectral analysis together with the high resolution Chandra observations of the diffuse emission from a few selected regions in the vicinity of the Westerlund 2 core in the 0.4 - 20 keV range. To fit the Westerlund 2 X-ray spectrum above a few keV one needs either a non-thermal power-law emission component, or a hot plasma with temperatures $\sim$ 5 keV. Our magnetohydrodynamic modeling of the plasma flows in Westerlund 2 shows substantially lower electron temperatures in the system and thus the presence of the non-thermal component is certainly preferable. A kinetic model of the particle acceleration demonstrated that the non-thermal component may originate from the synchrotron radiation of multi-TeV electrons and positrons produced in Westerlund 2 in accordance with the TeV photons detection from the source. However, the inverse Compton radiation of mildly relativistic electrons could explain the non-thermal emission as well.

A class of inflationary scenarios for primordial black hole (PBH) formation include a small barrier in the slope of the potential. There, the inflaton slows down, generating an enhancement of primordial perturbations. Moreoever, the background solution overcomes the barrier at a very low speed, and large backward quantum fluctuations can prevent certain regions from overshooting the barrier. This leads to localized bubbles where the field remains trapped behind the barrier. In such models, therefore, we have two distinct channels for PBH production: the standard adiabatic density perturbation channel and the bubble channel. Here, we perform numerical simulations of bubble formation, addressing the issues of initial conditions, critical amplitude and bubble expansion. Further, we explore the scaling behaviour of the co-moving size of bubbles with the initial amplitude of the field fluctuation. We find that for small to moderate non-Gaussianity $f_{\rm NL}\lesssim 2.6$, the threshold for the formation of vacuum bubbles agrees with previous analytical estimates arXiv:1908.11357 to $5\%$ accuracy or so. We also show that the mass distribution for the two channels is different, leading to a slightly broader range of PBH masses. The bubble channel is subdominant for small $f_{\rm NL}$, and becomes dominant for $f_{\rm NL}\gtrsim 2.6$. We find that the mass of PBHs in the bubble channel is determined by an adiabatic overdensity surrounding the bubble at the end of inflation. Remarkably, the profile of this overdensity turns out to be of type-II. This represents a first clear example showing that overdensities of type-II can be dominant in comparison with the standard type-I. We also comment on the fact that in models with local type non-Gaussianity (such as the one considered here), the occurrence of alternative channels can easily be inferred from unitarity considerations.

The alignment of striated intensity structures in thin neutral hydrogen (HI) spectroscopic channels with Galactic magnetic fields has been observed. However, the origin and nature of these striations are still debatable. Some studies suggest that the striations result solely from real cold-density filaments without considering the role of turbulent velocity fields, i.e., the velocity caustics effect in shaping the channel's intensity distribution. To determine the relative contribution of density and velocity in forming the striations in channel maps, we analyze synthetic observations of channel maps obtained with simulations that represent realistic magnetized multi-phase HI. We vary the thickness of the channel maps and apply the Velocity Decomposition Algorithm to separate the velocity and density contributions. In parallel, we analyze GALFA HI observations and compare the results. Our analysis shows that the thin channels are dominated by velocity contribution, and velocity caustics mainly generate the HI striations. We show that velocity caustics can cause a correlation between unsharp-masked HI structures and far-infrared emission. We demonstrate that the linear HI fibers revealed by the Rolling Hough Transform (RHT) in thin velocity channels originate from velocity caustics. As the thickness of channel maps increases, the relative contribution of density to fluctuations in channel maps also increases. As a result, more RHT-detected fibers tend to be perpendicular to the magnetic field. Conversely, the alignment with the magnetic field is the most prominent in thin channels. We conclude that similar to the Velocity Channel Gradients (VChGs) approach, RHT traces magnetic fields through the analysis of velocity caustics in thin channel maps.

Magnetic fields and their dynamical interplay with matter in galaxy clusters contribute to the physical properties and evolution of the intracluster medium. However, the current understanding of the origin and properties of cluster magnetic fields is still limited by observational challenges. In this article, we map the magnetic fields at hundreds-kpc scales of five clusters RXC J1314.4 -2515, Abell 2345, Abell 3376, MCXC J0352.4 -7401, and El Gordo using the innovative synchrotron intensity gradient technique in conjunction with high-resolution radio observations from JVLA and MeerKAT. We demonstrate that magnetic field orientation of radio relics derived from synchrotron intensity gradients is in very good agreement with that obtained with synchrotron polarization. Most important, synchrotron intensity gradients is not limited by Faraday depolarization in the cluster central regions and allows us to map magnetic fields in the radio halos of RXC J1314.4 -2515 and El Gordo. We find that magnetic fields in radio halos exihibit a preferential direction along the major merger axis and show turbulent structures at higher angular resolution. Results are consistent with expectations from numerical simulations which predict turbulent magnetic fields in cluster mergers that are stirred and amplified by matter motions.

We calculate the scalar-induced gravitational wave energy density in the theory of Ghost Inflation, assuming scale invariance and taking into account both the power spectrum- and trispectrum-induced contributions. For the latter we consider the leading cubic and quartic couplings of the comoving curvature perturbation in addition to two parity-violating quartic operators. In the parity-even case, we find the relative importance of the trispectrum-induced signal to be suppressed by the requirement of perturbativity, strengthening a no-go theorem recently put forth. The parity-odd signal, even though also bound to be small, is non-degenerate with the Gaussian contribution and may in principle be comparable to the parity-even non-Gaussian part, thus potentially serving as a probe of the Ghost Inflation scenario and of parity violating physics during inflation.

A dynamical system analysis is performed for a model of dissipative quintessential inflation realizing warm inflation at early primordial times and dissipative interations in the dark sector at late times. The construction makes use of a generalized exponential potential realizing both phases of accelerated expansion. A focus is given on the behavior of the dynamical system at late times and the analysis is exemplified by both analytical and numerical results. The results obtained demonstrate the viability of the model as a quintessential inflation model and in which stable solutions can be obtained.

We investigate the production of non-thermal dark matter (DM), $\chi$, during post-inflationary reheating era. For the inflation, we consider two slow roll single field inflationary scenarios - generalized version of Hilltop (GH) inflation, and Coleman-Weinberg (NM-CW) inflation with non-minimal coupling between inflaton and curvature scalar. Using a set of benchmark values that comply with the current constraints from Cosmic Microwave Background Radiation (CMBR) data for each inflationary model, we explored the parameter space involving mass of dark matter particles, $m_\chi$, and coupling between inflaton and $\chi$, $y_\chi$. For these benchmarks, we find that tensor-to-scalar ratio $r$ can be as small as $2.69\times 10^{-6}$ for GH and $3\times 10^{-3}$ for NM-CW inflation, both well inside $1-\sigma$ contour on scalar spectral index versus $r$ plane from Planck 2018+BICEP 3+Keck Array 2018 dataset, and testable by future cosmic microwave background (CMB) observations e.g. Simons Observatory. For the production of $\chi$ from the inflaton decay satisfying CMB and other cosmological bounds and successfully explaining total cold dark matter density of the present universe, we find that $y_\chi$ should be within this range ${\cal O}\left(10^{-8}\right) \gtrsim y_\chi\gtrsim {\cal O}\left(10^{-17}\right)$ for both inflationary scenarios. We also show that, even for the same inflationary scenario, the allowed parameter space on reheating temperature versus $m_\chi$ plane alters with inflationary parameters including scalar spectral index, tensor-to-scalar ratio, and energy scale of inflation.

In the string axiverse scenario, light primordial black holes may spin up due to the Hawking emission of a large number of light (sub-MeV) axions. We show that this may trigger superradiant instabilities associated with a heavier axion during the black holes' evolution, and study the coupled dynamics of superradiance and evaporation. We find, in particular, that the present black hole mass-spin distribution should follow the superradiance threshold condition for black hole masses below the value at which the superradiant cloud forms, for a given heavy axion mass. Furthermore, we show that the decay of the heavy axions within the superradiant cloud into photon pairs may lead to a distinctive line in the black hole's emission spectrum, superimposed on its electromagnetic Hawking emission.

We delve deeper into the potential composition of dark matter as stable scalar glueballs from a confining dark $SU(N)$ gauge theory, focusing on $N=\{3,4,5\}$. To predict the relic abundance of glueballs for the various gauge groups and scenarios of thermalization of the dark gluon gas, we employ a thermal effective theory that accounts for the strong-coupling dynamics in agreement with lattice simulations. We compare our methodology with previous works and find that our approach is more comprehensive and reliable. The results are encouraging and show that glueballs can account for the totality of dark matter in many unconstrained scenarios with a phase transition scale $20~{\rm MeV}\lesssim\Lambda\lesssim10^{10}~{\rm GeV}$, thus opening the possibility of exciting future studies.

We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon-induced fast neutrons and cosmogenic isotopes. A fiducial volume cut, as well as the pulse shape discrimination and the muon veto are applied to suppress the above backgrounds. It is shown that JUNO sensitivity to the thermally averaged dark matter annihilation rate in 10 years of exposure would be significantly better than the present-day best limit set by Super-Kamiokande and would be comparable to that expected by Hyper-Kamiokande.

Glitches, spin-up events in neutron stars, are of prime interest as they reveal properties of nuclear matter at subnuclear densities. We numerically investigate the glitch mechanism using analogies between neutron stars and magnetic dipolar gases in the supersolid phase. In rotating neutron stars, glitches are believed to occur when many superfluid vortices unpin from the interior, transferring angular momentum to the stellar surface. In the supersolid analogy, we show that a glitch happens when vortices pinned in the low-density inter-droplet region abruptly unpin. These supersolid glitches show remarkable parallels with neutron star glitches: they are characterized by a rapid spin-up followed by a long post-glitch spin-down due to relaxation towards a steady state. Dipolar supersolids offer an unprecedented possibility to test both the vortex and crystal dynamics during a glitch. Here, we explore the glitch dependence on the supersolid quality, finding strong suppression at the supersolid-to-solid transition. This provides a tool to study glitches originating from different radial depths of a neutron star. Benchmarking our theory against neutron star observations, our work will open a new avenue for the quantum simulation of stellar objects from Earth.

Human landing, exploration and settlement on Mars will require local compute resources at the Mars edge. Landing such resources on Mars is an expensive endeavor. Instead, in this paper we lay out how concepts from low-Earth orbit edge computing may be applied to Mars edge computing. This could lower launching costs of compute resources for Mars while also providing Mars-wide networking and compute coverage. We propose a possible Mars compute constellation, discuss applications, analyze feasibility, and raise research questions for future work.

We consider a type of k-inflation under the Hamilton-Jacobi approach. We calculate various observables such as the scalar power spectrum, the tensor-to-scalar ratio, the scalar spectra index for the case where the Hubble parameter is a power-law function of k-field. The model's parameters are constrained with Planck data and the concrete form of the potential is presented. The results show that the model can be in good agreement with observations.