Papers for Monday, Feb 12 2024

The transiting dust clouds that orbit the white dwarf J0328-1219 are devoid of small particles (< 0.1 micron). Observations show that fade amount doesn't depend on wavelength. This finding resembles a similar observation for white dwarf WD 1145+017, but the explanations for an absence of small particles in the two white dwarf systems may differ due to their different distances from the star.

Using JWST observations of a primary transit and two secondary eclipses for GJ 1214b, we determine an eccentricity that is more precise than a decade of HARPS data, which enables us to measure the stellar density to 2.62%. Coupled with a prior on the stellar mass from a dynamically calibrated K-$M_*$ relation, we determine $R_*$ to 1.13% -- 3 times more precise than any other published analysis of this system. Then, using the bolometric flux from an SED model, we determine $T_{\rm eff}$ to 1.39% -- 40% more precise than systematic floors from spectroscopy. Within the global model, these also improve the planetary radius and insolation. This is a proof of concept for a new method to determine accurate $R_*$ and $T_{\rm eff}$ to a precision currently achieved for only a small number of low-mass stars. By applying our method to all high signal-to-noise planetary transits and occultations, we can expand the sample of precisely measured stars without assuming tidal circularization and calibrate new relations to improve our understanding of all low-mass stars.

We study "buried" active galactic nuclei (AGNs) almost fully covered by circumnuclear material in ultra-/luminous infrared galaxies (U/LIRGs), which show weak ionized lines from narrow line regions. Employing an indicator of [O IV] 25.89-um or [Ne V] 14.32-um line to 12-um AGN luminosity ratio, we find 17 buried AGN candidates that are [O IV]-weak ($L_{\rm [O\,IV]}$/$L_{\rm 12,AGN} \leq -$3.0) or [Ne V]-weak ($L_{\rm [Ne\,V]}$/$L_{\rm 12,AGN} \leq -$3.4) among 30 AGNs in local U/LIRGs. For the [O IV]-weak AGNs, we estimate their covering fractions of Compton-thick (CT; $N_{\rm H} \geq 10^{24}$ cm$^{-2}$) material with an X-ray clumpy torus model to be $f^{\rm (spec)}_{\rm CT} = 0.55\pm0.19$ on average. This value is consistent with the fraction of CT AGNs ($f^{\rm (stat)}_{\rm CT} = 53\pm12$%) among the [O IV]-weak AGNs in U/LIRGs and much larger than that in Swift/BAT AGNs ($23\pm6$%). The fraction of [O IV]-weak AGNs increases from $27^{+13}_{-10}$% (early) to $66^{+10}_{-12}$% (late mergers). Similar results are obtained with the [Ne V] line. The [O IV] or [Ne V]-weak AGNs in late mergers show larger $N_{\rm H}$ and Eddington ratios ($\lambda_{\rm Edd}$) than those of the Swift/BAT AGNs, and the largest $N_{\rm H}$ is $\gtrsim$10$^{25}$ cm$^{-2}$ at ${\log}\lambda_{\rm Edd} \sim -$1, close to the effective Eddington limit for CT material. These suggest that (1) the circumnuclear material in buried AGNs is regulated by the radiation force from high-$\lambda_{\rm Edd}$ AGNs on the CT obscurers, and (2) their dense material with large $f^{\rm (spec)}_{\rm CT}$ ($\sim 0.5 \pm 0.1$) in U/LIRGs is a likely cause of a unique structure of buried AGNs, whose amount of material may be maintained through merger-induced supply from their host galaxies.

We present the first spectroscopic ALMA follow-up for a pilot sample of nine Radio-Selected NIRdark galaxies in the COSMOS field. These sources were initially selected as radio-detected sources (S(3GHz)>12.65 uJy), lacking an optical/NIR counterpart in the COSMOS2015 catalog (Ks>24.7 mag). Several studies highlighted how this selection could provide a population of highly dust-obscured, massive, and star-bursting galaxies. With these new ALMA observations, we assess the spectroscopic redshifts of this pilot sample of sources and improve the quality of the physical properties estimated through SED-fitting. Moreover, we measure the quantity of molecular gas present inside these galaxies and forecast their potential evolutionary path, finding that the RS-NIRdark galaxies could represent a likely population of high-z progenitors of the massive and passive galaxies discovered at z~3. Finally, we present some initial constraints on the kinematics of the ISM within the analyzed galaxies, reporting a high fraction (~55%) of double-peaked lines that can be interpreted as the signature of a rotating structure in our targets or with the presence of major mergers in our sample. Our results presented in this paper showcase the scientific potential of (sub)mm observations for this elusive population of galaxies and highlight the potential contribution of these sources in the evolution of the massive and passive galaxies at high-z.

A sizable fraction of the heavy elements synthesized by stars in galaxies condenses into sub-micron-sized solid-state particles, known as dust grains. Dust produces a wavelength-dependent attenuation, $A_\lambda$, of the galaxy emission, thereby significantly altering its observed properties. Locally, $A_\lambda$ is in general the sum of a power-law and a UV feature ('bump') produced by small, carbon-based grains. However, scant information exists regarding its evolution across cosmic time. Here, leveraging data from 173 galaxies observed by the James Webb Space Telescope in the redshift range z = 2 - 12, we report the most distant detection of the UV bump in a z ~ 7.55 galaxy (when the Universe was only ~ 700 Myr old), and show for the first time that the power-law slope and the bump strength decrease towards high redshifts. We propose that the flat $A_\lambda$ shape at early epochs is produced by large grains newly formed in supernova ejecta, which act as the main dust factories at such early epochs. Importantly, these grains have undergone minimal reprocessing in the interstellar medium due to the limited available cosmic time. This discovery opens new perspectives in the study of cosmic dust origin and evolution.

The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third $r$-process peak, but robustly produce a weak $r$-process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show that the PNS rotation rate significantly affects the thermodynamic conditions of the wind. We show that high entropy material is quasi-periodically ejected from the closed zone of the PNS magnetosphere with the required thermodynamic conditions to produce heavy elements. We show that maximum entropy $S$ of the material ejected depends systematically on the magnetar spin period $P_{\star}$ and scales as $S \propto P_{\star}^{-5/6}$ for sufficiently rapid rotation. We present results from simulations at a constant neutrino luminosity representative of $\sim 1-2$ s after the onset of cooling for $P_{\star}$ ranging from 5 ms to 200 ms and a few simulations with evolving neutrino luminosity where we follow the evolution of the magnetar wind until $10-14$ s after the onset of cooling. We estimate at magnetar polar magnetic field strength $B_0=3\times 10^{15}$ G, $10^{15}$ G, and $5\times 10^{14}$ G that neutron-rich magnetar winds can respectively produce at least $\sim 2-5\times 10^{-5}$ M$_{\odot}$, $\sim 3-4\times 10^{-6}$ M$_{\odot}$, and $\sim 2.5\times 10^{-8}$ M$_{\odot}$ of material with the required parameters for synthesis of the third $r-$process peak, within $1-2$ s, 10 s, and 14 s in that order after the onset of cooling. We show that proton-rich magnetar winds can have favorable conditions for production of $p-$nuclei, even at a modest $B_0=5\times 10^{14}$ G.

