Papers for Tuesday, Mar 11 2025

We present a dataset of high resolution spectra of the Sun and ten bright stars of the domains [6297-6303] and [6361-6367] Angtr{ö}m. Solar spectra were obtained in the quiet Sun at various distances from disk centre with the ground based Meudon Solar Tower and Themis telescope (12 mÅ resolution) and with the Solar Optical Telescope (SOT) onboard the Hinode satellite (21 mÅ resolution). Spectra of 10 bright stars (magnitude \< 2) were also got with Themis at 12 mÅ resolution. These spectral domains contain the faint and forbidden OI lines (6300.31 Å and 6363.79 Å) that are useful for the research of Oxygen abundance. The spectra shown here are freely available in FITS format to the research community.

The Immersion GRating INfrared Spectrometer (IGRINS) is a compact, high-resolution (R~45,000) near-infrared spectrograph spanning 1.45 to 2.45 um in a single exposure. We introduce the Raw and Reduced IGRINS Spectral Archive (RRISA), which provides public data access for all non-proprietary IGRINS data taken at McDonald Observatory's Harlan J. Smith Telescope, the Lowell Discovery Telescope (formerly Discovery Channel Telescope), and Gemini South. RRISA provides access to raw files, reduced data products, and cross-matched IGRINS targets with the SIMBAD, 2MASS, Gaia DR3, APOGEE2 DR17, and PASTEL catalogs. We also introduce version 3 of the IGRINS data reduction pipeline, IGRINS PLP v3, which implements an improved cosmic ray correction, pattern noise removal, and a new flexure correction that reduces telluric residuals. RRISA and supporting information can be found at this http URL.

