Papers for Monday, Dec 09 2024

votess is a library for computing parallel 3D Voronoi tessellations on heterogeneous platforms, from CPUs and GPUs, to future accelerator architectures. To do so, it leverages the SYCL abstraction layer to achieve portability and performance across these architectures. The core library is an implementation of a Voronoi cell-by-cell computation algorithm, producing the geometry of the cells and their neighbor connectivity information, rather than a full combinatorial mesh data structure. This simplifies the Voronoi tessellation and makes it more suitable to data parallel architectures than alternatives such as sequential insertion or the Bowyer-Watson algorithm. The library demonstrates significant performance improvements over established single-threaded programs and serves as a foundational tool for performance-critical applications, such as on-the-fly computations in hydrodynamical codes.

A chemical characterization of the Galactic Center is essential for understanding its formation and structural evolution. Trends of alpha-elements, such as Mg, Si, and Ca, serve as powerful diagnostic tools, offering insights into star-formation rates and gas-infall history. However, high extinction has previously hindered such studies. In this study, we present a detailed chemical abundance analysis of M giants in the Milky Way's Nuclear Star Cluster (NSC), focusing on alpha-element trends with metallicity. High-resolution, near-infrared spectra were obtained using the IGRINS spectrograph on the Gemini South telescope for nine M giants. Careful selection of spectral lines, based on a solar-neighborhood control sample of 50 M giants, was implemented to minimize systematic uncertainties. Our findings show enhanced alpha-element abundances in the predominantly metal-rich NSC stars, consistent with trends in the inner bulge. The NSC stars follow the high-[alpha/Fe] envelope seen in the solar vicinity's metal-rich population, indicating a high star-formation rate. The alpha-element trends decrease with increasing metallicity, also at the highest metallicities. Our results suggest the NSC population likely shares a similar evolutionary history with the inner bulge, challenging the idea of a recent dominant star formation burst. This connection between the NSC and the inner-disk sequence suggests that the chemical properties of extragalactic NSCs of Milky Way type galaxies could serve as a proxy for understanding the host galaxies' evolutionary processes.

In this study we present a detailed analysis of CN Lyn, an overlooked triple star system, by combining spectroscopic data from the literature, photometric \textit{TESS} data, and kinematic techniques. We updated the fundamental parameters of the known eclipsing components in the system with high precision. The chemical composition of both eclipsing components (Aab) and the third component (B) in the system were calculated with great accuracy. According to our analysis the mass, radius, and metallicity of the eclipsing components are $1.166_{-0.012}^{+0.013}\,M_\odot$, $1.786_{-0.014}^{+0.013}\,R_\odot$, and $-0.78_{-0.02}^{+0.02}$ dex for Aa and $1.143_{-0.012}^{+0.013}\,M_\odot$, $1.651_{-0.013}^{+0.014}\,R_\odot$, and $-0.55_{-0.02}^{+0.03}$ dex for Ab. The pair's age is $3.89_{-0.10}^{+0.10}$ Gyr. The mass, radius, metallicity, and age for B are $0.85_{-0.23}^{+0.23}\,M_\odot$, $1.436_{-0.023}^{+0.026}\,R_\odot$, $-1.83_{-0.11}^{+0.09}$ dex, and $12.5_{-2.5}^{+2.5}$ Gyr, respectively. It is also found that the triple system (AabB) satisfies the stability criteria for the hierarchical triple system. Kinematic and Galactic orbital parameters of CN Lyn were obtained from the astrometric and spectroscopic data of the system. Dynamical orbital analyses, taking into account the ages of the component stars in the central binary system (A) show that the CN Lyn originated at the metal-poor edge of the Galactic disk. The third component of the system was found to be a member of the halo population in terms of age, $\alpha$ elements and metal abundance. Given the different chemical abundances and age of B compared to A, this suggests that the third component was captured by the central system in a region with weak gravitational interactions far beyond the Galactic disc.

Using the Subaru/FOCAS IFU capability, we examine the spatially resolved relationships between gas-phase metallicity, stellar mass, and star-formation rate surface densities (Sigma_* and Sigma_SFR, respectively) in extremely metal-poor galaxies (EMPGs) in the local universe. Our analysis includes 24 EMPGs, comprising 9,177 spaxels, which span a unique parameter space of local metallicity (12+log(O/H) = 6.9 to 7.9) and stellar mass surface density (Sigma_* ~ 10^5 to 10^7 Msun/kpc^2), extending beyond the range of existing large integral-field spectroscopic surveys. Through spatially resolved emission line diagnostics based on the [NII] BPT-diagram, we verify the absence of evolved active galactic nuclei in these EMPGs. Our findings reveal that, while the resolved mass-metallicity relation exhibits significant scatter in the low-mass regime, this scatter is closely correlated with local star-formation surface density. Specifically, metallicity decreases as Sigma_SFR increases for a given Sigma_*. Notably, half of the EMPGs show a distinct metal-poor horizontal branch on the resolved mass-metallicity relation. This feature typically appears at the peak clump with the highest Sigma_* and Sigma_SFR and is surrounded by a relatively metal-enriched ambient region. These findings support a scenario in which metal-poor gas infall fuels episodic star formation in EMPGs, consistent with the kinematic properties observed in these systems. In addition, we identify four EMPGs with exceptionally low central metallicities (12+log(O/H) <~ 7.2), which display only a metal-poor clump without a surrounding metal-rich region. This suggests that such ultra-low metallicity EMPGs, at less than a few percent of the solar metallicity, may serve as valuable analogs for galaxies in the early stages of galaxy evolution.

Early JWST photometric studies discovered a population of UV faint ($\rm <L^{*}_{UV}$) $z \sim 6.5-8$ Lyman break galaxies with spectral energy distributions implying young ages ($\sim10$ Myr) yet relatively weak H$\beta$+[OIII] equivalent widths ($\rm EW_{H\beta+[OIII]} \approx 400$Å). These galaxies seemingly contradict the implicit understanding that young star-forming galaxies are ubiquitously strong H$\beta$+[OIII] emitters, i.e., extreme emission line galaxies (EW $\rm \gtrsim 750$Å). Low metallicities, high Lyman continuum escape fractions, and rapidly declining star-formation histories have been proposed as primary drivers behind low H$\beta$+[OIII] equivalent widths, but the blend of H$\beta$+[OIII] in photometric studies makes proving one of these scenarios difficult. We aim to characterize this peculiar population with deep spectroscopy from the JWST Advanced Deep Extragalactic Survey (JADES). We find that a significant subset of these galaxies at $z\gtrsim2$ with modest H$\beta$+[OIII] equivalent widths ($\rm \approx 300-600$Å) have high ionization efficiencies ($\rm \log \xi_{ion} \gtrsim 25.5~[Hz~erg^{-1}]$). Suppressed [OIII] EW values yet elevated H$\alpha$ and H$\beta$ EW values imply that the level of chemical enrichment is the primary culprit, supported by spectroscopic measurements of metallicities below 12+log(O/H)$\rm \approx 7.70~(10\%Z_{\odot})$. We demonstrate that integrated H$\beta$+[OIII] selections (e.g., H$\beta$+[OIII] EW $> 700$Å) exclude the most metal-poor efficient ionizers and favor 1) more chemically enriched systems with comparable extreme radiation fields and 2) older starbursting systems. In contrast, metallicity degeneracies are reduced in H$\alpha$ space, enabling the identification of these metal-poor efficient ionizers by their specific star-formation rate.

One-dimensional (1D) methods for simulating the common-envelope (CE) phase offer advantages over three-dimensional (3D) simulations regarding their computational speed and feasibility. We present the 1D CE method from Bronner et al. (2024), including the results of the CE simulations of an asymptotic giant branch star donor. We further test this method in the massive star regime by computing the CE event of a red supergiant with a neutron-star mass and a black-hole mass companion. The 1D model can reproduce the orbital evolution and the envelope ejection from 3D simulations when choosing suitable values for the free parameters in the model. The best-fitting values differ from the expectations based on the low mass simulations, indicating that the free parameters depend on the structure of the giant star. The released recombination energy from hydrogen and helium helps to expand the envelope, similar to the low-mass CE simulations.

According to modern cosmological models, galaxies are embedded within cosmic filaments, which supply a continuous flow of pristine gas, fueling star formation and driving their evolution. However, due to their low density, the direct detection of diffuse gas in cosmic filaments remains elusive. Here, we report the discovery of an extremely metal-poor ($[ X/H] \approx -3.7$), low-density ($\log_{10} n_{\rm H}/{\rm cm^{-3}} \approx -4$, corresponding to an overdensity of $\approx 5$) partial Lyman limit system (pLLS) at $z\approx3.577$ along the quasar sightline Q1317--0507, probing cosmic filaments. Additionally, two other low-metallicity ($[ X/H] \lesssim -2$) absorption systems are detected at similar redshifts, one of which is also a pLLS. VLT/MUSE observations reveal a significant overdensity of Ly-alpha emitters (LAEs) associated with the pLLS. The spatial distribution of the LAEs strongly suggests the presence of an underlying filamentary structure. This is further supported by the detection of a large Ly-alpha emitting nebula with a surface brightness of $\geq 10^{-19}~\rm erg~cm^{-2}~s^{-1}~arcsec^{-2}$, with a maximum projected linear size of $\approx 260$~pkpc extending along the LAEs.