Very recently, the Large High Altitude Air Shower Observatory (LHAASO) reported the observation of the very early TeV afterglow of the brightest-of-all-time GRB 221009A, recording the highest photon statistics in the TeV band ever from a gamma-ray burst. We use this unique observation to place stringent constraints on an energy dependence of the speed of light in vacuum, a manifestation of Lorentz invariance violation (LIV) predicted by some quantum gravity (QG) theories. Our results show that the 95% confidence level lower limits on the QG energy scales are $E_{\mathrm{QG},1}>10$ times of the Planck energy $E_\mathrm{Pl}$ for the linear, and $E_{\mathrm{QG},2}>6\times10^{-8}E_\mathrm{Pl}$ for the quadratic LIV effects, respectively. Our limits on the quadratic LIV case improve previous best bounds by factors of 5--7.

The evolution of the solar wind electron distribution function with heliocentric distance exhibits different features that are still unexplained, in particular, the increase of the Strahl pitch angle width. Wave-particle interactions between electrons and whistler waves are often proposed to explain these phenomena. We aim at quantifying the effect of whistler waves on suprathermal electrons as a function of heliocentric distance. We first perform a statistical analysis of whistler waves (occurrence and properties) observed by Solar Orbiter and Parker Solar Probe between 0.2 and 1 AU. The wave characteristics are then used to compute the diffusion coefficients in the framework of quasi-linear theory. These coefficients are integrated in order to deduce the overall effect of whistler waves on electrons along their propagation. 110,000 whistler wave packets are detected and characterized in the plasma frame. Most waves are aligned with the magnetic field and only about 0.5% of them have a propagation angle greater than 45{\deg}. Beyond 0.3 AU, almost exclusively anti-sunward waves (some of them are found sunward but are within switchbacks with a change of sign of the radial component of the background magnetic) are observed. These waves are therefore Strahl-aligned and not counter-streaming. At 0.2 AU we find both Strahl-aligned and counter-streaming waves. Beyond 0.3 AU, the integrated diffusion coefficients show that the observed waves can explain the measured Strahl pitch angle evolution and are effective in isotropizing the halo. Strahl diffusion is mainly due to whistler waves with an angle of propagation between 15{\deg} and 45{\deg}. Near 0.2 AU, counter-streaming whistler waves can diffuse the Strahl electrons more efficiently than the Strahl-aligned waves by two orders of magnitude.

High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the $(u)gri(z)$ bands, over a collective sky area of $\sim 38$ deg${\rm ^2}$. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependences on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first $\sim 3.5$ years of light curves during 2019A--2022A.

Recent observations of high-redshift galaxies ($z \lesssim 7$) reveal that a substantial fraction have turbulent, gas-rich disks with well-ordered rotation and elevated levels of star formation. In some instances, disks show evidence of spiral arms, with bar-like structures. These remarkable observations have encouraged us to explore a new class of dynamically self-consistent models using our hydrodynamic N-body simulation framework that mimic a plausible progenitor of the Milky Way at high redshift. We explore disk gas fractions of $f_{\rm gas} = 0, 20, 40, 60, 80, 100\%$ and track the creation of stars and metals. The high gas surface densities encourage vigorous star formation, which in turn couples with the gas to drive turbulence. We explore three distinct histories: (i) there is no ongoing accretion and the gas is used up by the star formation; (ii) the star-forming gas is replenished by cooling in the hot halo gas; (iii) in a companion paper, we revisit the models in the presence of a strong perturbing force. At low $f_{\rm disk}$ ($<0.3$), where $f_{\rm disk}$ is the mass fraction of stars relative to dark matter within 2.2 $R_{\rm disk}$, a bar does not form in a stellar disk; this remains true even when gas dominates the inner disk potential. For a dominant baryon disk ($f_{\rm disk} \gtrsim 0.5$) at all gas fractions, the turbulent gas forms a strong "radial shear flow" that leads to an intermittent star-forming bar within about 500 Myr; turbulent gas speeds up the formation of bars compared to gas-free models. For $f_{\rm gas} \lesssim 60\%$, all bars survive, but for higher gas fractions, the bar devolves into a central bulge after 1 Gyr. The star-forming bars are reminiscent of recent discoveries in high-redshift ALMA observations of gaseous disks.

JWST has recently sparked a new era of Lya spectroscopy, delivering the first measurements of the Lya escape fraction and velocity profile in typical galaxies at z~6-10. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust intrinsic models of Lya properties in galaxies at z~5-6 when the IGM is mostly ionized. Here we use new JWST observations from the JADES and FRESCO surveys to characterize statistical distributions of Lya velocity offsets, escape fractions, and EWs in z~5-6 galaxies that will be applicable to growing datasets at z>6. We find that galaxies with large Lya escape fractions (>0.2) are common at z~5-6, comprising 30% of Lyman break selected samples. Comparing to literature studies, our census suggests that Lya becomes more prevalent in the galaxy population toward higher redshift from z~3 to z~6, although we find that this evolution slows considerably between z~5 and z~6, consistent with modest attenuation from residual HI in the mostly ionized IGM at z~5-6. We find significant evolution in Lya velocity profiles between z~2-3 and z~5-6. At lower redshifts, the strongest Lya emitters often have line profiles peaking near the systemic redshift, reflecting escape through low HI density channels. At z~5-6, the strongest Lya emitters have profiles with flux emerging at typical redshifted velocities ~230km/s. The rarity of Lya emitters with peak flux near the systemic redshift at z~5-6 may reflect the influence of resonant scattering from residual HI in the IGM. This effect will make it challenging to use Lya peak offsets as a probe of Lyman continuum leakage at z~5-6. We use our z~5-6 Lya distributions to make predictions for typical Lya properties at z>8 and discuss implications of a recently-discovered Lya emitter at z=8.5 with a small peak velocity offset (156km/s).

