Papers for Thursday, Dec 19 2024

The Montage Image Mosaic Engine, first released in 2002, has found applicability across the electromagnetic spectrum to support data processing and visualization. This broad applicability has come about through its design as an Open Source ANSI-C toolkit (and Python binary extensions), with independent components to perform each step in the creation of a mosaic and with support for all WCS extensions. This design enables easy integration into custom environments, workflows and pipelines, and is the principal reason for its long lifetime. Here we emphasize the growing use of Montage in radio astronomy (37 peer-reviewed papers since 2020), and will focus on three high-profile applications: (1) Analysis of observations made with SKA precursor experiments, such as MeerKAT and the Murchison Wide-field Array, (2) Faraday tomography of LOFAR Two-Metre Sky Survey data (LoTSS-DR2), which explores the structure of the local interstellar medium, and (3) Identification of fast radio bursts.

The X-Ray Imaging and Spectroscopy Mission (XRISM) is the 7th Japanese X-ray observatory, whose development and operation are in collaboration with universities and research institutes in Japan, U.S., and Europe, including JAXA, NASA, and ESA. The telemetry data downlinked from the satellite are reduced to scientific products by the pre-pipeline (PPL) and pipeline (PL) software running on standard Linux virtual machines on the JAXA and NASA sides, respectively. We ported the PPL to the JAXA "TOKI-RURI" high-performance computing (HPC) system capable of completing $\simeq 160$ PPL processes within 24 hours by utilizing the container platform of Singularity and its "--bind" option. In this paper, we briefly show the data processing in XRISM and present our porting strategy of PPL to the HPC environment in detail.

The Gemini High-resolution Optical SpecTrograph (GHOST) at Gemini South started regular queue operations in early 2024, bringing a long-sought open-access capability to the astronomy community. This research note briefly describes an effort to provide easy-to-access reduced spectra for GHOST programs from all Gemini partner countries and encourage prompt data exploration and analysis. Since March 2024, over 4500 spectra have been reduced and made available to principal investigators (PIs). The aim is to increase demand for GHOST and expedite the publication of scientific results.

The Neptune desert is no longer empty. A handful of close-in planets with masses between those of Neptune and Saturn have now been discovered, and their puzzling properties have inspired a number of interesting theories on the formation and evolution of desert-dwellers. While some studies suggest that Neptune desert planets form and evolve similarly to longer-period Neptunes, others argue that they are products of rare collisions between smaller planets, or that they are the exposed interiors of giant planets (i.e., ``hot Jupiters gone wrong''). These origin stories make different predictions for the metallicities of Neptune desert host stars. In this paper, we use the homogeneous catalog of stellar metallicities from Gaia Data Release 3 to investigate the origins of Neptune desert-dwellers. We find that planets in the Neptune desert orbit stars that are significantly more metal-rich than the hosts of longer-period Neptunes ($p = 0.0016$) and smaller planets ($p = 0.00014$). In contrast, Neptune desert host star metallicities are statistically indistinguishable from those of hot Jupiter host stars ($p = 0.55$). Therefore, we find it relatively unlikely that Neptune desert planets formed and evolved similarly to longer-period Neptunes, or that they resulted from collisions between smaller planets, at least without another metallicity-selective process involved. A more straightforward explanation for this result is that planets in the desert truly are the exposed interiors of larger planets. Atmospheric spectroscopy of Neptune desert worlds may therefore provide a rare glimpse into the interiors of giant exoplanets.

We investigate primordial non-Gaussianity (NG) arising from the explicit $U(1)$ symmetry-breaking interactions during inflation involving a nearly massless axial component of a complex scalar field $P$. We analyze the induced NG parameter $f_{\mathrm{NL}}$ under scenarios where the axial field functions as either a curvaton or cold dark matter (CDM). In the curvaton framework, there is a conventional contribution to the local NG of $f_{\rm NL} \simeq -O(1)$. Additional positive local NG can result from either the self-interactions of axial field fluctuations, their interactions with a light radial partner, or kinetic mixing with the inflaton via $U(1)$ symmetry-breaking terms. We identify parameter regions where the interactions lead to cancellations, suppressing the overall local NG to $|f^{\rm loc}_{\mathrm{NL}}| \lesssim O(0.1)$, while leaving the trispectrum largely unaffected. In the CDM scenario, these interactions enhance the NG in the isocurvature fluctuations. Moreover, interactions between the axial field and another light scalar, such as a curvaton, can generate $O(1)$ curvature NG signals and significant mixed curvature-isocurvature NGs that are within the reach of future experiments with $\sigma(f^{\rm loc}_{\rm NL})\sim1$. We also explore the role of a heavy radial field in generating oscillating correlation signals, noting that such signals can dominate the shape of the mixed adiabatic-isocurvature bispectrum. In certain cases, an oscillatory isocurvature bispectrum signal may be observable in the future, aiding in distinguishing between certain types of the $U(1)$-breaking self-interactions of the axial field.

We present a spectroscopic analysis of MACS J0138$-$2155, at $z=0.336$, the first galaxy cluster hosting two strongly-lensed supernovae (SNe), Requiem and Encore, providing us with a chance to obtain a reliable $H_0$ measurement from the time delays between the multiple images. We take advantage of new data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope, covering a central $1 \rm \, arcmin^2$ of the lensing cluster, for a total depth of 3.7 hours, including 2.9 hours recently obtained by our Target of Opportunity programme. Our new spectroscopic catalogue contains reliable redshifts for 107 objects, including 50 galaxy cluster members with secure redshift values in the range $0.324 < z < 0.349$, and 13 lensed multiple images from four background sources between $0.767\leq z \leq 3.420$, including four images of the host galaxy of the two SNe. We exploit the MUSE data to study the stellar kinematics of 14 bright cluster members and two background galaxies, obtaining reliable measurements of their line-of-sight velocity dispersion. Finally, we combine these results with measurements of the total magnitude of the cluster members in the Hubble Space Telescope F160W band to calibrate the Faber-Jackson relation between luminosity and stellar velocity dispersion ($L \propto \sigma^{1/\alpha}$) for the early-type cluster member galaxies, measuring a slope $\alpha=0.25^{+0.05}_{-0.05}$. A pure and complete sample of cluster member galaxies and a reliable characterisation of their total mass structure are key to building accurate total mass maps of the cluster, mitigating the impact of parametric degeneracies, which is necessary for inferring the value of $H_0$ from the measured time delays between the lensed images of the two SNe.

In this paper we introduce the clus model, which has been newly implemented in the X-ray spectral fitting software package SPEX. Based on the 3D radial profiles of the gas density, temperature, metal abundance, turbulent, and inflow/outflow velocities, the clus model creates spectra for a chosen projected region on the sky. Additionally, it can also take into account the resonant scattering. We show a few applications of the clus model on simulated spectra of the massive elliptical galaxy NGC 4636, and galaxy clusters A383, A2029, A1795, A262, and the Perseus cluster. We quantify the effect of projection, as well as resonant scattering on inferred profiles of the iron abundance and temperature, assuming the resolution similar to Chandra ACIS-S and XRISM Resolve. Our results show that, depending on the mass of the object, as well as the projected distance from its core, neither a single-temperature, double-temperature, nor the Gaussian-shaped differential emission measure models can accurately describe the input emission measure distribution of these massive objects. The largest effect of projection as well as resonant scattering is seen for projected profiles of iron abundance of NGC 4636, where we are able to reproduce the observed iron abundance drop in its inner-most few kiloparsecs. Furthermore, we find that projection effects also influence the best-fit temperature, and the magnitude of this effect varies depending on the underlying hydrodynamical profiles of individual objects. In the core, the projection effects are the largest for A1795 and NGC 4636, while in the outskirts the largest difference between 2D and 3D temperature profiles are for Perseus and A1795, regardless of the instrumental resolution. These findings might potentially have an impact on cross-calibration studies between different instruments, as well as on the precision cosmology.

The field of exocomets has been built around the unmatched number of detections made in the circumstellar disc of the archetypal star Beta Pictoris. An exocomet detection in spectroscopy is identified by variable atomic absorption features in a stellar spectrum, associated with transiting gas in and trailing an exocomet coma. This paper presents the largest spectroscopic search for exocomet transits to date, which overcomes the limitations of biased samples of stars with debris discs, and instead looks through the $\approx$7500 stars in the HARPS archive for signs of exocomets in the CaII doublet (H:396.847nm and K:393.366nm). The search resulted in 155 candidate stars, which after filtering for false positives (e.g. binaries, stellar activity, etc.), were cut down to 22 stars. These 22 stars are classified into Tier1, 2, and 3 exocomet candidates, reflecting the confidence level of their exocomet detection. Our two best candidates (Tier1: Beta Pictoris, HD172555) and four lower confidence candidates (Tier2: Gl1, HIP5158, HD94771, HR1996) are discussed, yielding a detection rate of 0.03% (Tier1 only) and 0.1% (Tier1 & 2) in the HARPS sample. Both Tier1 stars are known exocomet host stars. These two young A-type stars correspond to 0.4% of all A-types in the sample, suggesting that detecting signs of exocomet transits using CaII is more likely around young A-type stars. Reanalysing a past HARPS study, we found no evidence to support the previously claimed four exocomet detections, indicating either that those detections are not robust or that we are only sensitive to the strongest signals.