The majority of most luminous quasars during the epoch of reionization accrete near or above the Eddington limit, marking the vigorous growth of primitive supermassive black holes (SMBHs). However, their subsequent evolution and environmental impact remain poorly characterized. We present JWST/NIRSpec prism IFU observations of HSC J2239+0207, a low-luminosity quasar at $z\sim6.25$ likely in a late stage of mass assembly with an overmassive SMBH relative to its host galaxy. Using H$\beta$ and H$\alpha$ broad emission lines, we estimate an SMBH mass $M_{\rm BH}\sim3\times10^8~M_{\odot}$ and confirm its sub-Eddington accretion at $\lambda_{\rm Edd}\sim0.4$. Strong FeII emission and a proximity zone of typical size suggest a metal-rich, highly evolved system. In the far-UV, this quasar presents strong broad-absorption-line features, indicative of high-velocity winds ($\nu\sim10^4~{\rm km/s}$). Meanwhile, minimal dust reddening is inferred from the quasar continuum and broad-line Balmer decrement, suggesting little dust along the polar direction. Most interestingly, we identify a gas companion $\sim$5 kpc from the quasar with a high [OIII]/H$\beta$ ratio ($\gtrsim10$), likely representing outflowing gas blown away by AGN feedback. These results highlight HSC J2239+0207 as a likely fading quasar in transition, providing rare insights into SMBH evolution, AGN feedback, and AGN-galaxy interactions in the early Universe.

We present Keck/KPIC phase II $K$-band observations of the non-transiting hot Jupiter HD 143105 b. Using a cross-correlation approach, we make the first detection of the planetary atmosphere at $K_p = 185^{+11}_{-13}\rm km\ s^{-1}$ and an inferior conjunction time 2.5 hours before the previously-published ephemeris. The retrieved $K_p$ value, in combination with orbital period, mass of the host star, and lack of transit detection, gives an orbital inclination of $78^{\circ+2}_{-12}$ and a true planet mass of 1.23$\pm0.10\rm\ M_J$. While the equilibrium temperature of HD 143105 b is in the transition regime between non-inverted and inverted atmospheres, our analysis strongly prefers a non-inverted atmosphere. Retrieval analysis indicates the atmosphere of HD 143105 b is cloud-free to approximately 1 bar and dominated by H$_2$O absorption ($\log \rm H_2O_{MMR} = -3.9^{+0.8}_{-0.5}$), placing only an upper limit on the CO abundance ($\log \rm CO_{MMR} < -3.7$ at 95% confidence). We place no constraints on the abundances of Fe, Mg, or $^{13}$CO. From these abundances, we place an upper limit on the carbon-to-oxygen ratio for HD 143105 b, $\rm C/O < 0.2$ at 95% confidence, and find the atmospheric metallicity is approximately $0.1\times$ solar. The low metallicity may be responsible for the lack of a thermal inversion, which at the temperature of HD 143105 b would likely require significant opacity from TiO and/or VO. With these results, HD 143105 b joins the small number of non-transiting hot Jupiters with detected atmospheres.

As star-forming galaxies approach or exceed a stellar mass around $10^{11} M_\odot$, they are increasingly likely to be quenched in a process generically called mass quenching. Central galaxies, which are quenched via mass rather than environmental quenching, therefore accumulate in a peak around this characteristic mass. While a number of processes may influence the shape of the quenched central stellar mass function (QCSMF), we find that its low-mass slope is strongly affected by the scatter in the mass of black holes at a given stellar mass, with higher scatters in the black hole population yielding shallower slopes. Higher scatters in the black hole mass spread out the stellar mass range over which quenching occurs, leading to shallower slopes. This trend holds across a variety of semi-analytic models and cosmological hydrodynamic simulations. A comparison with observations provides indirect evidence for a large scatter in black hole mass $\sigma(\log_{10}(M_\mathrm{BH})|M_*) \gtrsim 0.5$ dex, and a joint constraint on AGN feedback physics and the co-evolution of galaxies and black holes.

In the era of large-scale astronomical surveys, fast modeling of strong lens systems has become increasingly vital. While significant progress has been made for galaxy-scale lenses, the development of automated methods for modeling larger systems, such as groups and clusters, is not as extensive. Our study aims to extend the capabilities of the GIGA-Lens code, enhancing its efficiency in modeling multi-galaxy strong lens systems. We focus on demonstrating the potential of GPU-accelerated Bayesian inference in handling complex lensing scenarios with a high number of free parameters. We employ an improved inference approach that combines image position and pixelated data with an annealing sampling technique to obtain the posterior distribution of complex models. This method allows us to overcome the challenge of limited prior information, a high number of parameters, and memory usage. Our process is exemplified through the analysis of the compact group lens system DES J0248-3955, for which we present VLT/X-shooter spectra. We measure a redshift of $z = 0.69 \pm 0.04$ for the group, and $z = 1.2722 \pm 0.0005$ for one of the extended arcs. Our enhanced method successfully constrained a lens model with 29 free parameters and lax priors in a remarkably short time. The mass of the lens is well described by a single dark-matter halo with a velocity dispersion of $\sigma_v = (690 \pm 30) \, km \, s^{-1}$. The model predicts the presence of a second source at the same redshift and a third source at approximately $z \sim 2.7$. Our study demonstrates the effectiveness of our lens modeling technique for dealing with a complex system in a short time using ground-based data. This presents considerable potential within the context of large surveys such as LSST.

Numerical N-body simulations are commonly used to explore stability regions around exoplanets, offering insights into the possible existence of satellites and ring systems. This study aims to utilize Machine Learning (ML) techniques to generate predictive maps of stable regions surrounding a hypothetical planet. The approach can also be extended to planet-satellite systems, planetary ring systems, and other similar configurations. A dataset was generated using 10^5 numerical simulations, each incorporating nine orbital features for the planet and a test particle in a star-planet-test particle system. The simulations were classified as stable or unstable based on stability criteria, requiring particles to remain stable over a timespan equivalent to 10,000 orbital periods of the planet. Various ML algorithms were tested and fine-tuned through hyperparameter optimization to determine the most effective predictive model. Tree-based algorithms showed comparable accuracy in performance. The best-performing model, using the Extreme Gradient Boosting (XGBoost) algorithm, achieved an accuracy of 98.48%, with 94% recall and precision for stable particles and 99% for unstable particles. ML algorithms significantly reduce the computational time required for three-body simulations, operating approximately 100,000 times faster than traditional numerical methods. Predictive models can generate entire stability maps in less than a second, compared to the days required by numerical simulations. The results from the trained ML models will be made accessible through a public web interface, enabling broader scientific applications.

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present BNS_NURATES, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. BNS_NURATES targets an higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with M1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semi-leptonic charged-current processes, when $\mu^\pm$ are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of $\beta$-processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. The improved treatment for the reaction rates also modify the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

The modestly eccentric and non-coplanar orbits of the giant planets pose a challenge to solar system formation theories which generally indicate that the giant planets emerged from the protoplanetary disk in nearly perfectly circular and coplanar orbits. We demonstrate that a single encounter with a 2-50 Jupiter-mass object, passing through the solar system at a perihelion distance less than 20 AU and a hyperbolic excess velocity less than 6 km/s, can excite the giant planets' eccentricities and mutual inclinations to values comparable to those observed. We estimate that there is about a 1-in-100 chance that such a flyby produces a dynamical architecture similar to that of the solar system. We describe a metric to evaluate how closely a simulated system matches the eccentricity and inclination secular modes of the solar system. The scenario of a close encounter with a substellar object offers a plausible explanation for the origin of the moderate eccentricities and inclinations and the secular architecture of the planets.

We conducted a search for new ultracool companions to nearby white dwarfs using multiple methods, including the analysis of colors and examination of images in both the optical and the infrared. Through this process, we identified fifty-one previously unrecognized systems with candidate ultracool companions. Thirty-one of these systems are resolved in at least one catalog, and all but six are confirmed as co-moving companions via common proper motion and consistent parallax measurements (when available). We have followed up four co-moving companions with near-infrared spectroscopy and confirm their ultracool nature. The remaining twenty candidates are unresolved, but show clear signs of infrared excess which is most likely due to the presence of a cold, low-mass companion or a dusty circumstellar disk. Three of these unresolved systems have existing optical spectra that clearly show the presence of a cool stellar companion to the white dwarf primary via spectral decomposition. These new discoveries, along with our age estimates for the primary white dwarfs, will serve as valuable benchmark systems for future characterization of ultracool dwarfs.

As some of the most ancient materials in our Solar System, chondritic meteorites offer a valuable window into the early stages of planetary formation, particularly the accretion processes that built the most primitive asteroids. Until now, high energy shocks and collisions have been invoked to explain the deformation and fragmentation of chondrules, the main component of chondrites. However, simulating the cooling of chondrules using continuum mechanics and finite elements, we demonstrate that plastic deformation of chondrules can occur at low collision velocities of just a few meters per second and with kinetic energies less than tenths of a millijoule when temperatures exceed the glass transition temperature Tg ~ 1000 K. Conversely, below Tg, spontaneous chondrule cracking occurs due to differential thermal contraction between phases and is more pronounced in larger chondrules. Counterintuitively, our findings suggest that both ordinary and carbonaceous chondrites formed through similar low-energy processes, with varying degrees of ductility and brittleness depending on the amount of processed material. This implies that the environments where chondrites formed were likely less turbulent and more thermally active than previously thought.

The L9 dwarf SDSS J100711.74+193056.2 is situated 7$^{\circ}$.5 north of the nearest B-type star Regulus (d = 24.3+-0.2 pc), part of a stellar quadruplet. The object is at similar distance (d = 21.9+-1.0 pc) as Regulus, with a 3D separation of 3.9+0.6-0.5 pc ($\sim$1.6 tidal radii from Regulus), and shares tangential motion within 2 km/s, hinting at a physical connection. Near-infrared spectroscopy with Keck/NIRES finds that SDSS J100711.74+193056.2 also has a comparable radial velocity as Regulus A and B, a metallicity similar to Regulus B, and a spectral morphology consistent with the estimated 1--2 Gyr total age of Regulus's close pre-white dwarf companion. Taken together, these observations indicate that SDSS J100711.74+193056.2 is a very widely-separated and potentially physically-bound companion to Regulus, with a binding energy and mass ratio comparable to other wide star-brown dwarf systems. It joins a growing list of brown dwarfs at the L dwarf/T dwarf transition with independent constraints on physical properties such as age and metallicity.