The Gaia-ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted $10^5$ stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ($R\sim47,000$) and GIRAFFE medium ($R\sim20,000$) resolution spectra. The analysis of the Gaia-ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a bayesian Inference method developed specifically for the Gaia-ESO Survey by WG11. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6,231 with typical uncertainties for Teff, log g and [Fe/H] of 32~K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6,188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76,675 with typical uncertainties for Teff, log g and [Fe/H] of 64~K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64,177 of the stars.

Statistical analysis of samples of the orbits of celestial bodies is complicated by the fact that the Keplerian orbit is a multidimensional object, the coordinate representation of which nonlinearly depends on the choice of orbital elements. In this work, using the construction of the Fr\'echet mean, concepts of mean orbit and dispersion of the orbit family are introduced, consistent with the distance function on the orbit set. The introduced statistical characteristics serve as analogs of sample mean and variance of a one-dimensional random variable. Exact formulas for calculating the elements of mean orbits and dispersion quantities with respect to two metrics on the orbit space are derived. For a large sample of meteoroid orbits from the Geminid stream, numerical simulations of orbit evolution over 20,000 years in the past are conducted. By analyzing the dependency of statistical characteristics on time, estimates for the age of the stream and the gas outflow velocity are obtained under the assumption of the birth of the Geminids due to the rapid destruction of the cometary nucleus.

Abstract abridged. Eclipsing binary systems provide the opportunity to measure the fundamental parameters of their component stars in a stellar-model-independent way. This makes them ideal candidates for testing and calibrating theories of stellar structure and (tidal) evolution. Even without spectroscopic follow-up there is often enough information in their photometric time series to warrant analysis, especially if there is an added value present in the form of intrinsic variability, such as pulsations. Our goal is to implement and validate a framework for the homogeneous analysis of large numbers of eclipsing binary light curves, such as the numerous high-duty-cycle observations from space missions like TESS. The aim of this framework is to be quick and simple to run and to limit the user's time investment when obtaining, amongst other parameters, orbital eccentricities. We developed a new and fully automated methodology for the analysis of eclipsing binary light curves with or without additional intrinsic variability. Our method includes a fast iterative pre-whitening procedure. Orbital and stellar parameters are measured under the assumption of spherical stars of uniform brightness. We tested our methodology in two settings: a set of synthetic light curves with known input and the catalogue of Kepler eclipsing binaries. The synthetic tests show that we can reliably recover the frequencies and amplitudes of the sinusoids included in the signal as well as the input binary parameters. Recovery of the tangential component of eccentricity is the most accurate and precise. Kepler results confirm a robust determination of orbital periods, with 80.5% of periods matching the catalogued ones. We present the eccentricities for this analysis and show that they broadly follow the theoretically expected pattern as a function of the orbital period.

The circular orbits and elliptical orbits of moving objects in a gravitational field are essential information in astronomy. There have been many methods developed in the literature and textbooks to describe these orbits. In this report, I propose to use the vis-viva equation to construct a complex function to store the state of a moving object in elliptical orbits such that one can calculate its near future numerically. This state function is constructed by splitting its momentum into real and imaginary parts with one perpendicular to the radius of the mass center and the other parallel. This idea is inspired by the wavefunctions of electrons of hydrogen atoms in quantum mechanics. The equations are derived for one-body problems. Two-body problems can be constructed with subsets of one-body problems with the same center of mass, but different effective mass pinned there, which is significantly different from existing methods and provides the same results.

Since galaxy distribution reconstruction effectively reduces non-Gaussian terms in the power spectrum covariance matrix, it has attracted interest not only for Baryon Acoustic Oscillation (BAO) signals but also for various cosmological signal analyses. To this end, this paper presents a novel theoretical model that addresses infrared (IR) effects in the post-reconstruction galaxy power spectrum, including 1-loop corrections. In particular, we discuss the importance of incorporating non-perturbative effects arising from IR contributions into the displacement vector $\VEC{s}$ used for reconstruction. Consequently, post-reconstruction nonlinear damping of BAO can be described by a single two-dimensional Gaussian function. This is a phenomenon not observed when $\VEC{s}$ is considered to at a linear order in the Zel'dovich approximation. Furthermore, we confirm that the cross-power spectrum of the pre- and post-reconstruction density fluctuations lacks IR effect cancellations, and shows an exponential decay in both the cross-power spectrum and the associated shot-noise term.

We report sub-pc-scale observations of the 321-GHz H$_2$O emission line in the radio galaxy NGC 1052. The H$_2$O line emitter size is constrained in $< 0.6$ milliarcsec distributed on the continuum core component. The brightness temperature exceeding $10^6$ K and the intensity variation indicate certain evidence for maser emission. The maser spectrum consists of redshifted and blueshifted velocity components spanning $\sim 400$ km s$^{-1}$, separated by a local minimum around the systemic velocity of the galaxy. Spatial distribution of maser components show velocity gradient along the jet direction, implying that the population-inverted gas is driven by the jets interacting with the molecular torus. We identified significant change of the maser spectra between two sessions separated by 14 days. The maser profile showed a radial velocity drift of $127 \pm 13$ km s$^{-1}$ yr$^{-1}$ implying inward gravitational acceleration at 5000 Schwarzschild radii. The results demonstrate feasibility of future VLBI observations to resolve the jet-torus interacting region.

Clusters of galaxies have been found to host Mpc-scale diffuse, non-thermal radio emission in the form of central radio halos and peripheral relics. Turbulence and shock-related processes in the intra-cluster medium are generally considered responsible for the emission, though details of these processes are still not clear. The low surface brightness makes detection of the emission a challenge, but with recent surveys with high-sensitivity radio telescopes we are beginning to build large samples of these sources. The Evolutionary Map of the Universe (EMU) is a Southern Sky survey being performed by the Australian SKA Pathfinder (ASKAP) over the next few years and is well-suited to detect and characterise such emission. To assess prospects of the full survey, we have performed a pilot search of diffuse sources in 71 clusters from the Planck Sunyaev-Zeldovich (SZ) cluster catalogue (PSZ2) found in archival ASKAP observations. After re-imaging the archival data and performing both (u,v)-plane and image-plane angular scale filtering, we detect 21 radio halos (12 for the first time, excluding an additional six candidates), 11 relics (in seven clusters, and six for the first time, excluding a further five candidate relics), along with 12 other, unclassified diffuse radio sources. From these detections, we predict the full EMU survey will uncover up to ~254 radio halos and ~85 radio relics in the 858 PSZ2 clusters that will be covered by EMU. The percentage of clusters found to host diffuse emission in this work is similar to the number reported in recent cluster surveys with the LOw Frequency ARray (LOFAR) Two-metre Sky Survey (Botteon, et al. 2022a, A&A, 660, A78), suggesting EMU will complement similar searches being performed in the Northern Sky and provide us with statistically significant samples of halos and relics at the completion of the full survey.