The origin of the recombining plasma in several Galactic SNRs has been debated. A plausible mechanism would be a rapid cooling in the past, either by adiabatic or conductive process. A recent spectral study of W 49 B reported a possible charge exchange (CX) emission due to collisions between the shock-heated ejecta and external cold clouds, which could be a direct support for the conduction cooling scenario. However, a potentially large systematic uncertainty in the spectral analysis has not been examined. In this paper, we revisit the Suzaku spectrum of W 49 B with taking into account the systematic uncertainties in spectral codes and instrumental gain calibration. We find that the previously reported CX flux is fully attributable to dielectronic recombination satellite lines of high-shell transitions that are missing from the present version of the spectral codes. We also report refined Fe-group ejecta mass ratios, which, in comparison to those in the literatures, show a better agreement with theoretical expectations from nucleosynthesis models, either of Type Ia explosions or spherical core-collapse explosions.

Transforming the instrumental photometry of ground-based telescopes into a calibrated physical flux in a well-defined passband is a major challenge in astronomy. Along with the intrinsic instrumental difference between telescopes sharing the same filter, the effective transmission is continuously modified by the effects of the variable atmosphere of the Earth. We have developed a new approach to the absolute photometric calibration that simultaneously treats instrumental and atmospheric effects on an image-by-image basis by fitting the system transmission. This approach aims at breaking the 1% absolute photometric accuracy which limits current calibration methods for ground-based observatories. We fit the transmission, as a function of wavelength, for each image. The fit is done by comparing the instrumental fluxes of stars in the image to the synthetic photometry of the stars given their spectrum and the transmission function which have free parameters. A key element that enables this approach is the set of about 220 million low-resolution spectra measured by Gaia, which provides a large number of stellar calibrators in the image. We demonstrate the method using data from the Large Array Survey Telescope (LAST). We show that the residuals between observations and synthetic photometry of the Gaia spectra in the fitted transmission have a standard deviation $<$1% on an image-by-image basis, with no spatial and color dependencies. The median accuracy of the zero-point throughout the image is between 3-5 mmag, depending on the total image exposure. Furthermore we show that this method provides high stability over long temporal scales.

The aim of this work is to study the position of gas-rich and gas-poor galaxy clusters within the large-scale structure and, in particular, their distance to filaments. Our sample is built from 29 of the 34 clusters in the X-ray unbiased cluster sample (XUCS), a velocity-dispersion-selected sample for which various properties, including masses, gas fractions, and X-ray surface brightness were available in the literature. We compute the projected distance between each cluster and the spine of the nearest filament with the same redshift and investigate the link between this distance and the previously-mentioned properties of the clusters, in particular with their gas content. The average distance between clusters and filaments is larger for low X-ray surface brightness clusters than for those of high surface brightness, with intermediate brightness clusters being an intermediate case. Also the minimum distance follows a similar trend, with rare cases of low surface brightness clusters found at distances smaller than 2 Mpc from the spine of filaments. However, the Kolmogorov-Smirnov statistical test is not able to exclude the null hypothesis that the two distributions are coming from the same parent one. We speculate that the position of galaxy clusters within the cosmic web could have a direct impact in their gas mass fraction, hence on its X-ray surface brightness, since the presence of a filament can oppose resistance to the outward flow of gas induced by the central AGN and reduce the time required for this gas to fall inward after the AGN is shut. However, a larger sample of clusters is needed in order to derive a statistically-robust conclusion

Transitional millisecond pulsars (tMSPs) bridge the evolutionary gap between accreting neutron stars in low-mass X-ray binaries and millisecond radio pulsars. These systems exhibit a unique subluminous X-ray state characterized by the presence of an accretion disk and rapid switches between high and low X-ray emission modes. The high mode features coherent millisecond pulsations spanning from the X-ray to the optical band. We present multiwavelength polarimetric observations of the tMSP PSR J1023+0038 aimed at conclusively identifying the physical mechanism powering its emission in the subluminous X-ray state. During the high mode, we detect polarized emission in the 2-6 keV energy range, with a polarization degree of 12% +/- 3% and a polarization angle of -2deg +/- 9deg (1sigma) measured counterclockwise from the North celestial pole towards East. At optical wavelengths, we find a polarization degree of 1.41% +/- 0.04% and a polarization angle aligned with that in the soft X-rays, suggesting a common physical mechanism operating across these bands. Remarkably, the polarized flux spectrum matches the pulsed emission spectrum from optical to X-rays. The polarization properties differ markedly from those observed in other accreting neutron stars and isolated rotation-powered pulsars and are also inconsistent with an origin in a compact jet. Our results provide direct evidence that the polarized and pulsed emissions both originate from synchrotron radiation at the shock formed where the pulsar wind interacts with the inner regions of the accretion disk.

We report the discovery of an ultra-massive grand-design red spiral galaxy, named Zhúlóng (Torch Dragon), at $z_{\rm phot} = 5.2^{+0.3}_{-0.2}$ in the JWST PANORAMIC survey, identified as the most distant bulge+disk galaxy with spiral arms known to date. Zhúlóng displays an extraordinary combination of properties: 1) a classical bulge centered in a large, face-on exponential stellar disk (half-light radius of $R_{\rm e} = 3.7 \pm 0.1 \, \mathrm{kpc}$), with spiral arms extending across 19 kpc; 2) a clear transition from the red, quiescent core ($F150W-F444W=3.1$ mag) with high stellar mass surface density ($\log(\Sigma M_{\star}/M_{\odot} \, \mathrm{kpc}^{-2}) = 9.91_{-0.09}^{+0.11}$) to the star-forming outer regions, as revealed by spatially resolved SED analysis, which indicates significant inside-out galaxy growth; 3) an extremely high stellar mass at its redshift, with $\log (M_{\star}/M_{\odot})=11.03_{-0.08}^{+0.10}$ comparable to the Milky Way, and an implied baryon-to-star conversion efficiency ($\epsilon \sim 0.3$) that is 1.5 times higher than even the most efficient galaxies at later epochs; 4) despite an active disk, a relatively modest overall star formation rate ($\mathrm{SFR} =66_{-46}^{+89} ~M_{\odot} \, \mathrm{yr}^{-1}$), which is $>$0.5 dex below the star formation main sequence at $z \sim 5.2$ and $>$10 times lower than ultra-massive dusty galaxies at $z=5-6$. Altogether, Zhúlóng shows that mature galaxies emerged much earlier than expected in the first billion years after the Big Bang through rapid galaxy formation and morphological evolution. Our finding offers key constraints for models of massive galaxy formation and the origin of spiral structures in the early universe.

Fuzzy dark matter (FDM) granulations would drive orbital transport of stars in galactic disks, and in particular would produce roughly equal amounts of radial heating and radial migration. However, observations suggest that heating has been much less efficient than migration in our Galaxy. We argue that this decreases the amount of radial heating, $\mathcal{H}_\mathrm{FDM}$, that can safely be attributed to FDM. Consequently, lower bounds on the FDM particle mass $m$ derived through Galactic disk kinematics should be revised upwards; a rough estimate is $m \gtrsim 1.3\times 10^{-22} \mathrm{eV} \times [(\mathcal{H}_\mathrm{FDM}/\mathcal{H})/0.1]^{-1/2}$, where $\mathcal{H}$ is the total observed radial heating.

The results of the numerical investigations of the evolution of orbits of trans-Neptunian bodies at the 2 : 3 resonance with Neptune are presented. The gravitational influence of the four giant planets was taken into account. For identical initial values of the semimajor axes, eccentricities, and inclinations, but for different initial orbital orientations and initial positions in orbits, we obtained different types of variations in the difference {\Delta}{\Omega}={\Omega}-{\Omega}_N in the ascending-node longitudes of the body and Neptune, and in the perihelion argument {\omega}. When {\Delta}{\Omega} decreases and {\omega} increases during evolution, then most of the bodies leave the resonance in 20 Myr. In the case of an increase in {\Delta}{\Omega} and a decrease in {\omega}, the bodies stay in the resonance for a much longer time. Regions of eccentricities and inclinations, for which some bodies were in the {\eta}_18 secular resonance ({\Delta}{\Omega} is almost constant) and in the Kozai resonance ({\omega} is almost constant), were obtained to be larger than those predicted for small variations in the critical angle. Some bodies can at the same time be in both these resonances.

Galaxy clusters are important cosmological probes that have helped to establish the $\mathrm{\Lambda}$CDM paradigm as the standard model of cosmology. However, recent tensions between different types of high-accuracy data highlight the need for novel probes of the cosmological parameters. Such a probe is the turnaround density: the mass density on the scale where galaxies around a cluster join the Hubble flow. To measure the turnaround density, one must locate the distance from the cluster center where turnaround occurs. Earlier work has shown that a turnaround radius can be readily identified in simulations by analyzing the 3D dark matter velocity field. However, measurements using realistic data face challenges due to projection effects. This study aims to assess the feasibility of measuring the turnaround radius using machine learning techniques applied to simulated observations of galaxy clusters. We employ N-body simulations across various cosmologies to generate galaxy cluster projections. Utilizing Convolutional Neural Networks (CNNs), we assess the predictability of the turnaround radius based on galaxy line-of-sight velocity, number density, and mass profiles. We find a strong correlation between the turnaround radius and the central mass of a galaxy cluster, rendering the mass distribution outside the virial radius of little relevance to the model's predictive power. The velocity dispersion among galaxies also contributes valuable information concerning the turnaround radius. Importantly, the accuracy of a line-of-sight velocity model remains robust even when the data within $\mathrm{R_{200}}$ of the central overdensity are absent.