We increase the spectroscopic completeness of the 100 pc white dwarf sample in the SDSS footprint with 840 additional spectra. Our spectroscopy is 86% complete for white dwarfs hotter than $T_{\rm eff}= 5000$ K, where H$\alpha$ remains visible and provides reliable constraints on the atmospheric composition. We identify 2108 DA white dwarfs with pure hydrogen atmospheres, and show that ultramassive DA white dwarfs with $M\geq1.1~M_{\odot}$ are an order of magnitude less common below 10,000 K. This is consistent with a fraction of them getting stuck on the crystallization sequence due to $^{22}$Ne distillation. In addition, there are no ultramassive DA white dwarfs with $M\geq1.1~M_{\odot}$ and $T_{\rm eff}\leq6000$ K in our sample, likely because Debye cooling makes them rapidly fade away. We detect a significant trend in the fraction of He-atmosphere white dwarfs as a function of temperature; the fraction increases from 9% at 20,000 K to 32% at 6000 K. This provides direct evidence of convective mixing in cool DA white dwarfs. Finally, we detect a relatively tight sequence of low-mass DQ white dwarfs in color-magnitude diagrams for the first time. We discuss the implications of this tight DQ sequence, and conclude with a discussion of the future prospects from the upcoming ULTRASAT mission and the large-scale multi-fiber spectroscopic surveys.

The existence of one percent of lithium-rich giant stars among normal, lithium-poor giant stars continues to be poorly explained. By merging two catalogues, one containing 10,535 lithium-rich giant stars with lithium abundances ranging from 1.5 to 4.9 dex, and the other detecting infrared sources, we have found 421 clump giant stars and 196 first-ascending giant stars with infrared excesses indicating stellar mass losses. The clump stars are the most lithium-rich. Approximately 5.8 percent of these stars episodically lose mass in periods of approximately 10^4 years or less, while the remaining stars ceased their mass loss and maintained their lithium for nearly 10^7 years. We propose a scenario in which all giant stars with masses below two solar masses undergo prompt lithium enrichment with mass-ejection episodes. We suggest that mass loss results from internal angular-momentum transport. It is possible that a transitory instability, perhaps of magnetic origin, rapidly transports the nuclear material responsible for the lithium enrichment to the stellar surface and triggers shell ejections. Additionally, the strong mass loss in some lithium-rich stars during their evolution activates their chromospheres, as observed in ultraviolet spectra. Furthermore, intense episodical mass losses in these stages led to the observable formation of complex organic and inorganic particles, as detected in near-infrared spectra. In contrast to first-ascending giant stars, helium flashes during the clump can contribute to additional lithium enrichment alongside the aforementioned process. The combination of these two lithium sources may explain the much higher observed lithium abundances in clump stars, as well as their observed infrared excesses. If our scenario based on a universal and rapid lithium enrichment episode process is correct, it could explain the rarity of lithium-rich giant stars.

In this third paper in a series, we investigate the need of spectra denoising for the derivation of stellar parameters. We have used two distinct datasets for this work. The first one contains spectra in the range of 4450-5400 Å at a resolution of 42000 and the second in the range of 8400-8800 Å at a resolution of 11500. We constructed two denoising techniques, an autoencoder, and a Principal Component Analysis. Using random Gaussian noise added to synthetic spectra, we have trained a Neural Network to derive the stellar parameters Teff, log g, ve sin i, {\xi}t, and [M/H] of the denoised spectra. We find that, independently of the denoising technique, the stellar parameters accuracy values do not improve once we denoise the synthetic spectra. This is true with and without applying data augmentation to the stellar parameters Neural Network.

The Hubble constant $H_0$, the current expansion rate of the universe, is one of the most important parameters in cosmology. The cosmic expansion regulates the mutually approaching motion of a pair of celestial objects due to their gravity. Therefore, the mean pairwise peculiar velocity of celestial objects, which quantifies their relative motion, is sensitive to both $H_0$ and the dimensionless total matter density $\Omega_m$. Based on this, using the Cosmicflows-4 data, we measured $H_0$ for the first time via the galaxy pairwise velocity in the nonlinear and quasi-linear range. Our results yield $H_0=75.5\pm1.4$ km s$^{-1}$ Mpc$^{-1}$ and $\Omega_m=0.311^{+0.029}_{-0.028}$ . The uncertainties of $H_0$ and $\Omega_m$ can be improved to around 0.6% and 2%, respectively, if the statistical errors become negligible in the future.

About one third of Fermi Large Area Telescope (LAT) sources are unassociated. We perform multi-class classification of Fermi-LAT sources using machine learning with the goal of probabilistic classification of the unassociated sources. A particular attention is paid to the fact that the distributions of associated and unassociated sources are different as functions of source parameters. In this work, we address this problem in the framework of dataset shifts in machine learning.

Volatiles like $H_2O$ are present as ice in solids in the outer cold regions of protoplanetary disks and as vapor in the warm inner regions within the water snow line. Icy pebbles drifting inwards from the outer disk sublimate after crossing the snow line, enriching the inner disk with solid mass and water vapor. Meanwhile, proto-planets forming within the disk open gaps in the disk gas, creating traps against the inward drift of pebbles and in turn reducing water enrichment in the inner disk. Recent disk observations from millimeter interferometry and infrared spectroscopy have supported this broad picture by finding a correlation between the outer radial distribution of pebbles and the properties of inner water vapor spectra. In this work, we aim at further informing previous and future observations by building on previous models to explore pebble drift in disks with multiple gaps. We systematically explore multiple gap locations and their depths (equivalent to specific masses of planets forming within), and different particle sizes to study their impact on inner disk water enrichment. We find that the presence of close-in deep gaps carved by a Jupiter-mass planet is likely crucial for blocking icy pebble delivery into the inner disk, while planets with lower masses only provide leaky traps. We also find that disks with multiple gaps show lower vapor enrichment in the inner disk. Altogether, these model results support the idea that inner disk water delivery and planet formation are regulated by the mass and location of the most massive planets.

The Compton Spectrometer and Imager balloon payload (COSI-Balloon) is a wide-field-of-view Compton ${\gamma}$-ray telescope that operates in the 0.2 - 5 MeV bandpass. COSI-Balloon had a successful 46-day flight in 2016 during which the instrument observed the Crab Nebula, Cygnus X-1, and Centaurus A. Using the data collected by the COSI-Balloon instrument during this flight, we present the source flux extraction of signals from the variable balloon background environment and produce images of these background-dominated sources by performing Richardson-Lucy deconvolutions. We also present the spectra measured by the COSI-Balloon instrument, compare and combine them with measurements from other instruments, and fit the data. The Crab Nebula was observed by COSI-Balloon and we obtain a measured flux in the energy band 325 - 480 keV of (4.5 ${\pm}$ 1.6) ${\times}$ 10$^{-3}$ ph cm$^{-2}$ s$^{-1}$. The model that best fits the COSI-Balloon data combined with measurements from NuSTAR and Swift-BAT is a broken power law with a measured photon index ${\Gamma}$ = 2.20 ${\pm}$ 0.02 above the 43 keV break. Cygnus X-1 was also observed during this flight, and we obtain a measured flux of (1.4 ${\pm}$ 0.2) ${\times}$ 10$^{-3}$ ph cm$^{-2}$ s$^{-1}$ in the same energy band and a best-fit result (including data from NuSTAR, Swift-BAT, and INTEGRAL/ IBIS) was to a cutoff power law with a high-energy cutoff energy of 138.3 ${\pm}$ 1.0 keV and a photon index of ${\Gamma}$ = 1.358 ${\pm}$ 0.002. Lastly, we present the measured spectrum of Centaurus A and our best model fit to a power law with a photon index of ${\Gamma}$ = 1.73 ${\pm}$ 0.01.

Galactic bars can form via the internal bar instability or external tidal perturbations by other galaxies. We systematically compare the properties of bars formed through the two mechanisms with a series of controlled $N$-body simulations that form bars through internal or external mechanisms. We create three disk galaxy models with different dynamical ``hotness'' and evolve them in isolation and under flyby interactions. In the cold and warm disk models, where bars can form spontaneously in isolation, tidally-induced bars are promoted to a more ``advanced'' evolutionary stage. However, these bars have similar pattern speeds to those formed spontaneously within the same disk. Bars formed from both mechanisms have similar distributions in pattern speed--bar strength ($\Omega_p-A_2$) space and exhibit comparable ratios of co-rotation radius to bar length (${\cal R}={R_{\mathrm {CR}}}/{R_{\mathrm {bar}}}$). Dynamical analyses suggest that the inner stellar disk loses the same amount of angular momentum, irrespective of the presence or intensity of the perturbation, which possibly explains the resemblance between tidally and spontaneously formed bars. In the hot disk model, which avoids the internal bar instability in isolation, a bar forms only under perturbations and rotates more slowly than those in the cold and warm disks. Thus, if ``tidally-induced bars'' refer exclusively to those in galaxies that are otherwise stable against bar instability, they indeed rotate slower than internally-induced ones. However, the pattern speed difference is due to the difference in the internal properties of the bar host galaxies, not the different formation mechanisms.

The circumgalactic medium (CGM) is a key component needed to understand the physical processes governing the flows of gas around galaxies. Quantifying its evolution and its dependence on galaxy properties is particularly important for our understanding of accretion and feedback mechanisms. We select a volume-selected sample of 66 {\it isolated} star-forming galaxies (SFGs) at $0.4< z <1.5$ with $\log(M_\star/M_{\odot})> 9$ from the MusE GAs FLOw and Wind (MEGAFLOW) survey. Using MgII 2796,2803 absorptions in background quasars, we measure the covering fraction $f_c$ and quantify how the cool gas profile depends on galaxy properties (such as star-formation rate (SFR), stellar mass ($M_\star$) or azimuthal angle relative to the line of sight) and how these dependencies evolve with redshift. The MgII covering fraction of isolated galaxies is a strong function of impact parameter, and is steeper than previously reported. The impact parameter $b_{50}$ at which $f_c = $50\% is $b_{50}=50\pm7$kpc ($65\pm7$ kpc) for $W_r^{2796}>$0.5 Å($W_r^{2796}>0.1$ Å), respectively. It is weakly correlated with SFR ($\propto$ SFR$^{0.08\pm0.09}$) and decreases with cosmic time ($\propto (1+z)^{0.8 \pm 0.7}$), contrary to the expectation of increasingly larger halos with time. The covering fraction is also higher along the minor axis than along the major axis at the $\approx 2 \sigma$ level. The CGM traced by \MgII{} is similar across the isolated galaxy population. Indeed, among the isolated galaxies with an impact parameter below 55 kpc, all have associated MgII absorption with $W_r^{2796}>$0.3Å, resulting in a steep covering fraction $f_c(b)$.