The origin of X-ray emission from the resolved kiloparsec-scale jets of many active galactic nuclei (AGN) remains uncertain, particularly where the X-ray emission is separate from the radio-optical synchrotron component. Possible explanations include inverse-Compton upscattering of the Cosmic Microwave Background (IC-CMB) and synchrotron emission from a second electron population, alternatives which imply very different physical conditions within the jet. Until now, X-ray studies of resolved jets have been restricted to the soft X-rays (below $\sim$8 keV), often showing a hard spectral index indicating a spectral peak at higher energies. Here we present \emph{NuSTAR} and \emph{Chandra} observations of the nearby powerful radio galaxy Pictor A, in which we clearly detect the western hot spot at 4' from the host galaxy, the first detection of a kpc-scale jet above 10 keV. The \emph{NuSTAR} 3$-$20 keV spectrum is best fit by a powerlaw of index $\Gamma = 2.03\pm0.04$. Using a maximum likelihood model, we confirm previous reports of variations in the soft X-ray flux detected by \emph{Chandra} over the 2000 to 2015 period, at a 6.7$\sigma$ level, higher than the 3.4$\sigma$ significance level. This rises to $>8\sigma$ in the common 3-8 keV band using the full 22-year span of \emph{Chandra} and \emph{NuSTAR} observations. The variability of Pictor A rules out IC-CMB as the source of the X-ray emission, and the \emph{NuSTAR} observations are in keeping with a synchrotron origin for the X-ray emission from the western hotspot of Pictor A, with no sign of a cutoff to the flat spectrum.

By means of N-body simulations, we study early structure formation in the presence of a scaling distribution of cosmic string loops. Cosmic string loops dominate the high redshift halo mass function while the fluctuations seeded by the standard structure formation scenario dominate structure at low redshifts. In our study, the effects of the cosmic string loops are taken into account by displacing the dark matter particles and their velocities at the initial time of the simulation by amounts determined by the analytical analysis which makes use of the Zeldovich approximation. We find that the resulting halo mass function is to a good approximation given by the sum of the analytically determined cosmic string halo mass function and the halo mass function obtained from the standard $\Lambda$CDM model.

Stellar abundances for a large number of stars are key information for the study of Galactic formation history. Large spectroscopic surveys such as DESI and LAMOST take median-to-low resolution ($R\lesssim5000$) spectra in the full optical wavelength range for millions of stars. However, line blending effect in these spectra causes great challenges for the elemental abundances determination. Here we employ the DD-PAYNE, a data-driven method regularised by differential spectra from stellar physical models, to the DESI EDR spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperature $T_{\rm eff}$, surface gravity $\log g$, microturbulence velocity $v_{\rm mic}$, and abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, Ni. Given a spectral signal-to-noise ratio of 100 per pixel, internal precision of the label estimates are about 20 K for $T_{\rm eff}$, 0.05 dex for $\log~g$, and 0.05 dex for most elemental abundances. These results are agree with theoretical limits from the Cr\'amer-Rao bound calculation within a factor of two. The Gaia-Enceladus-Sausage that contributes the majority of the accreted halo stars are discernible from the disk and in-situ halo populations in the resultant [Mg/Fe]-[Fe/H] and [Al/Fe]-[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread a distance out to $\sim$100 kpc. The DESI sample has a significant higher fraction of distant (or metal-poor) stars than other existed spectroscopic surveys, making it a powerful data set to study the Galactic outskirts. The catalog is publicly available.

The work presented forms a part of the IASc-INSA-NASI Summer Research Fellowship (SRFP) project on the physical & chemical properties of planetary nebulae. Spectral observations of NGC3242 recorded using the VBT Observatory (IIA), Kavalur (TN) for orthogonal slit positions are used for the study. The spectral data were reduced and analyzed using the IRAF data analysis package following the standard procedure to obtain the one dimensional spectrum. Nebular Analysis is performed on the same using NEAT (Wesson et al, 2012)to obtain elemental abundances, densities and temperatures for the PNe and the results are compared to existing literature.

Minkowski tensors are comprehensive shape descriptors that robustly capture n-point information in complex random geometries and that have already been extensively applied in the Euclidean plane. Here, we devise a novel framework for Minkowski tensors on the sphere. We first advance the theory by introducing irreducible Minkowski tensors, which avoid the redundancies of previous representations. We, moreover, generalize Minkowski sky maps to the sphere, i.e., a concept of local anisotropy, which easily adjusts to masked data. We demonstrate the power of our new procedure by applying it to simulations and real data of the Cosmic Microwave Background, finding an anomalous region close to the well-known Cold Spot. The accompanying open-source software, litchi, used to generate these maps from data in the HEALPix-format is made publicly available to facilitate broader integration of Minkowski maps in other fields, such as fluid demixing, porous structures, or geosciences more generally.

HD 12098 is an roAp star pulsating in the most distorted dipole mode yet observed in this class of star. Using TESS Sector 58 observations we show that there are photometric spots at both the magnetic poles of this star. It pulsates obliquely primarily in a strongly distorted dipole mode with a period of $P_{\rm puls} = 7.85$ min ($\nu_{\rm puls} = 183.34905$ d$^{-1}$; 2.12210 mHz) that gives rise to an unusual quadruplet in the amplitude spectrum. Our magnetic pulsation model cannot account for the strong distortion of the pulsation in one hemisphere, although it is successful in the other hemisphere. There are high-overtone p~modes with frequencies separated by more than the large separation, a challenging problem in mode selection. The mode frequencies observed in the TESS data are in the same frequency range as those previously observed in ground-based Johnson $B$ data, but are not for the same modes. Hence the star has either changed modes, or observations at different atmospheric depth detect different modes. There is also a low-overtone p mode and possibly g modes that are not expected theoretically with the $> 1$ kG magnetic field observed in this star.