The LIGO-Virgo-KAGRA (LVK) collaboration has detected 90 gravitational wave events and will detect many more in its fourth observing run. Binary black hole (BBH) systems represent the overwhelming majority of these observations. We build a model for the population of the BBHs based on the observed distribution of metallicities in galaxies and state-of-the-art stellar evolution models implemented through the Stellar EVolution N-body (SEVN) code. We calculate the primary mass spectrum and merger rates of BBHs and find general agreement with the redshift evolution and mass ratio distribution inferred by the LVK collaboration. When comparing to the primary mass distribution, our results indicate that either the average IMF in dwarf galaxies must be top heavy, or most of the $30-40 \, M_\odot$ black holes must be formed through a dynamical capture mechanism. For masses greater than about $50 \, M_\odot$, the predicted number of BBH systems plummet to zero, revealing the well-known mass gap due to the pair instability mechanism and the mass loss in binary systems. We estimate the probability needed to fill part of the mass gap with mergers of dynamically-formed BBHs originating from the single black hole population.

In the context of Fuzzy Dark Matter (FDM) we study the core formation in the presence of an Ideal Gas (IG). Our analysis is based on the solution of the Schrödinger-Poisson-Euler system of equations that drives the evolution of FDM together with a compressible IG, both coupled through the gravitational potential they produce. Starting from random initial conditions for both FDM and IG, with dominant FDM, we study the evolution of the system until it forms a nearly relaxed, virialized and close to hydrostatic equilibrium core, surrounded by an envelope of the two components. We find that the core corresponds to Newtonian Fermion-Boson Stars (FBS). If the IG is used to model luminous matter, our results indicate that FBS behave as attractor core solutions of structure formation of FDM along with baryonic matter.

In astrophysical and cosmological analyses, the increasing quality and volume of astronomical data demand efficient and precise computational tools. This work introduces a novel adaptive algorithm for automatic knots (AutoKnots) allocation in spline interpolation, designed to meet user-defined precision requirements. Unlike traditional methods that rely on manually configured knot distributions with numerous parameters, the proposed technique automatically determines the optimal number and placement of knots based on interpolation error criteria. This simplifies configuration, often requiring only a single parameter. The algorithm progressively improves the interpolation by adaptively sampling the function-to-be-approximated, $f(x)$, in regions where the interpolation error exceeds the desired threshold. All function evaluations contribute directly to the final approximation, ensuring efficiency. While each resampling step involves recomputing the interpolation table, this process is highly optimized and usually computationally negligible compared to the cost of evaluating $f(x)$. We show the algorithm's efficacy through a series of precision tests on different functions. However, the study underscores the necessity for caution when dealing with certain function types, notably those featuring plateaus. To address this challenge, a heuristic enhancement is incorporated, improving accuracy in flat regions. This algorithm has been extensively used and tested over the years. NumCosmo includes a comprehensive set of unit tests that rigorously evaluate the algorithm both directly and indirectly, underscoring its robustness and reliability. As a practical application, we compute the surface mass density $\Sigma(R)$ and the average surface mass density $\overline{\Sigma}(<R)$ for Navarro-Frenk-White and Hernquist halo density profiles, which provide analytical benchmarks. (abridged)

Universal small-scale solar activity in quiet region are suggested to be a potential source of solar wind and the upper solar atmosphere. Here, with the high-resoltion 174 Å~imaging observations from the Solar Orbiter/Extreme Ultraviolet Imager (EUI), we investigate 59 EUV upflow-like events observed in the quiet Sun. Their average apparent (plane-of-sky) velocity, lifetime, and propagation distance are measured as 62 $\speed$, 68.6 s and 3.94 Mm, respectively. These upflow-like events exhibit dynamic characteristics but lack base brightening, featuring a hot front and subsequent cold plasma ejection. 39\% of the EUV upflow-like events exhibit recurrent characteristics. Unprecedented high-resolution 174 Å~observations reveal that some EUV upflow-like events exhibit blob-like fine structures and multi-strand evolutionary features, and some upflow-like events can cause localized haze-like plasma heating ahead of their spire region during the ejection process. A subset of the EUV upflow-like events covered by the Solar Dynamics Observatory reveals that they appear at the chromospheric networks. Through emission measure analysis, we found that these upflow-like events eject hot plasma of transient region or coronal temperature (an average of $\sim$10$^{5.5}$K). We suggest that EUV upflow-like events may be EUV counterparts of chromospheric spicules and/or transition region network jets, and play a role in heating localized corona above the network regions.

Fast type-I migration of (proto)planets poses a challenging problem for the core accretion formation scenario. We found that the dust-induced ``Streaming Torque (ST)'' may slow down or even reverse the planet migration in \cite{Hou2024}. But in realistic protoplanetary disks, dust diffusion induced by gas turbulence may have important influences on ST. We perform linear analysis to investigate the effects of dust diffusion on ST. The dependence of ST on the dust diffusion may provide better constraints on the turbulence strength and the stopping time $\tau$. We derive the dispersion relation for all the wave modes in the two-fluid system. The dust diffusion will smooth the short-wavelength structure of the the quasi-drift mode and split it into two predominant D-drift modes with opposite directions. The outgoing D-drift mode will contribute to a negative torque on planets, particularly when $\tau \sim 0.1$, which slightly shifts the zero-torque turning point. We explore how ST depends on the regimes of aerodynamic drag, dust mass fraction and disk scale height. We compare the radial wavenumbers of D-drift modes under different formulations of dust diffusion and find qualitative agreement. In all cases, $\tau$ at the zero-torque turning point, which determines the direction of planetary migration, consistently remains on the order of $\sim 0.1$, corresponding to large pebble-sized dust grains. This suggests that rapid dust coagulation can inhibit the inward migration of planets, implying that weak gas turbulence may enhance the survival of protoplanets.

Stellar halos around galaxies contain key information about their formation and assembly history. Using simulations, we can trace the origins of different stellar populations in these halos, contributing to our understanding of galaxy evolution. We aim to investigate the assembly of stellar halos and their chemical abundances in 28 galaxies from CIELO project with logMgal[9 and 11]Msun. Stellar halos were identified using the AM E method, focusing on the outer regions between the 1.5 optical radius and the virial radius. We divided the stellar populations based on their formation channel: exsitu, endodebris, and insitu, and analyzed their chemical abundances, ages, and spatial distributions. Additionally, we explored correlations between halo mass, metallicity, and alpha element enrichment. CIELO simulations reveal that stellar halos are predominantly composed of accreted material (exsitu and endodebris stars), in agreement with previous works. The mass fraction of these populations is independent of stellar halo mass, though their metallicities scale linearly with it. Exsitu stars tend to dominate the outskirts and be more alpha rich and older, while endodebris stars are more prevalent at lower radii and tend to be less alpha rich and slightly younger. Massive stellar halos require a median of five additional satellites to build 90 percent of their mass, compared to lower mass halos, which typically need fewer (median of 2.5) and lower-mass satellites and are assembled earlier. The diversity of accreted satellite histories results in well defined stellar halo mass metallicity and [alpha/Fe] [Fe/H] relations, offering a detailed view of the chemical evolution and assembly history of stellar halos. We find that the [alpha/Fe] [Fe/H] is more sensitive to the characteristics and star formation history of the contributing satellites than the stellar halo mass metallicity relationship

One of the main obstacles in exoplanet detection when using the radial velocity (RV) technique is the presence of stellar activity signal induced by magnetic regions. In this context, a realistic simulated dataset that can provide photometry and spectroscopic outputs is needed for method development. The goal of this paper is to describe two realistic simulations of solar activity obtained from SOAP-GPU and to compare them with real data obtained from the HARPS-N solar telescope. We describe two different methods of modeling solar activity using SOAP-GPU. The first models the evolution of active regions based on the spot number as a function of time. The second method relies on the extraction of active regions from the Solar Dynamics Observatory (SDO) data. The simulated spectral time series generated with the first method shows a long-term RV behavior similar to that seen in the HARPS-N solar observations. The effect of stellar activity induced by stellar rotation is also well modeled with prominent periodicities at the stellar rotation period and its first harmonic. The comparison between the simulated spectral time series generated using SDO images and the HARPS-N solar spectra shows that SOAP-GPU can precisely model the RV time series of the Sun to a precision better than 0.9 m/s. By studying the width and depth variations of each spectral line in the HARPS-N solar and SOAP-GPU data, we find a strong correlation between the observation and the simulation for strong spectral lines, therefore supporting the modeling of the stellar activity effect at the spectral level. These simulated solar spectral time series serve as a useful test bed for evaluating spectral-level stellar activity mitigation techniques.