In this work, we perform a detailed analysis to constrain the Hu-Sawicki $f(R)$ gravity model, using cosmic shear data from three prominent Stage-III weak lensing surveys: DES-Y3, KiDS-1000, and HSC-Y3. To accurately model the nonlinear matter clustering in the analysis of cosmic shear signals, we employ FREmu, a recently developed power spectrum emulator for the $f(R)$ gravity trained on the Quijote-MG simulations. This emulator achieves precise predictions, limiting the errors to 5% on scales of $0.009h\,{\rm Mpc}^{-1} < k < 0.5h\,{\rm Mpc}^{-1}$. Our findings reveal that cosmic shear data alone impose only weak constraints on the $f(R)$ parameter $\log_{10}|f_{R_0}|$. To improve these constraints, we incorporate state-of-the-art external observations, including data from the cosmic microwave background and baryon acoustic oscillations. The inclusion of these external datasets significantly enhances the constraints, yielding an upper limit of $\log_{10}|f_{R_0}| < -4.79$ at the 95% confidence level.

We investigate accretion onto an isolated black hole from uniform winds. If the winds are directed towards the black hole, then the accretion process can be well described by the classical Bondi-Hoyle Lyttleton or BHL accretion. If the wind is not directed towards the black hole and flows past it, then a smaller fraction of the flow can be attracted by the black hole, and this type of accretion cannot be described by the classical BHL, and we coin the second kind as the lateral BHL. We show that the classical BHL cannot form an accretion disk, while lateral BHL can form transient accretion disks. To describe the thermodynamics of the flow, we have used a variable adiabatic index equation of state which depends on the temperature of the flow as well as the composition of the gas. We show that the electron-proton gas forms an accretion disk, which disappears and forms a shock cone, only to form the disk again at a later time, while for flows with less protons, the accretion disk, once lost, does not reappear again. Only when the flow is pair-dominated does it form a persistent accretion disk. We also show that a shock cone is less luminous than the accretion disk.

Recently, the measurements of baryon acoustic oscillations (BAO) by the Dark Energy Spectroscopic Instrument (DESI) indicate a potential deviation from the standard $\Lambda$CDM model. Some studies suggest that the data points from the luminous red galaxies (LRG) survey in DESI BAO data may contribute to this discrepancy. In this work, our main goal is to investigate whether this deviation is caused by the parameterization of the equation of state (EoS) of dark energy (DE). Hence, we have examined four popular parameterized dark energy models in our analysis: the Chevallier-Polarski-Linder (CPL), Barboza-Alcaniz (BA), Jassal-Bagla-Padmanabhan (JBP), and Feng-Shen-Li-Li (FSLL) parameterizations. Considering that LRG1 and LRG2 data points may lead to deviation from the $\Lambda$CDM model, we use two versions of DESI BAO data, differing in whether these data points are included. Additionally, to break the parameter degeneracies and obtain robust constraint results, we introduce Type Ia supernovae (SNe Ia) and quasars (QSO) in our analysis. Our findings indicate that in these parameterizations, the deviation from ($w_0$,$w_{1}$)=(-1,0) becomes more pronounced when using the combined data from DESI BAO, SNe Ia, and QSO compilations. Here, $w_{0}$ and $w_{1}$ represent the EoS of DE. It suggests that the parameterizations of the EoS of DE have little impact on the deviation from the $\Lambda$CDM model. Besides, our analysis potentially hints that dark energy may have dynamic properties. In addition, the results obtained from different BAO datasets demonstrate that the LRG1 and LRG2 data points do indeed contribute to a deviation from the $\Lambda$CDM model. Finally, according to the statistical criteria, the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC), the joint constraints provide substantial observational support to the BA and FSLL models.

We report the discovery of a large ($\sim 10$ kpc diameter), massive ($\log(M_\star/M_\odot) = 10.15^{+0.01}_{-0.01}$), grand-design spiral galaxy with photometric redshift $z_{\text{phot}} = 4.03$ in the UNCOVER and Medium Band Mega Science surveys with JWST. This is the highest redshift spiral galaxy discovered with JWST so far. In the rest-frame near-UV and far-UV, we clearly see the beads-on-a-string pattern of star formation; in the rest-frame visible bands, each string appears as an arm. Spectral energy distribution modeling using the Bagpipes code is strongly constrained by detections and flux measurements in 21 JWST and HST filters. The stellar mass-weighted age is 228 Myr, implying that 50% of the stars in the galaxy formed after $z \sim 4.5$. This is a highly star-forming galaxy with a star formation rate (SFR) of $57.57^{+1.80}_{-1.90} \, M_\odot\, \text{yr}^{-1}$. We detect strong H-$\alpha$ + [NII] emission from the entire disk. The detection of a spiral galaxy at $z \sim 4$ indicates that massive and large spiral galaxies and disks were already in place merely 1.5 billion years after the Big Bang.

HD 60435 is a well-known rapidly oscillating (roAp) Ap star with a series of alternating even and odd degree modes, making it a prime asteroseismic target. It is also an oblique pulsator with rotational inclination, $i$, and magnetic/pulsation obliquity, $\beta$, such that both magnetic/pulsation poles are viewed over the rotation period, $P_{\rm rot} = 7.679696$ d, determined from rotational light variations. While some roAp stars have stable pulsation mode amplitudes over decades, HD 60435 is known to have amplitude variations on time scales as short as 1 d. We show from 5 yr of {\it TESS} observations that there is strong amplitude modulation on this short time scale with possible mode interactions. Most remarkably, HD 60435 stopped pulsating during the time span of the {\it TESS} observations. This is the first time that any pulsating star has been observed to cease pulsating entirely. That has implications for mode interaction, excitation and damping, and is relevant to the problem of why only some stars in many pulsation instability strips pulsate, while others do not. During a 24.45-d time span of the {\it TESS} data when there was mode stability for a dipole mode and a quadrupole mode, the oblique pulsator model constrained $i$ and $\beta$, which we used to model those modes with a magnetic pulsation model from which we determined a polar field strength of 4 kG, in good agreement with a known magnetic measurement. We modelled the frequency separations showing that they can constrain the global metallicity, something that is not possible from spectroscopy of the highly peculiar Ap atmosphere.

In recent years, the prospect of detecting gravitational waves sourced from a strongly first-order cosmological phase transition has emerged as one of the most exciting frontiers of gravitational wave astronomy. Cosmological phase transitions are an essential ingredient in the Standard Model of particle cosmology, and help explain the mechanism for creation of matter in the early Universe, provide insights into fundamental theories of physics, and shed light on the nature of dark matter. This underscores the significance of developing robust end-to-end tools for determining the resulting gravitational waves from these phase transitions. In this article we present PhaseTracer2, an improved version of the C++ software package PhaseTracer, designed for mapping cosmological phases and transitions in Standard Model extensions of multiple scalar fields. Building on the robust framework of its predecessor, PhaseTracer2 extends its capabilities by including new features crucial for a more comprehensive analysis of cosmological phase transitions. It can calculate more complex properties, such as the bounce action through the path deformation method or an interface with BubbleProfiler, thermodynamic parameters, and gravitational wave spectra. Its applicability has also been broadened via incorporating the dimensionally reduced effective potential for models obtained from DRalgo, as well as calculations in the MSbar and OS-like renormalisation schemes. This modular, flexible, and practical upgrade retains the speed and stability of the original PhaseTracer, while significantly expanding its utility.

The stability of weakly collisional plasmas is well represented by linear theory, and the generated waves play an essential role in the thermodynamics of these systems. The velocity distribution functions (VDF) characterizing kinetic particle behavior are commonly represented as a sum of anisotropic bi-Maxwellians. For the majority of in situ observations of solar wind plasmas enabled by heliospheric missions, a three bi-Maxwellian model is commonly applied for the ions, assuming that the VDF consists of a proton core, proton beam, and a single He ($\alpha$) particle population, each with their own density, bulk velocity, and anisotropic temperature. Resolving an $\alpha$-beam component was generally not possible due to instrumental limitations. The Solar Orbiter Solar Wind Analyser Proton and Alpha Sensor (SWA PAS) resolves velocity space with sufficient coverage and accuracy to routinely characterize secondary $\alpha$ populations consistently. This design makes the SWA PAS dataset ideal for examining effects of the $\alpha$-particle beam on the plasma's kinetic stability. We test the wave signatures observed in the magnetic field power spectrum at ion scales and compare them to the predictions from linear plasma theory, Doppler-shifted into the spacecraft reference frame. We find that taking into account the $\alpha$-particle beam component is necessary to predict the coherent wave signatures in the observed power spectra, emphasizing the importance of separating the $\alpha$-particle populations as is traditionally done for protons. Moreover, we demonstrate that the drifts of beam components are responsible for the majority of the modes that propagate in oblique direction to the magnetic field, while their temperature anisotropies are the primary source of parallel Fast Magnetosonic Modes in the solar wind.