We present our optical photometric observations of the 2022 eruption of the recurrent nova U Scorpii (U Sco) using 49,152 data points over 70 d following the optical peak. We have also analyzed its soft X-ray (0.3--1 keV) light curve by the Neil Gehrels Swift Observatory. During the 2022 eruption, the optical plateau stage started 13.8--15.0 d and ended 23.8--25.0 d after the optical peak. The soft X-ray stage started 14.6--15.3 d and ended 38.7--39.5 d after the optical peak. Both stages started later and had shorter durations, and the soft X-ray light curve peaked earlier and was less luminous compared to those during the U Sco 2010 eruption. These points suggest that there were differences in the envelope mass between the different cycles of the nova eruption. Furthermore, we have analyzed the optical eclipses during the 2022 eruption. The primary eclipse was first observed 10.4--11.6 d after the optical peak, earlier than the beginning of the optical plateau stage. This sequence of events can be explained by the receding ejecta photosphere associated with the expanding nova ejecta. We have determined the ingress and egress phases of the primary eclipses and estimated the outer radius of the optical light source centered at the white dwarf (WD). During the optical plateau stage, the source radius remained $\sim$1.2 times larger than the Roche volume radius of the primary WD, being close to the L1 point. When the optical plateau stage ended, the source radius drastically shrank to the tidal truncation radius within a few orbital periods. This previously unresolved phenomenon can be interpreted as a structural change in U Sco where the temporarily expanded accretion disk due to the nova wind returned to a steady state.

We investigate the long-term motion of Saturn's North-Pole Hexagon and the structure of its associated eastward jet, using Cassini ISS and ground-based images from 2008 to 2014. We show that both are persistent features that have survived the long polar night, the jet profile remaining essentially unchanged. During those years the hexagon vertices showed a steady rotation period of 10 hr 39 min 23.01 $\pm$ 0.01 s. Analysis of Voyager 1 and 2 (1980-1981) and HST and ground-based (1990-91) images shows a period shorter by 3.5s, due to the presence at the time of a large anticyclone. We interpret the hexagon as a manifestation of a vertically trapped Rossby wave on the polar jet and, because of their survival and unchanged properties under the strong seasonal variations in insolation, we propose that both hexagon and jet are deep-rooted atmospheric features that could reveal the true rotation of the planet Saturn.

Nearly a hundred progenitor-less, thin stellar streams have been discovered in the Milky Way, thanks to Gaia and related surveys. Most streams are believed to have formed from star clusters and it was recently proposed that extended star clusters -- rich in stellar-mass black holes (BHs) -- are efficient in creating streams. To understand the nature of stream progenitors better, we quantify the differences between streams originating from star clusters with and without BHs using direct $N$-body models and a new model for the density profiles of streams based on time-dependent escape rates from clusters. The QSG (Quantifying Stream Growth) model facilitates the rapid exploration of parameter space and provides an analytic framework to understand the impact of different star cluster properties and escape conditions on the structure of streams. Using these models it is found that, compared to streams from BH-free clusters on the same orbit, streams of BH-rich clusters: (1) are approximately five times more massive; (2) have a peak density three times closer to the cluster 1 Gyr post-dissolution (for orbits of Galactocentric radius > 10 kpc), and (3) have narrower peaks and more extended wings in their density profile. We discuss other observable stream properties that are affected by the presence of BHs in their progenitor cluster, namely the width of the stream, its radial offset from the orbit, and the properties of the gap at the progenitor's location. Our results provide a step towards using stellar streams to constrain the BH content of dissolved (globular) star clusters.

High-contrast eclipsing binaries with low mass M-dwarf secondaries are precise benchmark stars to build empirical mass-radius relationships for fully convective low-mass ($\rm M_{*} < 0.35\,M_{\rm sun}$) dwarf stars. The contributed light of the M-dwarf in such binaries is usually much less than one~per~cent at optical wavelengths. This enables the detection of circumbinary planets from precise radial velocity measurements. High-resolution cross-correlation techniques are typically used to detect exoplanet atmospheres. One key aspect of these techniques is the post-processing, which includes the removal of telluric and spectral lines of the host star. We introduce the application of such techniques to optical high-resolution spectra of the circumbinary planet-host TOI-1338/BEBOP-1, turning it effectively into a double-lined eclipsing binary. By using simulations, we further explore the impact of post-processing techniques for high-contrast systems. We detect the M-dwarf secondary with a significance of 11-$\sigma$ and measure absolute dynamical masses for both components. Compared to previous model-dependent mass measurements, we obtain a four times better precision. We further find that the post-processing results in negligible systematic impact on the radial velocity precision for TOI-1338/BEBOP-1 with more than $96.6\,$per~cent (1-$\sigma$) of the M-dwarf's signal being conserved. We show that these methods can be used to robustly measure dynamical masses of high-contrast single-lined binaries providing important benchmark stars for stellar evolution particularly near the bottom of the main sequence. We also demonstrate how to retrieve the phase curve of an exoplanet with high-resolution spectroscopy using our data.

We present new three-dimensional (3D) interstellar extinction maps in the $V$ and Gaia $G$ filters within 2 kpc of the Sun, a 3D differential extinction (dust spatial distribution density) map along the lines of sight in the same space, a 3D map of variations in the ratio of the extinctions in the $V$ and Gaia $G$ filters within 800 pc of the Sun, and a 2D map of total Galactic extinction through the entire dust half-layer from the Sun to extragalactic space for Galactic latitudes $|b|>13^{\circ}$. The 3D maps have a transverse resolution from 3.6 to 11.6 pc and a radial resolution of 50 pc. The 2D map has an angular resolution of 6.1 arcmin. We have produced these maps based on the Gaia DR3 parallaxes and Gaia, Pan-STARRS1, SkyMapper, 2MASS, and WISE photometry for nearly 100 million stars. We have paid special attention to the space within 200 pc of the Sun and high Galactic latitudes as regions where the extinction estimates have had a large relative uncertainty so far. Our maps estimate the extinction within the Galactic dust layer from the Sun to an extended object or through the entire dust half-layer from the Sun to extragalactic space with a precision $\sigma(A_\mathrm{V})=0.06$ mag. This gives a high relative precision of extinction estimates even at high Galactic latitudes, where, according to our estimates, the median total Galactic extinction through the entire dust half-layer from the Sun to extragalactic objects is $A_\mathrm{V}=0.12\pm0.06$ mag. We have shown that the presented maps are among the best ones in data amount, space size, resolution, precision, and other properties.