One way to constrain the evolutionary histories of galaxies is to analyse their stellar populations. In the local Universe, our understanding of the stellar population properties of galaxies has traditionally relied on the study of optical absorption and emission-line features. In order to overcome limitations intrinsic to this wavelength range, such as the age-metallicity degeneracy and the high sensitivity to dust reddening, we must use wavelength ranges beyond the optical. The near-infrared (NIR) offers a possibility to extract information on spectral signatures that are not as obvious in traditional optical bands. Moreover, with the current and forthcoming generation of instrumentation focusing on the NIR, it is mandatory to explore possibilities within this wavelength range for nearby-Universe galaxies. However, although the NIR shows great potential, we are only beginning to understand it. Widely used techniques such as a full spectral fitting and line strength indices need to be tested on systems that are as close to simple stellar populations as possible, and the result from the techniques need to be compared to the yields from a traditional optical analysis. We present a NIR spectral survey of extragalactic globular clusters (GCs). The set was composed of 21 GCs from NGC5128 that were obtained with SOAR/TripleSpec4, which covered the ~1.0-2.4$\mu$m range with a spectral resolution R=$\lambda/\Delta\lambda$ of 3500. These spectra cover H$\beta$ equivalent widths between 0.98$Å$ and 4.32$Å$, and [MgFe]' between 0.24$Å$ and 3.76$Å$. This set is ideal for performing absorption band measurements and a full spectral fitting, and it can be used for kinematic studies and age and abundance measurements. With this library, we expect to be able to probe the capabilities of NIR models, as well as to further improve stellar population estimates for the GCs around NGC5128.

The increasing population of objects in geostationary orbit has raised concerns about the potential risks posed by debris clouds resulting from fragmentation. The short-term evolution and associated hazards of debris generated by collisions in the geostationary region is investigated in this study. The initial distribution of two debris clouds is modeled using a single probability density this http URL combined distribution of the evolved clouds is determined by solving boundary value this http URL risks associated with these debris clouds are evaluated by calculating the instantaneous impact rate and cumulative collision this http URL probability of collisions with millimeter-sized fragments may increase to 1% within 36 hours, while the probability of collisions with fragments 5 cm or larger is approximately $10^{-5}$.These findings underscore the vulnerability of the geostationary region to space traffic accidents.

Many recent studies have pointed out significant discrepancies between observations and models of stellar populations in the near-infrared (NIR). With current and future observing facilities being focused in this wavelength range, properly assessing and solving these issues is of utmost importance. Here, we present the first application of the extragalactic globular cluster (GC) near-infrared spectroscopy survey, and present evidence that these GCs reveal an age zero-point problem of stellar population synthesis (SPS) models. This problem has already been identified in the optical range for the GCs of the Milky Way. Such an issue arises when derived GC spectroscopic ages appear older than the Universe itself. We extend this discussion for the first time to the NIR, specifically using the Pa$_{\rm\beta}$ line at 1.28~microns. We focus on the GCs of the nearby Centaurus~A galaxy using their NIR spectra. This work broadens our understanding of the age zero-point problem and emphasises the necessity to revisit and refine SPS models, especially in the NIR domain.

The origin of cool-core (CC) and non-cool-core (NCC) dichotomy of galaxy clusters remains uncertain. Previous simulations have found that cluster mergers are effective in destroying CCs but fail to prevent overcooling in cluster cores when radiative cooling is included. Feedback from active galactic nuclei (AGN) is a promising mechanism for balancing cooling in CCs; however, the role of AGN feedback in CC/NCC transitions remains elusive. In this work, we perform three-dimensional binary cluster merger simulations incorporating AGN feedback and radiative cooling, aiming to investigate the heating effects from mergers and AGN feedback on CC destruction. We vary the mass ratio and impact parameter to examine the entropy evolution of different merger scenarios. We find that AGN feedback is essential in regulating the merging clusters, and that CC destruction depends on the merger parameters. Our results suggest three scenarios regarding CC/NCC transitions: (1) CCs are preserved in minor mergers or mergers that do not trigger sufficient heating, in which cases AGN feedback is crucial for preventing the cooling catastrophe; (2) CCs are transformed into NCCs by major mergers during the first core passage, and AGN feedback is subdominant; (3) in major mergers with a large impact parameter, mergers and AGN feedback operate in concert to destroy the CCs.

We investigate cosmological constraints on local position invariance (LPI), a key aspect of the Einstein equivalence principle (EEP), through asymmetric galaxy clustering. The LPI asserts that the outcomes of the non-gravitational experiments are identical regardless of location in spacetime and has been tested through measurements of the gravitational redshift effect. Therefore, measuring the gravitational redshift effect encoded in galaxy clustering provides a powerful and novel cosmological probe of the LPI. Recent work by Saga et al. proposed its validation using the cross-correlation function between distinct galaxy samples, but their analysis focused solely on the dipole moment. In this paper, we extend their work by further analyzing a higher-order odd multipole moment, the octupole moment, in the constraints on the LPI-violating parameter, $\alpha$, expected from galaxy surveys such as Dark Energy Spectroscopic Instrument, Euclid space telescope, Subaru Prime Focus Spectrograph, and Square Kilometre Array. We demonstrate that combining the octupole and dipole moments significantly improves the constraints, particularly when the analysis is restricted to larger scales, characterized by a large minimum separation $s_{\rm min}$. For a conservative setup with $s_{\rm min}=15 {\rm Mpc}/h$, we find an average improvement of 11$\%$ compared to using the dipole moment alone. Our results highlight the importance of higher-order multipoles in constraining $\alpha$, providing a more robust approach to testing the EEP on cosmological scales.

We investigate the existence and properties of open superclusters within 500 pc of the Sun, with a particular focus on a newly identified OSC, designated HC8. Utilizing Gaia-derived advanced astrometric data and star cluster catalogs, we identify and analyze the member stars of various open clusters (OCs), deriving key parameters such as ages, distances, extinctions, and kinematic properties. Our findings reveal that HC8, a previously unrecognized OSC, encompasses several young clusters, including the newly discovered OC Duvia 1. We provide a detailed examination of HC8's star formation history and its spatial and kinematic characteristics. This study contributes to the growing body of evidence supporting the existence of primordial groups and enhances our understanding of the formation and evolution of star clusters in the Galaxy.

The pulsating variable star Mira (omikron Ceti) was observed by David Fabricius (Frisia) in 1596 and 1609. We review suggested previous detections (e.g. China, Hipparchos). We analyze all Mira records from Fabricius in their historical context. Fabricius measured the separation of Mira to other stars to \pm 1.6-1.7'. From his texts, we derive a brightness (slightly brighter than Hamal) of ca. 1.9 \pm 0.1 mag and a color index B-V \simeq 1.3-1.4 mag (`like Mars') for 1596 Aug 3 (jul.). Mira started to fainten 19 days later and was observed until mid/late Oct. We show why such a red star cannot be followed by the naked eye until ca. 6 mag: For Mira's color at disappearance and altitude from Frisia, the limit is reduced by ca. 1.0 mag. Since Fabricius connected the Mira brightening with the close-by prograde Jupiter, he re-detected it only 12 years later, probably shortly before a relatively bright maximum - discoveries are strongly affected by biases. A Mira period of 330.2 days is consistent with both the oldest data (from Fabricius 1596 to Hevelius 1660) and the most current data (VSX 2004-2023), so that we see no evidence for secular period or phase shifts. (We also present Fabricius' observations of P Cygni in 1602.)