From 2012 to 2023, the PRISMA-32 array was in operation at the Experimental Complex NEVOD (MEPhI, Moscow). The purpose of the array was to study extensive air showers by detecting the air-shower neutron and electron-photon components using unshielded neutron detectors. To expand the capabilities of this facility, including for the study of cosmic and geophysical phenomena with a neutron flux, its upgrade was carried out. During the upgrade, a dedicated measuring channel for studying variations of the neutron background and the processes affecting these variations was created. To achieve this, the photomultipliers, the integrating amplifiers, the digitalizing electronics and the high-voltage power supply system were replaced. The paper describes the structure of the upgraded array, which was named PRISMA-36, and presents the results of studying the characteristics of the main elements of its "variation" channel. A method for identifying signals caused by neutron capture and the determined criteria for their selection are discussed. An example of a Forbush decrease, caused by a X1.1-class flare and recorded with the variation channel of the PRISMA-36 array, is given.

The inference of stellar parameters (such as radius and mass) through asteroseismic forward modelling depends on the number, accuracy, and precision of seismic and atmospheric constraints. ESA's Gaia space mission is providing precise parallaxes which yield an additional constraint to be included in the model grid search. Using a handful of main-sequence benchmark stars, we perform a uniform characterisation of these stars. We assess the accuracy and precision of stellar parameters inferred from grid-based searches when a Gaia-based luminosity is combined with different stellar constraints. We also examine the precision needed for an interferometric radius (model-independent radius) to have a significant contribution towards the determination of stellar mass in the optimisation process. Our findings show that more precise stellar masses are inferred for some stars when seismic and spectroscopic constraints are complemented with a Gaia-based luminosity, with a scatter varying from 1.9 per cent to 0.8 per cent. However, the inferred stellar radii are underestimated when compared to the interferometric radii and yield a scatter of $\sim$1.9 per cent. In addition, we demonstrate that a precisely measured interferometric radius ($\lesssim$ 1 per cent) when applied in the optimisation process yields a mass with a precision $\lesssim$ 1.5 per cent. Finally, we find that when only $l=0$ mode oscillation frequencies are available, robust masses and radii are still attainable. However, this requires precise and numerous $l=0$ mode oscillations frequencies ($>$ 8) to be coupled with atmospheric constraints.

In the incoming years, cosmological surveys aim at measuring the sum of neutrino masses $\Sigma m_\nu$, complementing the determination of their mass ordering from laboratory experiments. In order to assess the full potential of large-scale structures (LSS), we employ state-of-the-art predictions from the effective field theory of LSS (EFTofLSS) at one loop to perform Fisher forecasts on the sensitivity (combining power spectrum and bispectrum) of ongoing and future surveys (DESI, MegaMapper) in combination with CMB measurements (Planck, Litebird and Stage-4). We find that the 1$\sigma$ sensitivity on $\Sigma m_\nu$ is expected to be 15 meV with Planck+DESI, and 7 meV with S4+MegaMapper, where $\sim 10\%$ and $30\%$ of the constraints are brought by the one-loop bispectrum respectively. To understand how robust are these bounds, we explore how they are relaxed when considering extensions to the standard model, dubbed `new physics'. We find that the shift induced on $\Sigma m_\nu$ by a $1\sigma$ shift on new physics parameters (we consider extra relativistic species, neutrino self-interactions, curvature or a time-evolving electron mass) could be $\mathcal O(10)$ meV for Planck+DESI, but it will be suppressed down to $\mathcal O(1)$ meV in S4+MegaMapper. Our study highlights the quantitative impact of including the bispectrum at one loop in the EFTofLSS, and the robustness of the sensitivity to $\Sigma m_\nu$ against potential new physics thanks to the synergy of cosmological probes.

The James Webb Space Telescope (JWST) has recently discovered a new population of objects at high redshift referred to as `Little Red Dots' (LRDs). Their nature currently remains elusive, despite their surprisingly high inferred number densities. This emerging population of red point-like sources is reshaping our view of the early Universe and may shed light on the formation of high-redshift supermassive black holes. Here we present a spectroscopically confirmed LRD CANUCS-LRD-z8.6 at $z_{\rm spec}=8.6319\pm 0.0005$ hosting an Active Galactic Nucleus (AGN), using JWST data. This source shows the typical spectral shape of an LRD (blue UV and red optical continuum, unresolved in JWST imaging), along with broad H$\beta$ line emission, detection of high-ionization emission lines (CIV, NIV]) and very high electron temperature indicative of the presence of AGN. This is also combined with a very low metallicity ($Z<0.1 Z_\odot$). The presence of all these diverse features in one source makes CANUCS-LRD-z8.6 unique. We show that the inferred black hole mass of CANUCS-LRD-z8.6 ($M_{\rm BH}=1.0^{+0.6}_{-0.4}\times 10^{8}\rm ~M_\odot$) strongly challenges current standard theoretical models and simulations of black hole formation, and forces us to adopt `ad hoc' prescriptions. Indeed if massive seeds, or light seeds with super-Eddington accretion, are considered, the observed BH mass of CANUCS-LRD-z8.6 at $z=8.6$ can be reproduced. Moreover, the black hole is over-massive compared to its host, relative to the local $M_{\rm BH}-M_*$ relations, pointing towards an earlier and faster evolution of the black hole compared to its host galaxy.

The motion of charged particles under the Lorentz force in the magnetosphere of neutron stars, represented by a dipole field in the Schwarzschild spacetime, can be determined by an effective potential, whose local extrema govern circular orbits both in and off the equatorial plane, which coincides with the symmetry plane of the dipole field. In this work, we provide a detailed description of the properties of these "conservative" circular orbits and, using the approximation represented by the Landau-Lifshitz equation, examine the role of the radiative back-reaction force that influences the motion of charged particles following both the in and off equatorial circular orbits, as well as the chaotic orbits confined to belts centered around the circular orbits. To provide clear insight into these dynamics, we compare particle motion with and without the back-reaction force. We demonstrate that, in the case of an attractive Lorentz force, the back-reaction leads to the charged particles falling onto the neutron star's surface in all scenarios considered. For the repulsive Lorentz force, in combination with the back-reaction force, we observe a widening of stable equatorial circular orbits; the off-equatorial orbits shift toward the equatorial plane and subsequently widen if they are sufficiently close to the plane. Otherwise, the off-equatorial orbits evolve toward the neutron star surface. The critical latitude, which separates orbital widening from falling onto the surface, is determined numerically as a function of the electromagnetic interaction's intensity.

In Lammer et al. 2024, we defined Earth-like Habitats (EH) as rocky planets in the habitable zone of complex life (HZCL) on which Earth-like N$_2$-O$_2$-dominated atmospheres with minor amounts of CO$_2$ can exist and derived a formula for estimating their maximum number in the Galaxy. Here, we apply this formula by considering only requirements that are already scientifically quantifiable. By implementing models for star formation rate, initial mass function, and galactic mass distribution, we calculate the spatial distribution of disk stars as functions of stellar mass and birth age. We apply models for the GHZ and evaluate the thermal stability of Earth-like atmospheres with various CO$_2$ mixing ratios by implementing the newest stellar evolution and upper atmosphere models. In addition, we include the rocky exoplanet frequency, the availability of oceans and subaerial land, and the potential large moon requirement by evaluating their importance and implementing these criteria from minima to maxima values. We also discuss factors that are not yet scientifically quantifiable but may be requirements for EHs to evolve. We find that EHs are rare by obtaining maximum numbers of $2.5^{+71.6}_{-2.4}\times10^{5}$ and $0.6^{+27.1}_{-0.59}\times10^{5}$ planets that can potentially host N$_2$-Earth-like atmospheres with maximum CO$_2$ mixing ratios of 10\% and 1\%, respectively, implying that a minimum of $\sim 10^3 - 10^6$ rocky HZCL planets are needed for 1 EH to evolve. Their actual number, however, may be substantially lower as several requirements are not included in our model; this also implies ETIs are significantly rarer still. Our results illustrate that neither every star can host EHs, nor that each rocky HZCL planet evolves such that it may be able to host complex animal-like life. The Copernican Principle therefore cannot be applied to infer that such life is common in the Galaxy.

In this hypothesis article, we discuss the basic requirements of planetary environments where aerobe organisms can grow and survive, including atmospheric limitations of millimeter-to-meter-sized biological animal life based on physical limits, and O$_2$, N$_2$, and CO$_2$ toxicity levels. By assuming that animal-like extraterrestrial organisms adhere to similar limits, we define Earth-like Habitats ($\eta_{\rm EH}$) as rocky exoplanets in the Habitable Zone of Complex Life that host N$_2$-O$_2$-dominated atmospheres with minor amounts of CO$_2$, at which advanced animal-like life can in principle evolve and exist. We then derive a new formula that can be used to estimate the maximum occurrence rate of such Earth-like Habitats in the Galaxy. This contains realistic probabilistic arguments that can be fine-tuned and constrained by atmospheric characterization with future space and ground-based telescopes. As an example, we briefly discuss two specific requirements feeding into our new formula that, although not quantifiable at present, will become scientifically quantifiable in the upcoming decades due to future observations of exoplanets and their atmospheres.

Star clusters are interesting laboratories to study star formation, single and binary stellar evolution, and stellar dynamics. We have used the exquisite data from $Gaia$'s data release 3 (DR3) to study 21 relatively rich and nearby open clusters with member numbers ($N_{\rm{cl}}$)$>500$. We have developed a non-parametric method to identify cluster members. Our method works well for clusters located in both sparse and crowded environments, hence, can be applied to a wide variety of star clusters. Since the member classification scheme does not make any assumptions on the expected distributions of potential cluster members, our method can identify members associated with clusters that are oddly shaped or have complex internal spatial or kinematic structures. In addition, since the membership determination does not depend on the proximity to any well-defined sequences on the color-magnitude diagram, this method easily identifies straggler members. Furthermore, for each of these clusters, we estimate essential cluster properties including age, metallicity, distance, and reddening using detailed Markov-Chain Monte Carlo parameter estimation. We report the full posteriors for these important cluster properties for all clusters in our study.

We analyzed the 7.92$\times 10^{11}$ cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from May 13, 2011, when the fully constructed experiment started to take data, to May 12, 2023. This dataset provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100-300 TeV energy range.