Whilst the formation of Solar system planets is constrained by meteoritic evidence, the geophysical history of low-mass exoplanets is much less clear. The bulk composition and climate states of rocky exoplanets may vary significantly based on the composition and properties of the planetesimals they form from. An important factor influenced by planetesimal composition is water content, where the desiccation of accreting planetesimals impacts the final water content of the resultant planets. While the inner planets of the Solar system are comparatively water-poor, recent observational evidence from exoplanet bulk densities and planetary formation models suggest that rocky exoplanets engulfed by substantial layers of high-pressure ices or massive steam atmospheres could be widespread. Here we quantify variations in planetesimal desiccation due to potential fractionation of the two short-lived radioisotopes 26Al and 60Fe relevant for internal heating on planetary formation timescales. We focus on how order of magnitude variations in 60Fe can affect the water content of planetesimals, and how this may alter the formation of extrasolar ocean worlds. We find that heating by 26Al is the dominant cause of planetesimal heating in any Solar system analogue scenario, thus validating previous works focussing only on this radioisotope. However, 60Fe can become the primary heating source in the case of high levels of supernova enrichment in massive star-forming regions. These diverging scenarios can affect the formation pathways, bulk volatile budget, and climate diversity of low-mass exoplanets.

Cool gaseous exoplanets ($1.75\ R_\oplus < R_\text{p} < 3\ R_\text{J}$, $200$ K $<T_\text{eq} < 1000$~K) are an as-yet understudied population, with great potential to expand our understanding of planetary atmospheres and formation mechanisms. In this paper, we outline the basis for a homogeneous survey of cool gaseous planets with Twinkle, a 0.45-m diameter space telescope with simultaneous spectral coverage from 0.5-4.5~$\mu$m, set to launch in 2025. We find that Twinkle has the potential to characterise the atmospheres of 36 known cool gaseous exoplanets (11~sub-Neptunian, 11~Neptunian, 14~Jovian) at an SNR $\geq$ 5 during its 3-year primary mission, with the capability of detecting most major molecules predicted by equilibrium chemistry to > $5\sigma$ significance. We find that an injected mass-metallicity trend is well-recovered, demonstrating Twinkle's ability to elucidate this fundamental relationship into cool regime. We also find Twinkle will be able to detect cloud layers at 3$\sigma$ or greater in all cool gaseous planets for clouds at $\leq$ 10 Pa pressure level, but will be insensitive to clouds deeper than $10^4$ Pa in all cases. With these results we demonstrate the capability of the Twinkle mission to greatly expand the current knowledge of cool gaseous planets, enabling key insights and constraints to be obtained for this poorly-charted region of exoplanet parameter space.

Context. Globular clusters are considered key objects for understanding the formation and evolution of the Milky Way. In this sense, their characterisation in terms of their chemical and orbital parameters can provide constraints to the chemical evolution models of the Galaxy. Aims. We use the heavy element abundances of globular clusters to trace their overall behaviour in the Galaxy, aiming to analyse potential relations between the hot H-burning and s-process elements. Methods. We measured the content of Cu I and s- and r-process elements (Y II, Ba II, La II, and Eu II) in a sample of 210 giant stars in 18 Galactic Globular Clusters from high-quality UVES spectra. The clusters span a large metallicity range, and the sample is the largest uniformly analysed for what concerns heavy elements in Globular Clusters. Results. Cu abundances did not show considerable spread in the sample nor correlation with Na, meaning that the Na nucleosynthesis process does not affect the Cu abundance. Most GCs closely follow the Cu, Y, Ba, La, and Eu field stars' distribution, revealing a similar chemical evolution. The Y abundances in mid-metallicity regime GCs (-1.10 dex <[Fe/H]<-1.80 dex) display a mildly significant correlation with the Na abundance, which should be further investigated. Finally, we did not find any significant difference between the n-capture abundances among GCs with Galactic and extragalactic origin.