Pre-stellar cores are the first steps in the process of star and planet formation. However, the dynamical and chemical evolution of pre-stellar cores is still not well understood. We aim at estimating the central density of the pre-stellar core IRAS16293E and at carrying out an inventory of molecular species towards the density peak of the core. We observed high-$J$ rotational transitions of N$_2$H$^+$ and N$_2$D$^+$, and several other molecular lines towards the dust emission peak using the Atacama Pathfinder EXperiment (APEX) telescope, and derived the density and temperature profiles of the core using far-infrared surface brightness maps from $Herschel$. The N$_2$H$^+$ and N$_2$D$^+$ lines were analysed by non-LTE radiative transfer modelling. Our best-fit core model consists in a static inner region, embedded in an infalling envelope with an inner radius of approximately 3000 au (21" at 141 pc). The observed high-J lines of N$_2$H$^+$ and N$_2$D$^+$ (with critical densities greater than 10$^6$ cm$^{-3}$) turn out to be very sensitive to depletion; the present single-dish observations are best explained with no depletion of N$_2$H$^+$ and N$_2$D$^+$ in the inner core. The N$_2$D$^+$/N$_2$H$^+$ ratio that best reproduces our observations is 0.44, one of the largest observed to date in pre-stellar cores. Additionally, half of the molecules that we observed are deuterated isotopologues, confirming the high-level of deuteration towards this source. Non-LTE radiative transfer modelling of N$_2$H$^+$ and N$_2$D$^+$ lines proved to be an excellent diagnostic of the chemical structure and dynamics of a pre-stellar core. Probing the physical conditions immediately before the protostellar collapse is a necessary reference for theoretical studies and simulations with the aim of understanding the earliest stages of star and planet formation and the time scale of this process.

Flares produced by certain classes of astrophysical objects may be sources of some ultra-high-energy particles, which, if they are photons, would group into clusters of events correlated in space and time. Identification of such clustering in cosmic-ray data would provide important evidence for possible existence of ultra-high-energy (UHE) photons and could potentially help identify their sources. We present an analysis method to search for space-time clustering of ultra-high-energy extensive air showers, namely the stacking method, which combines a time-clustering algorithm with an unbinned likelihood study. In addition, to enhance the capability to discriminate between signal (photon-initiated events) and background (hadron-initiated) events, we apply a photon tag. This involves using relevant probability distribution functions to classify each event as more likely to be either a photon or a hadron. We demonstrate that the stacking method can effectively distinguish between events initiated by photons and those initiated by hadrons (background). The number of photon events in a data sample, as well as the flare(s) duration can also be retrieved correctly. The stacking method with a photon tag requires only a few events to identify a photon flare. This method can be used to search for the cosmic ray sources and/or improve limits on the fluxes of UHE photons.

The stellar haloes of dwarf galaxies are becoming an object of interest in the extragalactic community due to their detection in some recent observations. Additionally, new cosmological simulations of very high resolution were performed, allowing their study. These stellar haloes could help shed light on our understanding of the assembly of dwarf galaxies and their evolution, and allow us to test the hierarchical model for the formation of structures at small scales. We aim to characterise the stellar haloes of simulated dwarf galaxies and analyse their evolution and accretion history. We use a sample of 17 simulated galaxies from the Auriga Project with a stellar mass range from 3.28x10^8 Msun to 2.08x10^10 Msun. We define the stellar halo as the stellar material located outside an ellipsoid with semi-major axes equal to 4 times the half light radius (Rh) of each galaxy. We find that the inner regions of the stellar halo (4 to 6 times the Rh) are dominated by in-situ material. For the less massive simulated dwarfs (M*<=4.54x10^8 Msun), this dominance extends to all radii. We find that this in-situ stellar halo is mostly formed in the inner regions of the galaxies and then ejected into the outskirts during interactions and merger events. In ~50% of the galaxies, the stripped gas from satellites contributed to the formation of this in-situ halo. The stellar haloes of the galaxies more massive than M*>=1x10^9 Msun are dominated by the accreted component beyond 6 Rh. We find that the more massive dwarf galaxies accrete stellar material until later times (t90~4.44 Gyr ago, being t90 the formation time) than the less massive ones (t90~8.17 Gyr ago), impacting on the formation time of the accreted stellar haloes. The galaxies have a range of 1 to 7 significant progenitors contributing to their accreted component but there is no correlation between this quantity and the galaxies' accreted mass.

We report an observed accretion rate of $\dot M_1 = (3.86\pm0.60)\times 10^{-11}$ $M_{\odot}$yr$^{-1}$ for the white dwarf in the short-period, intermediate polar EX Hya. This result is based upon the accretion-induced $4\pi$-averaged energy flux from 2.45 $\mu$m to 100 keV and the corresponding luminosity at the Gaia distance of 56.77 pc. Our result is in perfect agreement with the theoretical mass transfer rate from the secondary star induced by gravitational radiation (GR) and the spin-up of the white dwarf, $-\dot M_2 = (3.90\pm0.35)\times 10^{-11}$ $M_{\odot}$yr$^{-1}$; 24% of it is caused by the spin-up. The agreement indicates that mass transfer is conservative. The measured $\dot M_1$ obviates the need for angular momentum loss (AML) by any process other than GR. We complemented this result with an estimate of the mean secular mass transfer rate over $\sim 10^7$ yr by interpreting the non-equilibrium radius of the secondary star in EX Hya based on published evolutionary calculations. This suggests a time-averaged mass transfer rate enhanced over GR by a factor $f_{\mathrm{GR}} \gtrsim 2$. Combined with the present-day lack of such an excess, we suggest that an enhanced secular AML is due to an intermittently active process, such as the proposed frictional motion of the binary in the remnants of nova outbursts. We argue that EX Hya, despite its weakly magnetic nature, has evolved in a very similar way to non-magnetic CVs. We speculate that the discontinuous nature of an enhanced secular AML may similarly apply to the latter.

The strong gravitational potential of neutron stars (NSs) makes them ideal astrophysical objects for testing extreme gravity phenomena. We explore the potential of NS X-ray pulsed lightcurve observations to probe deviations from general relativity (GR) within the scalar-tensor theory (STT) of gravity framework. We compute the flux from a single, circular, finite-size hot spot, accounting for light bending, Shapiro time delay, and Doppler effect. We focus on the high-compactness regime, i.e., close to the critical GR value GM/(Rc^2) = 0.284, over which multiple images of the spot appear and impact crucially the lightcurve. Our investigation is motivated by the increased sensitivity of the pulse to the scalar charge of the spacetime in such high compactness regimes, making these systems exceptionally suitable for scrutinizing deviations from GR, notably phenomena such as spontaneous scalarization, as predicted by STT. We find significant differences in NS observables, e.g., the flux of a single spot can differ up to 80% with respect to GR. Additionally, reasonable choices for the STT parameters that satisfy astrophysical constraints lead to changes in the NS radius relative to GR of up to approximately 10%. Consequently, scalar parameters might be better constrained when uncertainties in NS radii decrease, where this could occur with the advent of next-generation gravitational wave detectors, such as the Einstein Telescope and LISA, as well as future electromagnetic missions like eXTP and ATHENA. Thus, our findings suggest that accurate X-ray data of the NS surface emission, jointly with refined theoretical models, could constrain STTs.

Galaxy cluster masses estimated from parametric modeling of weak lensing shear observations are known to be biased by inaccuracies in observationally determined centers. It has recently been shown that such systematic effects can be non-isotropic when centers are derived from X-ray or Compton-Y (Sunyaev-Zeldovich effect) observations, which is often the case in practice. This fact challenges current methods of accurately correcting for weak lensing mass biases using simulations paired with isotropic empirical miscentering distributions, in particular as the effect on determined masses is currently a dominant source of systematic uncertainty. We use hydrodanamical cosmological simulations taken from the Magneticum Pathfinder simulations to show that the non-isotropic component of the mass bias can be reduced to within one percent of the mass when considering the center of mass, rather than the bottom of the gravitational potential, as the reference center of a galaxy cluster.

The cosmic microwave background (CMB) and baryon acoustic oscillations (BAO) provide precise measurements of the cosmic expansion history through the comoving acoustic scale. The CMB angular scale measurement $\theta_*$ is particularly robust, constraining the ratio of the sound horizon to the angular diameter distance to last scattering independently of the late-time cosmological model. For models with standard early-universe physics, this measurement strongly constrains possible deviations from $\Lambda$CDM at late times. We show that the null energy condition imposes strict inequalities on the BAO observables $D_H(z)$, $D_M(z)$, $D_V(z)$ and $F_{\rm AP}(z)$ relative to $\Lambda$CDM predictions. These inequalities demonstrate that certain deviations from $\Lambda$CDM are impossible for any physical dark energy model that respects the null energy condition. We also identify the regions of parameter space in the CPL parameterization $w(a) = w_0 + w_a(1-a)$ that can give predictions consistent with both the null energy condition and the observed CMB scale. While current DESI DR1 BAO measurements exhibit slight joint-constraint parameter tensions with $\Lambda$CDM, this tension only arises in directions that are inconsistent with the null-energy condition, so $\Lambda$CDM is favoured by acoustic scale measurements unless the null-energy condition is violated.

The gradient and curvature of the Parker spiral interplanetary magnetic field give rise to curvature and gradient guiding centre drifts on cosmic rays. The plasma turbulence present in the interplanetary space is thought to suppress the drifts, however the extent to which they are reduced is not clear. We investigate the reduction of the drifts using a new analytic model of heliospheric turbulence where the dominant 2D component has both the wave vector and the magnetic field vector normal to the Parker spiral, thus fulfilling the main criterion of 2D turbulence. We use full-orbit test particle simulations of energetic protons in the modelled interplanetary turbulence, and analyse the mean drift velocity of the particles in heliolatitude. We release energetic proton populations of 10, 100 and 1000~MeV close to Sun and introduce a new method to assess their drift. We compare the drift in the turbulent heliosphere to drift in a configuration without turbulence, and to theoretical estimates of drift reduction. We find that drifts are reduced by a factor 0.2-0.9 of that expected for the heliospheric configuration without turbulence. This corresponds to a much less efficient suppression than what is predicted by theoretical estimates, particularly at low proton energies. We conclude that guiding centre drifts are a significant factor for the evolution of cosmic ray intensities in the heliosphere including the propagation of solar energetic particles in the inner heliosphere.

Different types of hot subdwarfs may have different origins, which will cause them to present different radial velocity (RV) variability properties. Only 6$\pm$4% of our single-lined He-rich hot subdwarfs that only show spectroscopic features of hot subdwarfs are found to be RV variable, which is lower than the fraction of single-lined He-poor sdB stars (31$\pm$3%). Single-lined sdB stars with effective temperatures ($T_{\rm eff}$) $\sim$ 25,000 $-$ 33,000 K show an RV-variability fraction of 34$\pm$5%, while lower RV-variability fractions are observed for single-lined sdB stars cooler than about 25,000 K (11$\pm$4%), single-lined sdB/OB stars with $T_{\rm eff}$ $\sim$ 33,000 $-$ 40,000 K and surface gravities about 5.7 $-$ 6.0 (13$\pm$3%), as well as single-lined sdO/B stars with $T_{\rm eff}$ $\sim$ 45,000 $-$ 70,000 K (10$\pm$7%). Single-lined hot subdwarfs with $T_{\rm eff}$ $\sim$ 35,000 $-$ 45,000 K located above the extreme horizontal branch (EHB) show a similar RV-variability fraction of 34$\pm$9% as single-lined sdB stars at about 25,000 $-$ 33,000 K. The largest RV-variability fraction of 51$\pm$8% is found in single-lined hot subdwarfs below the canonical EHB. The detected RV-variability fraction of our composite hot subdwarfs with an infrared excess in their spectral energy distributions is 9$\pm$3%, which is lower than that fraction of single-lined hot subdwarfs. Since the average RV uncertainty we measured in the LAMOST spectra is about 7.0 km/s, the lower detected RV-variability fraction for composite hot subdwarfs is expected because the RV amplitudes associated with long-period systems are lower.

Here we explore certain subtle features imprinted in data from the completed Sloan Digital Sky Survey IV (SDSS-IV) extended Baryon Oscillation Spectroscopic Survey (eBOSS) as a combined probe for the background and perturbed Universe. We reconstruct the baryon Acoustic Oscillation (BAO) and Redshift Space Distortion (RSD) observables as functions of redshift, using measurements from SDSS alone. We apply the Multi-Task Gaussian Process (MTGP) framework to model the interdependencies of cosmological observables $D_M(z)/r_d$, $D_H(z)/r_d$, and $f\sigma_8(z)$, and track their evolution across different redshifts. Subsequently, we obtain constrained three-dimensional phase space containing $D_M(z)/r_d$, $D_H(z)/r_d$, and $f\sigma_8(z)$ at different redshifts probed by the SDSS-IV eBOSS survey. Furthermore, assuming the $\Lambda$CDM model, we obtain constraints on model parameters $\Omega_{m}$, $H_{0}r_{d}$, $\sigma_{8}$ and $S_{8}$ at each redshift probed by SDSS-IV eBOSS. This indicates redshift-dependent trends in $H_0$, $\Omega_m$, $\sigma_8$ and $S_8$ in the $\Lambda$CDM model, suggesting a possible inconsistency in the $\Lambda$CDM model. Ours is a template for model-independent extraction of information for both background and perturbed Universe using a single galaxy survey taking into account all the existing correlations between background and perturbed observables and this can be easily extended to future DESI-3YR as well as Euclid results.

The propagation of disturbances in the solar atmosphere is inherently three dimensional (3D), yet comprehensive studies on the spatial structure and dynamics of 3D wavefronts are scarce. Here we conduct high resolution 3D numerical simulations to investigate filament eruptions, focusing particularly on the 3D structure and genesis of EUV waves. Our results demonstrate that the EUV wavefront forms a dome like configuration subdivided into three distinct zones. The foremost zone, preceding the flux rope, consists of fast-mode shock waves that heat the adjacent plasma. Adjacent to either side of the flux rope, the second zone contains expansion waves that cool the nearby plasma. The third zone, at the juncture of the first two, exhibits minimal disturbances. This anisotropic structure of the wavefront stems from the configuration and dynamics of the flux rope, which acts as a 3D piston during eruptions :compressing the plasma ahead to generate fast mode shocks and evacuating the plasma behind to induce expansion waves. This dynamic results in the observed anisotropic this http URL, with synthetic EUV images from simulation data, the EUV waves are observable in Atmospheric Imaging Assembly 193 and 211 angstrom, which are identified as the fast mode shocks. The detection of EUV waves varies with the observational perspective: the face on view reveals EUV waves from the lower to the higher corona, whereas an edge on view uncovers these waves only in the higher corona.

Motivated by observed discrepancies between ACT DR4 and Planck 2018 cosmic microwave background (CMB) anisotropy power spectra, particularly in the cross-correlation of temperature and E-mode polarization, we investigate challenges that may be encountered in the comparison of satellite and ground-based CMB data. In particular, we focus on the effects of Fourier-space filtering and masking involving bright point sources. We show that the filtering operation generates bright cross-shaped artifacts in the map, which stretch far outside typical point-source masks. If not corrected, these artifacts can add bias or additional variance to cross-spectra, skewing results. However we find that the effect of this systematic is not large enough to explain the ACT-Planck differences presented with ACT DR4.

High-precision astrometry offers a promising approach to detect low-frequency gravitational waves, complementing pulsar timing array (PTA) observations. We explore the response of astrometric measurements to a stochastic gravitational wave background (SGWB) in synergy with PTA data. Analytical, covariant expressions for this response are derived, accounting for the presence of a possible dipolar anisotropy in the SGWB. We identify the optimal estimator for extracting SGWB information from astrometric observations and examine how sensitivity to SGWB properties varies with the sky positions of stars and pulsars. Using representative examples of current PTA capabilities and near-future astrometric sensitivity, we demonstrate that cross-correlating astrometric and PTA data can improve constraints on SGWB properties, compared to PTA data alone. The improvement is quantified through Fisher forecasts for the SGWB amplitude, spectral tilt, and dipolar anisotropy amplitude. In the future, such joint constraints could play a crucial role in identifying the origin of SGWB signals detected by PTAs.

We investigate how the multiplicity of binary, triple and quadruple star systems changes as the systems evolve from the zero-age main-sequence to the Hubble time. We find the change in multiplicity fractions over time for each data set, identify the number of changes to the orbital configuration and the dominant underlying physical mechanism responsible for each configuration change. Finally, we identify key properties of the binaries which survive the evolution. We use the stellar evolution population synthesis code Multiple Stellar Evolution (MSE) to follow the evolution of $3 \times 10^4$ of each 1+1 binaries, 2+1 triples, 3+1 quadruples and 2+2 quadruples. The coupled stellar and orbital evolution are computed each iteration. The systems are assumed to be isolated and to have formed in situ. We generate data sets for two different black hole natal kick mean velocity distributions (sigma = 10 km/s and sigma = 50 km/s and with and without the inclusion of stellar fly-bys. Our fiducial model has a mean black hole natal kick velocity if sigma = 10 km/s and includes stellar fly-bys. Each system has at least one star with an initial mass larger than 10 solar masses. All data will be publicly available. We find that at the end of the evolution the large majority of systems are single stars in every data set (> 85%). As the number of objects in the initial system increases, so too does the final non-single system fraction. The single fractions of final systems in our fiducial model are 87.8 $\pm$ 0.2 % for the 2 + 2s, 88.8 $\pm$ 0.3 % for the 3 + 1s, 92.3 $\pm$ 0.2 % for the 2 + 1s and 98.9 $\pm$ 0.3 % for the 1 + 1s.

Due to the large dynamic ranges involved with separating the cosmological 21-cm signal from the Cosmic Dawn from galactic foregrounds, a well-calibrated instrument is essential to avoid biases from instrumental systematics. In this paper we present three methods for calibrating a global 21-cm cosmology experiment using the noise wave parameter formalisation to characterise a low noise amplifier including a careful consideration of how calibrator temperature noise and singularities will bias the result. The first method presented in this paper builds upon the existing conjugate priors method by weighting the calibrators by a physically motivated factor, thereby avoiding singularities and normalising the noise. The second method fits polynomials to the noise wave parameters by marginalising over the polynomial coefficients and sampling the polynomial orders as parameters. The third method introduces a physically motivated noise model to the marginalised polynomial method. Running these methods on a suite of simulated datasets based on the REACH receiver design and a lab dataset, we found that our methods produced a calibration solution which is equally as or more accurate than the existing conjugate priors method when compared with an analytic estimate of the calibrator's noise. We find in the case of the measured lab dataset the conjugate priors method is biased heavily by the large noise on the shorted load calibrator, resulting in incorrect noise wave parameter fits. This is mitigated by the methods introduced in this paper which calibrate the validation source spectra to within 5% of the noise floor.

Titan's abundant atmospheric N2 and CH4 gases are notable characteristics of the moon that may help constrain its origins and evolution. Previous work suggests that atmospheric CH4 is lost on geologically short timescales and may be replenished from an interior source. Isotopic and noble gas constraints indicate that N2 may derive from a mixture of NH3 ice and heating of organic matter. Here, we report experimental results from hydrothermal alteration of insoluble organic matter from the Murchison meteorite and analog insoluble organic matter at temperatures and pressures that are relevant to Titan's interior. Our results indicate both CH4 and CO2 are formed, with the ratio between the two depending on a multitude of factors, particularly temperature and, to a lesser degree, the dielectric constant of water and carbonyl abundance in the starting material. Sufficient CH4 is produced to source Titan's atmospheric reservoir if temperatures are greater than 250°C. Nitrogen is volatilized, primarily in the form of NH3, in sufficient abundances to source at least 50% of Titan's atmospheric N2. The isotopic characteristics of volatilized material relative to the starting organics are consistent with current constraints for the nature of the accreted complex organics and Titan's evolved atmosphere.

In this work, we investigate the generation of primordial black holes (PBHs) within the framework of single-field inflationary models and their compatibility with the cosmological history of the Universe. Our results suggest that, depending on the masses of the formed PBHs, single-field inflation models require more than fine-tuning a potential to induce ultra-slow roll; it necessitates a comprehensive understanding of the post-inflationary cosmological evolution. As an explicative example, we introduce a new model, based on a double inflection point and consistent with Cosmic Microwave Background observations, capable of generating sub-solar PBHs, whose merger could be potentially detectable by the LVK experiment.

Exoplanets on close-in orbit are subject to intense X-ray and ultraviolet (XUV) irradiation from their star. Their atmosphere therefore heats up, sometimes to the point where it thermally escapes from the gravitational potential of the planet. Nonetheless, XUV is not the only source of heating in such atmospheres. Indeed, close-in exoplanets are embedded in a medium (the stellar wind) with strong magnetic fields that can significantly vary along the orbit. The variations of this magnetic field can induce currents in the upper atmosphere, which dissipate and locally heat it up through Ohmic heating. The aim of this work is to quantify Ohmic heating in the upper atmosphere of hot exoplanets due to an external time-varying magnetic field, and to compare it to the XUV heating. Ohmic heating depends strongly on the conductivity properties of the upper atmosphere. A 1D formalism is developed to assess the level and the localization of Ohmic heating depending on the conductivity profile, and applied to the specific cases of Trappist-1 b and $\pi$ Men c. Ohmic heating can reach values up to 10$^{-3}$ erg s$^{-1}$ cm$^{-3}$ in the upper atmospheres of hot exoplanets. It is expected to be stronger the closer the planet is and the lower the central star mass is, as these conditions maximize the strength of the ambient magnetic field around the planet. We confirm that Ohmic heating can play an important role in setting the thermal budget of the upper atmosphere of hot exoplanets, and can even surpass the XUV heating in the most favorable cases. When it is strong, a corollary is that the upper atmosphere screens efficiently time-varying external magnetic fields, preventing them to penetrate deeper in the atmosphere or inside the planet itself. We find that both Trappist-1b and $\pi$ Men c are likely subject to intense Ohmic heating.

We study the non-Gaussian tail of the curvature fluctuation, $\zeta$, in an inflationary scenario with a transient ultra slow-roll phase that generates a localized large enhancement of the spectrum of $\zeta$. To do so, we implement a numerical procedure that provides the probability distribution of $\zeta$ order by order in perturbation theory. The non-Gaussianities of $\zeta$ can be shown to arise from its non-linear relation to the inflaton fluctuations and from the intrinsic non-Gaussianities of the latter, which stem from its self interactions. We find that intrinsic non-Gaussianities, which have often been ignored to estimate the abundance of primordial black holes in this kind of scenario, are important. The relevance of the intrinsic contribution depends on the rapidity with which the transient ultra slow-roll phase occurs, as well as on its duration. Our method cannot be used accurately when the perturbative in-in formalism fails to apply, highlighting the relevance of developing fully non-perturbative approaches to the problem.

There is increasing evidence that single-star evolutionary models are inadequate to reproduce all observational properties of massive stars. Binary interaction has emerged as a key factor in the evolution of a significant fraction of massive stars. In this study, we investigate the helium ($Y_{\mathrm He}$) and nitrogen ($\epsilon_{\mathrm N}$) surface abundances in a comprehensive sample of 180 Galactic O-type stars with projected rotational velocities $v\sin(i)\leq150{\mathrm km}\cdot{\mathrm s}^{-1}$. We found a subsample ($\sim20\%$ of the total, and $\sim80\%$ of the stars with $Y_{\mathrm He}\geq0.12$) with a $Y_{\mathrm He}$ and $\epsilon_{\mathrm N}$ combined pattern unexplainable by single-star evolution. We argue that the stars with anomalous surface abundance patterns are binary interaction products.

We investigate a cosmological model inspired by hybrid inflation, where two scalar fields representing dark energy (DE) and dark matter (DM) interact through a coupling that is proportional to the DE scalar field $1/\phi$. The strength of the coupling is governed solely by the initial condition of the scalar field, $\phi_i$, which parametrises deviations from the standard $\Lambda$CDM model. In this model, the scalar field tracks the behaviour of DM during matter-domination until it transitions to DE while the DM component decays quicker than standard CDM during matter-domination, and is therefore different from some interacting DM-DE models which behaves like phantom dark energy. Using \textit{Planck} 2018 CMB data, DESI BAO measurements and Pantheon+ supernova observations, we find that the model allows for an increase in $H_0$ that can help reduce the Hubble tension. In addition, we find that higher values of the coupling parameter are correlated with lower values of $\omega_m$, and a mild decrease of the weak-lensing parameter $S_8$, potentially relevant to address the $S_8$ tension. Bayesian model comparison, however, reveals inconclusive results for most datasets, unless S$H_0$ES data are included, in which case a moderate evidence in favour of the hybrid model is found.

Observations of strongly gravitationally lensed gravitational wave (GW) sources provide a unique opportunity for constraining their transverse motion, which otherwise is exceedingly hard for GW mergers in general. Strong lensing makes this possible when two or more images of the lensed GW source are observed, as each image essentially allows the observer to see the GW source from different directional lines-of-sight. If the GW source is moving relative to the lens and observer, the observed GW signal from one image will therefore generally appear blue- or redshifted compared to GW signal from the other image. This velocity induced differential Doppler shift gives rise to an observable GW phase shift between the GW signals from the different images, which provides a rare glimpse into the relative motion of GW sources and their host environment across redshift. We illustrate that detecting such GW phase shifts is within reach of next-generation ground-based detectors such as Einstein Telescope, that is expected to detect $\sim$hundreds of lensed GW mergers per year. This opens up completely new ways of inferring the environment of GW sources, as well as studying cosmological velocity flows across redshift.

Accurate modeling of how high-energy proton-proton collisions produce gamma rays through the decays of pions and other secondaries is needed to correctly interpret astrophysical observations with the Fermi-LAT telescope. In the existing literature on cosmic-ray collisions with gas, the focus is on the gamma-ray yield spectrum, $d N_\gamma/dE$. However, in some situations, the joint energy and angular distribution can be observed, so one needs instead $d^2 N_\gamma/dE \, d\Omega$. We provide calculations of this distribution over the energy range from the pion production threshold to $100~{\rm GeV}$, basing our results on FLUKA simulations. We provide the results in tabular form and provide a Python tool on GitHub to aid in utilization. We also provide an approximate analytic formula that illuminates the underlying physics. We discuss simplified examples where this angular dependence can be observed to illustrate the necessity of taking the joint distribution into account.

I derive a lower limit on the mass of an Unidentified Flying Object (UFO) based on measurements of its speed and acceleration, as well as the infrared luminosity of the airglow around it. If the object's radial velocity can be neglected, the mass limit is independent of distance. Measuring the distance and angular size of the object allows to infer its minimum mass density. The Galileo Project will be collecting the necessary data on millions of objects in the sky over the coming year.

Heavy axions that couple to both quantum electrodynamics and quantum chromodynamics with masses on the order of MeV - GeV and high-scale decay constants in excess of $\sim$$10^8$ GeV may arise generically in e.g. axiverse constructions. In this work we provide the most sensitive search to-date for the existence of such heavy axions using Fermi-LAT data towards four recent supernovae (SN): Cas A, SN1987A, SN2023ixf, and SN2024ggi. We account for heavy axion production in the proto-neutron-star cores through nuclear and electromagnetic processes and then the subsequent decay of the axions into photons. While previous works have searched for gamma-rays from SN1987A using the Solar Maximum Mission that observed SN1987A during the SN itself, we show that using Fermi Large Area Telescope data provides an approximately five orders of magnitude improvement in flux sensitivity for axions with lifetimes larger than around 10 years. We find no evidence for heavy axions and exclude large regions of previously-unexplored parameter space.

We explore a purely gravitational origin of observed baryon asymmetry and dark matter (DM) abundance from asymmetric Hawking radiation of light primordial black holes (PBH) in presence of a non-zero chemical potential, originating from the space-time curvature. Considering the PBHs are described by a Reissner-Nordström metric, and are produced in a radiation dominated Universe, we show, it is possible to simultaneously explain the matter-antimatter asymmetry along with right DM abundance satisfying bounds from big bang nucleosynthesis, cosmic microwave background and gravitational wave energy density due to PBH density fluctuation. We also obtain the parameter space beyond the semiclassical approximation, taking into account the quantum effects on charged PBH dynamics due to memory burden.

We conduct a systematic investigation of freeze-in during reheating while taking care to include both direct and indirect production of dark matter (DM) via gravitational portals and inflaton decay. Direct production of DM can occur via gravitational scattering of the inflaton, while indirect production occurs through scattering in the Standard Model radiation bath. We consider two main contributions to the radiation bath during reheating. The first, which may dominate at the onset of the reheating process, is produced via gravitational scattering of the inflaton. The second (and more standard contribution) comes from inflaton decay. We consider a broad class of DM production rates parameterized as $R_{\chi} \propto T^{n+6}/\Lambda^{n+2}$, and inflaton potentials with a power-law form $V(\phi) \propto \phi^{k}$ about the minimum. We find the relic density produced by freeze-in for each contribution to the Standard Model bath for arbitrary $k$ and $n$, and compare these with the DM density produced gravitationally by inflaton scattering. We find that freeze-in production from the gravitationally-produced radiation bath can exceed that of the conventional decay bath and account for the observed relic density provided that $m_{\chi} > T_{\rm RH}$, with additional $k$- and $n$-dependent constraints. For each freeze-in interaction considered, we also find $m_{\chi}$- and $T_{\rm RH}$-dependent limits on the BSM scale, $\Lambda$, for which gravitational production will exceed ordinary freeze-in production.

The current generation of Dark Matter Direct Detection Experiments has ruled out a large region of parameter space for dark matter, particularly in the ($10 - 1000$) GeV mass range. However, due to very low event rates, searching for dark matter in the heavy mass range, $\mathcal{O}$(TeV), is a daunting task requiring even larger volume detectors and long exposure times. We show that for a broad class of dark matter models of the type that these experiments are searching, including some of the most popular candidates, the heavy dark matter mass range can be ruled out in its entirety once we take into account the large corrections to Higgs mass imparted by such heavy dark matter. We show that such a limit is applicable to all types of dark matter i.e. scalar, vector, and fermionic, provided they couple directly with Higgs. By taking some simple and well studied dark matter models we show that the latest LZ limits can completely rule out such a dark matter except in a narrow range around $M_h/2$ mass.

In this paper, we investigate the geodesic motion of test particles in the spacetime surrounding a static, spherically symmetric black hole, which is described by an AdS-Schwarzschild-like metric and incorporates a quantum correction. This black hole also features phantom global monopoles, which modify the structure of the black hole space-time. We begin by deriving the effective potential governing the motion of test particles in this system and carefully analyze the impact of quantum correction in the presence of both phantom and ordinary global monopoles. Furthermore, we extend our study to include the spin-dependent Regge-Wheeler (RW) potential, which characterizes the dynamics of perturbations in this quantum-corrected black hole background. By examining this RW potential for various spin fields, we show how quantum corrections affect its form in the presence of both phantom and ordinary global monopoles. Our analysis demonstrate that quantum correction significantly alter the nature of the RW-potential, influencing the stability, and behavior of test particles and perturbations around the black hole.

In general the speed of Gravitational Waves (GWs) in Scalar-Tensor modifications of Einstein's gravity is different from the speed of Light. Nevertheless, it has been measured that their speeds are nearly the same. For the most general Scalar-Tensor theories classified to date that do propagate a graviton -- DHOST, including Horndeski and Beyond Horndeski (BH) theories -- we show that, remarkably, up to 5 self-consistent couplings of the scalar of Dark Energy (DE) to the Photon are enough to make their GWs luminal in a wide set of cases. There is at least a Luminal BH theory for which the GW decay into DE is suppressed in any cosmological background.

We analyze the effect of a Chaplygin dark fluid (CDF) core on neutron stars (NSs). To address this study, we focus on the relativistic structure of stellar configurations composed by a dark-energy core, described by a Chaplygin-like equation of state (EoS), and an ordinary-matter crust which is described by a polytropic EoS. We examine the impact of the rate of energy densities at the discontinuous surface, defined as $\alpha= \rho_{\rm dis}^-/\rho_{\rm dis}^+$, on the radius, total gravitational mass, oscillation spectrum and tidal deformability. Furthermore, we compare our theoretical predictions with several observational mass-radius measurements and tidal deformability constraints. These comparisons together with the radial stability analysis show that the existence of NSs with a dark-energy core is possible.

General Relativity suffers for two main problems which have not yet been overcome: it predicts spacetime singularities and cannot be formulated as a perturbative renormalizable theory. In particular, many attempts have been made for avoiding singularities, such as considering higher order or infinite derivative theories. The price to pay in both cases is to give up locality and therefore they are known altogether as non-local theories of gravity. In this paper, we investigate how to recognize the presence of non-local effects by exploiting the power emitted by gravitational waves in a binary system in presence of non-local corrections as $R\Box^{-1}R$ to the Hilbert-Einstein action. After solving the field equations in terms of the source stress-energy tensor $T_{\mu\nu}$ and obtaining the gravitational wave stress-energy pseudo-tensor, $\tau_{\mu\nu}$, we find that the General Relativity quadrupole formula is modified in a non-trivial way, making it feasible to find a possible signature of non-locality. Our final results on the gravitational wave stress-energy pseudo-tensor could also be applied to several astrophysical scenarios involving energy or momentum loss, potentially providing multiple tests for non-local deviations from General Relativity. We finally discuss the detectability of the massless transverse scalar mode, discovering that, although this radiation is extremely weak, in a small range around the model divergence, its amplitude could fall within the low-frequency Einstein Telescope sensitivity.

On April 1, 2471 bC an impressive, unpredictable phenomenon occurred over the Delta of the Nile: a total solar eclipse, with totality band almost centered on the sacred city of Buto, and with the capital Memphis very close (>95%) to totality. This date is compatible with existing chronologies for the reign of Pharaoh Shepseskaf, who adopted a clamorous symbolic break with respect to the tradition of solarized kings started by Khufu. Indeed his tomb was not built in view from Heliopolis and was not a pyramid, but a kind of unique monument resembling the symbolic shrine at Buto. The aim of the present paper is to investigate in a systematic way the possibility that the origin of this historical break, which marks the end of the 4th dynasty, can be identified precisely in the 2471 bC eclipse, therefore furnishing a new astronomical anchor for the chronology of the Old Kingdom.

In previous work we have developed a model-independent, effective description of quantum deformed, spherically symmetric and static black holes in four dimensions. The deformations of the metric are captured by two functions of the physical distance to the horizon, which are provided in the form of self-consistent Taylor series expansions. While this approach efficiently captures physical observables in the immediate vicinity of the horizon, it is expected to encounter problems of convergence at further distances. Therefore, we demonstrate in this paper how to use Padé approximants to extend the range of applicability of this framework. We provide explicit approximations of physical observables that depend on finitely many effective parameters of the deformed black hole geometry, depending on the order of the Padé approximant. By taking the asymptotic limit of this order, we in particular provide a closed-form expression for the black hole shadow of the (fully) deformed geometry, which captures the leading quantum corrections. We illustrate our results for a number of quantum black holes previously proposed in the literature and find that our effective approach provides excellent approximations in all cases.

We present an innovative charge detector with high resolution and wide dynamic range designed to fulfill the requirements of a monitoring system for a high energy ion beam. The detector prototype, constructed using Si photodiodes and a custom readout electronics, underwent extensive testing during HERD and AMS beam tests at CERN SPS facilities. Initial testing showcased the detector's exceptional performance, emphasizing both high resolution and a dynamic range capable of measuring nuclei with atomic numbers ranging from 1 to 80. The prototype's compatibility with fast, quasi real-time data analysis qualifies it as an ideal candidate for online applications. This article presents the results from the testing phase of the prototype, highlighting its capabilities and performance. Ongoing detector development, potential applications, and future developments aimed at enhancing the detector's functionality and versatility are also discussed.

The timeline of the expansion rate ultimately defines the interplay between high energy physics, astrophysics and cosmology. The guiding theme of this topical review is provided by the scrutiny of the early history of the space-time curvature through the diffuse backgrounds of gravitational radiation that are sensitive to all the stages of the evolution of the plasma. Due to their broad spectrum (extending from the aHz region to the THz domain) they bridge the macroworld described by general relativity and the microworld of the fundamental constituents of matter. It is argued that during the next score year the analysis of the relic gravitons may infirm or confirm the current paradigm where a radiation plasma is assumed to dominate the whole post-inflationary epoch. The role of high frequency and ultra-high frequency signals between the MHz and the THz is emphasized in the perspective of quantum sensing. The multiparticle final state of the relic gravitons and its macroscopic quantumness is also discussed with particular attention to the interplay between the entanglement entropy and the maximal frequency of the spectrum.

We discuss a model of neutrino mass based on the type I seesaw mechanism embedded in a spontaneously broken global lepton number framework with a $Z_2$ symmetry. We show that the resulting Majoron is a viable freeze-in dark matter candidate. Two right-handed neutrinos are assumed to have dominant off-diagonal masses suggesting resonant leptogenesis as the origin of baryon asymmetry of the Universe. Explicit higher dimensional lepton number violating operators, are shown to play a crucial role in simultaneously controlling both the Majoron production in the early Universe and the right handed neutrino mass splitting relevant for resonant leptogenesis. We perform a combined analysis of Majoron dark matter and leptogenesis, discussing the relative importance of self energy and vertex contributions to CP asymmetry, and explore the parameter space, leading to an intricate relation between neutrino mass, dark matter and baryon asymmetry.

We study the deflection of light rays in a cold, non-magnetized plasma using the worldline framework. Starting from Synge's Hamiltonian formalism, we construct a position-space action and use it perturbatively to calculate light bending angles. In the homogeneous case, the action reduces to that of a massive particle, allowing us to extract the bending angle of light in the presence of the medium using a well-known analogy. For the inhomogeneous case, we consider a power law model and construct Feynman rules in time to compute the purely plasma-induced corrections to the bending angle at Next-to-Leading-Order (NLO).