We present a machine learning method to assign stellar parameters (temperature, surface gravity, metallicity) to the photometric data of large photometric surveys such as SDSS and SKYMAPPER. The method makes use of our previous effort in homogenizing and recalibrating spectroscopic data from surveys like APOGEE, GALAH, or LAMOST into a single catalog, which is used to inform a neural network. We obtain spectroscopic-quality parameters for millions of stars that have only been observed photometrically. The typical uncertainties are of the order of 100K in temperature, 0.1 dex in surface gravity, and 0.1 dex in metallicity and the method performs well down to low metallicity, were obtaining reliable results is known to be difficult.

In this contribution we present the FAST, which is a comprehensive software suite that aims to streamline and automatically manage the forecast of atmospheric and astroclimatic parameters (provided respectively by Meso-Nh and Astro-Meso-Nh models) on large ground-based telescope installations. The forecast of the aforementioned parameters is becoming crucial for the operation of the large telescope installations which possess atmospheric-sensitive equipment equipped with Adaptive Optics (AO) systems. FAST performs automatically all the steps of an atmosphere forecast process: initialisation and forcing data, atmospheric simulation, postprocessing and managing of the this http URL role of such service is useful both in optimizing beforehand AO instruments to the next atmospheric conditions and in planning telescope observations (especially in "service mode") in order to maximize the scientific output. FAST was applied first to the ALTA Center project, which provides forecasts for the LBT telescope. Then it was extended to the more recent project FATE that is a similar forecast system applied to the VLT. Since its first version FAST evolved and it has has been modified to fit with the different technical specifications of the different projects gaining in modularity. It is now able to provide forecasts on different timescales (from days to hours before) and to provide forecast during night and day time. After several years of continuous development we can say that FAST reached full maturity and it is now ready for applications to other projects/sites.

Observed IR excesses indicate that protoplanetary discs evolve slowly for the majority of their lifetime before losing their near- and mid-IR excesses on short timescales. Photoevaporation models can explain this "two-timescale" nature of disc evolution through the removal of inner regions of discs after a few million years. However, they also predict the existence of a population of non-accreting discs with large cavities. Such discs are scarce within the observed population, suggesting the models are incomplete. We explore whether radiation-pressure-driven outflows are able to remove enough dust to fit observations. We simulate these outflows using cuDisc, including dust dynamics, growth/fragmentation, radiative transfer and a parameterisation of internal photoevaporation. We find that, in most cases, dust mass-loss rates are around 5-10 times too small to meet observational constraints. Particles are launched from the disc inner rim, however grains larger than around a micron do not escape in the outflow, meaning mass-loss rates are too low for the initial dust masses at gap-opening. Only systems that have smooth photoevaporation profiles with gas mass-loss rates $>\sim 5 \times 10^{-9}$ $M_\odot$ yr$^{-1}$ and disc dust masses $<\sim$1 $M_\oplus$ at the time of gap opening can meet observational constraints; in the current models these manifest as EUV winds driven by atypically large high-energy photon fluxes. We also find that the height of the disc's photosphere is controlled by small grains in the outflow as opposed to shadowing from a hot inner rim; the effect of this can be seen in synthetic scattered light observations.

Orphan gamma-ray burst afterglows are good candidates to learn more about the GRB physics and progenitors or for the development of multi-messenger analysis with gravitational waves. Our objective is to identify orphan afterglows in Rubin LSST data, by using the characteristic features of their light curves. In this work, we generated a population of short GRBs based on the Swift SBAT4 catalogue, and we simulated their off-axis afterglow light curves with afterglowpy. We then used the rubin_sim package to simulate observations of these orphan afterglows with Rubin LSST and proceeded with the characterisation of orphan light curves by extracting a number of parameters. The same parameters are computed for the ELAsTiCC (Extended LSST Astronomical Time-series Classification Challenge) data set, a simulated alert stream of the Rubin LSST data. We then started to develop a machine learning filter able to discriminate orphan-like events among all the variable objects. We present here the performance of our filter as implemented in the Fink broker and tested on the ELAsTiCC data set and our own Rubin pseudo-observation simulations.

The magnetic buoyancy (MBI) and Parker instabilities are strong and generic instabilities expected to occur in most astrophysical systems with sufficiently strong magnetic fields. In galactic and accretion discs, large-scale magnetic fields are thought to result from the mean-field dynamo action, in particular, the $\alpha^2\Omega$. Using non-ideal MHD equations, we model a section of the galactic disc in which the large-scale magnetic field is generated by an imposed $\alpha$-effect and differential rotation. We extend our earlier study of the interplay between magnetic buoyancy and the mean-field dynamo. We add differential rotation which enhances the dynamo and cosmic rays which enhance magnetic buoyancy. We construct a simple 1D model which replicates all significant features of the 3D simulations. We confirm that magnetic buoyancy can lead to oscillatory magnetic fields and discover that it can vary the magnetic field parity between quadrupolar and dipolar, and that inclusion of the differential rotation is responsible for the switch in field parity. Our results suggest that the large-scale magnetic field can have a dipolar parity within a few kiloparsecs of the galactic centre, provided the MBI is significantly stronger the the dynamo. Quadrupolar parity can remain predominant in the outer parts of a galactic disc. Cosmic rays accelerate both the dynamo and the MBI and support oscillatory non-linear states, a spatial magnetic field structure similar to the alternating magnetic field directions observed in some edge-on galaxies.

Although Cherenkov detectors of high-energy neutrinos in ice and water are often optimized to detect TeV-PeV neutrinos, they may also be sensitive to transient neutrino sources in the 1-100~GeV energy range. A wide variety of transient sources have been predicted to emit GeV neutrinos. In light of the upcoming IceCube-Upgrade, which will extend the IceCube detector's sensitivity down to a few GeV, as well as improve its angular resolution, we survey a variety of transient source models and compare their predicted neutrino fluences to detector sensitivities, in particular those of IceCube-DeepCore and the IceCube Upgrade. We consider the ranges of neutrino fluence from transients powered by non-relativistic shocks, such as novae, supernovae, fast blue optical transients, and tidal disruption events. We also consider fast radio bursts and relativistic outflows of high- and low-luminosity gamma-ray bursts. Our study sheds light on the prospects of observing GeV transients with existing and upcoming neutrino facilities.

Asteroseismic inferences of main-sequence solar-like oscillators often rely on best-fit models. However, these models cannot fully reproduce the observed mode frequencies, suggesting that the internal structure of the model does not fully match that of the star. Asteroseismic structure inversions provide a way to test the interior of our stellar models. Recently, structure inversion techniques were used to study 12 stars with radiative cores. In this work, we extend that analysis to 43 main-sequence stars with convective cores observed by Kepler to look for differences in the sound speed profiles in the inner 30% of the star by radius. For around half of our stars, the structure inversions show that our models reproduce the internal structure of the star, where the inversions are sensitive, within the observational uncertainties. For the stars where our inversions reveal significant differences, we find cases where our model sound speed is too high and cases where our model sound speed is too low. We use the star with the most significant differences to explore several changes to the physics of our model in an attempt to resolve the inferred differences. These changes include using a different overshoot prescription and including the effects of diffusion, gravitational settling, and radiative levitation. We find that the resulting changes to the model structure are too small to resolve the differences shown in our inversions.

Elemental abundances, particularly the C/O ratio, are seen as a way to connect the composition of planetary atmospheres with planet formation scenario and the disc chemical environment. We model the chemical composition of gas and ices in a self-gravitating disc on timescales of 0.5\,Myr since its formation to study the evolution of C/O ratio due to dust dynamics and growth, and phase transitions of the volatile species. We use the thin-disc hydrodynamic code FEOSAD, which includes disc self-gravity, thermal balance, dust evolution and turbulent diffusion, and treats dust as a dynamically different and evolving component interacting with the gas. It also describes freeze-out, sublimation and advection of four volatile species: H$_2$O, CO$_2$, CH$_4$ and CO. We demonstrate the effect of gas and dust substructures on the distribution of volatiles and C/O ratios, including the formation of multiple snowlines of one species, and point out the anticorrelation between dust-to-gas ratio and total C/O ratio emerging due to the contribution of oxygen-rich ice mantles. We identify time and spatial locations where two distinct trigger mechanisms for planet formation are operating and differentiate them by C/O ratio range: wide range of the C/O ratios of $0-1.4$ for streaming instability, and a much narrower range $0.3-0.6$ for gravitational instability (with the initial value of 0.34). This conclusion is corroborated by observations, showing that transiting exoplanets, which possibly experienced migration through a variety of disc conditions, have significantly larger spread of C/O in comparison with directly imaged exoplanets likely formed in gravitationally unstable outer disk regions. We show that the ice-phase C/O$\approx0.2-0.3$ between the CO, CO$_2$ and CH$_4$ snowlines corresponds to the composition of the Solar system comets, that represent primordial planetesimals.

In the era of precision cosmology, the ability to generate accurate and large-scale galaxy catalogs is crucial for advancing our understanding of the universe. With the flood of cosmological data from current and upcoming missions, generating theoretical predictions to compare with these observations is essential for constraining key cosmological parameters. While traditional methods, such as the Halo-Occupation Distribution (HOD), have provided foundational insights, they struggle to balance the need for both accuracy and computational efficiency. High-fidelity hydrodynamic simulations offer improved precision but are computationally expensive and resource-intensive. In this work, we introduce a novel machine learning approach that harnesses Convolutional Neural Networks (CNNs) and Diffusion Models, trained on the CAMELS simulation suite, to bridge the gap between computationally inexpensive dark matter simulations and the galaxy distributions of more costly hydrodynamic simulations. Our method not only outperforms traditional HOD techniques in accuracy but also significantly accelerates the simulation process, offering a scalable solution for next-generation cosmological surveys. This advancement has the potential to revolutionize galaxy catalog generation, enabling more precise, data-driven cosmological analyses.