Neutron star mergers are astrophysical `gold mines,' synthesizing over half of the elements heavier than iron through rapid neutron capture nucleosynthesis. The observation of the binary neutron star merger GW170817, detected both in gravitational waves and electromagnetic radiation, marked a breakthrough. One electromagnetic component of this event, the gamma ray burst GRB 170817A, has an unresolved aspect: the characteristics of its prompt gamma-ray emission spectrum. In this work, we propose that gamma-ray spectra in such GRBs may be influenced by de-excitations from isomeric transitions. Our study begins with a review of current knowledge on GRB structure and of r-process nucleosynthesis in neutron star collisions, focusing on the role of nuclear isomers in these settings. We then test our hypothesis by developing criteria to select representative isomers, based on known solar element abundances, for modeling GRB spectral characteristics. We integrate these criteria into an interactive web page, facilitating the construction and analysis of relevant gamma-ray spectra from isomeric transitions. Our analysis reveals that three isomers (zirconium, lead and yttrium) stand out for their potential to impact the prompt GRB spectrum due to their specific properties. This information allows us to incorporate nuclear isomer data into astrophysical simulations and calculate isomeric abundances generated by astrophysical r-processes in neutron star mergers and their imprint on the detected signal.

We fit various colour-magnitude diagrams (CMDs) of the Galactic globular clusters NGC\,6397 and NGC\,6809 (M55) by isochrones from the Dartmouth Stellar Evolution Database (DSED) and Bag of Stellar Tracks and Isochrones (BaSTI) for $\alpha$-enhanced [$\alpha$/Fe]$=+0.4$. For the CMDs, we use data sets from {\it HST}, {\it Gaia}, VISTA, and other sources utilizing 32 and 23 photometric filters for NGC\,6397 and NGC\,6809, respectively, from the ultraviolet to mid-infrared. We obtain the following characteristics for NGC\,6397 and NGC\,6809, respectively: metallicities [Fe/H]$=-1.84\pm0.02\pm0.1$ and $-1.78\pm0.02\pm0.1$ (statistic and systematic uncertainties); distances $2.45\pm0.02\pm0.06$ and $5.24\pm0.02\pm0.18$ kpc; ages $12.9\pm0.1\pm0.8$ and $13.0\pm0.1\pm0.8$ Gyr; reddenings $E(B-V)=0.178\pm0.006\pm0.01$ and $0.118\pm0.004\pm0.01$ mag; extinctions $A_\mathrm{V}=0.59\pm0.01\pm0.02$ and $0.37\pm0.01\pm0.04$ mag; extinction-to-reddening ratio $R_\mathrm{V}=3.32^{+0.32}_{-0.28}$ and $3.16^{+0.66}_{-0.56}$. Our estimates agree with most estimates from the literature. BaSTI gives systematically higher [Fe/H] and lower reddenings than DSED. Despite nearly the same metallicity, age, and helium enrichment, these clusters show a considerable horizontal branch (HB) morphology difference, which must therefore be described by another parameter. This parameter must predominantly explain why the least massive HB stars (0.58-0.63 solar masses) are only found within NGC 6809. Probably they have been lost by the core-collapse cluster NGC\,6397 during its dynamical evolution and mass segregation. In contrast, NGC\,6809 has a very low central concentration and, hence, did not undergo this process.

Light echoes of debris disks around active stars can reveal disk structure and composition even when disks are not spatially resolved. Unfortunately, distinguishing reflected light from quiescent starlight and unexpected post-peak flare structure is challenging, especially for edge-on geometries where the time delay between observed flare photons and light scattered from the near side of the disk is short. Here, we take advantage of the fact that scattered light from a dusty disk is polarized, depending on the location of the scattering site and the orientation of the disk relative to a distant observer. Filtering reflected light into its polarized components allows echoes to stand out in predictable ways. We test this idea with a simple model for a disk around an active M dwarf. Our results demonstrate that the use of polarimetric data of flaring stars can significantly enhance echo signals relative to starlight and yield more robust and accurate fits to disk parameters compared to analyses based on the total intensity alone.

Dust in the interstellar medium (ISM) is critical to the absorption and intensity of emission profiles used widely in astronomical observations, and necessary for star and planet formation. Supernovae (SNe) both produce and destroy ISM dust. In particular the destruction rate is difficult to assess. Theory and prior simulations of dust processing by SNe in a uniform ISM predict quite high rates of dust destruction, potentially higher than the supernova dust production rate in some cases. Here we show simulations of supernova-induced dust processing with realistic ISM dynamics including magnetic field effects and demonstrate how ISM inhomogeneity and magnetic fields inhibit dust destruction. Compared to the non-magnetic homogeneous case, the dust mass destroyed within 1 Myr per SNe is reduced by more than a factor of two, which can have a great impact on the ISM dust budget.

Since its discovery in 1992, GRS 1915+105 has been among the brightest sources in the X-ray sky. However, in early 2018, it dimmed significantly and has stayed in this faint state ever since. We report on AstroSat and NuSTAR observation of GRS 1915+105 in its unusual low/hard state during 2019 May. We performed time-resolved spectroscopy of the X-ray flares observed in this state and found that the spectra can be fitted well using highly ionised absorption models. We further show that the spectra can also be fitted using a highly relativistic reflection dominated model, where for the lamp post geometry, the X-ray emitting source is always very close to the central black hole. For both interpretations, the flare can be attributed to a change in the intrinsic flux, rather than dramatic variation in the absorption or geometry. These reflection dominated spectra are very similar to the reflection dominated spectra reported for Active Galactic Nuclei in their low flux states.

We study the possibility that primordial black holes (PBHs) can be formed from large curvature perturbations generated during the waterfall phase transition due to the effects of one-loop radiative corrections of Yukawa couplings between the inflaton and a dark fermion in a non-supersymmetric hybrid inflationary model. We obtain a spectral index, $n_s$ and a tensor-to-scalar ratio, $r$ consistent with the current Planck data. Our findings show that the abundance of PBHs are correlated to the dark fermion mass $m_N$ and peak in the GW spectrum. We identify parameter space where PBHs can be the entire dark matter (DM) candidate of the universe or a fraction of it. Our predictions are consistent with any existing constraints of PBH from microlensing, BBN, and CMB, etc. Moreover, the scenario is also testable via induced gravitational waves (GWs) from first-order scalar perturbations detectable in future observatories such as LISA and ET. For instance, with inflaton mass $m \sim 2\times 10^{12}$ GeV, $m_N \sim 5.4\times 10^{15}$ GeV, we obtain PBHs of around $10^{-13}\, M_\odot$ mass that can explain the entire abundance of DM and predict GWs with amplitude $\Omega_{\rm GW}h^2$ $\sim 10^{-9}$ with peak frequency $f$ $\sim$ $0.1$ Hz in LISA. By explicitly estimating fine-tuning we show the model has very mild tuning. We discuss successful reheating at the end of the inflationary phase via the conversion of the waterfall field into standard model (SM) particles. We also briefly speculate a scenario where the dark fermion can be the possible heavy right-handed neutrino (RHN) which is responsible for generating the SM neutrino masses via the seesaw mechanism. The RHN can be produced due to waterfall field decay and its subsequent decay may also explain the observed baryon asymmetry in the universe via leptogenesis.

Stellar population studies of massive early-type galaxies (ETGs) suggest that the stellar initial mass function may not be universal. In particular, the centres of ETGs seem to contain an excess of low-mass dwarf stars compared to our own Galaxy. Through high resolution MUSE IFU data, we carry out a detailed study of the stellar populations of eight massive ETGs. We use full spectrum fitting to determine ages, element abundances, and IMF slopes for spatially binned spectra. We measure flat gradients in age and [Mg/Fe] ratio, as well as negative gradients in metallicity and [Na/Fe]. We detect IMF gradients in some galaxies, with the centres hosting bottom-heavy IMFs and mass excess factors between 1.5-2.5 compared to a Kroupa IMF. The IMF slope below 0.5~M$_\odot$ varies for our galaxy sample between 1-2.8, with negative radial gradients, while the IMF slope between 0.5-1~M$_\odot$ has a steep value of $\sim$3 with mildly positive gradients for most galaxies. For M87, we find excellent agreement with the dynamical M/L as a function of radius. For the other galaxies, we find systematically higher M/L from stellar populations compared to orbit-based dynamical analysis of the same data. This discrepancy increases with NaI strength, suggesting a combination of calibration issues of this line and correlated uncertainties.

In this paper we investigate the vacuum stability of the non-minimally coupled Standard-Model Higgs during a phase of kinetic domination following the end of inflation. The non-minimal coupling to curvature stabilises the Higgs fluctuations during inflation while driving them towards the instability scale during kination, when they can classically overcome the potential barrier separating the false electroweak vacuum from the true one at super-Planckian field values. Avoiding the instability of the Standard-Model vacuum sets an upper bound on the inflationary scale that depends both on the strength of the non-minimal interaction and on the top quark Yukawa coupling. We find that for all combinations of parameters allowing for the existence of a barrier in the effective potential, classical vacuum stability is generically guaranteed. Interestingly enough, thanks to the explosive particle production in the tachyonic phase, the Higgs itself can be also appointed to the role of reheaton field responsible for the onset of the hot Big Bang era, setting an additional lower bound on the inflationary scale $\mathcal{H}_{\rm inf} \gtrsim 10^8 \text{ GeV}$. Overall, these constraints favour lower masses for the top quark, in agreement with the current measurements of the top quark pole mass. We perform our analysis semi-analytically in terms of the one-loop and three-loop running of the Standard-Model Higgs self-coupling and make use of lattice-based parametric formulas for studying the (re)heating phase derived in arXiv:2307.03774 . For a specific choice of $m_t=171.3 \text{ GeV}$ we perform also an extensive numerical scanning of the parameter space via classical lattice simulations, identifying stable/unstable regions and supporting the previous analytical arguments. For this fiducial value, the heating of the Universe is achieved at temperatures in the range $10^{-2} - 10^5 \text{ GeV}$.

Two way ranging data for spacecraft tracking of China's first Martian mission Tianwen-I is analysed. Shortly before the spacecraft entered the Mars parking orbit, the two way coherent microwave link between the spacecraft and the Earth resembles a long arm gravitational wave interferometer, with both the spacecraft and the Earth regarded as in an approximate free falling state. By carefully selecting and analysing data segments of the time series of the two way ranging data during this time span, a parametric statistical model is built for the data segments and an upper limit for the stochastic gravitational waves background (SGWB) is then estimated within the frequency window 0.1Hz to 0.1 mHz. The upper bound improves considerably on those obtained before. In particular, around the deci-Hz band, there is a three orders improvement on the bound obtained previously by the two way ranging data of the Chang e 3 mission. Scientific applications of the upper bound is then considered and a weak upper bound is worked out for axions which is a promising candidate for ultra light dark matter.

We compute GW170817 constraints on the ultralight isospin-violating mediator that has different couplings to protons and neutrons. Neutron stars, which primarily consist of neutrons, are the ideal places to probe the isospin-violating mediator. Such a mediator can significantly alter the dynamics of the binary neutron star mergers, through both the long-range Yukawa force and the new dipole radiation. We compute the gravitational waveform by taking into account the new physics effects due to the isospin-violating mediator and use the Bayesian inference to analyze the gravitational wave data in the GW170817 event. We find that in the parameter space where the isospin-violating force is screened by the Earth, the GW170817 data provide the leading constraints: the upper bound on the neutron coupling is $f_n \lesssim 10^{-19}$ in the mediator mass range of $\simeq(3\times10^{-16},\,5\times10^{-14})$ eV, surpassing constraints from other experiments, including MICROSCOPE and EW.

A galaxy cluster, such as RX J2129, sometimes produces two or more gravitationally lensed images of more distant galaxies. We attempt to regard pairs of these images as stereo pairs. While not successful due to the small disparity angles involved, we suggest that with the 1011 light amplification anticipated from the Solar Gravitational Lens (SGL), individual stars of the distant galaxy might be resolved, resulting in 3D images.

Third-generation dark matter detectors will be fully sensitive to the boron-8 solar neutrino flux. Because of this, the characterization of such a background has been the subject of extensive analyses over the last few years. In contrast, little is known about the impact of reactor neutrinos. In this letter we report on the implications of such a flux for dark matter direct detection searches. We consider five potential detector deployment sites envisioned by the recently established XLZD consortium: SURF, SNOLAB, Kamioka, LNGS and Boulby. By using public reactor data we construct five reactor clusters -- involving about 100 currently operating commercial nuclear reactors each -- and determine the net neutrino flux at each detector site. Assuming a xenon-based detector and a 50 tonne-year exposure, we show that in all cases the neutrino event rate may be sizable, depending on energy recoil thresholds. Of all possible detector sites, SURF and LNGS are those with the smallest reactor neutrino background. On the contrary, SNOLAB and Boulby are subject to the strongest reactor neutrino fluxes, with Kamioka being subject to a more moderate background. Our findings demonstrate that reactor neutrino fluxes should be taken into account in the next round of dark matter searches. We argue that this background may be particularly relevant for directional detectors, provided they meet the requirements we have employed in this analysis.

We describe an experimental setup developed aiming to irradiate samples under UV radiation for accelerated test for solar effects according to the relevant ECSS-ESA standards. This facility has been already used for projects belonging to large space programs (Cosmic Vision, Artes) for simulations up to 3500 equivalent sun hours. In particular, we detail the calculation of the UV dose delivered by Sun, the calibration of the detectors, the spatial distribution of the UV radiation on samples, the remote control of both samples temperature and lamp radiation, the samples heat dissipation and operation in a helium atmosphere.

Cosmological structure formation simulations of ultralight axion-like dark matter have shown that an axion star forms at the center of every dark matter halo in the Universe. These axion stars would then form in large numbers during the dark ages, $z \lesssim 70$. Axion stars would represent the densest axion environments in the Universe, and as such they can trigger collective processes that cannot otherwise occur for axions in vacuum. In particular, even though the lifetime of individual sub-eV axions decaying into a pair of photons is much larger than the age of the Universe, axion stars can decay into photons on very short time scales due to parametric resonance. In this talk, based on arXiv:2302.10206 and arXiv:2301.09769, I will discuss the cosmological implications of such decays. We show that massive enough axion stars will decay into a large number of radio photons which will in turn lead to heating and ionization during the dark ages which is strongly constrained by Planck. As a result, we find that couplings $10^{-14}\,{\rm GeV}^{-1} \lesssim g_{a\gamma\gamma} \lesssim 10^{-10}\,{\rm GeV}^{-1}$ are excluded by Planck for $10^{-14}\,{\rm eV}\lesssim m_a\lesssim 10^{-8}\,{\rm eV}$ within our benchmark model of axion star abundance. We also highlight that future measurements of the 21 cm line can have sensitivity to couplings at least one order of magnitude smaller.

We study the damping of density oscillations in the quark matter phase that might occur in compact stars. To this end we compute the bulk viscosity and the associated damping time in three-flavor quark matter, considering both nonleptonic and semileptonic electroweak processes. We use two different equations of state of quark matter, more precisely, the MIT bag model and perturbative QCD, including the leading order corrections in the strong coupling constant. We analyze the dependence of our results on the density, temperature and value of strange quark mass in each case. We then find that the maximum of the bulk viscosity is in the range of temperature from 0.01 to 0.1 MeV for frequencies around 1 kHz, while the associated minimal damping times of the density oscillations at those temperatures might be in the range of few to hundreds milliseconds. Our results suggest that bulk viscous damping might be relevant in the post-merger phase after the collision of two neutron stars if deconfined matter is achieved in the process.