The level of dust vertical settling and radial dust concentration in disks is of critical importance for understanding the efficiency of planet formation. We present the first uniform analysis of the vertical extent of millimeter dust for a representative sample of 33disks. We used radiative transfer modeling of archival high-angular-resolution (<=0.1") ALMA dust observations of inclined and ringed disks to estimate their vertical dust scale height, which was compared to estimated gas scale heights to characterize the level of vertical sedimentation. In all 23systems for which constraints could be obtained, we find that the outer parts of the disks are vertically settled. 5disks allow for the characterization of the dust scale height both within and outside approximately half the dust disk radius, showing a lower limit on their dust heights at smaller radii. This implies that the ratio between vertical turbulence and the Stokes number, $\alpha_z/\St$, decreases radially in these sources. For 21rings in 15disks, we also constrained the level of radial concentration of the dust, finding that about half of the rings are compatible with strong radial trapping. In most of these rings, vertical turbulence is found to be comparable to or weaker than radial turbulence, which is incompatible with the turbulence generated by the vertical shear instability at these locations. We further used our dust settling constraints to estimate the turbulence level under the assumption that the dust size is limited by fragmentation, finding typical upper limits around $\alpha_\text{frag}\leq10^{-3}$. In a few sources, we find that turbulence cannot be the main source of accretion. In the context of pebble accretion, we identify several disk regions that have upper limits on their dust concentration that would allow core formation to proceed efficiently, even at wide orbital distances outside of 50au.

Red supergiants (RSGs) are cool and evolved massive stars exhibiting enhanced mass loss compared to their main sequence phase, affecting their evolution and fate. However, the theory of the wind-driving mechanism is not well-established and the metallicity dependence has not been determined. We aim to uniformly measure the mass-loss rates of large samples of RSGs in different galaxies with $-0.7\lesssim[Z]\lesssim0$ to investigate whether there is a potential correlation with metallicity. We collected photometry from the ultraviolet to the mid-infrared for all our RSG candidates to construct their spectral energy distribution (SED). Our final sample includes 893 RSG candidates in the Small Magellanic Cloud (SMC), 396 in NGC 6822, 527 in the Milky Way, 1425 in M31, and 1854 in M33. Each SED was modelled using the radiative transfer code DUSTY under the same assumptions to derive the mass-loss rate. The mass-loss rates range from approximately $10^{-9} \ M_{\odot}$ yr$^{-1}$ to $10^{-5} \ M_{\odot}$ yr$^{-1}$ with an average value of $1.5\times10^{-7} \ M_{\odot}$ yr$^{-1}$. We provided a new mass-loss rate relation as a function of luminosity and effective temperature for both the SMC and Milky Way and compared our mass-loss rates with those derived in the Large Magellanic Cloud (LMC). The turning point in the mass-loss rate vs. luminosity relation differs by around 0.2 dex between the LMC and SMC. The mass-loss rates of the Galactic RSGs at $\log(L/L_\odot)<4.5$ were systematically lower than those determined in the other galaxies, possibly due to uncertainties in the interstellar extinction. We found 30-40% of the RSGs not to have any dust. The results for M31 and M33 are inconclusive because of source blending at distances above 0.5 Mpc given the resolution of Spitzer. Overall, we found similar mass-loss rates among the galaxies, indicating no strong correlation with metallicity.

Recent measurements of the Milky Way rotation curve found a sharp decline at around $15$-$20$ kpc from the center of the Galaxy, suggesting that the Galactic dark matter halo is much less massive than predicted by other dynamical tracers. To address this tension, we study the validity of the assumptions made in calculating the Milky Way's rotation curve. To do so, we apply Jeans' equation, the current standard approach of measuring rotation curves, to three cosmological zoom-in simulations of Milky Way-like galaxies from the FIRE-2 Latte suite. Using synthetic Gaia surveys, we replicate the sample selection process and calculation employed in measuring the Milky Way rotation curve. We examine four failure modes of this calculation and find that the measured curves deviate from the true curve by $5$-$20\%$ rather than below $5\%$, as estimated by previous works. Interestingly, there is a large galaxy-to-galaxy variance, and different systematics dominate different galaxies. We rederive the Milky Way's dark matter density profile with the rotation curve while incorporating systematics from the simulations. The posterior distribution of the density profiles is consistent with a fiducial NFW profile when assuming a gNFW profile for dark matter. We find that the virial mass, $7.32^{+1.98}_{-1.53}\times10^{11}~M_{\odot}$, consistent with other probes of the Milky Way's mass. However, we recommend that the field moves away from relying solely on the rotation curve when studying the dark matter profile, and adopts methods that incorporate additional probes and/or do not heavily depend on assumptions described in this study.

We present hydrodynamic simulations of a flavour-mixed two-component dark matter (2cDM) model that utilize IllustrisTNG baryonic physics. The model parameters are explored for two sets of power laws of the velocity-dependent cross sections, favoured on the basis of previous studies. The model is shown to suppress the formation of structures at scales $k\gtrsim 10^2\ h\text{ Mpc}^{-1}$ up to 40\% compared to cold dark matter (CDM) at redshifts $z\sim5-2$. We compare our results to structure enhancement and suppression due to cosmological and astrophysical parameters presented in the literature and find that 2cDM effects remain relevant at galactic and subgalactic scales. The results indicate the robustness of the role of nongravitational dark matter interactions in structure formation and the absence of putative degeneracies introduced by baryonic feedback at high $z$. The predictions made can be further tested with future Ly-$\alpha$ forest observations.

Type IIb supernovae (SNe IIb) often exhibit an early light curve excess (EE) preceding the main peak powered by radioactive nickel decay. The physical origin of this early emission remains an open question. Among the proposed scenarios, shock cooling emission-resulting from the interaction between the shockwave and extended envelopes-is considered the most plausible mechanism. The frequency of these events remains unconstrained. This study aims to quantify the frequency of EE in SNe IIb and investigate its physical origin by analyzing optical light curves from the Asteroid Terrestrial-impact Last Alert System (ATLAS) survey. We selected 74 SNe IIb from 153 spectroscopically classified events in the Transient Name Server (TNS) database, observed by ATLAS, with peak fluxes exceeding 150 {\mu}Jy and explosion epoch uncertainties lower than six days. Using light curve model fitting and outlier analysis, we identified SNe IIb exhibiting EE and analyzed their photometric properties. We found 21 SNe IIb with EE, corresponding to a frequency of approximately 28-40%, with the higher value obtained under the most stringent data cuts. The EE's duration and color evolution are consistent with shock cooling in extended hydrogen-rich envelopes. We also found that EE SNe IIb have longer rise times and faster post-peak decline rates than non-EE SNe IIb, while both groups share similar peak absolute magnitudes. Our findings suggest that EE and non-EE SNe IIb likely share similar initial progenitor masses but differ in ejecta mass properties, potentially due to varying degrees of binary interaction. This study provides constraints on the evolutionary pathways of SNe IIb progenitors as compact stars with and without extended hydrogen envelopes.

We present Gemini/GHOST high-resolution spectra of five stars observed in two low surface brightness Milky Way satellites, Sagittarius II (Sgr2) and Aquarius II (Aqu2). For Aqu2, the velocities and metallicities of the two stars are consistent with membership in a dark matter-dominated ultra faint dwarf galaxy (UFD). The chemical abundance ratios suggest inefficient star formation from only one or a few supernovae (e.g., low Na, Sr, Ba), and enriched potassium (K) from super-AGB stars. For Sgr2, the velocity and metallicity dispersions of its members are not clearly resolved and our detailed chemical abundances show typical ratios for metal-poor stars, with low dispersions. There is only one exception - we report the discovery of an r-process enhanced star (Sgr2584, [Eu/Fe]= $+0.7 \pm 0.2 $; thus, an r-I star). As r-I stars are found in both UFDs (Tuc III, Tuc IV, Grus II) and globular clusters (M15 and M92), then this does not help to further classify the nature of Sgr2. Our exploration of Sgr2 demonstrates the difficulty in classifying some of the faintest (ambiguous) satellites. We advocate for additional diagnostics in analysing the ambiguous systems, such as exploring radial segregation (by mass and/or chemistry), N-body simulations, and the need for dark matter to survive Galactic tidal effects. The spectra analysed in this paper were taken as part of the GHOST commissioning observations, testing faint observation limits (G < 18.8) and the single and double IFU observing modes.

In this paper, we present the Discovery simulations: a new pair of high-resolution N-body simulations motivated by the DESI Y1 BAO cosmological constraints on dark energy. The Discovery simulations were run with identical initial conditions, and differ only in their cosmological parameters. The first simulation is based on a flat $\Lambda\mathrm{CDM}$ cosmology, while the second is based on a $w_0 w_a\mathrm{CDM}$ cosmology, with particular parameter values chosen based on the DESI analysis which includes constraints from BAO with CMB priors. Both simulations evolve $6720^3$ particles in a box with a side length of $L_\mathrm{box} = 1.5$ Gpc, leading to a mass resolution of $\sim4 \times 10^8$ $\mathrm{M}_{\odot}$ in each simulation. In this work we demonstrate the impact of the $w_0 w_a\mathrm{CDM}$ cosmology on the matter power spectrum, halo mass function, and halo mass accretion rate. We also populate halos with galaxies using a novel forward model for in-situ star formation, and examine the way in which changes to cosmology manifest as changes in star formation history. The Discovery simulations provide a testbed for alternative cosmological probes that may offer additional constraining power beyond BAO, such as higher-order summary statistics and observables in the nonlinear regime. Halo catalogs from the Discovery simulations are publicly available and can be downloaded from the HACC Simulation Data Portal.

V1298 Tau is a very young and magnetically active star which hosts a benchmark multi-planetary system to study planet formation and evolutionary history at the earliest stages. We selected V1298 Tau for a first targeted follow-up at radio frequencies with the Karl G. Jansky Very Large Telescope (JVLA), the upgraded Giant Metrewave Radio Telescope (uGMRT), and the Sardinia Radio Telescope (SRT), to search for emission in the overall frequency range 0.55-7.2 GHz. Detecting radio emission from such a very active star is key to characterise its magnetosphere, allowing in principle to probe the strength of the coronal magnetic field and plasma density. Observations were carried out between Oct 2023 and Sept 2024: three epochs with uGMRT band-4 (0.55-0.75 GHz), 12 epochs with the JVLA using L (1-2 GHz) and C (4.5-6.5 GHz) bands, and three epochs with SRT using C-high band (6-7.2 GHz). We report the first detection of radio emission from V1298 Tau at different epochs using the JVLA. The emission has maximum peak flux densities of 91$\pm$10 and 177$\pm$6 $\mu$Jy/beam in the L- and C-band, respectively. From a comparison with contemporary optical photometry, we found that the detected emission with the highest fluxes are located around a phase of minimum of the photospheric light curve. Although the uGMRT and SRT observations could not detect the source, we measured 3$\sigma$ flux density upper limits in the range ~41-56 $\mu$Jy/beam using uGMRT, while with SRT we reached upper limits down to 13 mJy. The lack of a significant fraction of circular polarisation indicates that the observed flux is not due to electron cyclotron maser emission from star-planet interaction, and it is likely produced by gyrosynchroton/cyclotron emission from the corona triggered by stellar magnetic activity, although we cannot exclude thermal emission due to a lack of constraints on the brightness temperature.

Accurate neutrino transport is crucial for reliably modeling explosive astrophysical events like core-collapse supernovae (CCSNe) and neutron star mergers (NSMs). However, in these extremely neutrino-dense systems, flavor oscillations exhibit challenging nonlinear effects rooted in neutrino-neutrino forward scattering. Evidence is quickly accumulating that these collective phenomena can substantially affect explosion dynamics, neutrino signals, nucleosynthesis, and kilonova light curves. We review the progress made so far on the difficult and conceptually deep question of how to correctly include this physics in simulations of CCSNe and NSMs. Our aim is to take a broad view of where the problem stands, and so provide a critical assessment of where it is headed.

Widespread galactic winds emanate from the Large Magellanic Cloud (LMC), with the 30 Doradus starburst region generating the fastest and most concentrated gas flows. We report on the gas distribution, kinematics, and ionization conditions of the near-side outflow along 8 down-the-barrel sightlines using UV absorption-line observations from the HST's ULLYSES program for this region along with H I 21-cm observations from the GASS and GASKAP surveys. We find that within 1.7 degrees from the center of 30 Doradus, the wind reaches maximum speeds of $100-150\,\text{km}\,\text{s}^{-1}$ from the LMC's disk. The total integrated column densities of low-ions (O I, Si II, and Fe II) in the blueshifted wind, up to $v_{\rm LSR}=150\,\text{km}\,\text{s}^{-1}$, are highest near the center and decline radially outward. We estimate an outflow mass of $M_{\rm outflow,\,Si II}\approx(5.7-8.6)\,\times 10^{5} M_{\odot}$, outflow rate of $\dot{M}_{\rm outflow}\gtrsim0.02 M_{\odot}\,\text{yr}^{-1}$, and mass loading factor of $\eta\gtrsim0.10$ within 0.52 degrees from the center of 30 Doradus. The observed ion ratios$-$together with photoionization modeling$-$reveal that this wind is roughly $40-97\%$ photoionized. The metallicities and dust depletion patterns of the high-velocity absorbers at $v_{\rm LSR}\approx+120\,\text{km}\,\text{s}^{-1}$ can be explained by either a foreground Milky Way (MW) halo cloud or an outflow from the LMC. For the high-ions, Si IV and C IV are broader and kinematically offset from the low-ions, suggesting turbulent mixing layers (TMLs) existing in the wind. Finally, our hydrodynamical simulations of the Magellanic Clouds (MCs) and MW system suggest that the Magellanic Corona can protect the LMC winds from the ram-pressure forces exerted by the MW's halo.

How impulsive solar energetic particle (SEP) events are produced by magnetic-reconnection-driven processes during solar flares remains an outstanding question. Here we report a short-duration SEP event associated with an X-class eruptive flare on July 03, 2021, using a combination of remote sensing observations and in situ measurements. The in situ SEPs were recorded by multiple spacecraft including the Parker Solar Probe. The hard X-ray (HXR) light curve exhibits two impulsive periods. The first period is characterized by a single peak with a rapid rise and decay, while the second period features a more gradual HXR light curve with a harder spectrum. Such observation is consistent with in situ measurements: the energetic electrons were first released during the early impulsive phase when the eruption was initiated. The more energetic in situ electrons were released several minutes later during the second period of the impulsive phase when the eruption was well underway. This second period of energetic electron acceleration also coincides with the release of in situ energetic protons and the onset of an interplanetary type III radio burst. We conclude that these multi-messenger observations favor a two-phase particle acceleration scenario: the first, less energetic electron population was produced during the initial reconnection that triggers the flare eruption, and the second, more energetic electron population was accelerated in the above-the-looptop region below a well-developed, large-scale reconnection current sheet induced by the eruption.

Hickson compact groups (HCGs) offer an ideal environment for investigating galaxy transformation as a result of interactions. It has been established that the evolutionary sequence of HCGs is marked by an intermediate stage characterised by a substantial amount of HI in their intragroup medium (IGrM) in the form of tidal tails and bridges (Phase 2), rapidly followed by a final stage where no IGrM gas is found and where their member galaxies are highly HI-deficient (Phase 3). Despite numerous single-dish and interferometric HI studies on the HCGs, a clear HI picture of the groups within their large-scale environment still remains to be uncovered. Taking advantage of the MeerKAT's high column density sensitivity and large field-of-view, we aim to investigate the rapid transformation of HCGs from the intermediate to late phases, and establish a picture of their gas content variations in the context of their large-scale environments. We performed MeerKAT observations of six HCGs selected to represent the intermediate and late phases of the proposed evolutionary sequence. Combining the HI observations with data from recent wide-field optical surveys, we evaluated the HI deficiencies of galaxies in a ~30' radius of the HCGs. We find that galaxies surrounding both phases exhibit similar distributions in their gas content. Similarly, galaxies making up the cores of Phase 2 HCGs are comparable to their neighbours in terms of HI deficiencies. However, Phase 3 groups are over an order of magnitude more deficient than their surroundings, supporting previous findings that late-phase HCG galaxies are more evolved than their large-scale environments.

The Disk Detective project, a citizen science initiative, aims to identify circumstellar discs around stars by detecting objects with infrared (IR) excess using data from the Wide-field Infrared Survey Explorer (WISE). In this study, we investigate SIPS J2045-6332, a potential brown dwarf with significant IR excess in WISE and 2MASS bands, initially identified by project volunteers. Despite early indicators of a circumstellar disc, discrepancies between observed brightness and expected Spectral Energy Distribution (SED) models suggested unusual properties. To explore potential explanations, we created SED templates for spectral types M9 to L4 and compared them with SIPS J2045-6332's photometric data, revealing an excess brightness that points to either an unresolved low-mass companion or a young, inflated primary star. Further analysis of infrared spectral features and surface gravity indicators supports a youthful classification, estimating the object's age at 26-200 million years. Observations also suggest the presence of a mid L-type companion at a projected distance of 6.7 AU. This study highlights SIPS J2045-6332 as an intriguing system with unique IR characteristics and recommends follow-up observations with high-resolution telescopes to confirm the companion hypothesis and further characterize the system.

The Vera C. Rubin Observatory will conduct the Legacy Survey of Space and Time (LSST), promising to discover billions of galaxies out to redshift 7, using six photometric bands ($ugrizy$) spanning the near-ultraviolet to the near-infrared. The exact number of and quality of information about these galaxies will depend on survey depth in these six bands, which in turn depends on the LSST survey strategy: i.e., how often and how long to expose in each band. $u$-band depth is especially important for photometric redshift (photo-z) estimation and for detection of high-redshift Lyman-break galaxies (LBGs). In this paper we use a simulated galaxy catalog and an analytic model for the LBG population to study how recent updates and proposed changes to Rubin's $u$-band throughput and LSST survey strategy impact photo-z accuracy and LBG detection. We find that proposed variations in $u$-band strategy have a small impact on photo-z accuracy for $z < 1.5$ galaxies, but the outlier fraction, scatter, and bias for higher redshift galaxies varies by up to 50%, depending on the survey strategy considered. The number of $u$-band dropout LBGs at $z \sim 3$ is also highly sensitive to the $u$-band depth, varying by up to 500%, while the number of $griz$-band dropouts is only modestly affected. Under the new $u$-band strategy recommended by the Rubin Survey Cadence Optimization Committee, we predict $u$-band dropout number densities of $110$ deg$^{-2}$ (3200 deg$^{-2}$) in year 1 (10) of LSST. We discuss the implications of these results for LSST cosmology.

We propose a novel mechanism for angular momentum (AM) exchange between the crust and core of a neutron star (NS) via the gyromagnetic effect. Using extended hydrodynamics, we model the star by incorporating macroscopic AM and microscopic AM originating from neutron orbital and spin AM. We reveal that macroscopic dynamics in the crust can inform microscopic AM in the core leading to neutron spin polarization, and offer alternative scenario of (anti-)glitches. This work highlights the overlooked multi-scale AM interconversions in NS physics, paving the way for gyromagnetic astrophysics.

Constrained measurements of fundamental physical constants using astronomical observational data represent a powerful method for investigating potential new physics. In particular, the dispersion measure (DM) of fast radio bursts (FRBs), which probes the electron density along their propagation paths, may be influenced by the space-time variation of the fine-structure constant \(\alpha\). In this study, we analyze the cross-correlation signal between foreground galaxies and the DM of background FRBs to constrain the evolution of \(\alpha\). Assuming large-scale structure (LSS) galaxy surveys with the capabilities of the China Space Station Telescope (CSST) at \(z=0.15\) and { a mock FRB survey with \(N_{\text{FRB}}=10^5\) at \(z=0.4\), we test how well \(\alpha\) variation can be constrained}, with a standard deviation of \(\sigma(\Delta \alpha / \alpha) = 0.0007\) at \(z=0.15\). Furthermore, taking into account the nonminimal coupling between the scalar field and the electromagnetic field, the variation in \(\alpha\) can lead to the non-conservation of photon number along geodesics. This would result in a violation of the CDDR and affect the evolution of the Cosmic Microwave Background (CMB) temperature. In this work, we { obtain constraints results} on the CDDR parameter \(\eta\) and the parameter \(\beta\) governing CMB temperature evolution at \(z=0.15\), yielding \(\sigma(\eta) = 0.0004\) and \(\sigma(\beta) = 0.0006\), respectively. Finally, we relate the variation in \(\alpha\) to the time evolution of the proton-to-electron mass ratio, { reporting a standard deviation} of \(\sigma(\Delta \mu/\mu) = 0.002\) at $z=0.15$. Future FRB surveys hold significant potential for advancing our understanding of the evolution of fundamental physical constants.

This paper presents a new model-independent constraint on the Hubble constant ($H_0$) by anchoring relative distances from Type Ia supernovae (SNe Ia) observations to absolute distance measurements from time-delay strong Gravitational Lensing (SGL) systems. The approach only uses the validity of the cosmic distance duality relation (CDDR) to derive constraints on $H_0$. By using Gaussian Process (GP) regression to reconstruct the unanchored luminosity distance from the Pantheon$+$ compilation to match the time-delay angular diameter distance at the redshift of the lenses, one yields a value of \mbox{$H_0 = 75.57 \pm 4.415$ km/s/Mpc} at a 68\% confidence level. The result aligns well with the local estimate from Cepheid variables within the $1\sigma$ confidence region, indicating consistency with late-universe probes.

We examine mass outflows from a relativistic, viscous, advective, and magnetized accretion disk around a rotating black hole in presence of thermal conduction. We consider the disk is primarily threaded by the toroidal component of the magnetic field and an effective potential satisfactorily mimicked the spacetime geometry around the rotating black hole. With this, we self-consistently solve the coupled governing equations for inflow and outflow and compute the mass outflow rate $R_{\dot \rm m}$ (ratio of mass flux of inflow to outflow) in terms of the inflow parameters, namely energy ($\mathcal{E}$), angular momentum ($\lambda$), plasma-$\beta$ and conduction parameter ($\Upsilon_{\rm s}$) around weakly rotating ($a_{\rm k} \rightarrow 0$) as well as rapidly rotating ($a_{\rm k} =0.99$) black holes. Our findings reveal that the present formalism admits coupled inflow-outflow solutions across a wide range of inflow parameters yielding substantial mass loss. We observe that $R_{\dot{\rm m}}$ monotonically increases with $\Upsilon_{\rm s}$, irrespective of black hole spin. We also find that for a fixed $\Upsilon_{\rm s}$, when energy, angular momentum, and magnetic field strength of the inflowing matter is increased, $R_{\dot{\rm m}}$ is enhanced resulting the outflows even more pronounced. We further estimate the maximum outflow rate ($R^{\rm max}_{\dot{\rm m}}$) by varying the inflow parameters and find that thermal conduction leads to maximum mass outflow rate $R^{\rm max}_{\dot{\rm m}} \sim 24\%$ for rapidly rotating black hole of spin $a_{\rm k} = 0.99$. Finally, we employ our formalism to explain the kinetic jet power of $68$ radio-loud low-luminosity active galactic nuclei (LLAGNs), indicating that it is potentially promising to account for the observed jet power of substantial number of LLAGNs.

Context. Radical chemical reactions on cosmic dust grains play a crucial role in forming various chemical species. Among different radicals, the hydroxyl (OH) is one of the most important ones, with a rather specific chemistry. Aims. The goal of this work is to simulate the recombination dynamics of hydroxyl radicals and the subsequent formation of hydrogen peroxide (H$_{2}$O$_{2}$). Mthods. We employed neural network potentials trained on ONIOM(QM/QM) data, combining multi-reference (CASPT2) and density functional theory (DFT) calculations. This approach allows us to model the recombination of hydroxyl radicals on ice surfaces with high computational efficiency and accuracy. Results. Our simulations reveal that the initial position of the radicals plays a decisive role in determining recombination probability. We found that the formation of hydrogen-bond between radicals, competes with the formation of hydrogen peroxide, reducing the recombination efficiency, contrary to expectancy. This competition reduces the recombination probability for radicals that are initially formed approximately 3 Å apart. Recombination probabilities also depend on the kinetic energy of the added radicals, with values around 0.33 for thermal radicals and a wide range of values between 0.33 and 1.00 for suprathermal OH radicals. Conclusions. Based on our calculations we provide recommendations for introducing OH radical recombination into kinetic astrochemical models, differentiating between thermal and suprathermal radicals. The recombination behavior varies significantly between these two cases: while thermal radicals are sometimes trapped in hydrogen-bonded minima, the case of suprathermal radicals varies with the added energy. Our most important conclusion is that OH radical recombination probability cannot be assumed 1.0 for a wide variety of cases.

A number of polarization estimators have been developed for a variety of astrophysical applications to compensate measurements of linear polarization for a bias contributed by the instrumental noise. Most derivations of the estimators assume that the amplitude and orientation of the polarization vector are constant. This assumption generally is not valid for the radio emission from pulsars that fluctuates from pulse to pulse. The radio emission from pulsars, fast radio bursts, and magnetars can be elliptically polarized, and estimators of the total polarization and absolute value of the circular polarization are used in their observations. However, these estimators have not been formally developed to a level that is commensurate with those of linear polarization. Estimators are derived for circular, linear, and total polarization when the amplitude of the polarization vector is a constant or a random variable. Hybrid estimators are proposed for general application to pulsar polarization observations. They are shown to be more effective at removing instrumental noise than their commonly used counterparts.

The leaky-box model and the attendant concept of path-length distribution of cosmic rays were invented in the mid-1960's. Even though versatile computational packages such as GALPROP and DRAGON with the diffusion approach are now available for analyzing cosmic ray data, the concepts of the leaky-box and path-length distribution continue to be adopted extensively. We show here mathematically that there is a close correspondence between the two approaches: The path-length or resident-time of the leaky-box models are similar to 'impulse response functions' of complex dynamical systems and are intuitively transparent. The results provided by the leaky-box model are valid when used judiciously.

We present for the first time an analysis of high-frequency gravitational wave (GW) emission from proto-neutron stars (PNS) in core collapse supernovae (CCSNe) that combines spatial decomposition and modal decomposition to both source and characterize the emission. Our analysis is based on three-dimensional CCSN simulations initiated from two progenitors with differing mass and metallicity. We spatially decompose GW strains into five regions and show they are initially largest in the PNS surface layers from accretion and later largest from the Ledoux convective and convective overshoot regions within the PNS. We compute the fractional GW luminosity and observe that the majority of the luminosity moves from the same outer layers to deep within the PNS at comparable times. Using a self-consistent perturbative analysis, we investigate the evolution of the oscillation modes of the PNS. We find that the frequency of the evolving high-frequency component of the GW signal is well matched to the ${}^2g_2$-mode, the ${}^2g_1$-mode, and the ${}^2f$-mode over time. We show that the ${}^2g$-modes emit most of their power in GWs initially from the PNS surface region, but within a few 100 ms after bounce, it is the convective overshoot region of the PNS that emits the most GW power for the ${}^2g_1$-mode. Eventually, the ${}^2f$-mode is the dominant mode producing GWs, and they are emitted primarily from the convective overshoot region. Thus, with three interconnected analyses, we show that, while the GW emission is global, stemming from multiple regions in and around the PNS, we are able to source the dominant contributions to it. We find that high-frequency GW emission from the PNS in CCSNe is more complex than assessed by other methods, and dependent, first emitted mainly by ${}^2g$-modes driven by accretion onto the PNS and later emitted by the ${}^2f$-mode driven by sustained Ledoux convection.

We investigate the detectability of extended mid-infrared (MIR) emission associated with FU-Ori type objects (FUors) using the METIS coronagraphs on the 39-m Extremely Large Telescope (ELT). The imaging simulations were made for three representative filters ($\lambda$=3.8, 4.8, and 11.3 micron) of the METIS instrument. We demonstrate that the detectability of the extended MIR emission using these coronagraphs is highly dependent on the uncertain nature of the central FUor and its circumstellar environment in various contexts. These contexts are: (A) whether the central radiation source is either a flat self-luminous accretion disk or a star at near-infrared (NIR) wavelengths, (B) the size of the accretion disk for the bright central MIR emission at milliarcsecond scales, (C) whether the extended emission is due to either an optically thick disk or an optically thin envelope, and (D) dust grain models. Observations at $\lambda$=3.8 micron will allow us to detect the extended emission in many cases, while the number of cases with detection may significantly decrease toward longer wavelengths due to the fainter nature of the extended emission and high thermal background noise. In some cases, the presence of a binary companion can significantly hamper detections of the extended MIR emission. NIR and MIR imaging observations at existing 8-m class telescopes, prior to the METIS observations, will be useful for (1) reducing the many model uncertainties and (2) searching for binary companions associated with FUors, therefore determining the best observing strategy using METIS.

Magnetic fields significantly influence the structure of galaxies' interstellar media, but our understanding of magnetic field strengths and structures in external galaxies is severely limited. The Large Magellanic Cloud (LMC) offers a unique opportunity for improvement due to its proximity and large angular size, allowing for various detailed observations, particularly rich datasets of rotation measures and dispersion measures (RM and DM). However, interpreting these measurements is challenging due to the need for assumptions about the 3D structure for which we can only access line-of-sight integrated quantities. To address this, we conduct a suite of high-resolution magnetohydrodynamic (MHD) simulations of the LMC, incorporating star formation, star-by-star feedback, and ram pressure stripping by the Milky Way's circumgalactic medium (CGM), experienced as a circumgalactic wind in the frame of the LMC. Synthetic observations of these simulations allow us to identify parameters that closely match observed RM and DM values. Our best model, which is an excellent match to the real LMC, yields magnetic field strengths of $\sim 1.4~\mu{\rm G}$ (ordered) and $\sim 1.6~\mu{\rm G}$ (turbulent). In this model, Milky Way CGM wind experienced by the LMC plays a critical role in shaping the RM data, with the bulk of the RM signal arising not from the LMC's plane, but from warm, $\sim 10^4$ K, gas in a Reynolds layer region $\sim 1$ kpc off the plane where relatively dense material stripped from the LMC is partially ionised by hard extragalactic radiation fields. This finding suggests that we should be cautious about generalising inferences from the LMC to other galaxies that may not be shaped by similar interactions.

A new potential is presented for spherical galaxies. The technique of the construction of our model is similar to that given by An and Evans. In a special case, its mass density becomes a special one of the Hernquist model. Another special model is primarily discussed, and its intrinsic properties, such as velocity dispersions and surface densities, can be shown directly by numerical calculation. Its distribution functions in both isotropic and anisotropic cases can be expressed as some definite integrals which are suitable for numerical calculation.

Both simulations and observations indicate that the so-called missing baryons reside in the intergalactic medium (IGM) known as the warm-hot intergalactic medium (WHIM). In this article, we demonstrate that turbulence in the cosmic baryonic fluid is crucial for correctly understanding both the spatial distribution and the physical origins of the missing baryons in the universe. First, we find that dynamical effects cause the gas to be detained in low-density and intermediate-density regions, resulting in high baryon fractions, while prevent the inflow of the gas in high-density regions, leading to low baryon fractions. Second, turbulent energy is converted into thermal energy, and the injection and dissipation of turbulent energy have essentially reached a balance from $z=1$ to $0$. This indicates that the cosmic fluid is in a state of fully-developed turbulence within this redshift range. Due to turbulent heating, as redshift decreases, an increasing amount of warm gas is heated and transitions into the WHIM, and some even into hot gas.

Stars form from dense cores in turbulent molecular clouds. According to the standard scenario of star formation, dense cores are created by cloud fragmentation. However, the physical mechanisms driving this process are still not fully understood from an observational standpoint. However, the physical mechanisms driving this process are still not fully understood from an observational standpoint. Our goal is to investigate the process of cloud fragmentation using observational data from nearby clouds. Specifically, we aim to examine the role of self-gravity and turbulence, both of which are key to the dynamical evolution of clouds. We applied astrodendro to the Herschel H2 column density maps to identify dense cores and determine their mass and separation in two nearby low-mass clouds: the Polaris Flare and Lupus I clouds. We then compared the observed core masses and separations with predictions from models of gravitational and turbulent fragmentation. For turbulent fragmentation, the key scales are the cloud sonic scale and its corresponding mass. The average core masses are estimated to be 0.242 Msun for Lupus I and 0.276 Msun for the Polaris Flare. The core separations peak at 0.1 - 0.2 pc in both clouds. These separations are significantly smaller than the Jeans length but agree well with the cloud sonic scale. Additionally, the density probability distribution functions of the dense cores follow log-normal distributions, which is consistent with the predictions of turbulent fragmentation. These findings suggest that the primary process driving core formation in the observed low-mass star-forming regions is not gravitational fragmentation but rather turbulent fragmentation. We found no evidence that filament fragmentation plays a significant role in the formation of dense cores.

Astrometric measurements are significantly challenged by the relative motion between the point source and the telescope, primarily due to the difficulty in accurately determining the position of the point source at the mid-exposure moment. Especially when the trail is irregular in shape or results from nonuniform relative motion, determining the centroid of such a trail becomes significantly more challenging. To address this issue, a new centroiding algorithm for point-source trails has been developed. This algorithm employs a piecewise linear model to approximate the irregular trajectory of a point source. An estimated intensity distribution of the trail is constructed by integrating the point-spread function with the approximated trajectory. The cost function is defined as the difference between the estimated and observed trail intensity distributions, with an added smoothness constraint term. Optimizing this cost function yields a refined trajectory fit. A coarse-to-fine iterative approach is used to progressively converge on the true trajectory of the point source, ultimately determining both the trail's centroid and the trajectory of the point source. The efficacy of the algorithm is validated using synthetic images. Furthermore, this technique is applied to Cassini Imaging Science Subsystem images of several inner Saturnian satellites, successfully processing 267 astrometric observations. The results demonstrate the effectiveness of the algorithm in real astronomical applications.

We study the Kaiser formula for biased tracers in massive neutrino cosmologies by comparing its predictions with a large set of N-body simulations. In particular, we examine the ambiguity in the definition of the peculiar velocity contribution at linear level, whether it should be expressed in terms of the total matter velocity field or only the cold-matter component, comprised of cold dark matter and baryons. We revisit and extend previous qualitative studies on the topic with larger statistics and including a full fit of halo power spectrum and halo-matter cross-power spectrum measurements. We find clear evidence that the clustering of halos in redshift space is correctly described assuming that halo velocity is tracing the velocity field of cold-matter. The opposite assumption provides a worse fit to the simulated data and can lead to a spurious running of nuisance parameters in the perturbative galaxy power spectrum model.

We have developed an optical photon-counting imaging system, IMONY, as an instrument for short-scale time-domain astronomy. In this study, we utilized a Geiger avalanche photodiode array with a $4\times 4$ pixel configuration, with each pixel measuring \SI{100}{\micro m}. We developed a dedicated analog frontend board and constructed a data acquisition system with an FPGA to time-stamp each photon with a time resolution of \SI{100}{\ns}. We mounted a prototype model of the system on the 1.5-m Kanata telescope, intending to observe the Crab pulsar and conduct joint observations with Iitate and Usuda radio telescopes in Japan. We successfully demonstrated that IMONY could image the Crab pulsar as an expected point source and acquire the well-known pulse shape. We found that the time lag between the optical and radio main pulses was $304\pm$\SI{35}{\mu s}, consistent with previous studies.

In this proceeding, we provide a novel approach to study the General Relativistic Magnetohydrodynamic (GRMHD) accretion flows around rotating black holes (BHs). In doing so, we choose a sub-Keplerian distribution of angular momentum of the flow, which is necessary for the accreting matter to reach the event horizon of the BH. Further, we consider the convergent flow to be confined about the disk mid-plane and is threaded by both radial ($b^r$) and toroidal ($b^\phi$) magnetic field components. For simplicity, we neglect any motion along the vertical ($\theta$) direction, maintaining a vertical (hydrostatic) equilibrium about the midplane. With this, we describe the family of multi-trans-magnetosonic accretion solutions around rotating BHs and examine how the conserved magnetic flux ($\Phi$) and BH spin ($a_{\rm k}$) affect the accretion flow properties. Finally, we provide an insight into the thermal emissions from the magnetized disk for a specific set of accretion solutions.

The Dark Energy Survey (DES) recently released the final results of its two principal probes of the expansion history: Type Ia Supernovae (SNe) and Baryonic Acoustic Oscillations (BAO). In this paper, we explore the cosmological implications of these data in combination with external Cosmic Microwave Background (CMB), Big Bang Nucleosynthesis (BBN), and age-of-the-Universe information. The BAO measurement, which is $\sim2\sigma$ away from Planck's $\Lambda$CDM predictions, pushes for low values of $\Omega_{\rm m}$ compared to Planck, in contrast to SN which prefers a higher value than Planck. We identify several tensions among datasets in the $\Lambda$CDM model that cannot be resolved by including either curvature ($k\Lambda$CDM) or a constant dark energy equation of state ($w$CDM). By combining BAO+SN+CMB despite these mild tensions, we obtain $\Omega_k=-5.5^{+4.6}_{-4.2}\times10^{-3}$ in $k\Lambda$CDM, and $w=-0.948^{+0.028}_{-0.027}$ in $w$CDM. If we open the parameter space to $w_0$$w_a$CDM\$ (where the equation of state of dark energy varies as $w(a)=w_0+(1-a)w_a$), all the datasets are mutually more compatible, and we find concordance in the $[w_0>-1,w_a<0]$ quadrant. For DES BAO and SN in combination with Planck-CMB, we find a $3.2\sigma$ deviation from $\Lambda$CDM, with $w_0=-0.673^{+0.098}_{-0.097}$, $w_a = -1.37^{+0.51}_{-0.50}$, a Hubble constant of $H_0=67.81^{+0.96}_{-0.86}$km s$^{-1}$Mpc$^{-1}$, and an abundance of matter of $\Omega_{\rm m}=0.3109^{+0.0086}_{-0.0099}$. For the combination of all the background cosmological probes considered (including CMB $\theta_\star$), we still find a deviation of $2.8\sigma$ from $\Lambda$CDMin the $w_0-w_a$ plane. Assuming a minimal neutrino mass, this work provides further evidence for non-$\Lambda$CDM physics or systematics, which is consistent with recent claims in support of evolving dark energy.

The driving force behind outflows, often invoked to understand the correlation between the supermassive black holes powering active galactic nuclei (AGN) and their host galaxy properties, remains uncertain. We provide new insights into the mechanisms that trigger warm ionized outflows in AGN, based on findings from the MaNGA survey. Our sample comprises 538 AGN with strong [OIII]$\lambda$5007 emission lines, of which 197 are detected in radio and 341 are radio-undetected. We analyzed [OIII]$\lambda$5007 line in summed spectra, extracted over their central 500$\times$500 pc$^2$ region. The calculated Balmer 4000 $Å$ break, D$_n$4000 is larger than 1.45 for $\sim$95$\%$ of the sources indicating that the specific star-formation rate in their central regions is smaller than 10$^{-11.5}$ yr$^{-1}$, pointing to evidence of negative AGN feedback suppressing star-formation. Considering the whole sample, radio-detected sources show an increased outflow detection rate (56$\pm$7\%) compared to radio-undetected sources (25$\pm$3\%). They also show higher velocity, mass outflow rate, outflow power and outflow momentum rate. We noticed a strong correlation between outflow characteristics and bolometric luminosity in both samples, except that the correlation is steeper for the radio-detected sample. Our findings suggest (a) warm ionized outflows are prevalent in all types of AGN, (b) radiation from AGN is the primary driver of these outflows, (c) radio jets are likely to play a secondary role in enhancing the gas kinematics over and above that caused by radiation, and (d) very low star-formation in the central regions of the galaxies, possibly due to negative feedback of AGN activity.

The Soft X-ray Imager (SXI) is an X-ray CCD camera of the Xtend system onboard the X-Ray Imaging and Spectroscopy Mission (XRISM), which was successfully launched on September 7, 2023 (JST). During ground cooling tests of the CCDs in 2020/2021, using the flight-model detector housing, electronic boards, and a mechanical cooler, we encountered an unexpected issue. Anomalous charges appeared outside the imaging area of the CCDs and intruded into the imaging area, causing pulse heights to stick to the maximum value over a wide region. Although this issue has not occurred in subsequent tests or in orbit so far, it could seriously affect the imaging and spectroscopic performance of the SXI if it were to happen in the future. Through experiments with non-flight-model detector components, we successfully reproduced the issue and identified that the anomalous charges intrude via the potential structure created by the charge injection electrode at the top of the imaging area. To prevent anomalous charge intrusion and maintain imaging and spectroscopic performance that satisfies the requirements, even if this issue occurs in orbit, we developed a new CCD driving technique. This technique is different from the normal operation in terms of potential structure and its changes during imaging and charge injection. In this paper, we report an overview of the anomalous charge issue, the related potential structures, the development of the new CCD driving technique to prevent the issue, the imaging and spectroscopic performance of the new technique, and the results of experiments to investigate the cause of anomalous charges.

Observations of hierarchical triple star systems show that misalignments are common both between the angular momentum vector of the inner binary and the outer companion orbit, and between the outer binary orbit and a circumtriple gas disk. With analytic methods and n-body simulations we explore the dynamics of circumtriple orbits around a misaligned hierarchical triple star. Circumtriple test particle orbits nodally precess either about the outer binary angular momentum vector (circulating orbits) or about a stationary inclination that depends upon the binary properties (librating orbits). For a coplanar (or retrograde coplanar) triple star, the apsidal precession rate is maximal and the critical orbital radius outside of which all orbits are circulating is minimal. Polar alignment of a circumtriple gas disk requires nodal libration and therefore it can be more likely if there is a large misalignment between the inner and outer binary orbits. There are two values of the mutual misalignment, i_c and 180-i_c, for which the apsidal precession rate of the triple star is zero and polar alignment is possible at all orbital radii. For a circular inner binary orbit i_c=55, and it changes with eccentricity of the inner binary while being insensitive to other triple star parameters.

Rapid identification of candidates of high-value gamma-ray bursts (GRBs), including both high-$z$ and local events, is crucial for outlining subsequent observational strategy. In this paper, we present a model that enables an on-duty astronomer to rapidly identify candidates of local GRBs prior to spectroscopy, provided that these events have been localized at an arcseconds precision. After taking into account the mass distribution of the host galaxies of GRBs, the model calculates the two-dimensional cross-match probabilities between a localized GRB and its surrounding nearby galaxies, and then returns the best match with the highest probability. The model is evaluated not only by the observed GRB sample with redshifts up to $z=4$, but also through the simulated GRB samples. By using the recently published GLADE+ galaxies catalog with a completeness of 95\% up to 500Mpc, along with the NED-LVS catalog, the Precision and Recall of the model are determined to be 0.23-0.33 and 0.75, respectively, at the best performance. A dedicated web service, which will be integrated into the SVOM Science User Support System, has been developed to deploy the model.

The photon sphere defines the unstable circular orbit of photons in a black hole spacetime. Photons emitted by a source located inside the photon sphere can be gravitationally lensed by the black hole and have time delays when reaching the observer. These delays may lead to light echoes produced in the light curve if an accretion event in the vicinity of the horizon can be observed. In this work, we present fully analytical formulas with high accuracy to describe the change of the azimuthal angle and the travel time of those photons. By employing the analytical approaches, we find that the time delay between photons emitted from the interior of the photon sphere has a typical time scale of $2(\pi - \phi_\mathrm{S}) u_\mathrm{m}$ with $\phi_\mathrm{S}$ and $u_\mathrm{m}$ being respectively the azimuthal angle of the source and the impact parameter evaluated at the photon sphere, which can provide some clues on the future search for gravitational lensing signatures in the accretion inflow event.

This paper presents a comprehensive numerical framework for simulating radiation-plasma systems. The radiative transfer process spans multiple flow regimes due to varying fluid opacity across different regions, necessitating a robust numerical approach. We employ the multiscale unified gas-kinetic scheme (UGKS), which accurately captures photon transport phenomena from free streaming to diffusive wave propagation. The UGKS is also applied to the fluid model to address the significant mass disparity between electrons and ions, and their associated transport characteristics in both equilibrium continuum and non-equilibrium rarefied regimes. Our model explicitly incorporates momentum and energy exchanges between radiation and fluid fields in the coupled system, enabling detailed analysis of the complex interactions between electromagnetic and hydrodynamic phenomena. The developed algorithm successfully reproduces both optically thin and optically thick radiation limits while capturing the complex multiscale nonequilibrium dynamics of the coupled system. This unified treatment eliminates the need for separate numerical schemes in different regimes, providing a consistent and computationally effcient approach for the entire domain. The effectiveness and versatility of this approach are demonstrated through extensive numerical validation across a wide range of physical parameters and flow conditions.

The James Webb Space Telescope (JWST) has observed massive galaxies at high redshifts, which implies an earlier epoch of reionization (EoR) compared with the cosmic microwave background (CMB) results. In this paper, based on \texttt{Planck 2020} (NPIPE release), \texttt{ACT DR4} and \texttt{SPT-3G} data, if assumed a Harrison-Zel'dovich (HZ) primordial power spectrum in the standard cosmological model, we discover that the redshift or optical depth of reionization is larger than the case of a power-law (PL) primordial power spectrum. In HZ-$ \Lambda $CDM model, the redshift of reionization is $ z_\text{reio} = 9.11 \pm 0.61 $, which is consistent with the JWST result that $ z_\text{reio} \approx 8.9 $. Moreover, the cosmological tensions, i.e. Hubble ($H_0$) tension and $ S_8 $ tension are alleviated in HZ-$ \Lambda $CDM case. The Hubble constant is $ H_0 = 70.38 \pm 0.35 \, \text{km}/\text{s}/\text{Mpc}$ and the structure growth parameter is $ S_8 = 0.7645\pm 0.0094 $ in HZ-$ \Lambda $CDM model. We also consider two extensions of $ \Lambda $CDM, including $ \Lambda $CDM$ + A_\text{L} $ and $ \Lambda $CDM$ + \Omega_\text{k} $ models. But the extensions of $ \Lambda $CDM with a HZ spectrum meet more serious CMB anomalies, i.e. lensing anomaly and spatial curvature anomaly as compared with the extensions of $ \Lambda $CDM with a PL spectrum. We discuss that these two CMB anomalies may come from the degeneracy of cosmological parameters.

We derive circular velocity curves (CVCs) from stellar dynamical models for $\sim6{,}000$ nearby galaxies in the final data release of the SDSS-IV MaNGA survey with integral-field spectroscopy, exploring connections between the inner gravitational potential (traced by CVC amplitude/shape) and galaxy properties. The maximum circular velocity ($V_{\rm circ}^{\rm max}$) and circular velocity at the half-light radius ($V_{\rm circ}(R_{\rm e}^{\rm maj})$) both scale linearly with the stellar second velocity moment $\sigma_{\rm e}^2\equiv\langle V^2+\sigma^2\rangle$ within the half-light isophote, following $V_{\rm circ}^{\rm max} \approx 1.75\sigma_{\rm e}$ (8$\%$ error) and $V_{\rm circ}(R_{\rm e}^{\rm maj}) \approx 1.65\sigma_{\rm e}$ (9$\%$ error). CVC shapes (rising, flat, declining) correlate strongly with structural and stellar population properties: declining curves dominate in massive, early-type, bulge-dominated galaxies with old, metal-rich stars and early quenching, while rising CVCs prevail in disk-dominated systems with younger stellar populations and ongoing star formation. Using a unified bulge-disk-halo model, we predict CVC shapes with minimal bias, identifying three governing parameters: bulge-to-total mass ratio ($B/T$), dark matter fraction within $R_{\rm e}$, and bulge Sersic index. The distribution of CVC shapes across the mass-size plane reflects evolutionary pathways driven by (i) in situ star formation (spurring bulge growth) and (ii) dry mergers. This establishes CVC morphology as a diagnostic for galaxy evolution, linking dynamical signatures to structural and stellar population histories.

Europa, Jupiter's second Galilean moon, is believed to host a subsurface ocean in contact with a rocky mantle, where hydrothermal activity may drive the synthesis of organic molecules. Of these molecules, abiotic synthesis of aromatic amino acids is unlikely, and their detection on Europa could be considered a biosignature. Fluorescence from aromatic amino acids, with characteristic emissions in the 200-400 nanometer wavelength range, can be induced by a laser and may be detectable where ocean material has been relatively recently emplaced on Europa's surface, as indicated by geologically young terrain and surface features. However, surface bombardment by charged particles from the Jovian magnetosphere and solar ultraviolet (UV) radiation degrades organic molecules, limiting their longevity. We model radiolysis and photolysis of aromatic amino acids embedded in ice, showing dependencies on hemispheric and latitudinal patterns of charged particle bombardment and ice phase. We demonstrate that biosignatures contained within freshly deposited ice in high-latitude regions on the surface of Europa are detectable using laser-induced UV fluorescence, even from an orbiting spacecraft.

Context. Acetone (CH3COCH3) is a carbonyl-bearing complex organic molecule, yet interstellar observations of acetone remain limited. Studying the formation and distribution of CH3COCH3 in the interstellar medium can provide valuable insights into prebiotic chemistry and the evolution of interstellar molecules. Aims. We explore the spatial distribution of CH3COCH3 and its correlation with the O-bearing molecules acetaldehyde (CH3CHO) and methanol (CH3OH), as well as the N-bearing molecule ethyl cyanide (C2H5CN), in massive protostellar clumps. Methods. We observed 11 massive protostellar clumps using ALMA at 345 GHz, with an angular resolution of 0.7''-1.0''. Spectral line transitions were identified using the eXtended CASA Line Analysis Software Suite. We constructed integrated intensity maps of CH3COCH3, CH3CHO, CH3OH, and C2H5CN and derived their rotation temperatures, column densities, and abundances under the assumption of local thermodynamic equilibrium. Results. CH3COCH3 is detected in 16 line-rich cores from 9 massive protostellar clumps: 12 high-mass cores, 3 intermediate-mass cores, and 1 low-mass core. CH3CHO and CH3OH are also detected in all 16 cores, while C2H5CN is detected in 15. The integrated intensity maps reveal similar spatial distributions for CH3COCH3, CH3CHO, CH3OH, and C2H5CN. The line emission peaks of all four molecules coincide with the continuum emission peaks in regions without ultracompact HII regions. Significant correlations are observed in the abundances of these molecules, which also exhibit similar average temperatures. Conclusions. Our observational results, supported by chemical models, suggest that CH3COCH3, CH3CHO, and CH3OH originate from the same gas. The observed temperatures and abundances of CH3COCH3 are consistent with model predictions involving grain surface chemistry.

The spatial correlation between galaxies and the Ly$\alpha$ forest provides insights into how galaxies reionized the Universe. Here, we present initial results on the spatial cross-correlation between [OIII] emitters and Ly$\alpha$ forest at 5.4<z<6.5 from the JWST ASPIRE NIRCam/F356W Grism Spectroscopic Survey in z>6.5 QSO fields. Using data from five QSO fields, we find $2\sigma$ evidence for excess Ly$\alpha$ forest transmission at ~20-40 cMpc around [OIII] emitters at z=5.86, indicating that [OIII] emitters reside within a highly ionized IGM. At smaller scales, the Ly$\alpha$ forest is preferentially absorbed, suggesting gas overdensities around [OIII] emitters. Comparing with models including THESAN simulations, we interpret the observed cross-correlation as evidence for significant large-scale fluctuations of the IGM and the late end of reionization at z<6, characterized by ionized bubbles over 50 cMpc around [OIII] emitters. The required UV background necessitates an unseen population of faint galaxies around the [OIII] emitters. Furthermore, we find that the number of observed [OIII] emitters near individual transmission spikes is insufficient to sustain reionization in their surroundings, even assuming all [OIII] emitters harbour AGN with 100 % LyC escape fractions. Despite broad agreement, a careful analysis of ASPIRE and THESAN, using the observed host halo mass from the clustering of [OIII] emitters, suggests that the simulations underpredict the observed excess IGM transmission around [OIII] emitters, challenging our model of reionization. Potential solutions include larger ionized bubbles at z<6, more enhanced large-scale UV background or temperature fluctuations of the IGM, and possibly a patchy early onset of reionization at z>10. Current observational errors are dominated by cosmic variance, meaning future analyses of more QSO fields from JWST will improve the results.

We report the discovery and careful orbital determination of 64 new irregular moons of Saturn found in images taken using the Canada-France-Hawaii Telescope from 2019-2021, bringing the total number of saturnian irregulars to 122. By more than doubling the sample of saturnian irregular moon orbits, including pushing to smaller sizes, we can now see finer detail in their orbital distribution. We note the emergence of potential subgroups associated with each of Siarnaq and Kiviuq within the Inuit group. We find that in the inclination range 157-172 degrees the ratio of smaller moons (diameters less than 4 km) to larger moons (diameters greater than 4 km) is significantly larger than that of any other inclination range in the retrogrades. We denote this subset of the Norse group as the Mundilfari subgroup after its largest member. The incredibly steep slope of the Mundilfari subgroup's size distribution, with a differential power law index of q = 6, points to this subgroup being created by a recent catastrophic collision proposed in Ashton et al. (2021).

Hot cores, as a stage of massive star formation, exhibit abundant line emissions of COMs. We present a deep line survey of two isomers of C$_2$H$_6$O: ethanol (C$_2$H$_5$OH; EA), and dimethyl ether (CH$_3$OCH$_3$; DE) as well as their possible precursor CH$_3$OH towards 60 hot cores by using the ALMA 3 mm line observations. EA is detected in 40 hot cores and DE is detected in 59 hot cores. Of these, EA and DE are simultaneously detected in 39 hot cores. We calculate rotation temperatures and column densities of EA and DE by using the XCLASS software. The average rotation temperature of EA is higher than that of DE, whereas the average column density of EA is lower than that of DE. Combined with previous studies of hot cores and hot corinos, we find strong column density correlations among EA and DE ($\rho$ = 0.92), EA and CH$_3$OH ($\rho$ = 0.82), as well as DE and CH$_3$OH ($\rho$ = 0.80). The column density ratios of EA/DE versus the column densities of CH$_3$OH remain nearly constant with values within ~ 1 order of magnitude. These strong correlations and the stable ratios, suggest that EA, DE, and CH$_3$OH could be chemically linked, with CH$_3$OH potentially serving as a precursor for EA and DE. Compared with chemical models, the three different warm-up timescale models result in the systematic overproduction of EA and the systematic underproduction of DE. Therefore, our large sample observations can provide crucial constraints on chemical models.

The outer solar system is populated by a broad aggregate of minor bodies, which occupy orbits whose dynamical character ranges from long-term stable to rapidly diffusive. We investigate the chaotic properties of known distant trans-Neptunian objects (TNOs) by numerically integrating TNO clones and statistically analyzing their orbital diffusion. Comparing the measured diffusion with an analytical criterion yields a dynamically motivated separation into classes of stable, metastable and unstable objects. We then measure the level of clustering of the longitudes of perihelia and of the orbital poles, as functions of orbital distance and of their stability properties. Distant (meta)stable objects appear increasingly clustered in perihelion around $\varpi \sim 50^\circ$ for increasing semi-major axis, while the orbits of unstable objects are well described by two, roughly equally-populated groups of "clustered" and "anti-clustered" objects, with means around $\sim 25^\circ$ and $\sim 205^\circ$ respectively. We further find that, compared to the solar system's total angular momentum vector, the mean orbital poles of distant TNOs are significantly more misaligned for (meta)stable objects, while they remain roughly aligned for unstable objects. TNOs with intermediate orbital periods also appear to be misaligned with respect to the forced plane predicted by secular theory with the known planets. This gradation based on stability, if validated further by the upcoming VRO survey, necessitates a dynamical explanation.

Bilobed bodies represent a significant class of small extraterrestrial objects in the Solar System. We present a double harmonic-series approach to model the gravity of the lobes separately, thereby allowing their mass distributions to be constrained independently. We study an exemplary candidate contact binary, comet 67P/Churyumov-Gerasimenko, and establish two ellipsoidal harmonic series referred to the closest-bounding, non-overlapping ellipsoids of the respective lobes. We show that the model is more robust near the irregular body shape than a global spherical harmonic series, which is a customary product determined from spacecraft tracking. While they can be estimated directly from measurements, we demonstrate here that the double series can also be transformed from an existing spherical harmonic model by means of decomposition. We use simulations to analyze model solutions for a few heterogeneous mass distributions.

The ESA space mission Ariel requires bright sources that are stable at the level of 100ppm over 6 hours in order to accurately measure exoplanet atmospheres through transmission spectroscopy. To ensure this, in-flight instrument calibration can be performed by observing stellar calibrators. In this study, a stellar calibrator candidate list distributed over the sky is created and a flux variability analysis is performed to identify the best stellar calibrators for transit spectroscopy of exoplanet atmospheres with Ariel. A starting candidate sample of 1937 solar-type stars is created using the all-sky surveys Two Micron All Sky Survey and Gaia. Using stellar light curves from the Transit Exoplanet Survey Satellite (TESS), the flux variability of each star is characterised by computing its Lomb-Scargle periodogram and reduced chi-squared. This enables the elimination of stars with detectable variability from the sample. Approximately 22.2% of stars from the starting sample pass the selection as potential calibrators. These do not all necessarily meet Ariel's stability requirement, although some will. No correlation between flux stability and stellar properties is found, as long as the correct value ranges for the parameters are chosen, like a surface temperature between 5000 and 6300K. The only exception is stellar magnitude: Noise in TESS data increases as stars get dimmer, so, a high percentage of faint stars passes the selection since their variability is more likely hidden within the inherent TESS noise. Contrarily, stars brighter than 5mag cannot be used as calibrators. A list of 430 promising bright calibration targets distributed over the sky has been selected. These can potentially be used as stellar calibrators for the Ariel mission. Targets from this list will have to be further studied to determine which ones possess a flux stability better than 100ppm over 6 hours.

In this study, we used the 13.7m telescope at Qinghai Station to observe CO data in 48 galaxies in Virgo clusters, of which 41 sources observed CO signals. The properties of molecular gas are deduced by co-to-$H_2$ factor. We also collected and investigated the relationship between $M_{H_2}$ and other galactic properties ($M_B$, $L_K$, sfr, and $def_{HI}$). We found correlations between $M_{H_2}$ and $M_B$, $L_K$, sfr, and $def_{HI}$. It has the strongest correlation with Lk. There is a certain correlation between MB and MB, but the scattering is larger. For sfr, the more abundant the atomic and molecular gases, the stronger the sfr activity.

Recently, a powerful magnetic field was discovered in the hot helium star classified as a quasi-Wolf-Rayet star of ~2Msun, member of the HD45166 system. Upon its explosion as a core-collapse supernova, it is expected to produce a strongly magnetic neutron star, a magnetar. Among the key parameters governing the pre-supernova evolution is the amount of mass lost via stellar wind. However, the magnetic nature of this helium star is expected to affect its stellar wind making the estimation of the wind parameters uncertain. We report the first observations of HD45166 in X-rays with the XMM-Newton telescope and in radio with the VLA interferometer array. By placing the observation results in a theoretical framework, we aim to provide a reliable estimate of the wind strength of the magnetic qWR star. The X-ray properties are explained in the framework of the MCWS scenario, and the semi-analytic XADM model is applied to reproduce the X-ray emission. The thermal radio emission of the wind and its absorption effect on possible gyro-synchrotron emission from the underlying dipolar magnetosphere, sampled in 3D, are computed by integrating the radiative transfer equation. We did not detect radio emissions, this enabled us to set sensitive upper limits on the radio luminosity. The magnetic qWR star is a slow rotator, comparison with models reveals that the possible acceleration mechanisms occurring within its dynamical magnetosphere are not as efficient as in fast-rotating magnetic ApBp-type stars. In contrast, the system is detected in X-rays with log(L_X/L_bol)~ -5.6. Using suitable models, we constrain the mass lost from this magnetic quasi-Wolf-Rayet star as dot{M}~3e-10 Msun/yr. This novel empirical estimate of the mass-loss rate in a ~2Msun helium star confirms that it maintains super-Chandrasekhar mass till collapse and can produce a magnetar as its final evolutionary product.

Context: Several new phenomena have been surrounding the area of study of the repeating thermonuclear explosions called novae. For example, recurrent novae have been proven to be efficient cosmic ray hadronic accelerators thanks to the recent observations of RS Ophiuchi by different gamma-ray instruments. Novae have also demonstrated to have the ability to carve large cavities into the Interstellar Medium with parallelisms with the remnants of supernovae. We aim at calculating what is the effect of novae in their surrounding media and to which distances these effects dominate over the average quantities that are measured in the ISM. We calculate the filling factor of novae and their contribution to cosmic ray fluxes using cosmic ray propagation codes. To limit what is the atomic density of the Interstellar Medium (ISM) surrounding the region around RS Oph, we use Fermi-LAT observations of the region. The filling factor of novae in the Galaxy is not significant under all assumptions done in the paper. They do not dominate over the local cosmic ray fluxes, even at the lowest energies, for distances larger than a few parsec. The particle density of the ISM surrounding them is, however, very much modified, lowering it more than one order of magnitude with respect to galactic averages, confirming estimates done using other observatories. As a conclusion, even though at global galactic distances, novae do not seem to be dominating cosmic ray transport, they have the power to modify the conditions of their surrounding ISM over parsec distances.

The interstellar medium (ISM) has a number of tracers such as the Na I D 5890, 5896 AA absorption lines that are evident in the spectra of galaxies but also in those of individual astrophysical sources such as stars, novae or quasars. Here, we investigate narrow absorption features in the spectra of nearby supernovae (SNe) and compare them to local (< 0.5 kpc) and global host galaxy properties. With a large and heterogeneous sample of spectra, we are able to recover the known relations of ISM with galaxy properties: larger columns of ISM gas are found in environments that are more massive, more actively star-forming, younger and viewed from a more inclined angle. Most trends are stronger for local than global properties, and we find that the ISM column density decreases exponentially with the offset from the host galaxy centre, as expected for a gas distribution following an exponential radial profile. We also confirm trends for the velocity of galactic outflows increasing with radius. The current study demonstrates the capability of individual light sources to serve as ubiquitous tracers of ISM properties across various environments and galaxies.

Our goal is to find new candidate lambda Boo stars that belong to binary systems. A detailed abundance determination of some candidates could confirm their true lambda Boo nature, while the composition of eventual late-type companions could be used as a proxy for the initial composition of the lambda Boo stars. Results. We obtained a group of 19 newly identified binary systems which contain a candidate lambda Boo star, allowing to duplicate the number of lambda Boo stars currently known in multiple systems. This important group could be used in further studies of lambda Boo stars. For the first time, we performed a detailed abundance analysis of three of these binary systems which includes a candidate lambda Boo star and a late-type companion. We confirmed the true lambda Boo nature of the three early-type stars (HD 87304, HD 98069 and HD 153747), and obtained mostly a solar-like composition for their late-type components. Adopting as a proxy the late-type stars, we showed that the three lambda Boo stars were initially born with a solar-like composition. This is an important constraint for any scenario trying to explain the origin of lambda Boo stars. Also, the solar-like composition of the late-type stars supports the idea that lambda Boo stars are Population I objects, however, we caution that other explanations are also possible. Conclusions. We performed, for the first time, a detailed abundance analysis of binary systems including a lambda Boo star and a late-type companion. We obtained a solid indication that lambda Boo stars born from a solar-like composition and established an important constraint to test formation models of lambda Boo stars.[abridged]

Rayleigh-Lévy flights are simplified cosmological tools which capture certain essential statistical properties of the cosmic density field, including hierarchical structures in higher-order correlations, making them a valuable reference for studying the highly non-linear regime of structure formation. Unlike standard Markovian processes, they exhibit long-range correlations at all orders. Following on recent work on one dimensional flights, this study explores the one-point statistics and Minkowski functionals (density PDF, perimeter, Euler characteristic) of Rayleigh-Lévy flights in two dimensions. We derive the Euler characteristic in the mean field approximation and the density PDF and iso-field perimeter $W_{1}$ in beyond mean field calculations, and validate the results against simulations. The match is excellent throughout, even for fields with large variances, in particular when finite volume effects in the simulations are taken into account and when the calculation is extended beyond the mean field.

This study presents the first comprehensive Bayesian inference of neutron star matter, incorporating $\Delta$-resonances alongside hyperons and nucleons within a density-dependent relativistic hadron (DDRH) framework. Using constraints from nuclear saturation properties, chiral effective field theory ($\chi$EFT), NICER radius measurements, and tidal deformability data from GW170817, we systematically explore the impact of $\Delta$-resonances on the equation of state (EoS) of dense matter and neutron star observables. Our results demonstrate that the inclusion of $\Delta$-baryons softens the EoS at low densities while maintaining sufficient stiffness at high densities to support $2M_{\odot}$ neutron stars. This naturally reconciles neutron star radius constraints with the recent observation of the low-mass compact object in HESS J1731-347 while simultaneously exhibiting excellent agreement with GW170817 tidal deformability constraints, reinforcing the astrophysical viability of $\Delta$-admixed neutron stars. Additionally, $\Delta$-resonances are found to populate the outer layers of the neutron star core, which may have implications for neutron star mergers and their cooling. Furthermore, we show that the presence of $\Delta$-baryons might significantly influence neutron star cooling via the direct Urca process. We also investigate quasi-normal $f$-mode oscillations within a fully general relativistic framework, revealing strong correlations between the $f$-mode frequency, neutron star compactness, and tidal deformability. With the inclusion of $\Delta$-resonances and adherence to astrophysical constraints, we obtain $f_{1.4} = 1.97^{+0.17}_{-0.22}$ kHz and the damping time $\tau_{f_{1.4}} = 0.19^{+0.05}_{-0.03}$ s at the $1\sigma$ confidence level.

Effective field theory (EFT)-based full-shape analysis with simulation-based priors (SBPs) is a novel approach to galaxy clustering data analysis, which significantly boosts the constraining power by efficiently incorporating field-level simulation information from small scales. So far, SBPs have been mostly extracted from a large set of mock catalogs generated with the Halo Occupation Distribution (HOD) approach. We show that given a halo mass function model and assuming that the EFT parameters of halos depend only on the peak height, HOD-based priors can be computed analytically for the standard 7-parameter HOD model. We derive the relevant expressions for deterministic EFT parameters from the halo model. The halo model, however, fails to accurately describe stochastic EFT parameters, for which we use a physically motivated phenomenological prescription. We compare our analytic priors with simulations and find an excellent agreement. Our approach provides analytic insights into the physics behind the EFT parameters of galaxies, and allows one to reduce the computation time of SBPs to virtually nothing. As an application, we produce a set of 100,000 EFT parameters by sampling both HOD and cosmological models, and explicitly demonstrate that any cosmology dependence can be completely absorbed by relatively minor shifts of the HOD parameters. Our approach can be generalized to other variants of the HOD and cosmological models beyond $\Lambda$CDM.

In order to constrain ultra light dark matter models with current and near future weak lensing surveys we need the predictions for the non-linear dark matter power-spectrum. This is commonly extracted from numerical simulations or from using semi-analytical methods. For ultra light dark matter models such numerical simulations are often very expensive due to the need of having a very low force-resolution often limiting them to very small simulation boxes which do not contain very large scales. In this work we take a different approach by relying on fast, approximate $N$-body simulations. In these simulations, axion physics are only included in the initial conditions, allowing us to run a large number of simulations with varying axion and cosmological parameters. From our simulation suite we use machine learning tools to create an emulator for the ratio of the dark matter power-spectrum in mixed axion models - models where dark matter is a combination of CDM and axion - to that of $\Lambda$CDM. The resulting emulator only needs to be combined with existing emulators for $\Lambda$CDM to be able to be used in parameter constraints. We compare the emulator to semi-analytical methods, but a more thorough test to full simulations to verify the true accuracy of this approach is not possible at the present time and is left for future work.

Compact symmetric objects (CSOs) represent a key early stage in radio galaxy evolution, but their reliable identification remains challenging. We develop a method to identify CSOs by combining Gaia optical astrometry with VLBI radio imaging. We analyze 40 CSO candidates by overlaying Gaia DR3 positions on VLBI maps to locate their central engines. CSOs are confirmed when Gaia positions lie between symmetric radio lobes, while core-jet sources show optical positions coinciding with one end of the radio structure. We verify classifications using spectral indices, variability, and jet kinematics from multi-epoch VLBI observations. Our method identified 22 genuine CSOs and 10 core-jet sources, with 8 objects remaining ambiguous. Confirmed CSOs show kinematic ages from 20 to over 1000 years and hotspot speeds typically below 0.5c. Five nearby CSOs show optical-radio offsets despite strong CSO morphology, indicating host galaxy influence. The Gaia-VLBI method provides a reliable CSO identification tool. Our sample reveals diverse radio powers, suggesting multiple evolutionary paths. CSO evolution appears influenced by both intrinsic jet power and environmental factors, with high-power CSOs potentially evolving into large-scale radio galaxies while low-power CSOs often show confinement by their host environments.

We present mid-infrared transmission spectra from 2 to 23 microns of the 23 Atacama Desert chondrites of different types (carbonaceous Ornans and ordinary of H, L, and LL groups) as well as of some pure minerals (olivine and diopside). We focus on the characteristics of silicate at 10 and 20 microns, analyzing the influence of composition and grain size on peak strengths and spectral shapes. We present the first results of the Cosmic Dust Laboratory, a dedicated facility at the Universidad Diego Portales equipped with a VERTEX 80v vacuum Fourier transform infrared spectrometer. Through milling and sieving samples, we obtained different ranges of particle sizes to study the effect of grain size on the intensity and shape of the spectrum. The resulting spectral library can be compared with astronomical data of protoplanetary disks, debris disks, and even white dwarf disks obtained with instruments such as MIRI on board the James Webb Space Telescope and MATISSE on the Very Large Telescope Interferometer. We also present mass absorption coefficient values, which can be used for radiative transfer modeling of astronomical observations. This study aims to improve dust opacities for astronomical applications, with a focus on circumstellar disks.

Here we want to investigate X-shaped radio galaxy 3C 315, which is a FRII source, but it lies very close to the FRI/FRII borderline. We used publicly available data from Leahy's atlas of double radio-sources and NASA/IPAC Extragalactic Database (NED) in order to investigate its flux density, as well as the spectral index distribution. We obtained spectral index distributions between the frequencies: 1417 MHz - 2695 MHz and 1646 MHz - 2695 MHz. Our conclusion is that the synchrotron radiation is the dominant radiation mechanism over most of the area of 3C 315. Because of very poor hotspots (i.e. warmspots) we investigated which part of the source represents more active regions. The results of this study would be helpful to understand the evolutionary process of the galaxy 3C 315.

Ca {\sc i} 4227 Å\, line is a strong resonance line formed in the Solar chromosphere. At the limb, it produces the largest scattering polarization signal. So far, modeling the linear polarization in this line has been limited to the use of one-dimensional semi-empirical models of the solar atmosphere. In this paper, we use three-dimensional magnetohydrodynamical models of the solar atmosphere as well as 1.5D radiative transfer to understand the formation of linear polarization profiles due to resonance scattering in this line at a near limb position. Using three-dimensional magnetohydrodynamical models of the solar atmosphere, in this paper, we perform 1.5D radiative transfer calculations to understand the formation of linear polarization profiles due to resonance scattering in this line at a near limb position. We focus on studying the sensitivity of the resonance scattering polarization to the temperature and the density structures in the atmosphere. We do not include the effects of magnetic and velocity fields in this study. We use clustering analysis to identify linear polarization profiles with similar shape and group them accordingly for our study. We analyze the structure of the linear polarization profiles across 14 clusters, each representing different realizations of the solar atmosphere. Using source function ratio plots at various wing and core wavelength positions, we provide a qualitative explanation of linear polarization profiles in these clusters.

[Edited for arXiv] Source extraction in HI radio surveys is still often performed using visual inspection, but the efficacy of such procedures lacks rigorous quantitative assessment due to their laborious nature. Algorithmic methods are often preferred due to their repeatable results and speed. I here quantitatively assess visual source extraction using a large sample of artificial sources and a comparatively rapid source extraction tool, and compare the results with those from automatic techniques. I injected 4,232 sources into a total of 8,500 emission-free data cubes, with at most one source per cube. Sources covered a wide range of signal-to-noise and velocity widths. I blindly searched all cubes, measuring the completeness and reliability for pairs of signal-to-noise and line width values. Smaller control tests were performed to account for the possible biases in the search, which gave results in good agreement with the main experiment. I also searched cubes injected with artificial sources using algorithmic extractors, and compare these results with a set of catalogues independently reported from real observational data. I find that the results of visual extraction follow a tight relation between integrated signal-to-noise and completeness. Visual extraction compares favourably in efficacy with the algorithmic methods, tending to recover a higher fraction of fainter sources. Visual source extraction can be a surprisingly rapid procedure which gives higher completeness levels than automatic techniques, giving predictable, quantifiable results which are not strongly subject to the whims of the observer. For recovering the faintest features, algorithmic extractors can be competitive with visual inspection but cannot yet out-perform it, though their advantage in speed can be a significant compensating factor.

We explore the implications of finite-temperature quantum field theory effects on cosmological parameters within the framework of the $\Lambda$CDM model and its modification. By incorporating temperature-dependent corrections to the cosmological constant, we extend the standard cosmological model to include additional density parameters, $\Omega_{\Lambda_2}$ and $\Omega_{\Lambda_3}$, which arise from finite-T quantum gravitational effects. Using the Cosmic Linear Anisotropy Solving System, we analyze the impact of these corrections on the cosmic microwave background power spectrum and compare the results with the Planck 2018 data. Through brute-force parameter scans and advanced machine learning techniques, including quartic regression, we demonstrate that the inclusion of $\Omega_{\Lambda_2}$ and $\Omega_{\Lambda_3}$ improves the model's predictive accuracy, achieving high $R^2$ values and low mean squared error. The present work paves the way for future research into higher-order corrections and enhanced computational methods for cosmological parameter estimation.

We report on the validation and characterisation of two transiting planets around TOI-1453, a K-dwarf star in the TESS northern continuous viewing zone. In addition to the TESS data, we used ground-based photometric, spectroscopic, and high-resolution imaging follow-up observations to validate the two planets. We obtained 100 HARPS-N high-resolution spectra over two seasons and used them together with the TESS light curve to constrain the mass, radius, and orbit of each planet. TOI-1453 b is a super-Earth with an orbital period of $P_b$=4.314 days, a radius of $R_b$=1.17$\pm$0.06$R_{\oplus}$, and a mass lower than 2.32$M_{\oplus}$ (99$\%$). TOI-1453 c is a sub-Neptune with a period of $P_c$=6.589 days, radius of $R_c$=2.22$\pm$0.09$R_{\oplus}$, and mass of $M_c$=2.95$\pm$0.84$M_{\oplus}$. The two planets orbit TOI-1453 with a period ratio close to 3/2, although they are not in a mean motion resonance (MMR) state. We did not detect any transit timing variations in our attempt to further constrain the planet masses. TOI-1453 c has a very low bulk density and is one of the least massive sub-Neptunes discovered to date. It is compatible with having either a water-rich composition or a rocky core surrounded by a thick H/He atmosphere. However, we set constraints on the water mass fraction in the envelope according to either a water-rich or water-poor formation scenario. The star TOI-1453 belongs to the Galactic thin disc based on Gaia kinematics and has a sub-solar metallicity. This system is orbited by a fainter stellar companion at a projected distance of about 150 AU, classifying TOI-1453 b and c of S-type planets. These various planetary and stellar characteristics make TOI-1453 a valuable system for understanding the origin of super-Earths and sub-Neptunes.

Decades of progress have culminated in first light for high-energy neutrino astronomy: the identification of the first astrophysical sources of TeV-PeV neutrinos by the IceCube neutrino telescope, the active galactic nuclei NGC 1068 and TXS 0506+056. Today, the prospect of going beyond first light to build high-energy neutrino astronomy in earnest by discovering many more neutrino sources is hampered by the relatively low rate of neutrino detection and the limited view of the sky afforded by IceCube, the single cubic-kilometer-scale neutrino telescope in operation. Yet, this will not stand for much longer. Already today, and over the next 10-20 years, the combined observations of new neutrino telescopes, larger and distributed around the world, will have the potential for transformative progress. Together, they will increase the global rate of neutrino detection by up to 30 times and continuously monitor the entire sky. Within a new joint analysis network - the Planetary Neutrino Monitoring network (PLEnuM) - we make detailed forecasts for the discovery of steady-state astrophysical sources of high-energy neutrinos. We show that a combined analysis of global data will expedite source discovery - in some cases, by decades - and enable the detection of fainter sources anywhere in the sky, discovering up to tens of new neutrino sources.

Polarisation measurements of gamma-ray burst afterglows provide a powerful tool for probing the structure of relativistic jets. In this study, we revisit polarisation signals observed in gamma-ray burst afterglows, focusing on the effects of non-axisymmetric jet structures. To characterize these non-axisymmetric jets, we adopt a simple elliptical jet head model and investigate how deviations from axisymmetry influence the temporal evolution of polarisation properties, particularly around the jet break. Our results show that the polarisation degree curve typically exhibits two peaks for top-hat jets or a single peak for structured jets, even in the presence of an elliptical jet head. In non-axisymmetric jets, a complete drop in polarisation between peaks is generally absent, and the position angle rotation between the peaks can deviate significantly from 90 degrees. In single-peak cases, the polarisation position angle evolves gradually, contrasting with the constant position angle expected in axisymmetric jets. We also explore the implications of these findings for recent GRB events, including GRB 121024A, GRB 091018, GRB 020813, and GRB 210610B.

The James Webb Space Telescope (JWST) has transformed our understanding of early galaxy formation, providing an unprecedented view of the first billion years of cosmic history. These observations offer a crucial opportunity to probe the interplay between galaxy formation and reionization, placing stringent constraints on theoretical models. In this work, we build upon our previously developed semi-analytical framework that self-consistently models the evolving UV luminosity function (UVLF) of galaxies and the global reionization history while incorporating the effects of radiative feedback. Comparing our predictions with JWST and HST data, we identify a fundamental tension: models that match the UVLF fail to reproduce the observed evolution of galaxy clustering (bias) with redshift, and vice versa. To resolve this, we introduce a mass-dependent duty cycle linked to star formation duration. This duty cycle approaches unity at $z > 11$, requiring either enhanced UV radiation production or increased star formation efficiency to match the JWST UVLFs, while declining towards lower redshifts ($5 < z \leq 9$) to remain consistent with the bias measurements. Reconciling theory with observations requires that the characteristic star formation timescale increases from $\approx 80$ Myr at $z \approx 6$ to $\approx 120$ Myr at $z \approx 8$. Finally, our extended model, assuming a halo mass-independent escape fraction of $\approx 12\%$, produces a reionization history consistent with current constraints. These findings underscore the importance of jointly constraining high-redshift galaxy models using both UVLF and bias statistics to accurately interpret JWST data and refine our understanding of early cosmic evolution.

We show that it is possible to simulate realistic inhomogeneities during cosmological inflation with high precision using numerical relativity. Stochastic initial conditions are set in line with the Bunch-Davies vaccuum and satisfy the Hamiltonian and Momentum constraints of General Relativity to leading order in perturbation theory. The subsequent fully non-linear dynamical evolution is formulated within a family of geodesic gauges but can in principle be adapted to any choice of coordinates. We present 3 examples of inflationary dynamics: a simple quadratic potential, a potential with an inflection point and a strong resonance model. When perturbations are small, we recover standard predictions of cosmological perturbation theory, and we quantify strongly non-linear inhomogeneities when non-perturbative configurations emerge, such as in the strong resonance model. Our results pave the way towards the first realistic non-perturbative, and fully non-linear Numerical Relativity simulations of the early inflationary universe.

Metastable cosmic strings (MSCSs) are among the best-fitting explanations of the 2023 pulsar timing array (PTA) signal for gravitational waves at nanohertz frequencies. We propose the novel possibility that a network of MSCSs generating this signal originates from the multi-step spontaneous breaking of a gauged flavour symmetry. As a specific example, we construct a model of $SU(2)$ flavour symmetry in the context of $SU(5)$ grand unification, where the $SU(2)$ acts exclusively on the first two generations of the matter 10-plet, such that it is ``right for leptons'' and allows for large lepton mixing. The model explains the mass hierarchies of the Standard Model fermions, and predicts the string scale of the MSCSs in a range compatible with the 2023 PTA signal. Cosmic inflation is associated with the latter step of (two-step) family symmetry breaking, and the phase transition ending inflation generates the cosmic string network.

Nucleation in the supercooled Yukawa system is relevant for addressing current challenges in understanding a range of crystallizing systems including white dwarf (WD) stars. We use both brute force and seeded molecular dynamics simulations to study homogeneous nucleation of crystals from supercooled Yukawa liquids. With our improved approach to seeded simulations, we obtain quantitative predictions of the crystal nucleation rate and cluster size distributions as a function of temperature and screening length. These quantitative results show trends towards fast nucleation with short-ranged potentials. They also indicate that for temperatures $T > 0.9T_m$, where $T_m$ is the melt temperature, classical homogeneous nucleation is too slow to initiate crystallization but transient clusters of around 100 particles should be common. We apply these general results to a typical WD model and obtain a delay of approximately 0.6 Gyr in the onset of crystallization that may be observable.

A strongly self-interacting component of asymmetric dark matter can collapse and form compact objects, provided there is an efficient mechanism of energy evacuation. If the dark matter quantum number is not completely conserved but it is slightly violated due to some new physics e.g. at the Planck scale, dark matter particles can annihilate into Standard Model particles. Even tiny annihilation cross sections are sufficient to create observable luminosities. We demonstrate that these dark matter annihilations can trigger radial pulsations, causing a characteristic time modulation of the luminosities produced. We argue that in this scenario, the spectral features along with the properties of the oscillation can create a unique discovery signal for such objects in the sky.

Phase transitions of matter under changes of external environment such as temperature and magnetic field have attracted great interests to various quantum many-body systems. Several phase transitions must have occurred in neutron stars as well such as transitions from normal to superfluid/superconducting phases and crust formation. In this work, we extend the superfluid band theory, which has been formulated in our previous work [K. Yoshimura and K. Sekizawa, Phys. Rev. C 109, 065804 (2024)] based on the Kohn-Sham density functional theory (DFT) for superfluid systems, into the finite temperature and finite magnetic field systems. As a result of the finite temperature calculations, we find that the superfluidity of neutrons dissapears at around $k_\text{B}T=0.6$--$0.9\,$ MeV, and ``melting'' of nuclear slabs, that is, a structural change into the uniform matter, takes place at around $k_\text{B}T=2.5$--$4.5\,$ MeV. We also reveal that these transition temperatures exhibit a systematical dependence on the baryon densities. By turning on the magnetic field, we find that protons' spin gets polarized at around $B=10^{16}\,$G, whereas neutrons' spin is kept unpolarized on average up to around $B=10^{17}\,$G. Intriguingly, our microscopic calculations reveal that neutrons' spin is actually polarized locally inside and outside of the slab already at $B\sim10^{16}\,$G, while keeping the system unpolarized in total. As a conclusion, we have demonstrated validity and usefulness of the fully self-consistent superfluid nuclear band theory for describing neutron star matter under arbitrary temperature and magnetic field. Critical temperatures and magnetic fields have been predicted for 1) superfluid to normal transition, 2) crust formation, and 3) spin polarization, under conditions relevant to realistic neutron star environments.

We investigate the divergence-free parametric form of the deceleration parameter within the simplest non-minimal matter-geometry coupling in $f(R,T)$ gravity, where $R$ is the Ricci scalar and $T$ is the trace of the energy-momentum tensor. Specifically, we consider the linear model $f(R,T) = R + 2\lambda T$, where $\lambda$ governs the interaction between matter and geometry. Using this parametric form, we derive the Hubble parameter as a function of redshift $z$ and incorporate it into the modified Friedmann equations. Constraining the model with OHD and Pantheon data, we obtain precise estimates for $H_0$, the present deceleration parameter $q_0$, and its evolutionary component $q_1$, confirming a smooth transition between cosmic deceleration and acceleration. Further, we analyze the evolution of the energy density $\rho$ and total EoS parameter $\omega$ for different $\lambda$ values, highlighting deviations from $\Lambda$CDM and the role of $\lambda$ in shaping cosmic dynamics. In addition, we examine energy conditions, finding that the NEC and DEC are satisfied throughout evolution, while the SEC is violated at late times, supporting the observed acceleration. Our findings demonstrate that this divergence-free parameterization within $f(R,T)$ gravity offers a viable framework for explaining late-time cosmic acceleration while maintaining key observational and theoretical constraints.

We derive analytically some general features of the power-law sensitivity curve. They include an exact parametric equation, a formula for the peak sensitivity and a proof of convexity in log-log plot. A few conceptual points are also clarified.

We show that in the Starobinsky inflation model stochastic gravitational waves are produced when the scalaron - which is the massive scalar mode of the metric - decays into gravitons during reheating. This decay is accompanied by decay of scalaron into matter as well through a similar coupling, proving an efficient reheating stage. The stochastic gravitational waves thus produced have characteristic strain $h_c\sim 10^{-35}-10^{-34}$ in the frequency range $10^{5}-10^{12}\, {\rm Hz}$ which makes them accessible to resonant cavity searches for graviton to photon conversions. Their detection could conclusively validate the Starobinsky inflation model.

We present a new rotating black hole solution to the Einstein equations as an extension of the Kerr spacetime. To derive this solution, we use the Newman-Janis algorithm as a mathematical tool that reduces a general rotating metric into a tractable form by applying simple physical requirements. Interestingly, the solution we found may not be uniquely characterized by asymptotic parameters such as mass, angular momentum, and charge, thereby challenging the no-hair theorem. We also analyze in detail how this additional characteristics (``hair") affects the thermodynamic properties of the black hole.

The iconic, deep-MOND-limit (DML) relation between acceleration and mass, $a\sim (M\mathcal{A}_0)^{1/2}/r$, implies that, in MOND, accelerations cannot be linear in the mass distribution ($\mathcal{A}_0\equiv Ga_0$ is the DML constant, and $a_0$ the MOND acceleration). This leads to important idiosyncracies of MOND, such as a breakdown of the strong equivalence principle, and the resulting ``external-field effect''. I show that the DML axioms are, in themselves, consistent with a, possibly unique, nonrelativistic, action-based, linear formulation of the DML. This model suffers from important drawbacks, which may make it unacceptable as a basis for a full-fledged MOND theory. The model is unique among MOND theories propounded to date not only in being linear -- hence not exhibiting an external-field effect, for example -- but in constituting a modification of both Newtonian inertia and Newtonian gravity. This linear and time-local model inspires and begets several, one-parameter families of models. One family employs nonlinear, time-nonlocal kinetic terms, but still linear gravitational-field equations. Other families generalize the DMLs of AQUAL and QUMOND, modifying gravity as well as inertia. All families employ fractional time derivatives and possibly fractional Laplacians. At present, I cannot base some acceptable MOND theory on these models -- for example, I cannot offer a sensible umbrella theory that interpolates between these DML models and Newtonian dynamics. They are, however, quite useful in elucidating various matter-of-principle aspects of MOND; e.g., they help to understand which predictions follow from only the basic tenets of MOND -- so-called primary predictions -- and which are secondary, i.e., theory dependent. The models may also show the way to a wider class of MOND theories. (Abridged.)

We apply statistical analysis to search for processes responsible for turbulence in physical systems. In our previous studies, we have shown that solar wind turbulence in the inertial range of large magnetohydrodynamic scales exhibits Markov properties. We have recently extended this approach on much smaller kinetic scales. Here we are testing for the Markovian character of stochastic processes in a kinetic regime based on magnetic field and velocity fluctuations in the solar wind, measured onboard the Magnetospheric Multiscale (MMS) mission: behind the bow shock, inside the magnetosheath, and near the magnetopause. We have verified that the Chapman-Kolmogorov necessary conditions for Markov processes is satisfied for local transfer of energy between the magnetic and velocity fields also on kinetic scales. We have confirmed that for magnetic fluctuations, the first Kramers-Moyal coefficient is linear, while the second term is quadratic, corresponding to drift and diffusion processes in the resulting Fokker-Planck equation. It means that magnetic self-similar turbulence is described by generalized Ornstein-Uhlenbeck processes. We show that for the magnetic case, the Fokker-Planck equation leads to the probability density functions of the kappa distributions, which exhibit global universal scale invariance with a linear scaling and lack of intermittency. On the contrary, for velocity fluctuations, higher order Kramers-Moyal coefficients should be taken into account and hence scale invariance is not observed. However, the nonextensity parameter in Tsallis entropy provides a robust measure of the departure of the system from equilibrium. The obtained results are important for a better understanding of the physical mechanism governing turbulent systems in space and laboratory.

The current LAGEOS-LARES 2 experiment aims to accurately measure the general relativistic Lense-Thirring effect in the gravitomagnetic field of the spinning Earth generated by the latter's angular momentum $\boldsymbol{J}$. The key quantity to a priori analytically assess the overall systematic uncertainty is the ratio $\mathcal{R}^{J_2}$ of the sum of the classical precessions of the satellites' nodes $\Omega$ induced by the Earth's oblateness $J_2$ to the sum of their post-Newtonian counterparts. $In$ $principle$, if the sum of the inclinations $I$ of both satellites were $exactly$ $180^\circ$, the semimajor axes $a$ and the eccentricities $e$ being $identical$, $\mathcal{R}^{J_2}$ would $exactly$ vanish. Actually, it is $not$ so by a large amount because of the departures of the $real$ satellites' orbital configurations from their $ideal$ ones. Thus, $J_2$ impacts not only directly through its own uncertainty, but also $indirectly$ through the errors in all the other physical and orbital parameters entering $\mathcal{R}^{J_2}$. The consequences of this fact are examined in greater details than done so far in the literature. The Van Patten and Everitt's proposal in 1976 of looking at the sum of the node precessions of two counter-orbiting spacecraft in (low-altitude) circular polar orbits is revamped rebranding it POLAr RElativity Satellites (POLARES). (Abridged)