In this work, we use gas phase metallicities calculated from the Sloan Digital Sky Survey (SDSS) Mapping Nearby Galaxies at Apache Point (MaNGA) Data Release 17 (DR17) to assess the extent of potential biases in spaxels which are spatially adjacent to spaxels identified as non-star forming (non-SF) on a BPT diagram. We identify a sample of $\sim21,000$ such spaxels with calculable metallicities from the full metallicity catalogue ($\sim$1.57 million), representing a small fraction ($\sim1.3$ per cent) of the full metallicity sample. $\sim$23 per cent of all galaxies with at least one spaxel with a calculable metallicity also contain at least one spaxel with a calculated metallicity adjacent to a non-SF spaxel, with a typical galaxy hosting 9 non-SF-adjacent spaxels. From our suite of 6 different metallicity calibrations, we find that only the metallicity calibrations based entirely on the [NII]$_{6584}$/H$\alpha$ ratio are affected, showing systematic offsets to higher metallicities by up to $\sim$0.04 dex if they are located adjacent to a non-SF flagged spaxel, relative to a radially matched control sample. The inclusion of additional diagnostic diagrams (based on [OI]$_{6300}$~\&/or [SII]$_{6717+6731}$) is insufficient to remove the observed offset in the [NII]$_{6584}$/H$\alpha$ based calibrations. Using a stricter diagnostic line on the BPT diagram removes $\sim$94 per cent of identified bordering spaxels with metallicities for all metallicity calibrations, and removes the residual offset to higher metallicity values seen in [NII]$_{6584}$/H$\alpha$ calibrations. If science cases demand an exceptionally clean metallicity sample, we recommend either a stricter BPT cut, and/or a non-[NII]$_{6584}$/H$\alpha$ based metallicity calibration.

Recent James Webb Space Telescope observations of cool, rocky exoplanets reveal a probable lack of thick atmospheres, suggesting prevalent escape of the secondary atmospheres formed after losing primordial hydrogen. Yet, simulations indicate that hydrodynamic escape of secondary atmospheres, composed of nitrogen and carbon dioxide, requires intense fluxes of ionizing radiation (XUV) to overcome the effects of high molecular weight and efficient line cooling. This transonic outflow of hot, ionized metals (not hydrogen) presents a novel astrophysical regime ripe for exploration. We introduce an analytic framework to determine which planets retain or lose their atmospheres, positioning them on either side of the cosmic shoreline. We model the radial structure of escaping atmospheres as polytropic expansions - power-law relationships between density and temperature driven by local XUV heating. Our approach diagnoses line cooling with a three-level atom model and incorporates how ion-electron interactions reduce mean molecular weight. Crucially, hydrodynamic escape onsets for a threshold XUV flux dependent upon the atmosphere's gravitational binding. Ensuing escape rates either scale linearly with XUV flux when weakly ionized (energy-limited) or are controlled by a collisional-radiative thermostat when strongly ionized. Thus, airlessness is determined by whether the XUV flux surpasses the critical threshold during the star's active periods, accounting for expendable primordial hydrogen and revival by volcanism. We explore atmospheric escape from Young-Sun Mars and Earth, LHS-1140 b and c, and TRAPPIST-1 b. Our modeling characterizes the bottleneck of atmospheric loss on the occurrence of observable Earth-like habitats and offers analytic tools for future studies.

Plasmoids (or magnetic islands) are believed to play an important role in the onset of fast magnetic reconnection and particle acceleration during solar flares and eruptions. Direct imaging of flare current sheets and formation/ejection of multiple plasmoids in extreme ultraviolet (EUV) images, along with simultaneous X-ray and radio observations, offers significant insights into the mechanisms driving particle acceleration in solar flares. Here we present direct imaging of the formation and ejection of multiple plasmoids in flare plasma/current sheets and associated quasi-periodic pulsations (QPPs) observed in X-ray and radio wavelengths, using observations from SDO/AIA, RHESSI, and Fermi GBM. These plasmoids propagate bidirectionally upward and downward along the flare current sheet beneath the erupting flux rope during two successive flares associated with confined/failed eruptions. The flux rope exhibits evidence of helical kink instability with formation and ejection of multiple plasmoids in the flare current sheet, as predicted in an MHD simulation of a kink-unstable flux rope. RHESSI X-ray images show double coronal sources (``loop-top" and higher coronal sources) located at both ends of the flare current/plasma sheet. Moreover, we detected an additional transient faint X-ray source (6-12 keV) located between the double coronal sources, which was co-spatial with multiple plasmoids in the flare current sheet. X-ray (soft and hard) and radio (decimetric) observations unveil QPPs (periods$\approx$10-s and 100-s) associated with the ejection and coalescence of plasmoids. These observations suggest that energetic electrons are accelerated during the ejection and coalescence of multiple plasmoids in the flare current sheet.

We study the gamma-ray emission from millisecond pulsars within the Milky Way's globular cluster system in order to measure the luminosity function of this source population. We find that these pulsars have a mean luminosity of $\langle L_{\gamma}\rangle \sim (1-8)\times 10^{33}\, {\rm erg/s}$ (integrated between 0.1 and 100 GeV) and a log-normal width of $\sigma_L \sim 1.4-2.8$. If the Galactic Center Gamma-Ray Excess were produced by pulsars with similar characteristics, Fermi would have already detected $N \sim 17-37$ of these sources, whereas only three such pulsar candidates have been identified. We conclude that the excess gamma-ray emission can originate from pulsars only if they are significantly less bright, on average, than those observed within globular clusters or in the Galactic Plane. This poses a serious challenge for pulsar interpretation of the Galactic Center Gamma-Ray Excess.

Variability is a fundamental signature for active galactic nuclei (AGN) activity, and serve as an unbiased indicator for rapid instability happened near the center supermassive black hole (BH). Previous studies showed that AGN variability does not have strong redshift evolution, and scales with their bolometric luminosity and BH mass, making it a powerful probe to identify low-mass, low-luminosity AGNs at high redshift. JWST has discovered a new population of high-redshift galaxies likely hosting moderate accreting BHs ($10^6\,M_\odot$) -- the little red dots (LRDs, $z=4-10$). In this paper, we study the variability of a sample of 21 LRDs with V-shaped SEDs in three JWST deep fields that also have reliable HST observations in closely paired filters at 1-2 um (rest-frame UV), with the time difference between 6 and 11 years. This LRD sample covers a redshift range of $3<z<8$ with median -21.3$<M_\mathrm{UV}<$-18.4. Based on both photometry and imaging difference analyses, we find a mean magnitude difference of 0.12$\pm$0.24 mag, with none of the LRDs showing photometric variability at 3$\sigma$ significance. Extrapolation of SDSS quasar variability predicts a magnitude change of order 0.3 mag for our LRD sample. This suggests an upper limit of about 30\% AGN contribution to the total observed UV light in our sample of LRDs.

Hydrodynamic outflows, such as those observed escaping close-in gas giant planets, are not isothermal in structure. Their highly ionized nature allows them to cool adiabatically at distances beyond several planetary radii. The contrast between the hottest gas temperatures at around 10,000K and the coldest at around 1,000K triggers an excess population of the observable helium triplet. This excess is caused by the suppression of collisional de-excitation from the triplet state at cool temperatures. Using radiation-hydrodynamic simulations, we show that this helium triplet excess may explain the excess broadening seen in HD 189733b's observed transmission spectrum, demonstrating adiabatic cooling of its outflow, confirming its hydrodynamic nature on scales of several planetary radii. However, further observations are required to confirm this conclusion. Furthermore, we explore a range of electron transitions for neutral helium which were not considered in the previous literature. We find that the He$2^1$S state is unavailable as a potential reservoir for He$2^3$S electrons. Additionally, the de-excitation to the ground state must be considered for stellar spectra later than K2 in predicting the correct helium triplet population. Importantly, since triplet helium inherits momentum from ionized helium as it is generated by recombination, it is significantly less prone to fractionation than ground-state neutral helium. However at separations of $\gtrsim 0.05$~au, ionization at the flow base and drag on helium weaken, leading to significant fractionation of the then mostly neutral helium. This in turn, can cause a suppression of the Helium transit depth, even though the helium line width remains large.

In A&A 412, 35 (2003) Blanchard, Douspis, Rowan-Robinson, and Sarkar (BDRS) slightly modified the primordial fluctuation spectrum and produced an excellent fit to WMAP's CMB power spectrum for an Einstein-de Sitter (EdS) universe, bypassing dark energy. Curiously, they obtained a Hubble value of $H_0\approx46$, in sharp conflict with the canonical range $H_0\sim67-73$. However, we will demonstrate that the reduced value of $H_0\approx46$ achieved by BDRS is fully compatible with the use of variable speed of light in analyzing the late-time cosmic acceleration observed in Type Ia supernovae (SNeIa). In arXiv:2412.04257 [gr-qc] we uncovered a hidden aspect in a generic class of scale-invariant actions: the dynamics of the dilaton can induce a variation in the speed of light as $c\propto\chi^{1/2}$, causing $c$ to vary alongside $\chi$ across spacetime. For an EdS universe with varying $c$, besides the effects of cosmic expansion, light waves emitted from distant SNeIa are further subject to a refraction effect, which alters the Lemaitre redshift relation to $1+z=a^{-3/2}$. Based on this new formula, we achieve a fit to the SNeIa Pantheon Catalog exceeding the quality of the $\Lambda$CDM model. Crucially, our approach does not require dark energy and produces $H_0=47.2$ in strong alignment with the BDRS finding of $H_0\approx46$. Hence, BDRS's analysis of the (early-time) CMB power spectrum and our variable-$c$ analysis of the (late-time) Hubble diagram of SNeIa fully agree on two counts: (i) the dark energy hypothesis is avoided, and (ii) $H_0$ is reduced to $\sim47$, which also yields an age $t_0=2/(3H_0)=13.8$ Gy for an EdS universe, without requiring dark energy. Most importantly, we will demonstrate that the late-time acceleration can be attributed to the declining speed of light in an expanding EdS universe, rather than to a dark energy component.

This paper aims to explore the quasinormal modes (QNMs) and effective potential profiles of massless and rotating BTZ black holes within the frameworks of $f(\mathcal{R})$ and Ricci-Inverse ($\mathcal{RI}$) modified gravity theories, which, while producing similar space-time structures, exhibit variations due to distinct cosmological constants, $\Lambda_m$. We derive wave equations for these black hole perturbations and analyze the behavior of the effective potential $V_{\text{eff}}(r)$ under different values of mass $m$, cosmological constant $\Lambda_m$, and modified gravity parameters $\alpha_1$, $\alpha_2$, $\beta_1$, $\beta_2$, and $\gamma$. The findings indicate that increasing mass and parameter values results in a raised potential barrier, implying stronger confinement of perturbations and impacting black hole stability. Incorporating the generalized uncertainty principle, we also study its effect on the thermodynamics of rotating BTZ black holes, demonstrating how GUP modifies black hole radiation, potentially observable in QNM decay rates. Additionally, we investigate the motion of particles through null and timelike geodesics in static BTZ space-time, observing asymptotic behaviors for null geodesics and parameter-dependent shifts in potential for timelike paths. The study concludes that modified gravity parameters significantly influence QNM frequencies and effective potential profiles, offering insights into black hole stability and suggesting that these theoretical predictions may be tested through gravitational wave observations.

The existence of dark matter has long been extensively studied in the past few decades. In this study, we investigate the emission of neutrinos and photons from neutron stars (NSs) by employing the modified theory of gravity and the corresponding Tolman-Oppenheimer-Volkoff (TOV) system of equations. The extreme matter density and magnetic field inside the NSs provide a unique laboratory for studying fundamental physics, including the interplay between gravity and quantum field effects. The impact of a strong magnetic field has also been incorporated into the corresponding TOV equations. We here attempt to see how neutrinos and photons emissions from these compact objects are impacted by the modified TOV equations due to modified theory of gravity; f(R,T) gravity or scalar-tensor theory and strong magnetic fields. Our analysis focuses on how these modifications influence the structure, cooling, and photon/neutrino luminosities of NS. We computed the surface temperature of NSs for normal Einstein gravity and modified gravity theories with and without magnetic field for three EoSs; namely APR, FPS and SLY. On comparison of our predicted values of surface temperature with the observed surface temperature for three NSs, we find that modified gravity along with inside magnetic field-based predictions shows reasonable agreement with the corresponding observed values.

The quest to understand gravity's role in shaping the universe has led to the exploration of modified gravity theories. One such theory is Myrzakulov gravity, which incorporates both curvature and torsion. In this work, we investigate the effects of torsion within the framework of $f(R,T)$-gravity, a modification of General Relativity that includes both curvature and torsion. We present this theory in the Vielbein formalism, which offers a more flexible, geometric perspective on gravity. This formalism is particularly useful in Weitzenb"ock spacetime, where torsion significantly influences gravitational interactions. Our study aims to extend Myrzakulov's theory and explore its cosmological and astrophysical implications. We examine the effects of torsion on phenomena like black holes, gravitational waves, and neutron stars. These modifications could lead to observable deviations in black hole thermodynamics, gravitational wave propagation, and dense matter structure. This work provides new insights into the nature of spacetime and gravity, offering a fresh perspective on the fundamental forces governing the cosmos. It paves the way for future studies, both observational and theoretical, that could reveal new physics beyond General Relativity.

We consider a thermodynamically consistent approach for the computation of the masses, radii, and tidal deformabilities of compact stars consisting of two interacting fluids with separately conserved quantum numbers. We apply this interacting fluid approach to the case of compact stars of neutron star matter with the Higgs portal fermionic dark matter model for the first time in a thermodynamically consistent manner. The patterns for the mass-radius curves and the tidal deformability depend on the dark matter particle mass and are different from former studies. Compared to ordinary neutron star properties, we obtain smaller masses and radii for dark matter particle masses similar to the nucleon mass and, hence, smaller tidal deformabilities as a result of the softening of the equation of state due to the presence of dark matter. For dark matter particle masses below the nucleon mass and sizable chemical potentials with respect to the dark matter particle mass, there will be a dark halo instead of dark core. Our investigation provides the basis for studying mergers of compact stars where the two fluids of neutron star matter and dark matter are coupled kinetically to each other and are described by one combined energy-momentum tensor of the two interacting fluids but are chemically different with two separately conserved number currents.

The gravitational wave detectors used by the LIGO Scientific Collaboration, and the Virgo Collaboration are incredibly sensitive instruments which frequently detect non-stationary, non-Gaussian noise transients. iDQ is a statistical inference framework which leverages the use of auxiliary degrees of freedom monitored in the detectors to identify such transients. In this work, we describe the improvements to the iDQ pipeline made between the third and fourth observing run of the LIGO-Virgo-KAGRA (LVK) collaboration, and show the performance of these changes. We find that iDQ detects a total of 39,398 of the known 100,512 glitches identified by Omicron over the course of the second half of the third observing run. We construct a measure of the probability a glitch is present in the strain data of a given detector by combining information from iDQ and Omicron as well as extend the output of iDQ in a novel method which finds correlations between known glitch classifications, and auxiliary channels. We identify several channels over the course of O3b which frequently record instances of Scattered Light, Whistle, and Blip glitches and discuss use cases for this method in active observing runs.

After decades of experimental efforts, the DAMA/LIBRA(DL) annual modulation (AM) analysis on the ${\chi}$N (WIMP Dark Matter interactions on nucleus) channel remains the only one which can be interpreted as positive signatures. This has been refuted by numerous time-integrated (TI) and AM analysis. It has been shown that ${\chi}$e (WIMP interactions with electrons) alone is not compatible with the DL AM data. We expand the investigations by performing an AM analysis with the addition of ${\chi}$e long-range and short-range interactions to ${\chi}$N, derived using the frozen-core approximation method. Two scenarios are considered, where the ${\chi}$N and ${\chi}$e processes are due to a single ${\chi}$ (${\Gamma}^{1\chi}_{tot}$) or two different ${\chi}$s (${\Gamma}^{2\chi}_{tot}$). The combined fits with ${\chi}$N and ${\chi}$e provide stronger significance to the DL AM data which are compatible with the presence of additional physical effects beyond \c{hi}N alone. This is the first analysis which explores how ${\chi}$e AM can play a role in DL AM. The revised allowed regions as well as the exclusion contours from the other null AM experiments are presented. All DL AM allowed parameter spaces in ${\chi}$N and ${\chi}$e channels under both ${\Gamma}^{1\chi}_{tot}$ and ${\Gamma}^{2\chi}_{tot}$ are excluded at the 90\% confidence level by the combined null AM results. It can be projected that DL-allowed parameter spaces from generic models with interactions induced by two-WIMPs are ruled out.

The results of simultaneous measurements of noctilucent clouds (NLC) position in a number of ground-based locations are presented. Observational data of 14 bright NLC events over 5 years is used for building the altitude maps of cloud fields using triangulation technique updated for multi-location case. Statistical distribution of NLC altitude and its change during the summer season is considered. Mean NLC altitudes are compared with colorimetric technique based on the same data and simple radiation transfer model. This can be used to check the model and estimate the accuracy of single-camera technique of NLC altitude measurements. Results and methods are suggested for net ground-based survey of noctilucent clouds.

A crossover involving three-fermion clusters is relevant to the hadron-quark crossover, which, if occurring in a neutron star, could naturally reproduce the dense-matter equation of state recently deduced from simultaneous observations of neutron-star masses and radii. To understand the crossover mechanism, we examine the role of tripling fluctuations induced by the formation of three-fermion clusters. The phase-shift representation of fluctuations manifests an interplay of bound and scattering states, leading to non-monotonic momentum distributions of baryon-like clusters and peaked sound speed at finite densities. We demonstrate them by applying our approach to a nonrelativistic system of one-dimensional three-color fermions analogous to the hadron-quark matter.

We derive the differential age signal valid for cosmic chronometers (passively evolving galaxies) in any space-time that satisfies the following assumptions: (i) The space-time has a metric with Lorentzian signature and the connection is the Levi-Civita connection; (ii) the cosmic chronometers are collectively well approximated as a geodesic and irrotational congruence of time-like worldlines in the space-time; (iii) light travels on null geodesics and caustics on the observer's past light cone can be ignored; (iv) the space-time is cosmological, meaning that isotropic and positive expansion degrees-of-freedom dominate over anisotropic and negative expansion degrees-of-freedom when viewed on sufficiently large scales in the frame of the cosmic chronometers. The main result of the paper is an expression for the differential age signal that is written in terms of line-of-sight averages of the expansion rate along individual null lines, thus providing a kinematic interpretation of the differential age signal applicable to cosmological space-times satisfying (i)--(iv). We explain how this result indicates that the differential age signal is a robust probe of the volume-average expansion rate in very general statistically homogeneous and isotropic space-time scenarios where other probes of the volume-average expansion rate tend to yield biased results. We argue that this unique property of the differential age signal makes it an ideal measurement for constraining the expansion history model-independently.

We investigate the influence of dark matter on hybrid stars. Using a two-fluid approach, where normal and dark matter components interact only gravitationally, we explore how dark matter can trigger the appearance of quark matter in neutron stars for unprecedented low masses. Our findings reveal that dark matter increases the central pressure of neutron stars, potentially leading to the formation of hybrid stars with quark cores even at very low compact star masses. The critical mass for the appearance of quark matter decreases with increasing dark matter content. We introduce the concept of "masquerading hybrid stars", where dark matter admixed stars exhibit similar mass-radius relations to purely hadronic stars, making it challenging to distinguish between them based solely on these parameters. Additionally, we identify a unique class of objects termed "dark oysters", characterized by a large dark matter halo and a small normal matter core, highlighting the diverse structural possibilities for compact stars influenced by dark matter.